如何解决来自关节后方或边缘后方的RJAGS输出?
我已经使用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 举报,一经查实,本站将立刻删除。