问glmm包如何使用？

提问于 2024-11-21 13:01:11

回答 0关注 0查看 23

> c1<-makeCluster(4)
> binary_China_data$binary_anemia<-as.numeric(binary_China_data$binary_anemia)
> fit<-glmm(binary_anemia~0+PM2.5+SEX+NA13+NA14+nA4b+AA6+BG9+group_age,
+           random=list(binary_anemia~0+pro),varcomps.names=c("pro"),data=China_data,m=100,
+           family.glmm= bernoulli.glmm,doPQL=TRUE,debug=TRUE,cluster=c1)
> stopCluster(c1)
> summary(fit)

Fixed Effects:
           Estimate Std. Error z value Pr(>|z|)    
PM2.5             0          0  -1.326  0.18474    
SEX1              0          0   0.629  0.52932    
SEX2              0          0   0.359  0.71962    
NA131             0          0  -0.922  0.35644    
NA141             0          0   2.381  0.01728 *  
nA4b1             0          0  -1.887  0.05915 .  
AA61              0          0   2.750  0.00596 ** 
BG91              0          0   0.993  0.32090    
group_age1        0          0  -6.646 3.01e-11 ***
group_age2        0          0  -8.260  < 2e-16 ***
group_age3        0          0  -1.912  0.05584 .

> coef(fit)
       PM2.5         SEX1         SEX2        NA131        NA141        nA4b1         AA61         BG91 
 0.008826033 -0.669432525 -0.814801262  0.017406494  0.309346744 -0.386803903  0.397771561  0.092654976 
  group_age1   group_age2   group_age3   group_age4 
-0.687794855 -1.466765284 -2.022706397 -2.680593241

我是用glmm这个包进行混合效应模型分析，但是为什么我的summary（fit）里面固定效应的的系数都为0，但是运行这个包中的例子不会出现这个问题。可以通过coef（fit）查看系数，另外这个包运行很慢，还不如使用lme4？

lme4