plink --extract采坑问题记录

问题描述：使用–extract 参数后没能正常进行LD修剪

~/software/miniconda2/bin/plink --cow --allow-extra-chr -noweb --file /home/yanglv/soft_project/Target/target_result_test/result_5_3_clean_dropna2_filter_19sample/admixture_pairwise_50_10_0_2/test/admixture_pairwise_50_10_0_2 --geno 0.1 --make-bed --indep-pairwise 50 10 0.1# 322 variants removed due to missing genotype data (--geno).
# 70918 variants and 581 cattle pass filters and QC.
# Pruning complete.  47613 of 70918 variants removed.
# Marker lists written to
# /home/yanglv/soft_project/Target/target_result_test/result_5_3_clean_dropna2_filter_19sample/admixture_pairwise_50_10_0_2/test1_admixture.prune.in
# and
# /home/yanglv/soft_project/Target/target_result_test/result_5_3_clean_dropna2_filter_19sample/admixture_pairwise_50_10_0_2/test1_admixture.prune.out~/software/miniconda2/bin/plink --cow --allow-extra-chr -noweb --bfile plink --extract plink.prune.in --make-bed --out pruned_test1_admixtureNote: --noweb has no effect since no web check is implemented yet.
514998 MB RAM detected; reserving 257499 MB for main workspace.
Allocated 193124 MB successfully, after larger attempt(s) failed.
70918 variants loaded from .bim file.
581 cattle (0 males, 0 females, 581 ambiguous) loaded from .fam.
Ambiguous sex IDs written to pruned_test1_admixture.nosex .
--extract: 23836 variants remaining.
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 581 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.967899.
23836 variants and 581 cattle pass filters and QC.
Note: No phenotypes present.
--make-bed to pruned_test1_admixture.bed + pruned_test1_admixture.bim +
pruned_test1_admixture.fam ... done.

出现警告，大事不妙。
果然，惊奇的发现提取后的文件居然和原文件一模一样，即第二句完全没有起到作用。
猜测是我的map文件中SNP名称全用.来表示导致的。

1    .   0   100223803
1   .   0   100875122
1   .   0   102759472
1   .   0   102948333
1   .   0   103269622
1   .   0   106304192
1   .   0   109386355
1   .   0   112881004

随即修改第二列进行测试：第二列为SNP名称，需要与其他名称不同。整合第一列和第四列作为名称

1    1_100223803 0   100223803
1   1_100875122 0   100875122
1   1_102759472 0   102759472
1   1_102948333 0   102948333
1   1_103269622 0   103269622
1   1_106304192 0   106304192
1   1_109386355 0   109386355
1   1_112881004 0   112881004

plink --extract参数采坑问题记录相关推荐

MySQL单表数据量过千万，采坑优化记录，完美解决方案
MySQL单表数据量过千万,采坑优化记录,完美解决方案参考文章: (1)MySQL单表数据量过千万,采坑优化记录,完美解决方案 (2)https://www.cnblogs.com/ExMan/p/ ...
Centos7.9上利用cephadm安装Ceph Octopus 15.2的采坑记录，附带K8S挂载方法
Centos7.9上利用cephadm安装Ceph Octopus 15.2的采坑记录,附带K8S挂载方法 0.亮点 1 准备 1.1 修改历史记录 1.2 升级系统内核 1.3 配置免密登录问题1 ...
php给微信公众号接入聊天机器人程序+采坑记录
php给微信公众号接入聊天机器人程序今天逛了下我的公众号,突然心血来潮,想添加个自动聊天功能,于是-动手-!! 主要用到的api: 图灵机器人api 青云客智能聊天机器人API 茉莉机器人API 至 ...
微信支付采坑记录（java后端一：微信支付调用从哪一步开始）
普通商户模式微信支付之APP支付统一下单: 前言: 最近项目开发到微信支付的模块,经过一周的不懈努力,虽然微信支付流程跑通了,但当时的那些坑是真的烦,为了避免自己以后犯同样的错误,当然还有各位新接触微 ...
使用html2canvas,将页面转换成图片的采坑记录(Web/Taro h5)
使用html2canvas将页面转换成图片的采坑记录 "html2canvas": "^1.4.1","@tarojs/taro": &qu ...
MATLAB 非线性隐函数拟合采坑记录（使用 fsolve solve nlinfit lsqcurvefit函数）
MATLAB 非线性隐函数拟合采坑记录(使用 fsolve solve nlinfit lsqcurvefit函数) 问题描述解决思路错误示范1 代码思路原因解释模型更正更正模型1 更正模型 ...
iOS 微信SDK1.8.6后需要UniversalLink解决方案及采坑记录
项目最初因审核原因,一直使用iOS原生分享, 最近因项目需求要求, 接入微信分享, 以为和原来的没有区别, 但是接入时才发现改动的地方还是挺多的, 主要是需要配置UniversalLink和提包时的一 ...
H5拍照、预览、压缩、上传采坑记录
H5拍照.预览.压缩.上传采坑记录公司项目前段时间需要实现手机拍照上传的功能,本来以为用createObjectURL和canvas可以很轻松的实现,结果发现问题多多,特此记录下来. DEMO预览( ...
mysql8.0.19.0_分享MySql8.0.19 安装采坑记录
上篇文章给大家介绍了MySql8.0.19 安装过程,需要的朋友可以点击查看.https://www.jb51.net/article/178988.htm 1.ERROR 1820 (HY000): ...

plink --extract参数采坑问题记录

plink --extract采坑问题记录

问题描述：使用–extract 参数后没能正常进行LD修剪

plink --extract参数采坑问题记录相关推荐

最新文章

热门文章