第11章　生命表解析

11.1 生存率の計算方法

(1) 生存率と生命表解析

医学分野では、腫瘍の手術のように、ある治療を行うことによって患者がどの程度延命するか、つまり患者の生存率をどの程度上げることができるかによって、治療の効果を判定することがしばしばあります。 そもそも治療行為の最終的な目的は、病気を治すことによって、患者が健康で長生きできるようにすることですから、生存率は治療の究極的な評価指標ということになります。

臨床試験では、治療の評価指標のことを「エンドポイント(end point)」といいます。 これは生存率が究極の評価指標であり、患者の生死を観察することによって治療の評価を行うことに由来します。

しかしたいていの疾患は死亡までの期間が長く、生存率を評価指標にして臨床試験を行うことは困難です。 そこでそのような疾患では、生存率と深い関係があることが解明されていて、治療効果を短期間で評価することのできる臨床検査値、例えば血圧や血中脂質値などを暫定的な評価指標にします。 このような暫定的な評価指標のことを、「代用エンドポイント(surrogate end point)」といいます。 それに対して生存率のような究極の評価指標のことを、「真のエンドポイント(true end point)」と呼ぶことがあります。

この究極の評価指標である生存率を解析するための統計手法のことを、「生命表解析(life table analysis)」または「生存時間解析(survival time analysis)」といいます。 生命表解析は人口統計学の分野で古くから使われてきた手法であり、現在でも生命保険会社によって盛んに利用されています。 医学分野では主として腫瘍研究の分野で利用されていて、生存時間解析という名前で呼ばれることが多いようです。

生命表解析では、手術や投薬開始という明確な出発点から個体の追跡を開始し、死亡や疾患の再発といった非可逆的なイベント(event)が発現するまで観察を続けます。 ただし、途中で個体が偶発的な出来事で研究から脱落(drop out)したり、研究期間が終了して観察打ち切り(termination of the study)になることもあります。

原則としてイベントは非可逆的な反応でなければなりませんが、再手術、治療法の切り替え、効果や副作用の発現といった、緩い意味での非可逆的な反応にすることもできます。 また最初から脱落が発生することを想定しているため、脱落例が発生する可能性が高い臨床試験、例えば治癒率が非常に高い疾患の臨床試験などにも適しています。 これらのことから、生命表解析は医学分野で広く応用することが可能な手法といえるでしょう。

またこの手法ではイベントが発現するまで観察を続け、その発現率を問題にするため、原則として前向き研究で得られたデータに適用します。 ところが実際の医学研究では、後ろ向き研究で得られたデータにこの手法を間違って適用することがしばしばあります。 これは、医学研究では後ろ向き研究が多用されることに起因しているようです。 生命表解析を適用する場合は、このあたりのことに十分注意する必要があります。

今、22例の腫瘍患者について、12例にはA法という手術を施し、10例にはB法という手術を施して、予後を観察した結果が表11.1のようになったとします。

転帰項目の「脱落」は、研究期間が終了する前に偶発的な出来事で研究から脱落した症例であり、「打ち切り」は、研究期間が終了したため生存中であるにもかかわらず観察を打ち切った症例です。 これらの症例は、生命表解析ではどちらも同じように取り扱います。

(2) 古典的な生命表

人口統計学の分野では対象となる個体数が非常に多いため、表11.1のように個々の個体を正確に観察することは困難です。 そこで観察期間をある期間で区切り、その期間内に発生した死亡数と脱落または打ち切り数を観察し、表11.2のような「生命表(life table)」にまとめます。 そして、この生命表の中の累積生存率とその標準誤差の推移を折れ線グラフで描いたものが、図11.1の累積生存率曲線です。 なお表11.2では、脱落と打ち切りを区別せずに「脱落数」として集計しています。

基本的に、ある期間の死亡率は、期間中の死亡数を期間当初の生存数つまり対象数で割った値です。 例えば0〜12ヶ月におけるA群の死亡数は2例であり、期間当初の生存数は12例ですから、死亡率は2/12≒0.167になります。

