[论文解读] Weibull Generalized Exponential Distribution

Abdelfattah Mustafa, B. S. El-Desouky+1 more|arXiv (Cornell University)|Jun 2, 2016

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

一句话总结

该论文提出威布尔广义指数分布（WGED），一种具有三个参数的灵活寿命分布，通过威布尔-G族生成器扩展指数分布，以建模浴盆形风险率。在拟合真实玻璃纤维强度数据方面，WGED在AIC、BIC、HQIC和柯尔莫哥洛夫-斯米尔诺夫统计量方面均优于现有模型（如ED、GED、BED和BGED），表现出更优的拟合效果。

ABSTRACT

This paper introduces a new three-parameters model called the Weibull-G exponential distribution (WGED) distribution which exhibits bathtub-shaped hazard rate. Some of it's statistical properties are obtained including quantile, moments, generating functions, reliability and order statistics. The method of maximum likelihood is used for estimating the model parameters and the observed Fisher's information matrix is derived. We illustrate the usefulness of the proposed model by applications to real data.

研究动机与目标

• 开发一种新的三参数分布，通过提升对复杂风险率形状的建模灵活性，扩展指数分布。
• 实现对具有浴盆形或单峰故障率的数据的建模，而标准指数分布和威布尔分布无法有效捕捉这些特征。
• 通过在指数基分布中引入威布尔-G生成器，为现有广义指数族提供更具弹性的替代方案。
• 推导并分析关键统计特性，包括矩、分位数、顺序统计量和可靠性函数。
• 在真实世界数据上评估模型性能，并通过多种信息准则和拟合优度检验，与已建立的分布进行比较。

提出的方法

• WGED通过威布尔-G族生成器构建，其中基分布为尺度参数为λ的指数分布。
• 累积分布函数（CDF）定义为 $ F(x;a,b, heta) = 1 - \text{exp}\big(-a \big[\frac{G(x)}{1-G(x)}\big]^b \big) $，其中 $ G(x) = 1 - e^{-\theta x} $，且参数 $ a, b, \theta > 0 $。
• 概率密度函数（PDF）推导为 $ f(x) = ab\theta e^{-\theta x} \frac{(1 - e^{-\theta x})^{b-1}}{(e^{-\theta x})^{b+1}} \text{exp}\big(-a (1 - e^{-\theta x})^{-b} \big) $，可实现灵活的风险率形状。
• 通过级数展开和积分表示，解析推导出矩、生成函数、分位数和顺序统计量等统计特性。
• 采用最大似然估计（MLE）进行参数估计，并推导出观测信息矩阵，以计算标准误和置信区间。
• 通过真实玻璃纤维强度数据集（Smith和Naylor，1987年）上的AIC、BIC、HQIC、柯尔莫哥洛夫-斯米尔诺夫（K-S）检验和对数似然值，开展模型比较。

实验结果

研究问题

• RQ1 威布尔-G族生成器能否产生一种灵活的三参数扩展指数分布，以有效建模浴盆形风险率？
• RQ2 WGED的统计特性（如矩、分位数和顺序统计量）与基线分布及现有广义指数模型相比如何？
• RQ3 WGED在拟合真实世界可靠性数据方面是否显著优于现有模型（如ED、GED、BED和BGED）？
• RQ4 在参数估计的准确性与标准误估计方面，最大似然估计在WGED上的表现如何？
• RQ5 在经验数据上，信息准则（AIC、BIC、HQIC）和拟合优度检验（K-S）如何对WGED与竞争模型进行排序？

主要发现

• WGED模型表现出浴盆形风险率函数，适用于建模具有早期失效、稳定阶段和老化相关失效的系统。
• 在玻璃纤维强度数据集中，WGED在所有拟合模型中取得最低的AIC（35.656）、BIC（42.085）、HQIC（38.185）和K-S统计量（0.127）。
• WGED的对数似然值为-14.828，是所有模型中最高（最不负面）的，表明其在似然性方面拟合最佳。
• 参数 $ a $ 的95%置信区间为[0, 175.11]，$ b $ 为[3.989, 57.785]，$ \theta $ 为[0.146, 0.298]，方差-协方差矩阵表明参数具有稳定性。
• 轮廓对数似然图证实了所有三个参数的MLE解具有唯一性。
• Kaplan-Meier生存估计与拟合CDF图的可视化结果表明，WGED在所有竞争模型中与经验数据的拟合最为接近。