出版準備の整った回帰表を作成するためにRで利用できるツールはありますか? 私はいくつかの回帰モデルを比較する必要があるコース ペーパーに取り組んでおり、 Stata パッケージから、このような単一のテーブル内にそれらを入れ子にすることができれば非常に嬉しく思います。estout
を確認xtable
しましたが、同じ結果にはなりませんでした。ヒントをいただければ幸いです。
これが私が考えていることです:
R ウィキブックには、R での製品品質の出力に関する優れた情報源がいくつかあります。
ウィキブックにリストされているポール・ジョンソンのこの関数は、まさにあなたが探しているものだと思います:
http://pj.freefaculty.org/R/WorkingExamples/outreg-worked.R
booktabs 形式で機能するように関数を編集し、追加の属性を持つモデルを許可しました。
おそらく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}
xtable
これを行うことができますが、それはややハックです。
lm.x と lm.y という名前の 2 つの線形モデルを取り上げます。
次のコードを使用する場合:
myregtables <- rbind(xtable(summary(lm.x)), xtable(summary(lm.y)))
xtable
次に、両方の回帰モデルを含むテーブルを作成します。\hline
LaTeX に 1 つ (または 2 つ)追加すると、問題なく表示されるはずです。2 つのモデルのラベルとキャプションは 1 つだけです。私が言ったように、それはややハッキーな解決策です。