Криволинейные участки с 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).