> glmmY2<-glmmPQL(Y1~0+b1+b2+b3+b4+b5+b6+b7+b8+b9+b10+cy,random=~1|ay,
+ family=Gamma(link="log"),data=data1)
iteration 1
iteration 2
iteration 3
> summary(glmmY2)
Linear mixed-effects model fit by maximum likelihood
Data: data1
AIC BIC logLik
NA NA NA
Random effects:
Formula: ~1 | ay
(Intercept) Residual
StdDev: 0.1526869 0.477666
Variance function:
Structure: fixed weights
Formula: ~invwt
Fixed effects: Y1 ~ 0 + b1 + b2 + b3 + b4 + b5 + b6 + b7 + b8 + b9 + b10 + cy
Value Std.Error DF t-value p-value
b1 19.270985 0.2531549 35 76.12329 0.0000
b2 18.652922 0.2756090 35 67.67893 0.0000
b3 16.681047 0.3005444 35 55.50277 0.0000
b4 15.126077 0.3280902 35 46.10341 0.0000
b5 15.191968 0.3586030 35 42.36431 0.0000
b6 13.692755 0.3928639 35 34.85369 0.0000
b7 13.830753 0.4325537 35 31.97465 0.0000
b8 11.458363 0.4816133 35 23.79163 0.0000
b9 10.173265 0.5513203 35 18.45255 0.0000
b10 7.415907 0.6914696 35 10.72485 0.0000
cy 0.116504 0.0401557 35 2.90131 0.0064
我问一下random effects部分没有参数估计值,而是StdDev,在预测的时候怎么带入运算。