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)
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容