博文
用ASReml运行多点试验多性状的因子模型
||
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yzhlinscau/article/details/102239222
在VSNI官方论坛 link上，有人咨询能否运行多点试验多性状的因子模型？
问题是：多点试验多性状的因子模型，需要每个性状得到各自的因子模型结果。
答案是：可以运行。
疑问是：这种模型是否有意义？
简单示例如下：
run9<-asreml(cbind(yield,hh1)~trait+trait:Loc, random=~at(trait):us(trait):fa(Loc,1):Genotype, residual=~units:us(trait), maxit=10,data=MET)
运行结果如下：
Var(run9) at(trait, yield):trait:fa(Loc, 1):Genotype!trait_yield:yield 26.511314673 at(trait, yield):trait:fa(Loc, 1):Genotype!trait_hh1:yield 4.139143529 at(trait, yield):trait:fa(Loc, 1):Genotype!trait_hh1:hh1 6.208715293 at(trait, yield):trait:fa(Loc, 1):Genotype!1!var 0.000000000 at(trait, yield):trait:fa(Loc, 1):Genotype!2!var 0.087426566 at(trait, yield):trait:fa(Loc, 1):Genotype!3!var 0.090362905 at(trait, yield):trait:fa(Loc, 1):Genotype!4!var 0.261599691 at(trait, yield):trait:fa(Loc, 1):Genotype!5!var 0.000001600 at(trait, yield):trait:fa(Loc, 1):Genotype!6!var 0.000001600 at(trait, yield):trait:fa(Loc, 1):Genotype!1!fa1 0.407585698 at(trait, yield):trait:fa(Loc, 1):Genotype!2!fa1 0.394941929 at(trait, yield):trait:fa(Loc, 1):Genotype!3!fa1 0.255072194 at(trait, yield):trait:fa(Loc, 1):Genotype!4!fa1 0.375419515 at(trait, yield):trait:fa(Loc, 1):Genotype!5!fa1 0.301101708 at(trait, yield):trait:fa(Loc, 1):Genotype!6!fa1 0.104013083 at(trait, hh1):trait:fa(Loc, 1):Genotype!trait_yield:yield 6.208715293 at(trait, hh1):trait:fa(Loc, 1):Genotype!trait_hh1:yield 4.139143529 at(trait, hh1):trait:fa(Loc, 1):Genotype!trait_hh1:hh1 0.948916905 at(trait, hh1):trait:fa(Loc, 1):Genotype!1!var 7.445363318 at(trait, hh1):trait:fa(Loc, 1):Genotype!2!var 0.002515794 at(trait, hh1):trait:fa(Loc, 1):Genotype!3!var 0.000000000 at(trait, hh1):trait:fa(Loc, 1):Genotype!4!var 0.000000000 at(trait, hh1):trait:fa(Loc, 1):Genotype!5!var 0.000001600 at(trait, hh1):trait:fa(Loc, 1):Genotype!6!var 2.976343939 at(trait, hh1):trait:fa(Loc, 1):Genotype!1!fa1 -0.793744164 at(trait, hh1):trait:fa(Loc, 1):Genotype!2!fa1 3.816775300 at(trait, hh1):trait:fa(Loc, 1):Genotype!3!fa1 1.643122019 at(trait, hh1):trait:fa(Loc, 1):Genotype!4!fa1 -1.157603387 at(trait, hh1):trait:fa(Loc, 1):Genotype!5!fa1 0.002461397 at(trait, hh1):trait:fa(Loc, 1):Genotype!6!fa1 0.240477558
https://blog.sciencenet.cn/blog-1114360-1200915.html
上一篇：《动植物育种遗传数据分析》网购途径--更新
下一篇：[重磅推荐]免费版的动植物遗传评估软件Echidna
