全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 经管百科 爱问频道
1759 1
2016-03-23
悬赏 5 个论坛币 未解决
我现在做的就是一个广义非线性混合效应模型的ARMS抽样,全似然可以表示出来,比较复杂,现在想调用HI包里的arms函数:即arms(y.start, myldens, indFunc, n.sample, ...)来抽样,这个参数要怎么设置呢?帮助文件里给了例子,myldens表示目标密度函数,indFunc表示指示函数,具体怎么写呢?看的不太懂。是需要把参数的全条件后验分布的形式表示出来吗?带入那个参数呢?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2016-4-23 15:41:39
我现在对这个方法的理解是: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]
剩下的是其他的参数。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群