我需要在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
获取attributes(zph)
属性值(在本例中为0.57),但p
并未出现以执行zph$p
之类的操作。我还使用了plyr
包而没有结果。
有人可以帮助我吗?感谢。
答案 0 :(得分:1)
这是:
zph$table[ , "p" ]
您可以使用$
访问zph对象的成员,因为它是一个列表:
names( zph )
# returns:
# [1] "table" "x" "y" "var" "call" "transform"
然后查找zph$table
和presto。
在R中,通常可以通过$
(S3 OO框架)或@
(S4 OO框架)访问对象的不同“槽”。