我的研究是一个为期12周的研究,涉及两种饮食(饮食)。在时间(时间)= 0,6和12周的间隔进行结果测量,例如体重,腰围等。 (代码)代表主题。召回代表饮食召回问卷的数量;有3。Kcal是从给定召回中收集的热量摄入。
这就是数据框的样子
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时“权重”? >
感谢您的帮助,并抽出宝贵的时间阅读这篇文章。