ところがこの期間中に脱落例が1例あり、その症例を対象数から除外する必要があります。 そこで平均的に期間の中央で脱落が生じたと考え、脱落例の半分の例数を対象数から除外し、それを有効生存数(有効観察例数)とします。 その有効生存数を用いて死亡率と生存率を計算すると、次のようになります。

同様に、次の期間である12〜24ヶ月におけるA群の死亡率と生存率は次のようになります。

ここで、A群における24ヶ月後の生存率の計算方法を考えてみましょう。 24ヶ月後の生存者は、0〜12ヶ月間を生存し、さらに12〜24ヶ月間も生存した例です。 したがって、ある個体が24ヵ月後に生存する確率は、0〜12ヶ月の生存確率である0.826と、12〜24ヶ月の生存確率である0.882を掛けた値になります。 この掛け合わせた生存率のことを、「累積生存率(cumulative survival rate)」または「生命表生存率(life table survival rate)」といいます。

もし脱落例が無ければ、24ヵ月後の累積生存率は、単に24ヵ月後の生存数を開始時の対象例数で割った値になります。 しかし臨床試験などでは脱落例が生じるのが普通なので、表11.1のような生命表を作成して累積生存率を計算する必要があります。 これが古典的な生存率の計算方法であり、「カトラー・エデラー(Cutler-Ederer)法」または「生命保険数理法(acturial method、狭義の生命表法)」と呼ばれています。 _(注1)

(3) 正確な生命表

医学分野では対象となる個体数が比較的少ないため、表11.1のデータをそのまま用いて生命表を作成し、正確な生存率を計算することができます。 その正確な生命表が表11.3であり、その中の累積生存率の推移を折れ線グラフで描いたものが、図11.2の累積生存率曲線です。

表11.3の生存期間項目の(転帰)には、脱落を「(+)」、打ち切りを「+」で表記し、死亡は何も表記してありません。 ただし脱落と打ち切りの生存率計算上の取り扱いは同じですから、これらは単に表記上の区別にすぎません。

正確な生命表では、死亡または脱落が発生するたびに、その瞬間の死亡率と生存率、そして累積生存率を計算します。 例えばA群では、4ヵ月後に症例番号1番の症例が脱落しています。 この症例は4ヵ月後までは生存していましたので、この時点の瞬間死亡率は0で、瞬間生存率は12/12＝1になり、累積生存率も1になります。 表11.3で、1番の症例の「生存数／観察数」が括弧付きで(12/12)と表記されているのは、脱落の場合は生存数と観測数が同じになるため、瞬間生存率が1になって累積生存率が変わらないからです。

次に5ヵ月後に症例番号2番の症例が死亡しているので、この時点の瞬間死亡率と瞬間生存率、そして累積生存率は次のようになります。

同様に、その次は8ヵ月後に症例番号3番の症例が死亡しています。

もし脱落例が無ければ、古典的な生命表と同様に、ある時点の累積生存率は、単にその時点の生存数を開始時の対象例数で割った値になります。 これが正確な生存率の計算方法であり、「カプラン・マイヤー(Kaplan-Meier)法」と呼ばれています。 _(注2)

この計算法では死亡例が発生した時点で累積生存率がカタンと落ちるため、累積生存率曲線は図11.2のように階段状の折れ線グラフになります。 そして脱落または打ち切りの有無と時期がわかるように、それらが発生した時点に短い縦棒をプロットしてあります。

(注1)　カトラー・エデラー法による生存率の計算法を模式的に表すと、次のようになります。

s₁：観測開始時点t₀直後での最初の観察対象数
s_i：時点t_i-1直後での観察対象数
d_i：t_i-1からt_iの期間での死亡数
w_i：t_i-1からt_iの期間での脱落数
s_i+1=s_i-d_i-w_i：時点t_i直後での観察対象数

