data(s109)
id1 <- ideal(s109,
d=1,
meanzero=TRUE,
maxiter=1000, ## short run for demo purposes
burnin=100,
thin=10)
id1coda <- idealToMCMC(id1)
summary(id1coda)
Run the code above in your browser using DataLab