Вывод моей модели lmer (model1) в r следующий (время2 - квадрат времени, время - недели):
Фиксированные эффекты:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 13.592364 3.248658 165.110280 4.184 4.64e-05
time 2.077864 0.176357 266.696826 11.782 < 2e-16
time2 -0.027597 0.003191 240.196072 -8.648 7.59e-16
Drug -0.069304 4.283717 128.278880 -0.016 0.987
time:Drug 0.159494 0.096728 99.385047 1.649 0.102
Если бы я хотел рассчитать оценки для группы Наркотик = 0 против Наркотика = 1 для 0-44 недель (как смоделировано), я мог бы сделать:
13.592364 + 2.077864 * 44 - 0,027597 * 44 * 44 - 0,069304 * 0 + 0,159494 * 0 * 44 = 51,59
13,592364 + 2,077864 * 44 - 0,027597 * 44 * 44 - 0,069304 * 1 + 0,159494 * 1 * 44 = 58,54
Тем не менее, как я могу оценить / извлечь стандартные ошибки этих оценок продолжительности времени?
Спасибо!