bivariate model using bayesian by Luis A. Apiolaza‏ @zentree 14h

bi <- MCMCglmm(cbind(ht, su) ~ trait - 1,
random = ~ idh(trait):plot + us(trait):genotype,
rcov = ~ us(trait):units,
data = trial,
ginverse = list(genotype = Ainv),
prior = bivar,
family = c('gaussian', 'threshold'), burnin = 10000,
nitt = 110000, 
thin = 100)