始める無料で始める

Kaplan-Meier 解析

この演習では、共変量なし/ありの両方で Kaplan-Meier 解析を練習します。

survival パッケージはすでに読み込まれています。サバイバルオブジェクト survObj とデータ dataNextOrder も環境にあります。今回は、追加の共変量 voucher がデータに含まれています。これはカテゴリ変数で、顧客が最初の注文でバウチャーを使用したかどうかを表し、値は 0 または 1 です。

この演習はコースの一部です

Rで学ぶマーケティングアナリティクスのための機械学習

コースを見る

演習の手順

  • survfit() を使って(共変量なしの)Kaplan-Meier 解析を計算し、結果を fitKMSimple というオブジェクトに保存してください。従属変数(チルダ ~ の「左」側)は、今回もサバイバルオブジェクト survObj です。続いて、fitKMSimple を表示します。
  • 結果オブジェクト fitKMSimple をプロットし、軸ラベル(xlabylab 引数)を追加してください。
  • 次のステップとして、従属変数に survObj、共変量に voucher を用いた Kaplan-Meier 解析を計算してください。data 引数の指定を忘れないでください。
  • 新しいモデルの結果も同様にプロットし、軸ラベルを追加してください。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Compute and print fit
fitKMSimple <- ___(___ ~ ___)
___(fitKMSimple)

# Plot fit
plot(___,
     conf.int = FALSE, ___ = "Time since first purchase", ___ = "Survival function", main = "Survival function")

# Compute fit with categorical covariate
fitKMCov <- survfit(___ ~ ___, data = ___)

# Plot fit with covariate and add labels
plot(___, lty = 2:3,
     ___ = "Time since first purchase", ___ = "Survival function", main = "Survival function")
legend(90, .9, c("No", "Yes"), lty = 2:3)
コードを編集して実行