R 中生存包的 clogit 回归的稳健标准误差
Robust standard errors for clogit regression from survival package in R
我正在尝试从 R 中的生存包中获取 clogit 回归的稳健标准错误。这样做时，我正在尝试使用 [= 复制 Stata clogit 命令报告的标准错误13=] 选项。
我在 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, and here fails with the same error. Similarly, the clogit function 通过 sandwich 或 plm 包提取稳健标准误差的尝试都包含一个条件，即在使用 exact 方法（第 44 行）时停止尝试计算稳健标准误差。但是，conditional_logit$residuals 和 conditional_logit$score 存在于 clogit 回归对象中。
如果有人能帮助回答以下问题，我将不胜感激：
- 通常不可能或 "wrong" 计算 "exact" 条件逻辑回归的稳健标准误差？如果是这样，为什么 Stata 允许这样做？
- 如果不是：我如何计算 R 中 clogit 回归的稳健标准误差？
- 如果无法根据 clogit 回归对象中的数据计算稳健的标准误差：是否有另一个 R 包生成条件逻辑回归模型，该模型等同于 survival 包的 clogit 函数生成的模型，并且其中包括我计算稳健标准误差所需的数据？
更改原始 clogit 函数中处理结的方法。
方法命令的默认值为 method="exact"。例如，如果您使用 method="efron"，则 robust=TRUE 命令有效。
我正在尝试从 R 中的生存包中获取 clogit 回归的稳健标准错误。这样做时，我正在尝试使用 [= 复制 Stata clogit 命令报告的标准错误13=] 选项。
我在 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, and here fails with the same error. Similarly, the clogit function 通过 sandwich 或 plm 包提取稳健标准误差的尝试都包含一个条件，即在使用 exact 方法（第 44 行）时停止尝试计算稳健标准误差。但是，conditional_logit$residuals 和 conditional_logit$score 存在于 clogit 回归对象中。
如果有人能帮助回答以下问题，我将不胜感激：
- 通常不可能或 "wrong" 计算 "exact" 条件逻辑回归的稳健标准误差？如果是这样，为什么 Stata 允许这样做？
- 如果不是：我如何计算 R 中 clogit 回归的稳健标准误差？
- 如果无法根据 clogit 回归对象中的数据计算稳健的标准误差：是否有另一个 R 包生成条件逻辑回归模型，该模型等同于 survival 包的 clogit 函数生成的模型，并且其中包括我计算稳健标准误差所需的数据？
更改原始 clogit 函数中处理结的方法。
方法命令的默认值为 method="exact"。例如，如果您使用 method="efron"，则 robust=TRUE 命令有效。