if (FALSE) {
data(s109)
id1 <- ideal(s109,
d=1,
normalize=TRUE,
store.item=TRUE,
maxiter=500, ## short run for examples
burnin=100,
thin=10)
plot(id1)
id2 <- ideal(s109,
d=2,
store.item=TRUE,
maxiter=11e2,
burnin=1e2,
verbose=TRUE,
thin=25)
plot(id2,overlayCuttingPlanes=TRUE)
id2pp <- postProcess(id2,
constraints=list(BOXER=c(-1,0),
INHOFE=c(1,0),
CHAFEE=c(0,.25)))
plot(id2pp,overlayCuttingPlanes=TRUE)
}
Run the code above in your browser using DataLab