2015-04-20 37 views
0

在統計有很好的功能,以進行上lmglm對象殘餘&迴歸 缺失診斷,如:如何對spautolm/sarlm對象進行殘差和迴歸刪除診斷?

stats::plot.lmstats::influence.measures

實施例:

lm.object <- lm(sr ~ pop75 + dpi + ddpi, data = LifeCycleSavings) 
plot(lm.object) 
influence.measures(lm.object) 

我會喜歡在spdep::spautolm (類型「SA」)上執行這些(或類似的)診斷R「)對象,但plot.lminfluence.measures都不接受 spdep::spautolm/sarlm對象。

你能推薦一個解決方法嗎?

或者這些診斷是否不適合SAR類型迴歸?

回答

1

spautolm的對象有fittedresiduals檢索擬合值和殘差的方法。來自lm方法的一些輸出不能直接用於SAR模型,因爲這些是在不同的假設下(而不是獨立的誤差或觀察)計算的。

+0

好的,但是我怎麼做手工計算,比如槓桿,然後從'spdep :: spautolm'對象中檢索值?看到我的問題[這裏](https://stats.stackexchange.com/questions/148529/calculating-leverage-cooks-distance-for-a-weighted-spatial-simultaneous-autoreg) –

+0

取決於你是否想看看有條件或邊際殘差,見例如http://onlinelibrary.wiley.com/doi/10.1111/1467-9868.00119/abstract –

+0

恩,好問題。爲了檢查模型假設,我想我需要邊際殘差......我認爲'spautolmOBJECT $ fit $殘差是邊際殘差?!?鑑於哈斯特爾的論文,我不確定條件殘差在「spautolm」對象中會是什麼。響應 - SignalTrend?槓桿需要根據邊際殘差進行計算,對嗎? –