1

よくわかりません... これはそれほど難しいとは思いませんが、解決できません。実行する場合:

library(survival)    
leukemia.surv <- survfit(Surv(time, status) ~ 1, data = aml) 
plot(leukemia.surv, lty = 2:3)

生存曲線とその 95% 信頼区間が表示されます。95% CI の上限と下限を示す 2 本の線を表示する代わりに、95% の上限と下限の間の領域に陰影を付けます。

これは、polygon() などで行う必要がありますか? すべての座標は概要で見つけることができます...

> summary(leukemia.surv)
Call: survfit(formula = Surv(time, status) ~ 1, data = aml)

 time n.risk n.event survival std.err lower 95% CI upper 95% CI
5     23       2   0.9130  0.0588       0.8049        1.000
8     21       2   0.8261  0.0790       0.6848        0.996
9     19       1   0.7826  0.0860       0.6310        0.971
12     18       1   0.7391  0.0916       0.5798        0.942
13     17       1   0.6957  0.0959       0.5309        0.912
18     14       1   0.6460  0.1011       0.4753        0.878
23     13       2   0.5466  0.1073       0.3721        0.803
27     11       1   0.4969  0.1084       0.3240        0.762
30      9       1   0.4417  0.1095       0.2717        0.718
31      8       1   0.3865  0.1089       0.2225        0.671
33      7       1   0.3313  0.1064       0.1765        0.622
34      6       1   0.2761  0.1020       0.1338        0.569
43      5       1   0.2208  0.0954       0.0947        0.515
45      4       1   0.1656  0.0860       0.0598        0.458
48      2       1   0.0828  0.0727       0.0148        0.462

95% CI 領域をシェーディングする既存の関数はありますか?

4

2 に答える 2

4

からのデータを使用summary()して、信頼区間を多角形として独自のプロットを作成できます。

まず、 をsummary()オブジェクトとして保存します。プロット用のデータは、変数time、surv、upperおよびにありlowerます。

mod<-summary(leukemia.surv)

plot()関数を使用してプロット領域を定義できるようになりました。次に、polygon()プロットの信頼区間を使用します。ここでは、x 値と x 値を逆の順序で指定する必要があります。y 値にはlower値を使用し、値を尊重しupperます。lines()サバイバルライン追加機能付き。type="s"に引数を追加lines()すると、行がステップとして取得されます。

with(mod,plot(time,surv,type="n",xlim=c(5,50),ylim=c(0,1)))
with(mod,polygon(c(time,rev(time)),c(lower,rev(upper)),
                 col = "grey75", border = FALSE))
with(mod,lines(time,surv,type="s"))
于 2013-09-03T05:52:40.693 に答える
0

生存曲線に影付きの信頼区間をプロットする関数を開発しました。ここで見つけることができます：ggplot2を使用してRで生存曲線をプロットする

多分あなたはそれが便利だと思うかもしれません。

于 2014-04-13T12:42:01.163 に答える