2011-03-22 424 views
2

我試圖自動執行這個過程,通過這個過程,我可以將來自許多回歸模型的係數和標準錯誤輸出到單個表格中,以便通過xtable作爲html輸出。將回歸模型估計值合併到一張表中

一些類似的問題(like this one)已被定向到Paul Johnson的一個稱爲「outreg」的函數,但該頁面不再存在,我找不到代碼。

其他(like this one)使用似乎給我錯誤的解決方案,因爲我的模型並不都具有相同數量的變量。

爲了澄清我的任務......

  • 我有n個POLR(有序Logistic)從中我要輸出的係數和標準誤差模型。
  • 每個模型包含不同數量的預測變量。
  • 我需要一個data.frame(?)爲每個模型列或兩行的每個預測
  • 這不是關鍵的標準誤差是如何在相對於輸出到係數

模型有這樣的輸出具有連續多個預測:

>summary(model1)["coefficients"] 
$coefficients 

       Value Std. Error t value 
relGPA   0.8683499 0.04185389 20.74717 
mcAvgGPA  1.3885515 0.09688030 14.33265 
Deny|Waitlist 0.5707912 0.01553476 36.74283 
Waitlist|Accept 0.8398921 0.01618358 51.89779 
+1

我想這將適用於這裏:http://stats.stackexchange.com/questions/6856/aggregating-results-from-linear-model-runs-r – Chase 2011-03-22 23:25:00

+0

非常好,謝謝蔡斯!關於這個問題的圖表建議非常棒,但是mtable {memisc}解決方案正是我希望找到的解決方案。 – 2011-03-23 16:49:13

+0

PJ的網站上仍然提供outreg功能:http://pj.freefaculty.org/stat/ps706/ – ATMathew 2011-03-23 17:59:28

回答

3

您可以剪切和粘貼here。我已經用outreg.R修改了一些功能(t-統計,多頁面等),但它只在我的工作計算機上,所以我明天就會發布。


更新

Here's我扭捏版本,但想當初,它仍需要lmglm對象的列表。剪切&粘貼似乎太長了,所以這是一個鏈接到我的dropbox.com公用文件夾。

+0

感謝理查德,我很想看看你的修飾版本,因爲我沒有在PDF中的代碼運氣。我在函數的某處出現了有關使用「$」的錯誤。我認爲這可能與我的模型是polr()有關。或者我可能已經引入了從pdf轉移的錯誤。 – 2011-03-23 16:15:01

+0

我不熟悉'polr()',但'outreg.R'需要'lm'或'glm'對象。 – 2011-03-23 18:42:34

+0

感謝您貢獻您的outreg版本。我會將其標記爲已回答,因爲對於不執行序號邏輯模型的人來說,這是一個很好的解決方案。對於其他人來說,memsic中的mtable功能運行良好。 – 2011-03-23 20:32:57