操纵摘要(模型)的结果

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)$solutionsp值中,在第五列中.因此,您可以使用

summary(MC1)$solutions[,-5]