[論文レビュー] Bayesian Analysis of Modified Weibull distribution under progressively censored competing risk model

Arabin Kumar Dey, Abhilash Jha+1 more|arXiv (Cornell University)|May 21, 2016

Statistical Distribution Estimation and Applications参考文献 15被引用数 5

ひとこと要約

本稿では、2つの故障原因を伴う進行的型II打ち切り競合リスクモデル下での修正ワイブル分布のベイズ解析を提案する。参考事前分布のもとで、後方予測サンプルを生成するために決定的走査ギブス・サンプリングとスライス・サンプリングを組み合わせた手法を採用し、フォリクル細胞リンパ腫のシミュレーションおよび実データ解析を通じて、パラメータ λ₁, λ₂, α, β の中央値とHPD区間を計算することで、推定の頑健性を示している。

ABSTRACT

In this paper we study bayesian analysis of Modified Weibull distribution under progressively censored competing risk model. This study is made for progressively censored data. We use deterministic scan Gibbs sampling combined with slice sampling to generate from the posterior distribution. Posterior distribution is formed by taking prior distribution as reference prior. A real life data analysis is shown for illustrative purpose.

研究の動機と目的

• 進行的型II打ち切り下で修正ワイブル分布を用いたパラメトリックな競合リスク分析のためのベイズフレームワークの構築。
• 参考事前分布の非可解性に起因する後方分布の非可解性を克服するため、ギブス・サンプリングとスライス・サンプリングを組み合わせた後方分布推定手法の開発。
• シミュレーションスタディおよびフォリクル細胞リンパ腫の実データ応用を通じて、提案手法の性能を評価。
• 競合リスク下での修正ワイブル分布パラメータのベイズ推定量および最高後方密度（HPD）区間の推定。

提案手法

• 2つの競合原因を伴う進行的型II打ち切りデータの尤度関数を、各原因に対して独立な修正ワイブル分布を仮定して導出する。
• パラメータ (λ₁, λ₂, α, β) に対して参考事前分布を用いることで、客観的ベイズ推論を実現するが、これにより後方分布の表現が非可解となる。
• 非可解な後方分布を扱うために、決定的走査ギブス・サンプリングとスライス・サンプリングを組み合わせたハイブリッドアルゴリズムを用いて後方分布のシミュレーションを実施する。
• 修正ワイブル分布は、生存関数 $ S(t) = \exp\left( -\lambda \alpha \left(1 - \exp\left( \left(\frac{t}{\alpha}\right)^\beta \right) \right) \right) $ で定義され、パラメータ λ, α, β を持つ。
• スライス・サンプリングでは、ステップアウト手順を用い、混合性と収束性を向上させるために幅の適応的選択を実施する。
• 後方要約統計量には中央値、HPD区間、およびモデル妥当性の検証に用いるトレースプロットが含まれる。

実験結果

リサーチクエスチョン

• RQ1 提案されたベイズ手法は、進行的型II打ち切りと競合リスク下での修正ワイブル分布のパラメータをどの程度正確に推定できるか？
• RQ2 ギブス・サンプリングとスライス・サンプリングの組み合わせは、参考事前分布に起因する非可解な後方分布を効果的に扱えるか？
• RQ3 さまざまな打ち切りスキームと標本サイズ下で、後方推定値と信用区間はどのように振る舞うか？
• RQ4 本モデルは、再発と死亡という2つの原因を伴う実世界の生存データ（例：フォリクル細胞リンパ腫）に対して適切なフィットを示すか？

主な発見

• シミュレーションスタディ1では、λ₁, λ₂, α, β の後方中央値はそれぞれ 0.7918, 0.6261, 0.2063, 0.0833 であり、収束性と精度が良好であることが示された。
• シミュレーションスタディ2では、λ₁ = 0.6300, λ₂ = 0.6987, α = 0.2155, β = 0.0978 であり、さまざまなシナリオでも一貫した性能を示した。
• シミュレーションスタディ3では、λ₁ = 0.9181, λ₂ = 0.8261, α = 0.3063, β = 0.0833 であり、安定したトレースプロットとヒストограмが収束を確認した。
• フォリクル細胞リンパ腫データ（ケース2）では、後方中央値は λ₁ = 0.034, λ₂ = 0.004, α = 30.228, β = 0.536 であり、α のHPD区間が広く、高い不確実性を示した。
• 原因別データのヒストограмと適合した修正ワイブル分布の比較により、分布仮定の妥当性が確認された。
• 全パラメータのトレースプロットとヒストограмは、混合性と収束性が良好であり、ギブス・スライス・サンプリング手法の信頼性を裏付けた。