我想将 Cox 比例风险模型的模型拟合与似然比检验（或等价物）进行比较。这通常通过将拟合与anova()函数进行比较来实现。但是，当我有以集群形式组织的数据并因此必须以这种方式定义时，这不起作用。见下面的例子（肺数据集包含在生存包中）：
library(survival)
data(lung)
lung_cox <- coxph(Surv(time, event = status) ~ age + cluster(sex), data = lung)
lung_cox1 <- coxph(Surv(time, event = status) ~ 1 + cluster(sex), data = lung)
anova(lung_cox, lung_cox1)
我在函数后得到的错误anova()是：
anova.coxphlist(c(list(object), dotargs), test = test) 中的错误：不能做具有稳健方差的 anova 表
有没有办法解决这个问题？我对R中的生存分析很陌生，所以也许我不知道提供这个的某个包？
非常感谢！