Криволинейные участки с sjplot

Question

У меня есть смешанная модель (с использованием 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)

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).