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,nit