4

R中是否有可用的工具来生成可发布的回归表?我正在写一篇课程论文,我需要在其中比较几个回归模型,如果我可以让它们嵌套在像这样的一个表中,我会很高兴,来自estoutStata 包。

我已经检查过xtable,但无法达到相同的结果。任何提示将不胜感激。

这是我的想法:

多重嵌套回归

4

3 回答 3

3

您可能想要mtable“memisc”包中的功能。它具有关联的 LaTeX 输出参数:

==========================================================================
                                              Model 1   Model 2   Model 3 
--------------------------------------------------------------------------
Constant                                     30.628***  6.360*** 28.566***
                                             (7.409)   (1.252)   (7.355)  
Percentage of population under 15            -0.471**            -0.461** 
                                             (0.147)             (0.145)  
Percentage of population over 75             -1.934              -1.691   
                                             (1.041)             (1.084)  
Real per-capita disposable income                       0.001    -0.000   
                                                       (0.001)   (0.001)  
Growth rate of real per-capita disp. income             0.529*    0.410*  
                                                       (0.210)   (0.196)  
--------------------------------------------------------------------------
sigma                                          3.931     4.189     3.803  
R-squared                                      0.262     0.162     0.338  
F                                              8.332     4.528     5.756  
p                                              0.001     0.016     0.001  
N                                             50        50        50      
==========================================================================

这是你得到的 LaTeX 代码:

texfile123 <- "mtable123.tex"
write.mtable(mtable123,forLaTeX=TRUE,file=texfile123)
file.show(texfile123)
#------------------------
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Calls:
% Model 1:  lm(formula = sr ~ pop15 + pop75, data = LifeCycleSavings) 
% Model 2:  lm(formula = sr ~ dpi + ddpi, data = LifeCycleSavings) 
% Model 3:  lm(formula = sr ~ pop15 + pop75 + dpi + ddpi, data = LifeCycleSavings) 
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{tabular}{lcD{.}{.}{7}cD{.}{.}{7}cD{.}{.}{7}}
\toprule
&&\multicolumn{1}{c}{Model 1} && \multicolumn{1}{c}{Model 2} && \multicolumn{1}{c}{Model 3}\\
\midrule
Constant                                    &  & 30.628^{***} &&  6.360^{***} && 28.566^{***}\\
                                            &  &  (7.409)     &&  (1.252)     &&  (7.355)    \\
Percentage of population under 15           &  & -0.471^{**}  &&              && -0.461^{**} \\
                                            &  &  (0.147)     &&              &&  (0.145)    \\
Percentage of population over 75            &  &  -1.934      &&              &&  -1.691     \\
                                            &  &  (1.041)     &&              &&  (1.084)    \\
Real per-capita disposable income           &  &              &&   0.001      &&  -0.000     \\
                                            &  &              &&  (0.001)     &&  (0.001)    \\
Growth rate of real per-capita disp. income &  &              &&  0.529^{*}   &&  0.410^{*}  \\
                                            &  &              &&  (0.210)     &&  (0.196)    \\
\midrule
sigma                                       &  &     3.931    &&     4.189    &&     3.803   \\
R-squared                                   &  &     0.262    &&     0.162    &&     0.338   \\
F                                           &  &     8.332    &&     4.528    &&     5.756   \\
p                                           &  &     0.001    &&     0.016    &&     0.001   \\
N                                           &  &    50        &&    50        &&    50       \\
\bottomrule
\end{tabular}
于 2012-05-15T16:49:05.320 回答
3

R wikibook 有一些关于 R 中生产质量输出的良好资源。

我认为 wikibook 中列出的 Paul Johnson 的这个功能正是您正在寻找的:

http://pj.freefaculty.org/R/WorkingExamples/outreg-worked.R

我编辑了该函数以供自己使用,以使用 booktabs 格式并允许具有额外属性的模型:

http://chandlerlutz.com/R/outregBkTabs.r

于 2012-05-15T18:02:37.623 回答
1

xtable可以做到这一点,但它有点黑客。

取两个线性模型,命名为 lm.x 和 lm.y。

如果您使用以下代码:

myregtables <- rbind(xtable(summary(lm.x)), xtable(summary(lm.y)))

xtable然后将生成一个包含两个回归模型的表。如果您\hline在 LaTeX 中添加一个(或两个),那么它应该看起来不错。这两个模型仍然只有一个标签和标题。正如我所说,它有点像一个 hacky 解决方案。

于 2012-05-15T15:09:17.860 回答