ei.MD.bayes.  Since we are 
drawing each parameter one at a time, target acceptance rates are between 0.4 to 0.6.tuneMD(formula, covariate = NULL, data, ntunes = 10,  totaldraws = 10000, ...)cbind(col1, col2, ...) ~ cbind(row1, row2, ...) with rows as the predictor and columns 
as the response~ xformula and
covariateei.MD.bayesei.MD.bayes