||
在ASReml-R V4的referrence manual第60页中,终于列出了千呼万唤的计算标准误函数vpredict(),据该函数介绍其可用于遗传力、遗传相关等指标的标准误计算。然我在运行自己的例子时,发现计算结果有误。
首先按手册运行其例子,发现运行结果也与手册的不同,不知是何故?应该是手册里的结果有误。
实际运行结果如下:
虽然运行结果与手册里不同,但看起来vpredict函数似乎可以使用。但是运行我自己的一份数据,却发现结果貌似不行。
结果如下:
不论是亲本模型还是个体模型,遗传力的标准误在100左右,这个很夸张!用我自己编码的pin3()函数,得到的标准误是ok的。
不知道是我使用不当,还是作者验证vpredit函数时所用数据集偏少。我也期待其他用户是否会有类似问题。当然目前我所用的ASReml-R V4仍是beta版,如存在一些问题也不奇怪。
V4,在很多方面已经与V3完全不同了,这对于习惯V3的用户来说,不是件好事,至少对我而言,挺痛苦的。手册过于简单,所以要使用V4的新功能,也非易事,我倒希望其他用户能出本V4手册,让我也偷偷懒。^_^
ps: 因VSNi官方在处理ASReml-R V4的license方面存在问题,所以连beta版都下线了,具体何时上线,时间未知。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2025-1-10 17:28
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社