死亡率・生存分析 計算ツール｜患者数・追跡期間・死亡数から粗死亡率と5年生存率を即算出、KM法・打切り・年齢調整まで
臨床研究・疫学調査の要である粗死亡率(%/年)と5年生存率(%)を、患者数・追跡年数・死亡数から瞬時に算出。Kaplan-Meier法、Cox比例ハザードモデル、打切り(censoring)、年齢調整死亡率、SMR(標準化死亡比)、コホート研究と症例対照研究の設計選択、生命表法との違いまで、厚生労働省人口動態統計・WHO国際疾病分類ICDに基づき解説します。臨床試験、大学病院・地域基幹病院の医療の質評価、公衆衛生・疾病対策の実務にご活用ください。
死亡率・生存分析ツール
計算式と単位系
基本式
粗死亡率(%/年) = 死亡数 ÷ 患者数 ÷ 追跡年数 × 100
5年生存率(%) = (1 − 年死亡率/100)⁵ × 100 (指数モデル近似)
本計算は、追跡期間中の死亡確率が一定と仮定する指数モデル(constant hazard)による近似です。1,000人を5年追跡し50人が死亡した場合、年死亡率は50÷1000÷5×100=1.0%、5年生存率は(1-0.01)⁵×100=約95.1%となります。より正確な生存曲線はKaplan-Meier法で描きます。
密度(rate)と確率(probability)
「死亡率」には2つの意味があります。(1)率(rate)は単位時間当たりの発生密度で、通常「人年(person-year)」を分母に取り、%/年で表現します。(2)確率(risk)は一定期間内の発生確率で、%(単位なし)で表現します。1,000人年で10死亡なら死亡率=10/1000=1.0/100人年、1,000人5年追跡で50死亡なら死亡確率=50/1000=5.0%(累積)です。
人年(person-year)ベース
粗死亡率(/100人年) = 死亡数 ÷ 総観察人年 × 100
各患者の観察期間が異なる場合(途中脱落・研究終了)は、単純に「患者数」ではなく「総観察人年」を分母にします。1,000人を平均4年追跡(=4,000人年)で50死亡なら、死亡率は50/4000×100=1.25/100人年。この方法はKM法と同じ考え方に基づく標準的な疫学指標です。
打切り(censoring)の扱い
臨床研究では全患者を最終まで観察できるとは限りません。研究終了時点・追跡不能・別死因での死亡等で観察が終了する状態を「打切り(censoring)」と呼びます。
| 打切りの種類 | 発生原因 | 統計処理 |
|---|---|---|
| 右打切り(最も一般的) | 研究終了時に生存、追跡脱落 | KM法の標準対象 |
| 左打切り | 登録前に既にイベント発生 | 登録日を起点に再設定 |
| 区間打切り | 2つの受診日の間で発生 | 区間中点を仮定または区間モデル |
| 競合リスク | 目的以外の死因(例:がん研究中の心疾患死) | Fine-Grey/累積発生関数で処理 |
打切りを無視して単純に「死亡数÷患者数」で計算すると、追跡不十分な例が「生存扱い」になり死亡率を過小評価します。逆に打切り例を除外(complete case)すると母集団を歪めます。KM法・生命表法はこれを適切に処理します。
Kaplan-Meier法の基礎
Kaplan-Meier(KM)法は、打切りを適切に扱いながら生存曲線を推定するノンパラメトリック手法で、臨床研究の標準です。各死亡時点で以下を計算します。
S(tₖ) = S(tₖ₋₁) × (1 − dₖ/nₖ) │ dₖ=時点tₖの死亡、nₖ=直前の生存者数
時点tₖ直前の生存者数nₖ人のうちdₖ人が死亡した場合、その時点の条件付き生存率は(1-dₖ/nₖ)、累積生存率は前時点までの生存率と掛け合わせます。打切り例はその時点でnₖから除外されますが「死亡」ではないため、S(t)は下がりません。生存曲線は階段状に降下し、95%信頼区間はGreenwoodの公式で算出します。
実務ではSPSS・SAS・R(survivalパッケージ)・EZR(自治医大)等で描画・検定します。2群比較はログランク検定(log-rank test)、リスク比較はハザード比(HR)を用います。
Cox比例ハザードモデル
Cox比例ハザードモデルは、複数の説明変数(年齢・性別・治療群・病期等)の影響を同時に評価する半パラメトリック回帰モデルです。KM法が「群の生存曲線を比較」するのに対し、Cox回帰は「共変量を調整したハザード比」を算出できます。
h(t|X) = h₀(t) × exp(β₁X₁ + β₂X₂ + … + βₚXₚ)
ハザード比(HR)=exp(β)は、他条件を固定した時のリスク倍率を示します。HR=2.0なら死亡リスクが2倍、HR=0.5なら半減。「比例ハザード仮定(HRが時間に依存しない)」が前提のため、Schoenfeld残差プロットで仮定確認が必要です。
| HR | 解釈 | 臨床的意義 |
|---|---|---|
| 1.0 | 差なし | 基準群と同等 |
| 0.5-0.8 | リスク低下 | 治療効果あり |
| 1.2-1.5 | リスク軽度上昇 | 要注意因子 |
| 2.0以上 | リスク大幅上昇 | 強い予後不良因子 |
年齢調整・SMR・PMR
年齢調整死亡率
集団間の粗死亡率比較は、年齢構成の違いで歪みが生じます。直接法(基準人口を使用)と間接法(SMR)で年齢調整します。日本の統計では長らく1985年昭和60年モデル人口が使われていましたが、2020年以降は「2015年平成27年モデル人口」に更新されました。
SMR(標準化死亡比)
SMR(%) = 観察死亡数 ÷ 期待死亡数 × 100
「基準集団の年齢別死亡率を対象集団の年齢構成に適用した期待死亡数」に対し、実際の観察死亡数の比。SMR=100が基準集団と同等、150で1.5倍のリスク、80で2割少ないと解釈します。がん登録・職業病調査・地域保健で標準的に使われます。
PMR(比例死亡比)
特定死因の死亡が全死亡に占める比率。特定原因の重要性評価に使いますが、他死因の水準に依存する相対指標のため、SMRより解釈に注意が必要です。
研究デザイン(コホート/症例対照)
| デザイン | 方向 | 適用場面 | 算出指標 |
|---|---|---|---|
| 前向きコホート研究 | 曝露→疾病 | 予後研究・介入試験・一次予防 | 発生率・相対危険度・NNT |
| 後向きコホート研究 | 過去記録から追跡 | 職業病・環境曝露 | SMR・相対危険度 |
| 症例対照研究 | 疾病→曝露 | 稀な疾病・多要因 | オッズ比(OR) |
| ネステッドケースコントロール | コホート内でOR算出 | 大規模コホートの効率化 | OR(条件付き) |
| 横断研究 | 時点断面 | 有病率調査 | 有病率(prevalence) |
| 介入試験(RCT) | ランダム割付 | 治療効果評価 | HR・RR・NNT |
臨床研究のエビデンスレベルはRCT>コホート>症例対照>症例集積の順で、メタアナリシスがRCTを統合するとGRADE基準で最高評価になります。ただし稀な疾病・長期曝露評価では症例対照・コホートしか実施できない領域も多く、目的に応じた選択が必要です。
疾患別5年生存率の目安
国立がん研究センター院内がん登録・患者調査・全国がん罹患モニタリング集計(MCIJ)による、代表的疾患の5年相対生存率の目安です。医療機関・病期(ステージ)・年齢層で大きく変動します。
| 疾患 | 5年相対生存率 | 備考 |
|---|---|---|
| 甲状腺がん | 92-95% | 予後良好の代表 |
| 前立腺がん | 93-99% | 早期発見で極めて良好 |
| 乳がん | 92-93% | ステージⅠ99%、Ⅳ35% |
| 大腸がん | 72-75% | ステージⅠ97%、Ⅳ20% |
| 胃がん | 66-70% | 早期胃がん97% |
| 肺がん | 34-40% | ステージⅠ80%、Ⅳ5% |
| 肝がん | 35-40% | 肝硬変合併で低下 |
| 膵がん | 10-13% | 予後不良の代表 |
| 心筋梗塞(急性) | 85-90%(1年) | PCI・二次予防で改善 |
| 脳梗塞 | 70-80%(1年) | rt-PA・血栓回収療法 |
| 透析導入(慢性腎不全) | 55-60% | 年齢・合併症で大きく変動 |
症例数設計と検出力(Power)
臨床研究では事前の症例数設計(sample size calculation)が必須です。効果量・有意水準・検出力・脱落率から必要症例数を算出し、生物統計家がレビューします。
| 要素 | 標準設定 | 意味 |
|---|---|---|
| 有意水準α | 0.05(両側) | 第1種過誤の許容確率 |
| 検出力1-β | 0.80-0.90 | 差があるとき正しく検出する確率 |
| 効果量(HR等) | 試験デザインで指定 | 臨床的に意義のある差 |
| 追跡期間 | 1-10年 | イベント発生に必要な期間 |
| 脱落率 | 10-30% | 予測される脱落を上乗せ |
| アロケーション比 | 1:1が標準 | 群間の割付比率 |
Schoenfeld公式・Freedman公式が生存分析の標準式。R(powerSurvEpi)・PASS・nQueryなどのソフトで算出します。中間解析(interim analysis)を予定する場合、多重比較補正でαを分割する必要があります。
典型的な症例数の目安
5年生存率80%vs90%の比較(HR=0.53、両側α0.05、検出力80%、追跡5年、脱落10%)では片群約200例が必要。差が小さくなる(80%vs85%、HR=0.72)と片群1,000例超が必要になります。稀な疾病では多施設共同研究・国際共同研究で症例を集約するのが定石です。
業界・現場での使い方
臨床研究・治験
新薬・新治療法の有効性評価にKM法・Cox回帰を必須で使用。主要評価項目(全生存期間OS・無病生存期間DFS)の中央値、HR、95%CIをICH-E9に基づいて報告。PMDA・FDA申請の基礎データ。
公衆衛生・疫学
都道府県・保健所単位の疾病負担計算、生活習慣病対策の効果測定、がん検診の死亡率減少効果検証。年齢調整死亡率・SMRで地域比較・時系列比較を行う。
医療の質評価
DPCデータ・NDBデータを用いた病院間パフォーマンス比較。標準化死亡比(SMR)・在院日数調整で急性期医療の質を評価。医療機能情報提供制度の指標にも活用。
がん登録
全国がん登録・院内がん登録の5年相対生存率算出は国立がん研究センターが集計・公表。地域がん診療連携拠点病院の指定要件・機能評価にも生存率が使われる。
製薬企業・医療機器メーカー
市販後調査(PMS)・製造販売後臨床試験での長期安全性・有効性の評価。KM法による無増悪生存期間(PFS)・全生存期間(OS)は主要エンドポイント。
保険数理・生命保険
生命表・生存曲線・死亡率テーブルは生命保険料算定の基礎。標準生命表(日本アクチュアリー会)・簡易生命表(厚労省)を職務上使用。予定利率・予定死亡率で保険設計。
よくある間違い・注意点
- 打切りを無視して単純除算 — 追跡不十分な症例を「生存」扱いすると死亡率を過小評価。KM法・生命表法で打切りを適切に処理する必要があります。単純除算は追跡完了率が90%以上でないと大きく歪みます。
- 粗死亡率だけで集団比較 — 高齢者の多い地域の粗死亡率は当然高い。年齢調整死亡率・SMRで比較しないと集団の実質的リスクを見誤ります。老健地域・工業地域・都市部の比較で頻発。
- 比例ハザード仮定の未検証 — Cox回帰の前提「HRが時間非依存」が破れると結果に大きなバイアス。Schoenfeld残差プロットまたは対数-対数プロットで必ず検証します。
- 競合リスクの無視 — がん研究で心疾患死を「打切り」扱いすると累積発生を過大評価。Fine-Grey競合リスク回帰・累積発生関数(CIF)で対処します。
- 多重比較の未調整 — 多数のサブグループ解析でp値を出すと偶然の有意差(第1種の過誤)が増加。事前登録・Bonferroni・Holm補正・FDR制御で対応します。
- KM曲線の視覚的比較で結論 — 交差する曲線・末期の少数例による曲線の乱れは統計検定なしで判断すべきでない。ログランク検定・95%信頼区間・数字リスク表(number at risk)を必ず併記。
- 相対生存率と実測生存率の混同 — がん統計の「5年相対生存率」は同年齢一般集団と比べた相対値。実測は他死因も含む。用途に応じた指標選択が必要です。
- 症例対照研究でRRを算出 — 症例対照研究で相対危険度(RR)は原則算出できず、オッズ比(OR)が正しい指標。稀な疾病ではORとRRが近似します。
関連する規格・法規
- ICH-E9(統計原則) ― 臨床試験のための統計的原則。生存分析・多重比較・欠測値・中間解析の国際ガイドライン。
- ICH-E9(R1)(推定量) ― Estimand framework。目的と推定量の関係を明確化する新枠組み。
- CONSORT声明 ― ランダム化比較試験の報告基準。生存曲線・HRの記載要件。
- STROBE声明 ― 観察研究の報告基準。コホート・症例対照・横断研究の必須項目。
- ICD-10/ICD-11 ― WHO国際疾病分類。死因統計の標準分類。日本は2022年ICD-11採用移行中。
- がん登録等の推進に関する法律 ― 全国がん登録の法的根拠。全がん患者の生存率算出の基礎。
- 統計法 ― 人口動態統計・患者調査等の基幹統計の実施根拠。
- GCP省令 ― 医薬品の臨床試験の実施基準。データマネジメント・生存分析の品質基準。
- 人を対象とする生命科学・医学系研究に関する倫理指針 ― 研究倫理・IRB審査の統合ガイドライン。
参考文献・公的資料
統計手法・法制度・データベースの詳細は、下記の公的機関・学会の資料を参照してください。
- 厚生労働省 人口動態統計 ― 出生・死亡・婚姻・離婚・死産の月報・年報。死因別・年齢別の公式統計。
- 厚生労働省 簡易生命表・完全生命表 ― 平均余命・生存率テーブルの公式データ。
- 国立がん研究センター がん統計 ― 全国がん罹患・死亡・5年生存率の集計。院内がん登録・全国がん登録の連携公開。
- WHO ICD-10/ICD-11 ― 国際疾病分類の公式サイト。死因コード・診断分類。
- 医薬品医療機器総合機構(PMDA) ― 医薬品・医療機器の承認審査。ICH-E9の日本版適用ガイド。
- 日本疫学会 ― 疫学研究の学術団体。ガイドライン・教育資料。
- 日本計量生物学会 ― 臨床研究統計・生存分析の学術団体。
- 日本アクチュアリー会 ― 標準生命表・保険数理の公的団体。
- 厚生労働省 医学研究倫理 ― 人を対象とする研究の倫理指針の公式解説。
- EQUATOR Network ― CONSORT・STROBE等の報告ガイドラインの国際ハブ。
よくある質問（FAQ）
1,000人・5年追跡・50死亡の年死亡率は?
50÷1000÷5×100=1.0%/年。5年生存率は(1-0.01)⁵×100=約95.1%(指数モデル近似)。より正確にはKM法で描画し、95%信頼区間・数字リスク表を併記します。
KM法と生命表法の違いは?
KM法は「死亡が起きた各時点」で生存率を更新する連続時間モデル。生命表法(actuarial)は「一定区間(1年ごと等)」で計算する離散時間モデル。KM法の方が精度が高く、現代の臨床研究では標準です。生命表法は保険数理・大規模集団統計で使用。
打切り(censoring)って何?
患者が「イベント(死亡)を経験せずに観察が終了する状態」。研究終了時に生存、追跡不能、別死因の死亡等が該当。単純除算では「生存」扱いになり死亡率を過小評価するため、KM法・生命表法・Cox回帰で適切に処理します。
ハザード比(HR)って何?
Cox比例ハザードモデルで推定されるリスク比。HR=1.0で差なし、0.5でリスク半減、2.0で2倍リスク。「比例ハザード仮定(HRが時間非依存)」が前提で、95%信頼区間とp値、そしてSchoenfeld残差での仮定検証が必要です。
年齢調整死亡率はなぜ必要?
高齢者の多い地域は粗死亡率が当然高いため、年齢構成が異なる集団を比較すると誤解を招きます。年齢調整で「同じ年齢構成なら」というシナリオに揃えて比較します。日本は2015年モデル人口を使用。
SMR(標準化死亡比)って何?
期待死亡数に対する観察死亡数の比。基準集団の年齢別死亡率を対象集団の年齢構成に適用して期待値を計算。SMR=100が基準と同等、150で1.5倍リスク、80で2割少ない。がん登録・職業病調査・地域保健で標準的。
コホート研究と症例対照研究の違いは?
コホートは「曝露→疾病」の前向き追跡で相対危険度(RR)を算出。症例対照は「疾病→曝露」の後向き調査でオッズ比(OR)を算出。稀な疾病はコホートが非効率で症例対照が有利、一般的疾病はコホートが強力です。
5年相対生存率と実測生存率の違いは?
相対生存率は「同年齢一般集団と比較した相対値」で、他死因を除外した「がん自体の予後」を評価。実測生存率は他死因を含む実際の生存率。高齢者・長期予後の議論では相対生存率、若年・急性期は実測生存率を使い分けます。
ログランク検定って何?
2群以上の生存曲線が同一かを検定するノンパラメトリック検定。各死亡時点の期待死亡数と観察死亡数の差からカイ二乗統計量を算出。KM曲線比較の標準検定で、Wilcoxon-Gehan検定は初期の差に感度が高い代替です。
競合リスク(competing risk)って何?
目的とするイベント以外の原因でイベント発生が妨げられる状態。がん研究で心疾患死が競合リスク。KM法で「打切り」扱いすると累積発生を過大評価するため、Fine-Grey競合リスク回帰または累積発生関数(CIF)で対処します。
相対危険度(RR)とオッズ比(OR)の違いは?
RR=リスクの比(コホート研究向け)、OR=オッズの比(症例対照研究向け)。稀な疾病(発生率10%以下)ではRR≈OR。頻度が高いとORはRRから乖離しRRを過大評価する傾向。臨床解釈は原則RRが直感的です。
臨床研究に必要な症例数は?
効果量(HR)、有意水準α(通常0.05)、検出力(通常0.80)、追跡期間、脱落率から計算。5年生存率80%vs90%の比較で片群200例前後が典型。生物統計家がSchoenfeld公式・シミュレーションで事前算出します。
統計解析ソフトの選び方は?
臨床研究の主流はR(survivalパッケージ、無料、透明性高)とSAS(PROC LIFETEST/PHREG、規制申請向け、有償)。SPSSは初心者向け、EZR(自治医大)は日本語でR操作が可能。PMDA申請ではSASが多い一方、学術論文はRが増加中です。
症例数はどう決める?
効果量・有意水準α(0.05)・検出力(0.80-0.90)・追跡期間・脱落率(10-30%)から算出。5年生存率80%vs90%の比較で片群200例、80%vs85%で1,000例超が必要。Schoenfeld公式・Freedman公式が標準で、生物統計家のレビューが必須です。
中間解析の注意点は?
研究途中に効果判定を行う中間解析は多重比較補正でαを分割(O'Brien-Fleming・Pocock法等)。事前に解析計画書で明記し、独立データモニタリング委員会(IDMC)による独立審査が国際基準です。
臨床研究の統計チェックリスト
論文投稿・PMDA申請・研究計画審査で必ず確認される項目:
- 研究デザインの明示 ― 前向き/後向き、単施設/多施設、無作為化の有無、盲検化(単盲検/二重盲検)
- 主要評価項目と副次評価項目 ― 事前登録(UMIN・ClinicalTrials.gov)との整合
- 症例数計算根拠 ― 効果量・α・検出力・脱落率の設定根拠
- 統計解析計画書(SAP) ― 解析対象集団(ITT/FAS/PPS)、欠測値処理、多重比較補正
- KM曲線と数字リスク表 ― 曲線図・95%CI・打切りマーク・時点別リスク集団数
- 比例ハザード仮定の検証 ― Schoenfeld残差プロット・対数-対数プロット
- 感度分析 ― 欠測値補完・per-protocol解析・サブグループ解析
- 臨床的意義との対応 ― 統計的有意差と臨床的有意差の区別