[論文レビュー] A novel high-power test for continuous outcomes truncated by death

Andreas Kryger Jensen, Theis Lange|arXiv (Cornell University)|Oct 27, 2019

Statistical Methods and Inference参考文献 3被引用数 7

ひとこと要約

本論文は、死亡によって打ち切りられる連続的アウトカム（例：ランダム化比較試験における生活機能（QoL）スコア）に対して、著者らが提案する新しい高パワーな統計的検定を提示する。死亡した被験者には最低スコアが割り当てられるが、この方法は生存とQoLアウトカムを同時にモデル化する半パラメトリックな実現尤度アプローチを用いる。この手法により、完全な統計的パワーを維持し、信頼区間を伴う効果推定値が得られ、従来のウィルコクソン検定よりも優れている。ウィルコクソン検定は、死亡被験者に割り当てられる同一スコアによるタイ（同一値）と、ずれた帰無仮説により、深刻なパワー損失を受ける。

ABSTRACT

Patient reported outcomes including quality of life (QoL) assessments are increasingly being included as either primary or secondary outcomes in randomized controlled trials. While making the outcomes more relevant for patients it entails a challenge in cases where death or a similar event makes the outcome of interest undefined. A pragmatic - and much used - solution is to assign diseased patient with the lowest possible QoL score. This makes medical sense, but creates a statistical problem since traditional tests such as t-tests or Wilcox tests potentially looses large amounts of statistical power. In this paper we propose a novel test that can keep the medical relevant composite outcome, but preserve full statistical power. The test is also applicable in other situations where a specific value (say 0 days alive outside hospitals) encodes a special meaning. The test is implemented in an R package which is available for download.

研究の動機と目的

• ランダム化比較試験（RCTs）において、死亡被験者に最低のQoLスコアを割り当てる標準的手法が引き起こす統計的パワー損失を是正すること。
• 生存とQoLスコアの複合アウトカム（例：死亡を最小スコアとしてコード化）の臨床的解釈可能性を保ちながら、完全なパワーを維持する手法を開発すること。
• 生存と連続的アウトカムの組み合わせに対する治療効果の単一のp値と、信頼区間を伴う効果推定値を提供すること。
• ベースライン共変量を組み込むことで、RCTを超えて観察研究への適用可能性を拡大すること。
• ウィルコクソン型検定の限界を克服すること。これらの検定は、アウトカム分布に特異的成分が存在するため無効であり、帰無仮説がずれている。

提案手法

• アウトカムを二成分混合分布としてモデル化：生存者（Y|A=1）には連続分布、非生存者（A=0）には固定値εにおけるデルタ関数（点質量）を割り当てる。
• 生存確率と生存者における平均QoLに対する治療効果のない帰無仮説を同時に検定する。帰無仮説は H₀: E[Y|A=1,R,X] = E[Y|A=1,R=0,X] および P(A=1|R,X) = P(A=1|R=0,X) として定式化される。
• 連続的アウトカム部にパラメトリックな形を仮定せず、半パラメトリック実現尤度を用いて結合分布を推定する。
• ベースライン共変量Xを条件付き推定方程式を用いて組み込み、交絡要因の調整を可能にする。
• 実現尤度に基づく尤度比検定を用いて、検定統計量とp値を導出。同時に、効果対比の信頼領域を提供する。
• パラメトリックおよび半パラメトリックな両バージョンを用いて、生存およびQoLの両成分に対する効果推定値と信頼区間を提供する。

実験結果

リサーチクエスチョン

• RQ1 高死亡率を示すRCTにおいて、死亡被験者に最低のQoLスコアを割り当てる標準的手法は、統計的パワーを著しく損なうか？
• RQ2 特異的成分（固定値ε）を含むアウトカム分布に対しても、生存とQoLを同時にモデル化する統合的検定が、完全なパワーを維持できるか？
• RQ3 さまざまな分布仮定と治療効果パターン下で、提案手法の半パラメトリック実現尤度法は、ウィルコクソン検定およびt検定に比べてパワーで優位か？
• RQ4 ベースライン共変量の調整を伴いながら、解釈可能な効果推定値と信頼区間を提供できるか？
• RQ5 生存者におけるQoLスコアの重い尾を持つ（裾が長い）分布に対して、本手法は頑健性を示すか？

主な発見

• 設定1（死亡率に変化なし、生存者におけるQoLが向上）では、本手法がウィルコクソン検定よりも顕著に高いパワーを達成した。ウィルコクソン検定は、死亡被験者に割り当てられた固定値によるタイの影響を強く受けていた。
• 設定2（QoLに変化なし、治療群で死亡率が低下）では、ウィルコクソン検定が1自由度の検定であるためわずかに高いパワーを示したが、本手法は依然として妥当な推論と効果推定値を提供した。
• 設定3および4（生存者におけるQoL向上と死亡率上昇という相反する効果）では、ウィルコクソン検定は明確な治療効果にもかかわらず非常に低いパワーを示したが、本手法は高いパワーを維持した。
• 重い尾を持つQoL分布（設定4）では、半パラメトリックアプローチがパラメトリックバージョンを著しく上回り、優れた頑健性を示した。
• 本手法は、治療効果のない帰無仮説を強く棄却した（p < 0.001）。一方、ウィルコクソン検定およびt検定は有意性を検出できず、本手法の感度の高さが浮き彫りになった。
• 治療対比の同時信頼領域は原点を含まず、有意な統合的効果を確認した。信頼区間は情報量が多く、解釈可能であった。