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

在VSNI官方论坛 link上，有人咨询能否运行多点试验多性状的因子模型？

问题是：多点试验多性状的因子模型，需要每个性状得到各自的因子模型结果。
答案是：可以运行。
疑问是：这种模型是否有意义？

简单示例如下：

asreml.options(ai.sing=T)
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

yzhlinscau

发布了15 篇原创文章 · 获赞 3 · 访问量 1803

私信关注

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

猜你喜欢