[英]How can I extract data from mcmc RJAGS
我正在運行以下模型,其中每個投票周都有一個度量值。 每一方產生100多個alpha向量。 可能的問題是,我該如何以繪制折線圖的方式來存放它們?
model{
## measurement
for(i in 1:NPOLLS){
p1[i] ~ dnorm(alpha1[WoY[i]] + pollster1[org[i]], prec1[i])
p2[i] ~ dnorm(alpha2[WoY[i]] + pollster2[org[i]], prec2[i])
p3[i] ~ dnorm(alpha3[WoY[i]] + pollster3[org[i]], prec3[i])
}
## transition
for(i in 2:NPERIODS){
alpha1[i] ~ dnorm(alpha1[i-1],phi1)
alpha2[i] ~ dnorm(alpha2[i-1],phi2)
alpha3[i] ~ dnorm(alpha3[i-1],phi3)
}
pollster1[1] <- -sum(pollster1[2:NPOLLSTERS])
pollster2[1] <- -sum(pollster2[2:NPOLLSTERS])
pollster3[1] <- -sum(pollster3[2:NPOLLSTERS])
## priors
phi1 ~ dgamma(5000,1)
alpha1[1] ~ dunif(.2,.4)
phi2 ~ dgamma(5000,1)
alpha2[1] ~ dunif(.2,.4)
phi3 ~ dgamma(10000,1)
alpha3[1] ~ dunif(.1,.3)
for(i in 2:NPOLLSTERS){
pollster1[i] ~ dnorm(0,.01);
pollster2[i] ~ dnorm(0,.01);
pollster3[i] ~ dnorm(0,.01);
} }
jags函數為mmcm.list對象提供此類對象的summary
通用函數。 在jags
函數返回的新對象上使用str
和name
函數,您需要找到類似mcmc
東西。 mcmc
子對象是mcmc.list
類的矩陣, mcmc.list
包含mcmc.list
其中n row
表示迭代次數, n col
表示估計參數。 這樣您就可以提取所需的內容。 mcmc
對象的長度等於鏈數。 它們還具有通用的plot.mcmmc
來繪制后密度。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.