患者の痛み治療を比較する

新しい痛み治療の実験があり、新治療を受けた患者と旧治療を受けた患者を比較したいとします。各グループに対して Kaplan–Meier 推定量を当てはめ、サバイバル曲線を並べて可視化します。

DataFrame は recur で、以下の列が含まれます。

• age: 患者の年齢
• treat: 受けた治療（新治療なら 0、旧治療なら 1）
• duration: 治療後の経過時間（時間）
• censor: 事象が打ち切りかどうか

pandas は pd として読み込み済み、lifelines から KaplanMeierFitter クラスがインポートされ、matplotlib から pyplot モジュールが plt としてインポートされています。

さらに、KaplanMeierFitter のインスタンス kmf と、サブプロット図 ax が用意されています。