我的数据子集有一些问题(子集仍然有 600 个值)。对于实验,我有两个时间点,每个时间点嵌套了三个处理 (TT),其中嵌套了 5 个重复培养物(A 到 E)。在这些文化中的每一种中,都有 20 个个体有机体的值。对于这个子集,我想查看治疗之间的差异以及同一治疗中时间点之间的差异。
我正在使用 R 3.1.2 和 nlme 包
我的代码如下:
model4a <-lme(Velocity~WeekTT, random=~1|Week/TT/Culture, method = "REML", data=body2, weights = varIdent(form=~1|WeekCulture), control=lmeControl(opt="optim"))
STR如下:
data.frame': 600 obs. of 4 variables:
$ Culture : Factor w/ 5 levels "A","B","C","D",..: 1 2 3 4 5 1 2 3 4 5
$ Treatment: Factor w/ 3 levels "T20","T25","T25F": 1 1 1 1 1 2 2 2 2 2
$ Week : num 10 10 10 10 10 10 10 10 10 10 ...
$ Velocity : num 259 279 265 275 256 ...
这是我得到的结果和底部错误的屏幕截图(正在发布一些可重现的数据,但认为这可能是一个简单的代码错误)。
在过去的 3 周里,我一直在研究不同的模型,我认为这是一个简单的问题,只是我的大脑已经疲惫不堪,而且我想多了。