1

我有这个生存数据，描述了十年时间研究中三种服务（沙龙、餐厅和快递）的死亡率。

数据包含三个变量：服务类型（1=saloon、2=restaurant 和 3=express）、年数（1 到 11 的整数，其中 11 表示大于 10 年）和审查员。

我有两个问题：

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

2）如何在 R 中拟合时间相关的 Cox 模型？

这是我的代码：

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

1 回答 1

0

pkg 生存中的 cox.zph 函数检查相称性。请注意，使用“服务类型”作为列名提供了编程麻烦，可以通过允许用句点替换空格来轻松避免这些麻烦，就像 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
于 2016-03-08T00:31:21.763 回答