曲线图与sjplot

时间:2016-09-20 23:54:37

标签: r plot mixed-models quadratic sjplot

我有一个混合模型(使用nlme),时间和时间是二次的(time2)。我怎样才能让sjp.int(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)

1 个答案:

答案 0 :(得分:0)

您的二次项应具有相同的名称,即您应该使用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)

enter image description here

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