data(cardiac)
cardiac <- Meth(cardiac)
summary(cardiac)
# Visually check exchangeability
plot( cardiac )
plot( perm.repl( cardiac ) )
BA.est(cardiac)
# Run MCmcmc using BRugs for an insufficient amount of iterations
if (FALSE) card.mi.ir <- MCmcmc( cardiac,
beta=FALSE, random=c("mi","ir"),
n.iter=100, trace=T )
print( card.mi.ir )
Run the code above in your browser using DataLab