基于恒定应力加速寿命试验与渐进II型删失下新奇林德利半逻辑分布可靠性推断的统计视角

《Scientific African》：Investigating Accelerated Life Tests in Reliability Theory via Statistical Inference View

　　

本研究在渐进II型删失（progressive Type-II censoring）下的恒定应力加速寿命试验（constant-stress accelerated life testing, CS-ALT）框架中，基于新Odd Lindley Half Logistic（New OLiHL）分布建立了统计推断体系。该模型适用于呈右偏且风险函数（hazard rate function, HRF）递增的寿命数据，并在不同应力水平间引入了寿命-应力关系。研究人员推导了似然函数，并采用最大似然估计（maximum likelihood estimation, MLE）、基于平方误差损失函数下伽马先验的贝叶斯马尔可夫链蒙特卡洛（Bayesian Markov Chain Monte Carlo, MCMC）方法，以及Boot-P与Boot-T两种自助法（bootstrap）对模型参数进行估计；同时给出了Fisher信息矩阵、渐近置信区间、贝叶斯可信区间与自助置信区间。蒙特卡洛模拟研究表明，均方误差与区间长度总体随样本量增大而减小，估计性能得以提升。将该方法应用于真实回火钢寿命数据，证实所提模型在各应力水平下均拟合良好，MLE与贝叶斯估计结果相近，且贝叶斯可信区间通常给出更精确的区间估计。

研究背景与动机

工业系统中人机可靠性直接影响质量与安全，但在正常工作条件下高可靠性产品失效数据难以获取，加速寿命试验（ALT）通过对样品施加高于额定水平的应力诱发早期失效，再将数据外推至正常使用条件。现有模型常难以适配复杂风险函数形态，且既有的应力-寿命关系（如Arrhenius模型、逆幂模型）可能无法涵盖全部制造异常。为此，研究人员在恒定应力ALT（CS-ALT）结合渐进II型删失（PT-II censoring）场景下，引入新Odd Lindley Half Logistic（New OLiHL）分布，以兼顾右偏、单峰及风险函数递增的寿命数据特征，构建一套完整的统计推断方案。

论文发表于《Scientific African》，其意义在于为可靠性工程提供一类可处理复杂删失结构与多应力水平的新建模工具，并系统比较频率学派与贝叶斯学派估计在有限样本下的表现。

主要关键技术方法

研究人员采用五类核心技术：其一，基于OLiHL族导出New OLiHL分布的密度、累积与风险函数，并通过分位函数结合Lambert-W函数生成随机数；其二，在k水平CS-ALT与PT-II删失下构造含加速因子α_hi的似然函数，应力-形状参数关系取ln λ_i=a+b φ_i；其三，用最大似然估计（MLE）求解非线性方程并借Fisher信息矩阵给渐近区间；其四，设λ₀、α、θ服从伽马先验，以平方误差损失函数运行MCMC（Gibbs抽样）得后验与可信区间；其五，实施Boot-P（参数自助）与Boot-T（枢轴量自助）重复抽样。模拟研究按既定应力φ_i与删失方案R_ij生成New OLiHL样本，实证部分采用公开回火钢寿命数据集。

研究结果

Statistical methodology

给出OLiHL基线密度f(x)=λ²/(4(1+λ))(1+e^x)e^{?λ/2(ex?1)+x}及New OLiHL经尺度变换t=xθ后的pdf、CDF与HRF，证明λ≈0.81近对称，分布可左右偏且峰态多变；借分位函数Q(u)=θ ln(1?2/λ[λ+1+W_?1{(1+λ)(u?1)e^?(1+λ)}])用Uniform(0,1)随机数造样。列表呈示均值、标准差、偏度、峰度与变异系数随λ、θ的变化规律：λ升则均值与标准差降、CV升；θ升则均值与标准差升、形状系数不变。

Test Assumptions

设s₀< />₁<…< />_k为有序应力，第i级投n_i件，预置删失方案R_i1,…,R_imi满足ΣR_ij+m_i=n_i；首失效抽R_i1件、次失效抽R_i2件直至第m_i次失效。假定各应力下失效时间服从OLiHL，寿命-应力联律取ln λ_i=a+b φ_i（φ_i为1/(?z)、ln z或z），加速因子α_hi=λ_i/λ₀=exp(b(φ_i?φ₀))，h_i=(φ_i?φ₀)/(φ₁?φ₀)且h_k>…>h₁=1。

Maximum Likelihood Estimation (MLE)

推导PT-II删失下联合似然L(λ₀,α,θ)=∏_i=1^kC_i∏_j=1^m_if_Ti(t_ij)[1?F_Ti(t_ij)]^R_ij，对数似然对λ₀、α、θ求偏导得三组非线性方程，无闭式解须数值解；由Fisher信息矩阵逆阵给渐近协方差与置信区间。

Bayesian Estimation and MCMC

取λ₀、α、θ独立伽马先验，平方误差损失下后验期望即贝叶斯估计；MCMC用Gibbs从全条件分布迭代，获后验均值与最高后验密度（HPD）可信区间，避免大样本近似偏差。

Bootstrap Methods

Boot-P从拟合New OLiHL重抽样本再估参数；Boot-T以(估?真)/标准误为枢轴量建区间，二者供小样本下与MLE、贝叶斯互校。

Simulation Study

按n_i、m_i、R_ij、φ_i多情景造数，重复千次：MSE与区间长随n增大而降，θ估计偏倚最小，λ₀与α受应力间距影响；贝叶斯区间覆盖率最接近名义水平，Boot-T略窄于Boot-P。

Real Data Application

回火钢寿命数据分三应力水平拟合New OLiHL，AIC、KS检验优于Weibull与广义半正态对照；MLE与MCMC点估计差<3%，贝叶斯可信区间宽度平均比MLE渐近区间窄12%，验证模型实用性与估计一致性。

讨论与结论翻译

讨论指出，New OLiHL在PT-II删失CS-ALT中能灵活刻画递增风险与右偏寿命，MLE、MCMC、Boot-P、Boot-T四类估计互补：大样本下MLE渐近有效，小样本下MCMC可信区间覆盖更稳，自助法提供非参数校验。局限是应力-寿命仍依赖预设φ(·)形式，若制造异常致非线性漂移则外推存偏。结论为：所建推断框架在模拟与回火钢实证中均表现良好，贝叶斯途径在区间精度上略优，该分布可作为可靠性寿命试验的一类候选基线模型。