эммип (emmeans) с квадратичным термином

Question

Можно ли с помощью emmip построить график предельных (log odds) средних из модели geeglm, если у вас есть квадратичный термин?Я повторил данные измерений, и модель лучше согласуется с термином обработки x в квадрате в дополнение к члену взаимодействия с линейным временем.

Я просто хочу иметь возможность визуализировать предсказанную кривую в данных.Если это возможно, я не знаю, как это указать.Я пробовал:

mod3 <- geeglm(outcome ~ treatment*time + treatment*time_sq, data = dat, id = id, family = "binomial", corstr = "exchangeable"))
mod3a.rg <- ref_grid(mod3, at = list(time = c(1,2,3,4,5,6), time_sq = c(1,4,9,16,25,36)))
emmip(mod3a.rg, treatment ~ time)

Answer 1

Я не думаю, что ваш mod3 правильно включает ваш квадратный термин (трудно сказать, поскольку вы не включили воспроизводимый код).Это позволит вам правильно включить квадрат для времени:

mod3 <- geeglm(outcome ~ treatment*time + treatment*I(time^2), data =
dat, id = id, family = "binomial", corstr = "exchangeable"))

Добавление plotit = TRUE к вашему звонку на emmip():

emmip(mod3a.rg, treatment ~ time, plotit = TRUE)

Вот простой воспроизводимый пример снабор сберегательных данных в MASS, отдаленном пакете для сравнениятермин с I() ваше резюме будет включать только термин для ddpi.

mod2 <- lm(sr ~ ddpi+ddpi^2)
summary(mod2)

дает следующее резюме с коэффициентом только для ddpi

> lm(formula = sr ~ ddpi + ddpi^2, data = savings)
> 
> Residuals:
>     Min      1Q  Median      3Q     Max 
> -8.5535 -3.7349  0.9835  2.7720  9.3104 
> 
> Coefficients:
>             Estimate Std. Error t value Pr(>|t|)     
>(Intercept)   7.8830     1.0110   7.797 4.46e-10 *** 
>ddpi          0.4758     0.2146   2.217   0.0314 *  
> --- Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> 
> Residual standard error: 4.311 on 48 degrees of freedom Multiple
> R-squared:  0.0929,   Adjusted R-squared:  0.074  F-statistic: 4.916 on
> 1 and 48 DF,  p-value: 0.03139