来自关节后方或边缘后方的RJAGS输出？

编程问答 2022-06-07

问题描述

我已经使用R中的RJAGS包创建了MCMC。RJAGS要求您在文本字符串中定义模型。定义我模型的先验分布和似然函数的文本字符串如下：

modelString="# BEGIN

model{

for(l in 1:25){
  for(j in 1:25){
      cov.matrix[l,j]<- theta^abs(l-j)
    }
  }

#------------- Likelihood
  for(i in 1:tau){
  K[i]~ dbin(p[i],N[i]) 
  p[i]<- pnorm(((qnorm(lambda,1) - sqrt(rho) *ss[i])/sqrt(1-rho)),1
               )
  }
  
#------------- Priors

  lambda ~ dunif(0,0.25)
  rho ~ dbeta(12.4,50.4)
  theta ~ dunif(0,1)
  ss ~ dmnorm(rep(0,tau),cov.matrix)  
}

" # END

嵌入在其余代码中的此模型字符串针对lambda，rho，theta和s1，...，s25生成观察值列表。我的问题是，这些观察样本列表是否来自λ，rho，θ，S1-25的联合后验分布？还是从每个参数的边际后验分布中取样？

解决方法

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

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

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

bayesianjagsmcmcrrrjags