新分布的吉布斯抽样

编程问答 2022-05-27

问题描述

我想运行以下密度的 Gibbs 采样：

其中 $beta$ 是一个正参数，它取小值，Z 是归一化常数。 为此，我首先使用以下代码生成了一个新的绝对连续分布

f <- Abscontdistribution(d = function(x) 1/(sqrt(2*pi))*exp(-0.5*x^2)*cosh(beta*x),withStand = TRUE)

然后我使用库 jags 来运行 Gibbs Sampling：

require(rjags)
library(coda)
writeLines("model {
           y ~ f(x)
           x ~ dnorm (0,1)
           }",con="jags01.txt")
jagsModel <- jags.model('jags01.txt',data=list ('y'=1))
postSamples <- coda.samples(jagsModel,'x',n.iter = 100)
plot(postSamples,density=TRUE,main="Prova",auto.layout=FALSE)

不幸的是它不起作用。有人可以解释我为什么吗？ 提前致谢。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

jagsmcmcrr