library( MethComp )
sessionInfo()

data( ox )
ox <- Meth( ox )
summary( ox )

plot( ox )

DA.reg( ox )

par( mar=c(3,3,1,3), mgp=c(3,1,0)/1.6 )
BA.plot( ox )
par( mar=c(3,3,1,1), mgp=c(3,1,0)/1.6 )
BA.plot( ox, pl.type="conv", repl.conn=TRUE )
BA.est( ox )

( AR.ox <- AltReg(ox,linked=TRUE,trace=TRUE) )

MCmcmc( ox, code.only=TRUE )
MC.ox <- MCmcmc( ox, n.iter=500 )

print( MC.ox )
plot( MC.ox )
trace.MCmcmc( MC.ox )
post.MCmcmc( MC.ox )