计算xpclr

1.conda安装xpclr

首先安装流程很轻松

复制代码
conda create -n xpclr  -c bioconda xpclr
conda activate xpclr
xpclr -h

2.按照要求准备文件

XPCLR - 简书 (jianshu.com)

根据教程准备文件,vcf,计算好的map,以及样本文件txt

其实官网也有介绍

3.使用软件计算xpclr

复制代码
for k in {1..3}; do   
xpclr 
--out chr${k} 
--format vcf 
--input gwas.phased${k}.vcf 
--samplesA g4.txt 
--samplesB g8.txt 
--map g2_new${k}.map 
--chr ${k} 
--size 500 
--step 500; done

在这里遇到了第一个报错:

这个报错很好解决,网上有解决教程,Python版的xpclr的艰辛debug之旅 - 简书 (jianshu.com)

根据这篇文章,在这个网站里面GitHub - xuzhougeng/xpclr: Code to compute the XP-CLR statistic to infer natural selection

把这个的地方给改了,再次运行。

它一直在输出。需要等待一会,中途会有log文件报错,应该不影响计算结果

结果可以看到文件

补充一点:

生信分析18:如何鉴定选择清除区域 - 简书 (jianshu.com)

选择信号分析XP-CLR原版Python版使用简介,特别提示窗口是遗传距离,需要根据不同物种修改 (qq.com)

根据这两篇文章。可以看这个,这个参数可以加上去,比计算xpehh的map文件多了两个00

还有一个问题是,但是我觉得两个群体只需要比较就行,没有强制要求位置

4.结果解读

有一些数据缺失,比起用原版来说

有一些可以参考和理解的文章

群体遗传分析必备教程:使用vcftools等工具批量计算Pi、Fst、TajimaD、XP-CLR的方法与脚本 (qq.com)

相关推荐
love530love11 小时前
命令行创建 UV 环境及本地化实战演示—— 基于《Python 多版本与开发环境治理架构设计》的最佳实践
开发语言·人工智能·windows·python·conda·uv
Ann's Blog1 天前
pycharm配conda环境
ide·pycharm·conda
饭来_2 天前
Python 使用环境下编译 FFmpeg 及 PyAV 源码(英特尔篇)
python·ubuntu·ffmpeg·conda
love530love4 天前
使用 Conda 工具链创建 UV 本地虚拟环境全记录——基于《Python 多版本与开发环境治理架构设计》
开发语言·人工智能·windows·python·机器学习·conda
Kaydeon6 天前
【Anaconda】Conda 虚拟环境打包迁移教程
人工智能·pytorch·python·conda
love530love6 天前
《Anaconda 精简路径治理》系列 · 番外篇Conda 虚拟环境路径结构方案全解——六种路径布局对比、优劣与治理建议
运维·人工智能·windows·python·conda
匆匆整棹还8 天前
服务器清理空间--主要是conda环境清理和删除
运维·服务器·conda
红藕香残玉簟秋9 天前
【python学习】windows使用conda管理python虚拟环境
windows·python·conda
思念殇千寻9 天前
使用conda隔离编译环境
conda
LiQiang3310 天前
conda activate 时报错: CondaError: Run ‘conda init‘ before ‘conda activate‘
conda