総仕上げ:Cox モデル
コースの締めくくりとして、直近の2章で少しだけ使った肺がんデータセットを振り返りましょう。おさらいすると、このデータセットには North Central Cancer Treatment Group による進行肺がん患者の生存に関する情報が含まれています。イベントは status 変数に格納されており、個人が生存しなかった場合は値が 2 です。パフォーマンススコア(変数 performance)は、医師が評価した、患者が通常の生活活動をどの程度行えるか(悪い=0、良い=100)を表します。特定のパフォーマンススコアと生存時間の関連を知りたいと考えています。
この演習はコースの一部です
Rで学ぶサバイバル分析
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Compute Cox model and survival curves
cxmod <- ___(___(___, ___) ~ ___, data = lung)
new_lung <- data.frame(performance = ___)
cxsf <- ___(cxmod, data = ___, newdata = ___, conf.type = "none")