GWAS分析-说人话(16)- "一键提取“所有染色体需要的表型
2020-08-14 本文已影响0人
医学小蛋散
前言
这个其实是“GWAS分析-说人话(7)keep提取我们需要的表型数据 III”的续编IV
没错,只要你的需求(daoshi)不是什么变态,
一般都是要所有染色体都提取数据的。
留意哦,简短粗暴:
Terminal中发动for循环!~
for ch in `seq 22`; do /Users/xxx/plink --noweb --file /Volumes/xxx/chr$ch --keep /xxx/你想要提取的文件.txt --recode --make-bed --out 你想要的文件名$ch; done
是的,以上~
进一步了解--keep功能,请到:
1.GWAS分析-说人话(5)keep提取我们需要的表型数据 I