现在在做混合模型:
model.fit1 <- lme(TOTAL ~ DBH , random = ~ DBH | region, data=rd)
> summary(model.fit1)
Fixed effects: TOTAL ~ DBH
Value Std.Error DF t-value p-value
(Intercept) -2.110174 0.20635881 242 -10.22575 <2e-16 ***
DBH 2.419176 0.07943456 242 30.45495 <2e-16 ***
> random.effects(model.fit1)
(Intercept) DBH
1 0.2580542 -0.1199950
2 0.2626176 -0.1408667
3 -0.1931359 0.1010073
4 -0.5674468 0.2393904
5 0.6746693 -0.2188693
6 -0.4347583 0.1393333
现在的问题是怎么计算出random.effects中每一个 (Intercept) 和 DBH的P值 (p-value)