新分布的吉布斯抽样

如何解决新分布的吉布斯抽样

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

$\frac{1}{Z} \frac{1}{\sqrt{2\pi}e^{-0.5x^2}}cosh(\beta x)$

其中 $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)

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