Question

我在R中使用nlme包运行混合效应模型。我的数据包括82只动物（重复），这82只动物按3个品种（定义为分类变量）分组，我的连续变量是时间平方，回应是我的。混合模型的形式是：

model<-lme(MY~time + breed*time + time-squared,random=1|Animal, data=mydata)

结果如下：

Linear mixed-effects model fit by REML
 Data: na.omit(centred_phuong) 
       AIC      BIC    logLik
  93698.27 93769.46 -46840.13
Random effects:
 Formula: ~1 | Cow_code
        (Intercept) Residual
StdDev:    1.283306 2.453689
Fixed effects: MY ~ time + time * Breed + time_squared 
                         Value Std.Error    DF   t-value p-value
(Intercept)          1.3398578 0.2495665 20048   5.36874  0.0000
DFC                  0.0523516 0.0003675 20048 142.43468  0.0000
Breed2               0.4998856 0.3521235    78   1.41963  0.1597
Breed3              -0.3683371 0.3520760    78  -1.04619  0.2987
Time_squared         0.0001213 0.0000025 20048  48.14845  0.0000
Time:Breed2          0.0084160 0.0005011 20048  16.79463  0.0000
Time:Breed3         -0.0086297 0.0005028 20048 -17.16272  0.0000

正如我从这些结果中所理解的，该模型将品种2和3的截距和斜率与品种1的截距和斜率进行比较，并且品种2和品种2之间没有比较？无论如何要指明这个吗？

我希望能得到你的建议！感谢

与nlme成对比较

0 个答案: