yzhlinscau的个人博客分享 http://blog.sciencenet.cn/u/yzhlinscau

博文

用ASReml运行多点试验多性状的因子模型

已有 1919 次阅读 2019-10-7 09:32 |个人分类: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
收藏 IP: 202.105.47.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 02:22

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部