我如何從mcmc中提取數據RJAGS

我正在運行以下模型，其中我爲每個投票周進行了測量。每個派對產生超過100個阿爾法矢量。可能的問題是，我怎樣才能以一種我可以繪製折線圖的方式儲存它們？我如何從mcmc中提取數據RJAGS

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); 
} }

來源

2014-03-31 user45367

對mcmc對象的總結？ –

V形凹口功能provid一個mmcm.list對象與這個類對象的summary通用函數。使用str和name功能在您的新對象返回jags功能，你需要找到像mcmc。 mcmc子對象是包含n row迭代次數和n col估計參數的類別mcmc.list的矩陣。所以你可以提取你想要的東西。 mcmc對象的長度等於鏈的數量。他們還有一個通用的功能來繪製後部密度。

來源

2014-05-26 08:03:26 Philippe

我如何從mcmc中提取數據RJAGS

回答

相關問題