[论文解读] 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

一句话总结

本文提出了一种在具有两种失效原因的逐步右删失竞争风险模型下，对修正威布尔分布进行贝叶斯分析的方法。该方法结合确定性扫描吉布斯抽样与切片抽样，在参考先验下生成后验样本，通过模拟研究和滤泡性细胞淋巴瘤的真实数据分析展示了稳健的估计性能，其中对参数 λ₁、λ₂、α 和 β 计算了后验中位数和最高后验密度（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.

研究动机与目标

• 开发一种基于修正威布尔分布的逐步右删失竞争风险参数化贝叶斯分析框架。
• 通过结合吉布斯抽样与切片抽样来处理参考先验导致的不可解析后验分布，以应对贝叶斯推断中非可解析参考先验的挑战。
• 通过模拟研究和滤泡性细胞淋巴瘤的真实数据应用，评估所提方法的性能。
• 在竞争风险下估计贝叶斯估计量和最高后验密度（HPD）区间，针对修正威布尔分布的参数。

提出的方法

• 推导了具有两种竞争失效原因的逐步右删失数据的似然函数，假设每种原因均服从独立的修正威布尔分布。
• 为参数 (λ₁, λ₂, α, β) 使用参考先验以确保客观贝叶斯推断，尽管这导致后验表达式不可解析。
• 通过结合确定性扫描吉布斯抽样与切片抽样的混合算法进行后验模拟，以处理不可解析的后验分布。
• 修正威布尔分布由生存函数 $ S(t) = \exp\left( -\lambda \alpha \left(1 - \exp\left( \left(\frac{t}{\alpha}\right)^\beta \right) \right) \right) $ 定义，参数为 λ、α、β。
• 切片抽样使用逐步外推程序，并采用自适应宽度选择以改善混合效果与收敛性。
• 后验汇总包括中位数、HPD区间以及用于模型验证的轨迹图。

实验结果

研究问题

• RQ1 所提出的贝叶斯方法在逐步右删失竞争风险下，对修正威布尔分布参数的估计效果如何？
• RQ2 吉布斯抽样与切片抽样结合是否能有效处理由参考先验引起的不可解析后验分布？
• RQ3 在不同删失方案与样本大小下，后验估计与可信区间的表现如何？
• RQ4 该模型是否能很好地拟合现实世界生存数据，例如具有复发与死亡竞争风险的滤泡性细胞淋巴瘤数据？

主要发现

• 在模拟研究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 区间较宽，表明不确定性较高。
• 通过将特定原因的生存数据直方图与拟合的修正威布尔分布进行比较，验证了模型拟合效果，确认了分布假设的合理性。
• 所有参数的轨迹图与直方图均显示良好的混合性与收敛性，支持了吉布斯-切片抽样方法的可靠性。