R中是否有可用的工具来生成可发布的回归表?我正在写一篇课程论文,我需要在其中比较几个回归模型,如果我可以让它们嵌套在像这样的一个表中,我会很高兴,来自estout
Stata 包。
我已经检查过xtable
,但无法达到相同的结果。任何提示将不胜感激。
这是我的想法:
您可能想要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}
R wikibook 有一些关于 R 中生产质量输出的良好资源。
我认为 wikibook 中列出的 Paul Johnson 的这个功能正是您正在寻找的:
http://pj.freefaculty.org/R/WorkingExamples/outreg-worked.R
我编辑了该函数以供自己使用,以使用 booktabs 格式并允许具有额外属性的模型:
xtable
可以做到这一点,但它有点黑客。
取两个线性模型,命名为 lm.x 和 lm.y。
如果您使用以下代码:
myregtables <- rbind(xtable(summary(lm.x)), xtable(summary(lm.y)))
xtable
然后将生成一个包含两个回归模型的表。如果您\hline
在 LaTeX 中添加一个(或两个),那么它应该看起来不错。这两个模型仍然只有一个标签和标题。正如我所说,它有点像一个 hacky 解决方案。