Question

我有survival data描述了十年时间研究中三种服务（沙龙，餐馆和快递）的死亡率。

数据包含三个变量：服务类型（1 =轿车，2 =餐厅，3 =快递），年份（1到11之间的整数，其中11表示大于10年）和审查。

我有两个问题：

1）我已经安装了Cox比例风险模型，但检查比例风险假设的方法是什么。也就是说，我们假设每个人的危险比和基线危险与时间无关。

2）如何在R？

中拟合时间相关的Cox模型

以下是我的代码：

   #Cox Proportional Hazards
   cox <- coxph(Surv(Years, Censor) ~ data$`Service Type` )
   summary(cox)

Answer 1

pkg生存中的cox.zph函数检查比例性。请注意使用＆＃39;服务类型＆＃34;作为列名提供编程麻烦，通过允许句点替换空格可以轻松避免，就像read.table的默认操作一样：

data <- read.table(url("http://www.stat.ufl.edu/~winner/data/bizmort.dat"), col.names=c("Service Type","Years",  "Censor")
# Also note that the censoring indicatior is reversed so will use 1-Censor
require(survival)

cox <- coxph(Surv(Years, 1-Censor) ~ factor(Service.Type), data=data )
# The test fro proportionality:

> cox.zph(cox)
                         rho chisq     p
factor(Service.Type)2 0.0306  0.98 0.322
factor(Service.Type)3 0.0429  1.91 0.167
GLOBAL                    NA  2.33 0.312

R中的Cox比例风险模型和时间依赖的Cox模型

1 个答案: