0

死亡率を予測する R の次のコードを使用して設定された Cox 比例ハザード モデルがあります。共変量 A、B、および C は単に交絡 (年齢、性別、人種) を避けるために追加されていますが、私たちは予測子 X に本当に関心があります。ここで、X は死亡リスクを下げる薬です (私たちの仮説)。

そこで、他の共変量 (A、B、C) を調整して 2 つのグループ (X=0 と X=1) を比較するコックス生存プロットを作成したいと思います。

私は最初のステップから始めました

cox.model <- coxph(Surv(time, dead) ~ A + B + C + X, data = df)

次に、このコードを使用してコックスプロットを作成しました

ggsurvplot(survfit(cox.model), data = df, palette = "#2E9FDF",
       ggtheme = theme_minimal(),
       risk.table = TRUE)

(X = 1とX = 0)を比較するのではなく、1つの(累積)行だけでプロットを作成しました。リスクテーブルは1つのグループにすぎませんでした。

私の質問は、そのプロットで比較をどのように表示できるかです。また、リスク テーブルに 2 つのグループを表示する方法は?

注: KM を実行したとき、共変量 A+ B+ C を調整せずに比較するのは簡単でした。

4

0 に答える 0