R 中生存包的 clogit 回归的稳健标准误差

Question

我正在尝试从 R 中的生存包中获取 clogit 回归的可靠标准错误。在此过程中，我尝试使用该选项复制 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 包提取稳健标准错误（如此处、此处、此处和此处建议）的尝试都会失败并出现相同的错误。类似地，clogit 函数包含一个在使用该方法时停止尝试计算稳健标准误差的条件exact（第 44 行）。但是，conditional_logit$residuals 和conditional_logit$score 存在于clogit 回归对象中。

如果有人能帮助回答以下问题，我将不胜感激：

• 计算“精确”条件逻辑回归的稳健标准误差通常是不可能或“错误”的吗？如果是这样，为什么 Stata 允许这样做？
• 如果不是：我如何计算 R 中 clogit 回归的稳健标准误差？
• 如果无法根据 clogit 回归对象中的数据计算稳健的标准误差：是否有另一个 R 包可以生成条件逻辑回归模型，该模型与 Survival 包的 clogit 函数生成的模型等效，并且包含以下数据：我需要计算稳健的标准误差？

Answer 1

更改原始 clogit 函数中处理关系的方法。

方法命令的默认值为method="exact"。method="efron"例如，如果您使用，则该robust=TRUE命令有效。