2013-05-29 57 views
0

我有公司財務信息組織的各公司的財務信息數據。我希望在保持公司不變的情況下,將其中一列的價值觀與其他人的價值觀相比較。有沒有一種簡單的方法可以用lm()表示法寫出來?按數據組分組的迴歸係數R

我已經嘗試使用:

reg <- lmList(lead2.dDA ~ paudit1 + abs.d.GINDEX + logcapx + logmkvalt + 
       logmkvalt2|pp, data=reg.df) 

其中pp是公司名稱的載體,但好像我回歸的所有數據一次(並沒有按公司名稱分開)這將返回係數。

+2

'lmList'調用看起來應該已經工作。請提供一個可重複的示例,說明它無法正常工作。 –

回答

2

lm()中用於估計單獨的迴歸係數組的方便且明顯鮮爲人知的語法包括使用嵌套運算符/。在這種情況下,它看起來像:

reg <- lm(lead2.dDA ~ 0 + pp/(paudit1 + abs.d.GINDEX + logcapx + 
      logmkvalt + logmkvalt2), data=reg.df) 

確保pp因素而不是數字。還要注意,爲了使其工作,必須抑制整體攔截;在新的公式中,我們對每個組都有不同的「截距」。

一對夫婦的評論:

  • 儘管得到這樣將匹配那些lmList()給出的迴歸係數,應該指出的是,與lm()我們估計只有在所有組的單個剩餘方差,而lmList()會估計每個組的單獨殘差。
  • 就像我在我之前的評論中提到的那樣,你給出的lmList()語法看起來應該是可行的。既然你說沒有,那麼這讓我期望這個問題真的是別的東西(儘管很難說沒有一個可重複的例子),所以我發佈的解決方案似乎也可能會失敗,原因不明。如果您需要更詳細的指導,請提供更多信息;幫助我們幫助你。
+0

謝謝傑克!我不關心殘差,所以lm()方法很好。 – user2303635