带有sjplot的曲线图

Question

我有一个混合模型(使用nlme)和时间和时间平方(time2)。我怎样才能让sjp.int (of sjplot)要求绘制曲线呢？

   model = lme(MR ~ Age* time + sex*time + ed* time+ time* HMD+time2*HMD, random=~1|subject, na.action="na.omit",data=mydata);summary(model)

我将sjplot调用为：

  sjp.int(model, type="eff", swap.pred=T, mdrt.values="quart", show.ci=T)

Answer 1

您的二次项应该具有相同的名称，即您应该使用poly()或I()^2。请参见以下示例：

library(sjPlot)
library(sjmisc)
data(efc)
fit <- lm(tot_sc_e ~ c161sex * e17age + c161sex * I(e17age^2), data = efc)
sjp.int(fit)

​

​

请注意，交互作用项的预测器的顺序很重要，并且在公式中应该是相同的(即在上面的情况下，先是c161sex，然后是e17age)。