require(ltable)
data(tdata, package="ltable")
## For better illustration You should increase draw and burnin pars
pres<-MCPower(Counts~smoker +contraceptive +tromb +
contraceptive*tromb, scale_max=1.5, effect="contraceptive*tromb",
data=tdata, draw=1000, burnin=300)
plot(pres)
plot(pres, stencil=3)
Run the code above in your browser using DataLab