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 |
本文来自博客园,作者:橙子牛奶糖(陈文燕),转载请注明原文链接:https://www.cnblogs.com/chenwenyan/p/5855278.html
标签:
haploview
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 从 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生成工具