haploview出现“more than two alleles”的解决方法

弹出“more than two alleles”的错误是因为ped文件中一个SNP位点上存在两个以上的等位基因,haploview连锁分析时默认为只有两个等位基因,因此我们要去掉超过两位等位基因的SNP才能做连锁分析。

解决方案一:针对vcf格式的文件;

用到命令:--min-alleles 2 --max-alleles 2

先在linux上下载安装vcftools软件,下载地址“https://vcftools.github.io/man_latest.html”

具体命令如下 :

1
2
/vcftools --gzvcf  your.vcf.gz --remove-indels --recode --recode-INFO-all --out yourvcf_snps
/vcftools --vcf /to/your/pathway/yourvcf_snps.vcf --min-alleles 2 --max-alleles 2 --plink --out /to/your/pathway/file

  

“/vcftools”指vcftools软件在终端上的路径

“file”是指你想保存的文件名字。

VCFtools的所有命令见此网页:https://vcftools.github.io/man_latest.html

解决方案二:如果原始文件是plink格式,也可以用以下格式修改:

1
/software/plink --noweb --file filetext --snps-only --recode --out filetext_snponly

  

 

posted @   橙子牛奶糖  阅读(2629)  评论(28编辑  收藏  举报
编辑推荐:
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
· 开发者必知的日志记录最佳实践
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 上周热点回顾(3.3-3.9)
· AI 智能体引爆开源社区「GitHub 热点速览」
· 写一个简单的SQL生成工具
点击右上角即可分享
微信分享提示