R anova.coxph Cox 模型的偏差分析。

说明

基于对数偏似然，计算一个或多个 Cox 模型拟合的偏差表分析。

用法

## S3 method for class 'coxph'
anova(object, ...,  test = 'Chisq')

参数

object	coxph 类的对象
...	更多coxph 对象
test	一个字符串。适当的测试是卡方，所有其他选择都会导致不进行测试。

细节

指定单个对象可以对该拟合的偏差表进行顺序分析。也就是说，随着公式的每一项依次相加，模型 Cox log-partial-likelihood 中的减少量以表的行加上对数似然本身的形式给出。稳健的方差估计通常用于模型可能是mis-specified的情况，例如每个受试者有多个事件。在这种情况下，似然值的比较没有意义(差异不再具有卡方分布)，并且 anova 将拒绝打印结果。

如果指定了多个对象，则表中的一行表示每个模型的自由度和对数似然。对于除第一个模型之外的所有模型，还给出了自由度和 loglik 的变化。 (这仅在模型嵌套时才具有统计意义。)按照惯例，从最小到最大列出模型，但这取决于用户。

该表可以选择包含比较每行 loglik 减少情况的测试统计数据(和 P 值)。

值

类 "anova" 的对象继承自类 "data.frame" 。

警告

仅当两个或多个模型拟合到同一数据集时，anova 之间的比较才有效。如果存在缺失值，这可能会出现问题。

例子

fit <- coxph(Surv(futime, fustat) ~ resid.ds *rx + ecog.ps, data = ovarian) 
anova(fit)
fit2 <- coxph(Surv(futime, fustat) ~ resid.ds +rx + ecog.ps, data=ovarian)
anova(fit2,fit)
 

也可以看看

coxph、anova。