Question

我正在尝试在JAGS中运行分层Dirichlet模型，但是我没有更新，因此必须做错什么。我尝试使用伽玛分布对其进行近似：

#Creating some data
set.seed(555)
cat1=rbeta(15,20,60)
cat2=rbeta(15,20,80)
cat3=rbeta(15,20,160)
cat4=1-cat1-cat2-cat3

dat_dirich=list(
  dirF=cbind(cat1, cat2, cat3, cat4),
  AL=c(1,1,1,1),
  K=4
)

require(runjags)

{
  catch_mod="model{

    for(y in 1:15){
      dirF[y,1:K]~ddirch(alpha_dirF1B[y,1:K])

      #Approximation with Gamma
      for(a in 1:K){
          alpha_dirF1B[y,a]<-P2[y,a]/sum(P2[y,1:K])
          P2[y,a]~dgamma(alpha_dirF1[a,y],1)#1 or kappa_dirF
      }
      #hierarchical structure   
      alpha_dirF1[1:K,y]~ddirch(AL*kappa_dirF)
    }
    kappa_dirF~dunif(0.1,5000)   #kappa_dirF~dlnorm(0,0.01)   
  }"
}

Answer 1

您必须在代码中的某个地方调用run.jags()。看起来应该像

results <- run.jags(catch_mod, data = dat_dirich, n.chains = XXX, ...)

这在runjags（https://cran.r-project.org/web/packages/runjags/vignettes/quickjags.html）的插图中有描述

Dirichlet参数不会在JAGS中更新

1 个答案: