0

我正在运行以下模型,在该模型中我对每个投票周进行了测量。它为每一方产生超过 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);
 } }
4

1 回答 1

0

jags 函数为 mmcm.list 对象提供此类对象的summary通用函数。新对象上的使用str和函数由函数返回,您需要找到类似. 子对象是一个包含具有迭代次数和估计参数的后验的类矩阵。所以你可以提取你想要的。对象的长度等于链的数量。它们还具有绘制后验密度的通用函数。namejagsmcmcmcmcmcmc.listn rown colmcmcplot.mcmmc

于 2014-05-26T08:03:26.213 回答