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。