如何在 plink 中按家庭 ID 拆分我的 bed bim 和 fam 文件？答案

【问题标题】：How do I split my bed bim and fam files by family ID in plink?如何在 plink 中按家庭 ID 拆分我的 bed bim 和 fam 文件？
【发布时间】：2015-04-07 05:17:24
【问题描述】：

我想在 plink 中按家庭 ID 拆分我的 bed、bim 和 fam 文件。例如，我可以通过输入染色体来执行此操作：

p-link --bfile datafile1 --chr 1 --make-bed --out datafile2

这就为 1 号染色体创建了一个单独的 bed、bim 和 fam 文件。

但是，在我的 fam 文件中，我有几个不同的家庭 ID，我想以此来分割文件。例如

Family 1 TS11339 0 0 2 -9
Family 1 TS11233 0 0 2 -9
Family 2 TF99288 0 0 2 -9
Family 2 YJ22997 0 0 2 -9
Family 3 YF00277 0 0 2 -9
Family 3 YY92779 0 0 2 -9

我想要 Family1.bed -.bim 和 -.bam，Family2.bed -.bim 和 -.bam 等。

我该怎么做？（抱歉知识贫乏 - 刚开始使用 plink）。

【问题讨论】：

标签： split plink

【解决方案1】：

以下命令应在 UNIX shell 上拆分 plink 数据集：

for fam in $(awk '{print $1}' datafile1.fam | sort | uniq); do echo $fam | plink --bfile datafile1 --keep-fam /dev/stdin --make-bed --out $fam; done

【讨论】：

如果您想将文件拆分为每个只有一个系列的新文件，这非常有效，特别是如果您有很多系列并且不想手动提取每个系列。非常感谢！

【解决方案2】：

简单一点：plink --keep-fam。

【讨论】：

【解决方案3】：

我想通了，你可以通过plink如下：

p-link --bfile datafile1 --keep LIST1.txt --make-bed --out

其中 LIST1.txt 是家庭 ID 和个人的列表。所以说如果你想提取家庭 1，LIST1.txt 将包括：

家庭 1 TS11339

家庭 1 TS11233

更多信息可以在 plink 网站上找到。

【讨论】：