PopSizeABC--评估历史有效群体大小Ne

文件 - PopSizeABC - Forge DGA (inra.fr)

1、打开上述链接下载cattle_data.tar.gzcomp_stat1.zipestim.zip,这三个文件,并解压到到PopSizeABC的目录下,并在目录下建一个res目录,方便结果文件的保存

2.用python运行stat_from_vcf_ex1.py,生成一个stat后缀文件,保存在res里。pop为物种品种,n为物种单倍体数量,比如10头牛,牛是二倍体生物,这里n设为20,mac看情况,一般设置为0.2*n,macld同。

3.用python运行simul_data_ex1.py,生成两个s1后缀的文件,保存在res里,outfile设置为ex1_物种名,其他同上设值。

4.用R运行abc_cv_ex1.R,需要安装abc这个R包,输入文件为上一步的s1后缀文件,同上设置参数

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容