我需要在 R 函数p
的输出中提取特定行的属性。cox.zph
为了让您有疑问,我将我的问题描述如下,并通过一个示例逐步说明:
require('survival')
# I create the simplest test data set
test1 <- list(time=c(4,3,1,1,2,2,3),
status=c(1,1,1,0,1,1,0),
x=c(0,2,1,1,1,0,0),
sex=c(0,0,0,0,1,1,1))
# Fit a stratified model
coxmodel <- coxph(Surv(time, status) ~ x + strata(sex), test1)
然后,我使用cox.zph
功能:
zph <- cox.zph(coxmodel)
使用以下输出:
rho chisq p
x 0.354 0.322 0.57
我试图p
通过使用来获取属性值(在本例中为 0.57),attributes(zph)
但p
为了执行类似zph$p
. 我也使用plyr
了没有结果的包。
有人可以帮助我吗?谢谢。