gamma分布_求倒Gamma分布的可信区间迭代实例

在Bayes统计中,可信区间的求解,往往囿于所用分布的不对称,多数需要通过迭代求得,比如倒Gamma的可信区间的求解就需要这种迭代,下面是一个迭代求解实例。

function Ch3P72Eg326(b,tMid,p) %名称释义:Ch3P72Eg326.m. %实现功能:实施例题3.2.6的具体迭代计算，求倒Gamma分布的可信区间HPD. %输入参数:函数含3输入参数:含义如下： %　　　　:(1),b,IGA的尺度参数. %　　　　:(2),tMid,参数theta的初估值，一般使用后验中位数tMid设置. %　　　　:(3),p,小概率临界值,缺省默认0.1. %输出参数:函数默认无输出参数. %函数调用:实现函数功能不需要调用子函数,但含内部匿名函数两个: %　　　　:(1),bpi,用来处理后验概率分布 %　　　　:(2),Ft,用来处理后验分布函数 %参考文献:实现函数算法,参阅了以下文献资料: %　　　　:(1),马寨璞,MATLAB语言编程[M],北京:电子工业出版社,2017. %原始作者:马寨璞,wdwsjlxx@163.com. %创建时间:2020.11.07,09:36:24. %原始版本:1.0 %验证说明:本函数在MATLAB2020a等版本运行通过. %使用样例:直接运行:Ch3P72Eg326 % if nargin<1               % 缺省设置     b=42868;             % 尺度值     tMid=14289;         % 众数 end if nargin<3||isempty(p)  %设置缺省的临界概率(显著性水平)     p=0.1; end %两个匿名函数，分别计算概率密度值和累积概率值