Question

我有大约100个来自MCMCglmm的模型，它们的输出类似于：

> MC1 <- MCMCglmm(...)
> summary(MC1)

 Iterations = 20001:99991
 Thinning interval  = 10
 Sample size  = 8000 

 DIC: 10924.52 

 G-structure:  ~school

         post.mean l-95% CI u-95% CI eff.samp
school    0.1753   0.1059   0.2554     1529

 R-structure:  ~units

  post.mean l-95% CI u-95% CI eff.samp
units         1        1        1        0

 Location effects: bl ~ +bm1 + bm2 + bm3 + bm4 + bm5 + bm6 + bm7 + bm8

                   post.mean  l-95% CI  u-95% CI eff.samp  pMCMC    
(Intercept)        -7.381791 -8.105984 -6.657302     1212 <1e-04 ***
bm1                 0.062922  0.063024  0.078611     1028 <1e-04 ***
bm2                -0.015807 -0.016998 -0.019064     1732 <1e-04 ***
bm3                 0.005978  0.003845  0.000207     2124 <1e-04 ***
bm4                 0.223856  0.105453  0.342821     1999 <1e-04 ***
bm5                 0.044622 -0.394179  0.523072     1758   0.88    
bm6                 3.899881  3.672857  3.997223     2976 <1e-04 ***
bm7                 0.547813  0.341128  0.749568     2916 <1e-04 ***
bm8                 0.658511  0.541424  0.783192     2196 <1e-04 ***
---

现在我需要将截取和bm1-bm8的固定效果表放到数据帧中（pvalues除外 - 我对那些不感兴趣）。任何人都可以帮助我这样做，以便可以使用更多相同类型的模型轻松复制它吗？

Answer 1

通过检查str(summary(MC1))，我们发现此信息包含在summary(MC1)$solutions中，其中p值位于第五列。因此，您可以使用

summary(MC1)$solutions[,-5]

操纵摘要（模型）的结果

1 个答案: