從R生存包中clogit迴歸的魯棒性標準錯誤

Question

我試圖從R中的生存包中獲得可靠的標準錯誤以進行迴歸分析。在此過程中，我嘗試複製由Stata clogit報告的標準錯誤命令與vce(robust)選項。

我中的R式是

conditional_logit <- clogit(dependent_variable ~ independent_variable + some_controls + strata(year), method= "exact", data = data_frame) 

添加robust = TRUE函數參數失敗，錯誤：

Error in residuals.coxph(fit2, type = "dfbeta", weighted = TRUE) : 
score residuals are not available for the exact method 

任何嘗試來提取通過夾心或PLM包穩健標準誤差的建議here，here,here和here失敗，並出現相同的錯誤。同樣，使用exact方法（第44行）時，clogit function包含停止嘗試計算強健標準錯誤的條件。但是，clogit迴歸對象中存在conditional_logit $ residuals和conditional_logit $分數。

我會感激如果有人可以幫助回答下列問題：

• 它是一般不可能或「錯」來計算「精確」條件Logistic迴歸穩健的標準差？如果是這樣，爲什麼Stata允許這樣做？
• 如果不是的話：我怎樣才能計算R中阻塞回歸的穩健標準誤？
• 如果不可能基於clogit迴歸對象中的數據計算可靠的標準錯誤：是否有另一個R包產生的條件邏輯迴歸模型等同於生存包的clogit函數產生的條件邏輯迴歸模型，並且包括我需要計算強大的標準錯誤的數據？

Answer 1

更改處理原始clurit函數中關係的方法。

method命令的默認值是method="exact"。例如，如果使用method="efron"，則robust=TRUE命令有效。