mcmc.list에서 plot.bugs 및 plot.jags로 생성 한 것과 비슷한 플롯을 어떻게 생성 할 수 있습니까? [닫은]


9

R은 R2WinBUGS :: bugsR2jags : jags 함수에 의해 생성 된 bugsjags객체 에서 멋진 요약 플롯을 출력 할 수있는 것으로 보입니다 .

그러나 rjags패키지를 사용하고 있습니다. 나는 함수의 결과 플롯 할 때 rjags::coda.samples사용 R2WinBUGS::plot.mcmc.list결과를 각 매개 변수에 대한 진단 플롯 (매개 변수 밀도, 체인 시계열 자기 상관)입니다.

아래는 Andrew Gelman의 자습서 "Running WinBuugs and OpenBugs from R"에서 생성하고 싶은 플롯 유형입니다 . 이들은를 사용하여 제작되었습니다 plot.pugs.

문제는 즉 plot.bugsbugs하면서, 인수로서 오브젝트 plot.mcmc.list의 출력을 얻어 coda.samples.

다음은 예제입니다 coda.samples.

 library(rjags)
 data(LINE)
 LINE$recompile()
 LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
 plot(LINE.out)

내가 필요한 것은

  • 에 의해 생성 된 것과 유사한 정보가 풍부한 유사한 한 페이지 요약 도표를 생성하는 방법 plot.bugs
  • LINE.outbugs 객체 로 변환하는 함수 또는

여기에 이미지 설명을 입력하십시오

답변:


2

답이 없기 때문에 적어도 지금까지 얻은 것을 게시 할 것입니다.

패키지 의 as.bugs.array기능은 R2WinBUGS이 목적으로 만들어졌습니다. 설명서 ( ?as.bugs.array) 에 따르면 :

버그가 아닌 Markov chain 시뮬레이션 결과를 버그 객체로 변환하는 함수. 주로 plot.bugs로 결과를 표시하는 데 사용됩니다.

따라서 LINE.out올바른 변수를 플롯하지는 않지만 예제에서 플롯을 얻을 수 있습니다 .

plot(as.bugs.array(sims.array = as.array(LINE.out)))        

를 변환하는 올바른 방법을 결정하려면 약간의 작업이 더 필요 LINE.out하며, LINE.samples객체 example(jags.samples)를 시작하기가 더 쉬울 수 있습니다.


0

다음은 저에게 효과적입니다.

require(R2jags)
m <-jags(data=d,inits=i,pars,n.iter=1000,n.chains=3,model.file="foo.txt",DIC=F)
m <- autojags(m)
plot(m)

다음은 재현 가능한 예입니다.

example(jags)
plot(jagsfit)

1
그것은 유용한 단서이지만 mcmc.list(내가 말할 수있는 한)로 시작하는 문제를 해결하지 못합니다 .
David LeBauer
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.