data(s109)
id1 <- ideal(s109,
d=1,
meanzero=TRUE,
store.item=TRUE,
maxiter=1000,
burnin=100,
thin=10)
plot(id1) ## same as plot1d(id1) for d=1
plot(id1,showAllNames=TRUE)
## long run, too long for examples
idLong <- ideal(s109,
d=1,
meanzero=TRUE,
store.item=TRUE,
maxiter=251e3,
burnin=1000,
thin=1e3)
pdf(file="s109.pdf",
paper="special",
height=16,width=8.5) ## big graphing area
plot(idLong,showAllNames=TRUE)
dev.off()
Run the code above in your browser using DataLab