library(MethComp) data(ox) ox <- Meth(ox) system.time( MCox <- MCmcmc( ox, n.iter=10000 ) ) system.time( Jox <- MCmcmc( ox, n.iter=10000, program="jags" ) ) system.time( tJox <- MCmcmc( ox, n.iter=10000, program="jags", Transform="pctlogit" ) ) save( MCox, Jox, tJox, file="oxex.Rdata" ) MCox Jox MethComp(MCox)$VarComp MethComp( Jox)$VarComp MethComp(tJox)$VarComp # load( file="oxex.Rdata" ) pdf( "JAGS-conv.pdf" ) par( mar=c(3,3,1,1), mgp=c(3,1,0)/1.6 ) plot( MethComp( Jox), axlim=c(0,100), points=TRUE, repl.conn=TRUE, col.points="blue", grid=0:10*10, xaxs="i", yaxs="i" ) dev.off() pdf( "JAGS-BA.pdf" ) par( mar=c(3,3,1,1), mgp=c(3,1,0)/1.6 ) plot( MethComp( Jox), axlim=c(0,100), pl.type="BA", points=TRUE, repl.conn=TRUE, col.points="blue", grid=0:10*10, xaxs="i", yaxs="i" ) dev.off() pdf( "JAGS-tconv.pdf" ) par( mar=c(3,3,1,1), mgp=c(3,1,0)/1.6 ) plot( MethComp(tJox), axlim=c(0,100), points=TRUE, repl.conn=TRUE, col.points="blue", grid=0:10*10, xaxs="i", yaxs="i" ) dev.off() pdf( "JAGS-tBA.pdf" ) par( mar=c(3,3,1,1), mgp=c(3,1,0)/1.6 ) plot( MethComp(tJox), axlim=c(0,100), pl.type="BA", points=TRUE, repl.conn=TRUE, col.points="blue", grid=0:10*10, xaxs="i", yaxs="i" ) dev.off()