r - 带有 sjplot 的曲线图

Question

我有一个时间和时间二次方（time2）的混合模型（使用 nlme）。我怎样才能让 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()or 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）。