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
## Not run: card.mi.ir <- MCmcmc( cardiac,
# beta=FALSE, random=c("mi","ir"),
# n.iter=100, trace=T )
# print( card.mi.ir )## End(Not run)
Run the code above in your browser using DataLab