観察期間は等間隔にするのが普通ですが、等間隔にしなければならないというわけではありません。 ある左開右閉区間(ti-1,ti]についての死亡率、生存率、累積生存率、そして累積生存率の標準誤差は次のようになります。

有効生存数：
死亡率：
生存率：

累積生存率：
　P₀=1≧P₁≧…≧P_i




　　※共分散C(p_k,p_k')=0のため

　　※グリーンウッド(Greenwood)の近似式

表11.2のA群について実際に計算すると、次のようになります。

・t₁：12ヵ月後



・t₂：24ヵ月後




・t₃：36ヵ月後

48ヵ月後と60ヵ月後は死亡例がないため、累積生存率とその標準誤差は変化しません。

(注2)　カプラン・マイヤー法による生存率の計算法は次のようになります。

n：観測開始時点直後での最初の観察対象数
s_i：i番目の症例が死亡または脱落する直前の観察対象数

・i番目の症例が死亡または脱落した時の瞬間生存率：
死亡例の場合：
脱落例の場合：
累積生存率：…死亡例だけ掛け合わせる
…死亡例だけ足し合わせる

表11.3のB群について実際に計算すると、次のようになります。

最後の症例は35ヵ月後に死亡例しているため、35ヵ月後の累積生存率は0になります。

この方法で累積生存率を計算した時、累積生存率が0.5になる時点つまり「50％生存時間(Median Sarvival Time、MST)」と、その95％信頼区間を求める方法が考案されています。 累積生存率がπになる時点tを挟む前後の観察時期をu、vとすると、100π％生存時間とその100(1-α)％信頼区間を次のようにして求めることができます。

100π％生存時間：
100α％信頼区間は、カプラン・マイヤー法で計算した各時点のP_iとSE(P_i)、そして正規分布の100α％点の値t(∞,α)から、各時点の100α％信頼区間=P_i±t(∞,α)・SE(P_i)を計算し、
　　下限 τ_L：時点t以前で、π≦P_i-t(∞,α)・SE(P_i)となる最初の時点(脱落例発生時点を除く)
　　上限 τ_U：時点t以後で、π≧P_i+t(∞,α)・SE(P_i)となる最初の時点

この方法は、図11.2の累積生存率曲線の上下に100α％信頼区間つまりt(∞,α)×SEを描き、その下限がπ以上になるπの直近の時点を下限にし、その上限がπ以下になる直近の時点を上限にする、という方法に相当します。

この式にπ=0.5、α=0.05を代入すれば、MSTとその95％信頼区間を求めることができます。 表11.3のA群について、50％生存時間とその95％信頼区間を計算すると次のようになります。

u=35　　v=36

時点 2= 5ヶ月：0.909-1.960×0.087≒0.744＞0.5
時点 3= 8ヶ月：0.818-1.960×0.116≒0.597＞0.5
時点 4=13ヶ月：0.727-1.960×0.134≒0.472＜0.5
時点 6=27ヶ月：0.623-1.960×0.15≒0.338＜0.5
時点 7=28ヶ月：0.519-1.960×0.157≒0.221＜0.5
MST=35.11ヶ月：0.5
時点10=36ヶ月：0.346+1.960×0.176≒0.680＞0.5

この結果から、MSTは約35.11ヶ月であり、これ以前で95％信頼区間の下限が最初に0.5以上になるのは時点3つまり8ヵ月後であることと、MST以後では最後の時点10でも95％信頼区間の上限が0.5より大きいことがわかります。 したがってMSTの95％信頼区間の下限は8ヶ月、上限は56ヶ月以上ということになります。

ちなみにB群のMSTは15ヶ月、その95％信頼区間の下限は4ヶ月、上限は25ヶ月になります。

最終更新日：2009年8月1日 　第10章第5節へ 第2節へ

　webmaster@snap-tck.com 　Copyleft (C) 2000-2012 SNAP(Sugimoto Norio Art Production)