在使用plink1.9版本时遇到core dumped报错,运行命令如下:
/share/work/biosoft/plink/plink_v1.9/plink --vcf /share/nas1/wangq/project/pop_test/pop_cucu_demo/work/ann/1.qc_filter/filtered_snp.vcf.gz --make-bed --extract ld.prune.in --out LDfiltered --recode vcf-iid --keep-allele-order --allow-extra-chr --set-missing-var-ids @:#
报错信息如下:
PLINK v1.90b3.32 64-bit (24 Feb 2016) https://www.cog-genomics.org/plink2
(C) 2005-2016 Shaun Purcell, Christopher Chang GNU General Public License v3
Logging to LDfiltered.log.
Options in effect:
--allow-extra-chr
--extract ld.prune.in
--keep-allele-order
--make-bed
--out LDfiltered
--recode vcf-iid
--set-missing-var-ids @:#
--vcf /share/nas1/wangq/project/pop_test/pop_cucu_demo/work/ann/1.qc_filter/filtered_snp.vcf.gz
128701 MB RAM detected; reserving 64350 MB for main workspace.
--vcf: LDfiltered-temporary.bed + LDfiltered-temporary.bim +
LDfiltered-temporary.fam written.
9696 variants loaded from .bim file.
9696 missing IDs set.
114 people (0 males, 0 females, 114 ambiguous) loaded from .fam.
Ambiguous sex IDs written to LDfiltered.nosex .
--extract: 435 variants remaining.
Using 1 thread (no multithreaded calculations invoked.
Before main variant filters, 114 founders and 0 nonfounders present.
Calculating allele frequencies... done.
Total genotyping rate is 0.952531.
435 variants and 114 people pass filters and QC.
Note: No phenotypes present.
--make-bed to LDfiltered.bed + LDfiltered.bim + LDfiltered.fam ... done.
--recode vcf-iid to LDfiltered.vcf ... done.
Segmentation fault (core dumped)
尝试各种方法后使用2.0版本plink后没有再遇到该报错。
此外,我们在网易云课堂上有各种教学视频,有兴趣可以了解一下:
1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4分文章简单快速,学习链接:基因家族分析实操课程
2. 转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接:转录组(有参)结果解读;转录组(无参)结果解读
3. 转录组数据深入挖掘技能-WGCNA,提升你的文章档次,学习链接:WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘?学习链接:转录组标准分析后的数据挖掘
6. 更多学习内容:linux、perl、R语言画图,更多免费课程请点击以下链接:
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!