4

出版準備の整った回帰表を作成するためにRで利用できるツールはありますか? 私はいくつかの回帰モデルを比較する必要があるコース ペーパーに取り組んでおり、 Stata パッケージから、このような単一のテーブル内にそれらを入れ子にすることができれば非常に嬉しく思います。estout

を確認xtableしましたが、同じ結果にはなりませんでした。ヒントをいただければ幸いです。

これが私が考えていることです:

複数のネストされた回帰

4

3 に答える 3

3

R ウィキブックには、R での製品品質の出力に関する優れた情報源がいくつかあります。

ウィキブックにリストされているポール・ジョンソンのこの関数は、まさにあなたが探しているものだと思います:

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

booktabs 形式で機能するように関数を編集し、追加の属性を持つモデルを許可しました。

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

于 2012-05-15T18:02:37.623 に答える
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 に答える
1

xtableこれを行うことができますが、それはややハックです。

lm.x と lm.y という名前の 2 つの線形モデルを取り上げます。

次のコードを使用する場合:

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

xtable次に、両方の回帰モデルを含むテーブルを作成します。\hlineLaTeX に 1 つ (または 2 つ)追加すると、問題なく表示されるはずです。2 つのモデルのラベルとキャプションは 1 つだけです。私が言ったように、それはややハッキーな解決策です。

于 2012-05-15T15:09:17.860 に答える