# Variance terms
coefplot2(mpl1,ptype="vcov",intercept=TRUE,main="Random effect variance")

# Fixed effects
coefplot2(mpl1,intercept=TRUE,main="Fixed effect coefficient")