Rp
の関数の出力で特定の行の属性を抽出する必要があります.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
パッケージも使用しましたが、結果はありませんでした。
誰でも私を助けることができますか?ありがとう。