3

我需要在 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了没有结果的包。

有人可以帮助我吗?谢谢。

4

1 回答 1

3

这里是:

zph$table[ , "p" ]

您可以使用 访问 zph 对象的成员$,因为它是一个列表:

names( zph )
# returns:
# [1] "table"     "x"         "y"         "var"       "call"      "transform"

然后抬头一看zph$table

在 R 中,通常可以通过$(S3 OO 框架)或@ (S4 OO 框架)访问对象的不同“槽”。

于 2013-07-17T09:00:16.560 回答