我想将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）中的错误： 不能创建方差很大的方差分析表
是否有解决此问题的方法？我对R中的生存分析非常陌生，所以也许我不知道有提供此功能的特定软件包？
非常感谢！