> 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,在预测的时候怎么带入运算。