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

Question

我正在尝试从 R 中的生存包中获取 clogit 回归的稳健标准错误。在这样做时，我正在尝试使用 vce(robust) 选项复制 Stata clogit 命令报告的标准错误。

我在 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

按照here、here、here 和here 的建议，任何通过 Sand 或 plm 包提取稳健标准错误的尝试都会失败，并出现相同的错误。同样，clogit function 包含一个条件，用于在使用exact 方法（第 44 行）时停止尝试计算稳健标准误差。但是，clogit 回归对象中存在 conditional_logit$residuals 和 conditional_logit$score。

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

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

Answer 1

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

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