使用ReporteRs输出lmer汇总信息

时间:2015-04-27 05:37:07

标签: r ms-word lmer reporters

我正在使用ReporteRs中的R个包将统计测试结果和图表输出到Word 2010

我在lmer包中使用lme4进行一些分析,我想将summary(foo.lmer)的结果输出到Word。我希望捕获lmerMod对象中的所有内容并将其放入Word中。例如,这里有一个summary(foo.lmer)结果:

Linear mixed model fit by REML ['lmerMod']
Formula: value ~ Drug * Time + (Time + 1 | ID)
Data: DWeight.Female

REML criterion at convergence: 2296.9

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.3165 -0.6132  0.0343  0.6476  2.3088 

Random effects:
Groups   Name        Variance  Std.Dev. Corr
    ID   (Intercept) 573.98186 23.9579      
         Time          0.07738  0.2782  0.29
Residual              15.61168  3.9512      
Number of obs: 360, groups:  ID, 40

Fixed effects:
                 Estimate Std. Error t value
(Intercept)     270.44441    5.83694   46.33
Drug              1.24479    4.11450    0.30
Time              1.30481    0.07641   17.08
Drug:Time        -0.11055    0.05386   -2.05

Correlation of Fixed Effects:
        (Intr)   Drug  Time  
Drug    -0.758              
Time     0.211 -0.160       
Drug:Tm -0.160  0.211 -0.758

如何使用固定宽度和硬回车存储此信息,以便我可以使用ReporteRs输出此信息?我可以复制和粘贴,但我正在尝试自动化输出,我经常重写输出Word文档。

固定效果输出是最重要的部分,但将所有输出复制到Word中非常有用。

更新和解决方案 感谢在回复中使用capture.output的建议,这段代码对我有用,仍然使用ReporteRs:

 D1male.Results <- capture.output(summary(D1Male.lmer))
 myBWdoc = docx()
 myBWdoc = addTitle( myBWdoc, "Drug 1 Results", level = 1 )
 myBWdoc = addTitle( myBWdoc, "Repeated Measures Models", level = 2 )
 myBWdoc = addParagraph( myBWdoc, D1male.Results)

如果你想把你的回答作为答案我可以接受并提出它。

1 个答案:

答案 0 :(得分:2)

尝试使用capture.output()函数。像这样:

Results <- capture.output(summary(your.lmer.model))

它将摘要信息存储在新变量中。之后,您可以在例如ReporteRs段落功能。