與R中

一個新手的問題，可靠的集羣標準誤差Logistic迴歸：沒有人知道如何在執行中的R集羣標準誤差的迴歸？在Stata它只是logit Y X1 X2 X3, vce(cluster Z)，但不幸的是我還沒有想出如何在R做同樣的分析。與R中

來源

2013-05-11 danilofreire

'sandwich'包中的'vcovHC（）'函數也可能是有用的（不確定它是否適用於邏輯迴歸估計） – 2013-05-11 21:34:57

如果您從Stata遷移，您可能會發現名爲「plm」的包有用。此外，有一個名爲包''pcse''通過操縱估計 – 2013-05-12 06:36:38

後方差協方差矩陣實現面板修正的標準差非常感謝您的答覆，本和休伯特。我也會測試你所建議的軟件包，看看它們是否適用於物流估算。再次感謝！ – danilofreire 2013-05-13 22:25:46

你可能想看看rms（迴歸建模策略）封裝。所以，lrm是邏輯迴歸模型，如果fit是你輸出的名字，你有這樣的事情：

fit=lrm(disease ~ age + study + rcs(bmi,3), x=T, y=T, data=dataf) 

fit 

robcov(fit, cluster=dataf$id) 

bootcov(fit,cluster=dataf$id)

您在模型語句指定x=T，y=T。 rcs表示三節有限的三次樣條曲線。

來源

2013-05-11 18:00:52

非常感謝！它創造了奇蹟！我會更仔細地閱讀rms的手冊，看看是否有按國家和年份聚集係數的方法。再一次，謝謝你！ – danilofreire 2013-05-13 22:27:12

這個答案已經非常好，但如果它完全可複製的話，它可以得到改進。我不知道變量來自哪裏，輸出是什麼，爲什麼'rcs（bmi，3）'是必須的。 – MERose 2018-01-15 14:46:26

我一直在敲我的頭靠在這個問題在過去兩天;我神奇地發現看起來似乎是一個很好的東西的新軟件包 - 例如，我也在分析一些集羣健壯的Tobit模型，並且該軟件包內置了該功能。且不說語法比我見過的所有其他解決方案更清潔多（我們談論的清潔近Stata的水平）。

因此，對於你的玩具比如，我會運行：

library(Zelig) 
logit<-zelig(Y~X1+X2+X3,data=data,model="logit",robust=T,cluster="Z")

的Et瞧！

來源

2015-04-01 23:50:43 MichaelChirico

哇，這似乎是「只是工作」的方式，我的R代碼似乎從來沒有。這是新功能嗎？如果沒有，爲什麼澤利格不是在R中解決這個問題的標準方法？ – Philip 2015-05-05 03:35:31

不知道，但我希望它變得如此。 [The project]（http://zeligproject.org）的確顯得雄心勃勃！ [Google小組]（https://groups.google.com/forum/m/#!forum/zelig-statistical-software）雖然看起來並不那麼活躍，但不知道進展有多快。 – MichaelChirico 2015-05-05 04:11:26

不幸的是，我認爲這個命令在最新版本的「Zelig」（在CRAN上）不起作用。我剛剛運行了幾個模型，並且沒有'cluster'參數，標準錯誤也完全一樣。我相信從4.0版本開始就是這樣，我最後一次使用這個軟件包。 – danilofreire 2015-07-01 05:07:39

回答

相關問題