医療統計と疫学

Kaplan-Meierの生存曲線をRで理解する

生存解析をする上で欠かせない”Kaplan-Meier(カプランマイヤー)”。もはやその方法(Kaplan-Meier estimator)を用いるのに、元文献を提示しなくても良いほどメジャーなものになりました。今回は、そのKaplan-Meierの理屈をRとともに解説したいと思います。

Censoringのない場合

観察期間中のdrop-out(”loss to follow-up”)や観察期間中のeventなし(”administrative censoring”)といったcensoringがない場合は、すべての患者において観察期間中にeventが発生したことになります。

求めたいゴールであるS(t)は、eventまでの期間がtより長くなる確率(ex. tより長く生存する確率)です。S(8)であれば、8 weeksより後にeventが起こる確率であす。今回の例では計21人中、青字で示した8人となる確率、すなわち、P(T>8) = 8/21 = 0.381となります。Censoringがなければ、どのタイミングであっても分母はtotal数の21であり、すべてのS(t)を簡単に求められます。

(すなわち、Kaplan-Meierを用いるまでもありません)

手計算で簡単に求められますが、折角なのでRを使ってみましょう。survivalというpackageを用います。

とすると、アウトプットは以下のようになります。

t=8の時のsurvival probabilityは0.381となっていますね。

Censoringのある場合

まずは用語の説明をします。生存解析においては、censoringとなった人には、その時点の時間の横に”+”を付けて表します。10週でloss to follow-upとなれば、”10+”と記します。

そして、tjを時点j, djを時点jにおけるevent発生数, cjを時点jにおけるcensorとなった人数として、それぞれ表します。

censoringがある場合、S(t)の計算は上記のように単純にはいきません。なぜなら、censoringとなった人が、どのタイミングでeventが起こったかわからないから、すなわち、分母が必ずしもトータル人数になりません。

そこで、Kaplan-Meierという方法を用いてS(t)を推測することになります。

Kaplan-Meier Estimator

tが連続ではなく時点・時点で区切られている(Day1, Day2…など)discrete caseを考えてみましょう。

詳細は上記のような数式から導かれますが、S(t) = P(T > t)は結局は、

それぞれの時点での[1 – {eventが起こった人数(dj)/リスクのある人数(rj)}]を、時点jまで掛け合わせたもの

になります。といってもわかりにくいですよね。Kaplan-Meierは、一度自分で表を作ることをお勧めします。それぞれの時点tjにおける、dj, cj, rjを書き、そこから1-dj/rj、そしてS(tj)を計算します。

上記の例を表にすると、以下のようになります。

この一番右列にあるS(t)が、Kaplan-Meier estimatorを用いて推測したS(t)です。Kaplan-Meierも、ちゃんと手計算できますね。

こちらも、Rを用いて算出してみましょう。

アウトプットは、以下のようになります。

いかがでしょうか。手計算と同じS(t)が得られたでしょうか。

計算してみればわかりますが、censoringがない場合もKaplan-Meierを使うことができ、上記と同じ結果になります。

Kaplan-Meierの生存曲線

最後に、Kaplan-Meier Estimatorを用いた生存曲線を描いてみましょう。縦軸にS(t)、横軸に時間をとります。上の手計算で作成した表を元に描くことができますが、ここではRを使って描いてみましょう。先ほど定義したkmを用います。

ちゃんとKaplan-Meierの生存曲線が描けました。ちなみに、生存曲線と直交している短い線は、censoringの患者を表しています。

まとめ

Kaplan-Meierを理解するためには、上記のような表を一度でいいから自分で作ることがとても重要です。Event発生だけでなく、それが起こるまでの時間にも興味がある、途中でfollow-upできなくなりoutcomeがわからない人がいる、といった場合には、視覚的にも有効な方法ですので、是非理解しておいてください。

 

医療従事者に必要な統計学と疫学(目次)へ

Reference

Andrea Bellavia. BST 223: Applied Survival Analysis. Harvard T.H. Chan School of Public Health

ABOUT ME
木村聡
福岡県の研修病院で初期研修修了後、大学に入局。米国オハイオ州に臨床研究で留学するも、知識の欠如を痛感。ハーバード公衆衛生大学院に進学し、MPHを取得。マサチューセッツ工科大学メンバーとの共同研究などに関与。 日本で麻酔・集中治療医として働いた後、オーストラリアで臨床留学も経験。 書籍『絶対あきらめない医学留学』著者。 乗り越える壁を見つけ続けることは、なかなか簡単ではありませんよね。アラフォー目前、様々な壁にぶち当たり、それなりに多くの経験をしてきました。私の挑戦や経験・知識、失敗談などが、他の誰かの刺激になり、役に立つことを切に願っています。 プロフィールをもっと詳しく見る

POSTED COMMENT

  1. 歩海 より:

    私はブラウンでMPH、ワシントン大学の疫学博士課程にいます。(現在アラフォー独身女性w、私の専門分野であるインドの伝統医学の研究のために進学しました。KーMとてもわかりやすかった上にRのコードまでついていてとてもありがたかったです!これからもこのような情報とてもお待ちしています!!

COMMENT

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

error: