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)
bivariate model using bayesian by Luis A. Apiolaza @zentree 14h
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- Import Most Recent(导入最新)按钮让应用程序忽略来源导入最新的传递函数。左下角的下拉菜单列表用于...
- AIDAPro可用于鉴定转换域(线性)的模型。AIDAPro拥有一个可执行各种变换的计算引擎。然而,AIDAPro...
- 这些元素可被分为三类:1.变量模块,包括:•三种类型的输入模块:MV(操作变量),DV(干扰变量)和UNM(不可测...
- 2017/10/1 荷香第九堂,灵性的音乐像阳光一样照耀你 要义:好的音乐,会让你不自觉放松,心里的门会一扇一扇的...