0

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

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 之间没有比较？无论如何要指定这个吗？

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

4

0 回答 0