死亡率を予測する 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 を調整せずに比較するのは簡単でした。