NTDBとTQPにおける生理学的データの欠落：完全症例分析は実際にどの程度外傷死亡率に偏りをもたらすのか

掲載 2026 年 9 月 20 日

完全症例分析は、外傷登録研究における暗黙のデフォルトです。GCS、血圧、または呼吸数が欠落している行を削除し、残りのデータでモデルを適合させ、結果を報告します。2018年に行われた539件の外傷登録論文の系統的レビューでは、半数以上がデータの完全性について全く言及していないか、欠落データについて言及しているものの、その処理方法については述べていないことがわかりました。つまり、文書化されていない完全症例分析がデフォルトとして行われていることは、例外というよりむしろ標準に近いということです。私たちは、このデフォルトが実際にどれだけのコストがかかるかを正確に把握するために、2つの統計エンジンで実行した独立したシミュレーションを構築しました。

なぜこの欠落は偶然ではないのか

救急外来のGCS、収縮期血圧、呼吸数は、臨床医が患者を十分に評価できなかったために欠落しており、評価が最も困難な患者は、気管挿管、鎮静、または重篤な状態で到着した患者に偏っている。ビクトリア州外傷登録簿の2,520症例の分析により、このパターンが直接確認され、呼吸数、GCS、気管挿管の修飾フラグ、収縮期血圧が、標準的な外傷リスク調整に必要な変数の中で最も多くの欠落観測値を占めており、どの患者に欠落があるかは偶然ではなく臨床的論理によって予測されることが判明した。統計学者はこれを「ランダムではない欠落」と呼ぶ。なぜなら、値が欠落する確率は値自体に依存し、重度の外傷が鎮静と気管挿管を引き起こすのとまったく同じ割合で、低いGCSと不安定なバイタルサインがカルテから削除されるからである。完全症例分析は暗黙のうちにその逆を仮定している。

シミュレーション

我々は、複製ごとに 10,000 人の外傷入院患者からなる合成コホートを生成した。このコホートには、部分的に相関する 2 つの潜在的重症度コンポーネントがあり、1 つは解剖学的損傷の重症度と血行動態の状態 (ISS、収縮期血圧、呼吸数) を、もう 1 つは神経学的状態 (GCS) を決定づけている。これは、実際の多発外傷で見られる解剖学的損傷と神経学的損傷の間の相関が完全ではなく中程度であることと一致する。真の死亡率は、年齢、メカニズム、ISS、GCS、収縮期血圧、および呼吸数を組み込んだロジスティック モデルから生成された。GCS、収縮期血圧、および呼吸数の欠損は、文献で説明されているとおり、患者自身の真の生理学的値と挿管指標の関数として課せられ、欠損率は現実的な範囲、すなわち GCS で約 30%、収縮期血圧で約 24%、呼吸数で約 16% に収まった。これは、公表されている外傷登録監査で報告されているこれら 3 つの変数の欠損率と桁が似ている。

