[論文レビュー] Generalized inverse xgamma distribution: A non-monotone hazard rate model

Harsh Tripathi, Abhimanyu Singh Yadav+2 more|arXiv (Cornell University)|Nov 3, 2018

Statistical Distribution Estimation and Applications参考文献 10被引用数 4

ひとこと要約

本稿では、逆xガンマ分布のパワー変換により導かれる、非単調なハザードレートを示す柔軟な2パラメータの寿命モデルである一般化逆xガンマ分布（GIXGD）を提案する。このモデルは、キビグマの生存データに適用され、AIC、BIC、コルモゴロフ＝スミルノフ統計量を用いて、他の競合分布よりも優れた適合度を示した。パラメータおよび信頼性関数の推定には最尤推定法が用いられた。

ABSTRACT

In this article, a generalized inverse xgamma distribution (GIXGD) has been introduced as the generalized version of the inverse xgamma distribution. The proposed model exhibits the pattern of non-monotone hazard rate and belongs to family of positively skewed models. The explicit expressions of some distributional properties, such as, moments, inverse moments, conditional moments, mean deviation, quantile function have been derived. The maximum likelihood estimation procedure has been used to estimate the unknown model parameters as well as survival characteristics of GIXGD. The practical applicability of the proposed model has been illustrated through a survival data of guinea pigs.

研究の動機と目的

• 逆xガンマ分布をパワー変換を用いて一般化することで、より柔軟な寿命分布の構築を目的とする。
• 信頼性および臨床研究でよく見られる非単調ハザードレートをモデリングすることを目的とする。
• モーメント、分位数関数、信頼性特性などの主要統計的特性の明示的表現を導出することを目的とする。
• 実データへの適合と既存の寿命分布との比較を通じて、モデルの性能を評価することを目的とする。
• 生存分析における1および2パラメータの逆分布族に対する強力な代替モデルを提供することを目的とする。

提案手法

• X ∼ IXGD(θ) である逆xガンマ分布を基盤とし、Y = X^{1/α} によるパワー変換を用いて一般化逆xガンマ分布（GIXGD）を提案し、形状母数αを導入する。
• GIXGDの確率密度関数（PDF）および累積分布関数（CDF）を導出し、α = 1 のときにはIXGDに帰着することを示す。
• 生存関数 S(y;α,θ) = 1 - F(y;α,θ) およびハザードレート関数 H(y;α,θ) = f(y;α,θ)/S(y;α,θ) を定義する。
• モーメント、逆モーメント、条件付きモーメント、平均偏差、分位数関数の解析的表現を導出する。
• パラメータ α および θ の推定に最尤推定法（MLE）を用い、生存関数およびハザード関数の推定には不変性の性質を応用する。
• キビグマの生存データセットを用いて、コルモゴロフ＝スミルノフ（K-S）検定および情報量基準（AIC、BIC、CAIC、HQIC）を用いてモデルの性能を検証する。

実験結果

リサーチクエスチョン

• RQ1 既存の逆分布族と比較して、非単調ハザードレートをより柔軟にモデリングできる一般化逆xガンマ分布（GIXGD）を構築できるか？
• RQ2 GIXGDの統計的特性（モーメント、分位数関数、信頼性指標など）は、基盤となる逆xガンマ分布と比較してどのように異なるか？
• RQ3 GIXGDは、逆ワイブル分布、逆リンドレー分布、一般化指数分布といった競合分布よりも、実生存データに対してより良い適合度を示すか？
• RQ4 さまざまな打ち切り条件および標本サイズ下で、MLEがGIXGDのパラメータおよび信頼性関数を推定する際に示す性能特性は何か？
• RQ5 非単調ハザードパターンが観察される状況において、GIXGDは信頼性および生存分析における実用的代替モデルとして有効であるか？

主な発見

• GIXGDは非単調なハザードレートを示し、初期には増加し、その後低下する傾向を示す。この特性により、信頼性および臨床研究における複雑な故障パターンのモデリングに適している。
• キビグマの生存データに対して、6つの競合分布と比較してGIXGDが優れた適合度を示した。AIC（787.98）、BIC（792.54）、CAIC（793.54）、HQIC（789.79）、K-S統計量（0.1432）のすべてにおいて、他のモデルより最小値を記録した。
• 72件の観察データに基づき、形状母数および尺度母数の最尤推定値はそれぞれ \hat{\alpha} = 1.624 および \hat{\theta} = 641.75 であった。
• y = 54, 70, 99, 112日における生存確率の推定値は、それぞれ0.626、0.475、0.308、0.260であり、対応するハザードレートは0.0177、0.0165、0.0136、0.0124であった。
• 一般化指数分布（GED）、逆ワイブル分布（IWD）、ガンマ分布（GD）と比較して、GIXGDは優れた性能を示した。特に、ガンマ分布は非常に高いK-S統計量（0.998）を示し、適合度が著しく低いことが判明した。
• 本研究は、GIXGDが正の歪度を示し、非単調ハザードレートを示すデータを生存分析でモデリングする強力な候補であることを確認した。特に、基盤となるIXGDでは不十分に柔軟な場合に有効である。