我现在对这个方法的理解是:HI包里的arms函数:即arms(y.start, myldens, indFunc, n.sample, ...)中:myldens表示目标密度函数,一般是我们要抽样的参数的对数满条件分布,是自定义函数;indFunc表示指示函数,给的应该是待估计参数的取值范围。
新的疑问是:我在估计服从先验分布为wishart分布的协方差矩阵时老是估计不粗来值,老说我没有提供 indFunc函数:
这里我来估计2*2协方差矩阵的var1,var2和相关系数rou:
para_cov.new=arms(para_cov.kesi,myldens,function(x...){(x[1]>0)*(x[2]<=1)*(x[2]>=-1)*(x[3]>0)},n.sample=Iter,param)
para_cov.kesi:为待估计参数的初始值
myldens为我自己定义的log condition full distribution
取值范围:协方差矩阵的var1,var2要求大于等于0 ,而相关系数rou取值在[-1,1]
剩下的是其他的参数。