plink 软件中依据LD信息对数据进行去冗余

 

001、

plink --file file --sheep --allow-extra-chr --indep-pairwise 50 5 0.2 --set-missing-var-ids @:# --keep-allele-order  --out tmp_result       ## 依据ld过滤

plink --file tmp_result --sheep --allow-extra-chr --extract tmp_result.prune.in --keep-allele-order --set-missing-var-ids @:# --recode tab --out result_LD_pruned     ## 提取过滤后的位点

 

 

reference:

01、https://zhuanlan.zhihu.com/p/9441695073

02、Lv FH, Cao YH, Liu GJ, et al. Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity, Introgression, and Agronomically Important Loci. Mol Biol Evol. 2022 Feb 3;39(2)

posted @   小鲨鱼2018  阅读(10)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律
历史上的今天:
2024-02-19 linux系统中的tr -dc选项
2024-02-19 linux 中生成随机数
2024-02-19 linux 中利用awk进行数值计算
2024-02-19 linux 中 grep匹配多个可能的字符
2024-02-19 linux中grep匹配连续出现指定次数的字符组合
2024-02-19 linux 中 grep命令 $?返回值 0、1、2的意义
2024-02-19 while循环与until循环的区别
点击右上角即可分享
微信分享提示