mcmc RJAGS에서 데이터를 추출하는 방법

Question

각 폴링 주에 대해 측정 한 다음 모델을 실행 중입니다. 각 파티마다 100 개 이상의 알파 벡터를 생성합니다. 질문 할 수 있습니다. 선 그래프를 그릴 수있는 방법으로 어떻게 비축 할 수 있습니까?

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

Answer 1

들쭉날쭉 기능은 객체의 클래스에 대한 summary 일반적인 기능으로 mmcm.list 객체를 제공하기. str 및 name 함수를 새 객체에 사용하면 jags 함수가 반환되고 mcmc과 같은 것을 찾아야합니다. mcmc 하위 오브젝트는 mcmc.list의 매트릭스이며, 예를 들어 n row 이터 레이션 수와 n col 추정 매개 변수를 포함합니다. 그래서 당신은 당신이 원하는 것을 추출 할 수 있습니다. mcmc 개체의 길이는 체인 수와 같습니다. 그들은 또한 후방 밀도를 그리기위한 일반 plot.mcmmc 기능을 가지고 있습니다.