大家好,我是邓飞。
最近有小伙伴问起如何检查系谱数据是否正常,这里,分享一个我之前写的一个R包,里面有函数check_pedigree,可以说是非常方便了。
命令行安装:
library(devtools)
install_github("dengfei2013/learnasreml")
本地安装包:(公众号回复:learnasreml)
功能介绍:
模拟一个数据:
运行程序检测代码:
library(learnasreml)
check_pedigree(ped)
结果:
可以看到,系谱是有问题的,父母本交叉个数1个,编号为2。
立一个Falg吧,这个包好久没有维护了,而自己又写了很多实用的函数,都可以打包进去,后面更新吧。
PS:话说函数输出中文的汇总结果是不是很666……