只有map文件和ped文件怎么通过plink进行确实基因筛选
时间: 2024-10-10 08:06:04 浏览: 39
libplinkio:一个小的C和Python库,用于读取PLINK基因型文件
在使用PLINK进行基因筛选时,如果你已经有map(遗传图)文件和ped(个体信息)文件,你可以执行以下操作:
1. **导入和预处理**[^1]:
```bash
plink --file output # 假设你的ped和map文件名为output.ped和output.map
```
这会加载ped和map文件,对数据进行初步的预处理,如删除缺失值和非单核苷酸多态性(SNP)。
2. **选择特定的标记(variants)**:
如果你知道要筛选的特定基因位置(例如基于基因ID或染色体/位置),可以使用`--extract`选项指定这些位置:
```bash
plink --file output --extract "chr1:10000-20000" # 选择chr1上的10000-20000bp区域内的SNPs
```
3. **保存筛选后的数据**:
完成筛选后,你可以使用`--make-bed`命令创建一个新的bed文件,其中只包含你感兴趣的SNPs:
```bash
plink --make-bed --biallelic-only --out selected_snps # 新的bed文件命名为selected_snps
```
现在,`selected_snps.bed`文件就包含了你所筛选的确切基因SNPs。
阅读全文