次に、各複製に対して 4 つの分析アプローチを実行しました。 1 つは、欠損値がない場合に研究者が実行する分析の代わりとなる、真の、完全に観測された生理学的値に適合するオラクル モデルです。 2 つは、3 つの変数のいずれかが欠損しているすべての患者を除外する完全ケース モデルで、文書化されていないデフォルトの慣行を反映しています。 3 つは、欠損値をその変数の全体平均で置き換える単純な単一平均代入モデルで、投稿された原稿でまだ見られるショートカットです。 4 つは、NTDB 固有の文献で推奨されているアプローチに従う適切な多重代入モデルで、死亡結果自体が代入モデルに含まれ、結果は Rubin のルールによって代入間でプールされます。 パイプライン全体は、異なる代入実装 (scikit-learn の反復代入と R の mice パッケージ）を使用し、Pythonで300回、Rで100回繰り返し実行します。以下の各効果の方向と概算サイズの両方について、両者の一致が結果を保証するものであり、単一のスクリプトの特異性ではありません。

発見：完全症例分析により報告された死亡率がほぼ半減

シミュレーションされたコホート全体の実際の死亡率は平均5.50%でした。GCS、収縮期血圧、呼吸数がすべて揃っている患者に分析を限定すると（これは記録されていない完全症例分析で行われていることと全く同じです）、観察された死亡率は2.68%に低下し、半分以下に減少しました。この差は統計的なものではなく、機械的なプロセスによって生じています。生理学的評価を完全に記録するには状態が不安定すぎるという理由で患者を分析から除外することは、死亡する可能性が最も高い患者を除外することと全く同じであり、実際に治療された集団よりも体系的に健康なサブグループに基づいて計算された死亡率を、どれだけ高度な下流モデリングを用いても修正することはできません。この証拠に基づけば、完全症例外傷コホートからの生の、あるいはわずかに調整された死亡率を引用するベンチマーク研究や品質レポートは、丸め誤差レベルのずれではなく、実際の死亡率を大幅に過小評価した数値を報告していることになります。

我々の結果は、68の病院にわたる63,020人の実際のNTDB患者に基づいて構築されたモンテカルロシミュレーションが、関連するものの異なる疑問、すなわち欠損データ処理方法の選択によって、観察死亡率と期待死亡率の差に基づいてどの病院が品質外れ値と判定されるかが変わるかどうかという疑問に対して得た結果と一致している。その先行研究では、基礎となるコホートを固定したまま分析方法を変更していた。我々の研究では、基礎となる欠損メカニズム自体を合成的に変更することで、推定される真の値が常に既知となるようにしている。これにより、方法間の順位の不一致だけでなく、正確なバイアスを述べることが可能になる。

発見その２：調整後の効果は、完全データ制約下でも予想以上に良好な結果を示した

モデルの対数オッズ係数として表される、GCSの死亡率に対するリスク調整効果は、オラクルモデルでは平均-0.300でした。3つの変数すべてが観察された患者の約半数に限定された完全症例モデルでは、係数は平均-0.301となり、オラクル値に十分近いため、損失は方向性バイアスというよりも、より小さく均質なサンプルによる精度の低下として説明する方が適切です。この結果の背後には、特定の有用な非対称性があります。この状況での完全症例分析による害は、ほとんどの外傷論文が実際に主要な発見として報告している条件付きリスク調整係数ではなく、報告された死亡率と症例構成という周辺量に集中しています。論文のGCSの調整オッズ比は、記載されているコホート死亡率が大きく間違っている場合でも、ほぼ正しい可能性があります。これは、抄録の主要な統計だけをざっと見る査読者が見逃してしまう種類の選択的損傷です。

発見３：単純な平均値補完は、完全ケース分析が保持する調整効果を弱める

欠損している GCS、収縮期血圧、呼吸数の値を各変数の全体平均で置き換えるという、少なくともすべての患者を分析対象に残すための近道は、平均 GCS 係数が-0.165となり、オラクル値と比較して約半分の減衰となった。その理由は、微妙なものではなく、直接的なメカニズムによって説明できる。欠損しているすべての GCS を同じ定数で埋めると、回帰分析が死亡率への影響を推定するために必要な GCS の変動が正確に破壊され、しかも、真の GCS が最も極端だった患者ほど、そもそも欠損している可能性が最も高いため、その影響が最も大きくなる。平均値による補完は、完全症例分析の死亡率バイアスを、研究が実際に推定しようとしている調整済み効果の、異なる、そしてこのシミュレーションではより大きなバイアスと交換する。

発見4：結果を含めた多重代入により、損失の大部分が回復する

多重代入法を代入モデルに死亡結果を含めて適用し、ルービンのルールに従って代入間でプールしたところ、GCS係数の平均値は-0.263となり、単純な平均代入法で残されたギャップの大部分を回復するとともに、完全症例分析のように最も重症の患者を除外するのではなく、すべての患者を死亡率の計算に含めました。これは、NTDB固有の方法論文献が少なくとも2009年以来推奨してきた、結果を含む多重代入法という設計選択と同じであり、本シミュレーションでは、実際のレジストリの特異性ではなく、メカニズムを分離するために構築された完全に合成されたデータを使用して、単一のデータセットとは独立して、その推奨が成り立つ理由を再現しています。

これを自分のレジストリ抽出に適用する

4つのチェック項目は、このシミュレーションを実際の参加者使用ファイルに適用するものです。GCS、収縮期血圧、呼吸数の欠損率を個別に定量化し、問題の原因となっている変数が隠蔽される単一の統合された完全性数値を報告しないようにします。各変数の欠損が、挿管フラグ、ヘリコプターによる到着、ISS自体など、完全に観察された他の重症度マーカーと関連しているかどうかをテストします。明確な関連性は、このシミュレーションでモデル化された情報のある欠損の実際的な特徴だからです。死亡結果と、研究で報告されているその他の結果を、多重代入モデルの予測因子として含めて多重代入を実行します。結果は分析段階のみに属するという仮定に基づいて除外しません。また、レビュー担当者が修正の大きさを直接確認できるように、多重代入結果の代わりに完全ケースの結果を報告するのではなく、完全ケースの結果と並べて報告します。

この種のチェックに使える無料ツール

• 臨床レジストリデータをクリーニングするための6つの手法欠損値監査は、補完決定に先立って実施されるべきである。
• STROBE報告チェックリストチュートリアル欠損データフローを原稿自体にどのように記載すべきかについては、
• NTDBとTQPデータベースの説明TQP研究が対処する必要のある、欠損値以外のより広範な設計上の問題に対応するため

原稿にとってこれは何を意味するのか

3つの実用的な結論が直接導き出されます。外傷登録コホートから得られた死亡率または症例致死率については、分母がどのように構築されたかを特に精査する必要があります。なぜなら、完全症例制限によって、出力にエラーメッセージや明らかな警告が表示されることなく、報告された率が半分に減少する可能性があるからです。生理学的共変量の調整オッズ比は、生の率よりも完全症例処理に対して比較的頑健ですが、サンプルサイズが減少すると精度が低下し、信頼区間が広がります。また、不完全な症例を削除するか、サンプル平均で欠損値を埋めるという2つの一般的な近道のうち、結果変数を含めた多重代入法がどちらよりも優れており、これはEditverse独自のNTDBおよびTQP外傷登録研究方法論が、実施するすべてのプロジェクトで採用している方法です。

Python 300回とR 100回の複製における結果概要

完全症例解析のみが報告死亡率を変更する。なぜなら、この方法は、患者に関する欠落情報を補完するのではなく、解析対象となる患者を変更する唯一の方法だからである。2つの補完法のみが死亡率計算に全コホートを保持し、そのうち、結果包含多重補完法は、リスク調整済み効果量も保護する。

このシミュレーションの限界

これは、1 つのメカニズムを明確に分離するために作成された合成図であり、実際の TQP データの再分析ではありません。Editverse は TQP 参加者の使用ファイル契約を保持しておらず、クライアントに代わってデータセットを取得していないためです。欠損率、相関構造、および死亡率モデルは、特定の登録年に合わせて調整されたものではなく、全体を通して引用されている公開文献に基づいて現実的な範囲に収まるように選択されています。クライアントとの契約では、クライアント自身の抽出データに基づいて作業を行い、そのファイルの実際の欠損パターンを直接特徴付け、NTDB および TQP の説明で説明されているとおり、主要な多重代入結果とともに完全ケース感度分析を報告します。

Professional Services

Editverseは、NTDBおよびTQPのトラウマ転帰研究を構築し、生理学的欠損は転帰を含む多重代入法によって処理され、完全症例感度分析が併せて報告されます。これは、ここで検証されたとおりです。

参考情報

1. Shivasabesan G、Mitra B、O'Reilly GM。外傷登録における欠損データ：系統的レビュー。 損傷. 2018;49(9):1641-1647. doi: 10.1016/j.injury.2018.03.035
2. O'Reilly GM、Cameron PA、Jolley DJ。どの患者に欠損データがあるか？外傷登録における欠損データの分析。 損傷. 2013;43(11):1917-1923. doi: 10.1016/j.injury.2012.07.185
3. Moore L、Hanley JA、Turgeon AF、Lavoie A、Emond M. 国立外傷データバンクにおける欠損生理学的データの補完のための多重代入モデル。 アメリカ外科学会誌. 2009;209(5):572-579. doi: 10.1016/j.jamcollsurg.2009.07.004
4. Glance LG、Osler TM、Mukamel DB、Meredith W、Dick AW。欠損データの処理における統計的手法が外傷センターの質に与える影響。 外科の年報. 2009;249(1):143-148. doi: 10.1097/SLA.0b013e31818e544b
5. Newgard CD. 州の外傷登録における院外データの欠損に対する多重代入法の妥当性。 アカデミック救急医学. 2006;13(3):314-324. doi: 10.1197/j.aem.2005.09.011
6. O'Reilly GM、Jolley DJ、Cameron PA、Gabbe B. 行方不明：外傷システムベンチマークへの欠損データ処理方法の適用に関するケーススタディ。 アカデミック救急医学. 2011;17(10):1122-1129. doi: 10.1111 / j.1553-2712.2010.00887.x