Моим исследованием было 12-недельное исследование, включающее две диеты (диета). Измерения результата, такие как вес, окружность талии и т. Д., Проводились во время (время) = 0,6 и 12-недельные интервалы. (Код) обозначает предметы. Напомним обозначает номер диетического отзыва анкеты; Есть 3. Калорий - калорийность рациона, полученная из данного отзыва
Вот как выглядит фрейм данных
code treatment time weight recall kcal
1 a 1 80 1 1500
1 a 2 78 2 1400
1 a 3 75 3 1450
2 b 1 82 1 1600
2 b 2 81 2 1700
2 b 3 81 3 1650
3 a 1 95 1 1200
3 a 2 90 2 1100
3 a 3 85 3 1150
4 b 1 97 1 1900
4 b 2 95 2 1850
4 b 3 96 3 1875
Я хотел бы создать модель, в которой я использую вес в момент времени = 1 в качестве фиксированного эффекта. И вес во время = 2 и время = 3 в качестве результата.
Я узнал здесь вчера, как делать вес в момент времени = 1 как фиксированный эффект. Я использовал следующий код.
mydata <- mydata[order(mydata$code, mydata$time), ]
mydata$weight1 <- with(my data, ave(weight, code, FUN = function (x) x[1]))
Это моя модель
lmer(weight ~ treatment * time + weight1 + weight1 * time + (time | code), REML = FALSE, data = mydata)
Однако в этой модели «вес» результата - это результат во время = 1,2 и 3. Как ограничить этот результат только «весом» при исходных значениях времени = 2 и времени = 3?
Спасибо за вашу помощь и нашли время, чтобы прочитать этот пост.