1. 引言
逆伽马分布在统计学、物理学、药理学等领域都有广泛的运用,例如在雷达工程中逆伽马分布作为统计特性严重拖尾的杂波建模的复合高斯过程的纹理分量 [1] ;在药理学中灭绝条件下随机抗生素耐药性模型中的患者携带对药物敏感的细菌弱收敛于逆伽马分布 [2] ;在误差分析和随机控制中也会经常运用该分布。所以对逆伽马分布进行深入的研究,具有重要的学术意义。
统计决策中统计量的优劣由损失函数的选择和参数估计问题的精确程度决定。目前,研究对称损失下一般分布参数的Bayes估计的文献有很多,如王亚楠等 [3] 讨论了复合LINEX对称损失下逆高斯分布参数倒数Bayes估计;岑泰林等 [4] 讨论了复合LINEX对称损失下广义Pareto分布形状参数
的Bayes估计;吕佳等 [5] 讨论了复合LINEX损失下艾拉姆咖分布参数的贝叶斯估计;李俊华 [6] 讨论了复合LINEX对称损失函数下Burr分布参数的估计;但逆伽马分布有关的Bayes估计的文章很少,如丁新月等 [7] 讨论了Mlinex损失函数下逆伽马分布尺度参数的Bayes估计。
本文在复合Linex对称损失函数下对双参数逆伽马分布在形状参数已知的情况下,研究了其尺度参数的估计。
设X为随机变量。假如它的密度函数为
(1)
称X服从参数为
的逆伽马分布,记作
,其中
为形状参数,
为尺度参数。
设
为来自总体X独立同分布的样本,则此样本的似然方程为
(2)
其中
。
复合Linex对称损失函数由张睿 [8] 提出,其形式为
(3)
2. θ的Bayes估计
引理1 对任何的先验分布
,在复合Linex对称损失函数下,
的Bayes估计为
(4)
定理1 设
是
的一组观察值,则在损失函数(3)下,取
为
的先验分布,
的Beyes估计为
其中
。
证明 由
为逆伽马分布
的先验分布,可得
的先验密度函数为
(5)
根据样本的联合密度函数(2)和式(5),得到
的后验密度函数为
(6)
其中
,从而
的后验分布为
。那么
(7)
(8)
将(7)和(8)带入式(4)可得
3. θ的E-Bayes估计
对
中的超参数a、b进行估计。设超参数a、b的先验分布为
和
,若
的联合分布为均匀分布,即
,此时
的多层先验密度函数为
(9)
定理2 在复合Linex对称损失函数下,对于逆伽马分布(1),若参数
的先验分布为
,超参数a和b的先验分布为D上的均匀分布,则参数
的E-Bayes估计为
(10)
其中
。
证明 根据E-Bayes的定义得
4. θ的多层Bayes估计
定理3 在复合Linex对称损失函数下,对于逆伽马分布(1),若参数的先验分布为
,则参数
的多层Bayes估计为
(11)
其中
。
证明 根据Bayes定理,
的多层后验密度为
从而有
故在复合Linex对称损失函数下
的多层Bayes估计为
5. 数值模拟
利用R语言进行数值模拟,生成一组
,形状参数
、尺度参数
真值为1的逆伽马分布随机样本,并根据定理1中参数
的Bayes估计
,在参数
时,对参数a和c取不同的值时
的估计值,模拟结果如表1。
利用相同的随机样本,根据定理2中参数
的E-Bayes估计
和定理3中参数
的多层Bayes估计的表达式
在参数
时求估计值,结果分别为表2与表3。
稳健性:
由表1,表2,表3可以看出当
且适当选择参数
以及d的值时,
,
,
的极差都非常小。从统计决策中稳健性角度考虑,参数
的三种Bayes估计都很稳健,其中Bayes估计最稳健,符合统计决策中稳健性。
精确度:
由表1,表2,表3可以看出当
时,容易求出
,
,
的偏差(
,其中
为参数
的估计量,
为参数
的真值),区间分别为[0.000303, 0.011622],[0.018552, 0.025928],[0.005407, 0.12818]。可见偏差非常小,所以其精确度也很高,其中Bayes估计的精确度最高,符合统计决策的要求。
基金项目
国家自然科学基金(11801488);新疆师范大学教学研究与改革(SDJG2020-30);新疆师范大学科研发展专项(XJNUZX202001)。
NOTES
*通讯作者。