1. 原版
安裝比較簡(jiǎn)單耐朴。
wget https://reich.hms.harvard.edu/sites/reich.hms.harvard.edu/files/inline-files/XPCLR.tar
tar xvf XPCLR.tar
直接運(yùn)行bin下的XPCLR即可衡楞，若不能運(yùn)行吱雏，則編譯下：
cd src
make
make install
原版分析時(shí)準(zhǔn)備文件過(guò)程較為繁瑣，因此更建議使用Python版。
$ /project/biosoft/XPCLR/bin/XPCLR -h
Usage:
XPCLR -xpclr hapmapInput1 hapmapInput2 mapInput outFile -w gWin(Morgan) snpWin gridSize(bp) chrN -p corrLevel
-w1: gWin sets the size of a sliding window(units: 100cM),sWin sets # of SNPs in a window. otherwise, no sliding window
-p1:the input genotpe is already phased. -p0: the input genotype is not phased
corrLevel: the value is on (0,1], set corrLevel equal to 0 if no correction is needed
2. Python版本
安裝
conda create -n xpclr -c bioconda xpclr
conda activate xpclr
報(bào)錯(cuò)：
$ xpclr -h
Traceback (most recent call last):
File "/home/miniconda3/envs/xpclr/bin/xpclr", line 5, in <module>
import xpclr
File "/home/miniconda3/envs/xpclr/lib/python2.7/site-packages/xpclr/__init__.py", line 3, in <module>
from xpclr import methods
File "/home/miniconda3/envs/xpclr/lib/python2.7/site-packages/xpclr/methods.py", line 11, in <module>
from functools import lru_cache
ImportError: cannot import name lru_cache
重新安裝pip install lru_cache仍然報(bào)錯(cuò)歧杏。
由于python版本引起镰惦，直接修改methods.py中from functools import lru_cache的為：
try:
from functools import lru_cache
except ImportError:
from backports.functools_lru_cache import lru_cache
運(yùn)行
xpclr --format vcf --input /project/04.sweep/sample750_miss0.6_impute/meanDP3.miss0.6.maf0.01.impute.rename.vcf \
--samplesA /project/04.sweep/sample750_miss0.6_impute/List/w-l-c/Cultivar.list \
--samplesB /project/04.sweep/sample750_miss0.6_impute/List/w-l-c/Wild.list \
--chr 1 --maxsnps 600 --size 1000 --step 1000 --out test_out
報(bào)錯(cuò)：
Traceback (most recent call last):
File "/home/miniconda3/envs/xpclr/bin/xpclr", line 195, in <module>
main()
File "/home/miniconda3/envs/xpclr/bin/xpclr", line 88, in main
"No permission to write in the specified directory: {0}".format(outdir)
AssertionError: No permission to write in the specified directory:
xpclr第84行fn = args.out修改為：
fn = os.path.abspath(args.out)
運(yùn)行繼續(xù)報(bào)錯(cuò)：
2023-03-17 19:05:48 : INFO : running xpclr v1.1.0
2023-03-17 19:05:48 : INFO : Loading VCF
Traceback (most recent call last):
File "/home/miniconda3/envs/xpclr/bin/xpclr", line 196, in <module>
main()
File "/home/miniconda3/envs/xpclr/bin/xpclr", line 103, in main
gdistkey=args.gdistkey)
File "/home/miniconda3/envs/xpclr/lib/python2.7/site-packages/xpclr/util.py", line 112, inload_vcf_format_data
pos1, geno1 = load_vcf_wrapper(vcf_fn, chrom, samples1)
File "/home/miniconda3/envs/xpclr/lib/python2.7/site-packages/xpclr/util.py", line 94, in load_vcf_wrapper
callset = allel.read_vcf(
AttributeError: 'module' object has no attribute 'read_vcf'
網(wǎng)上沒(méi)有我的同類(lèi)安裝相關(guān)錯(cuò)誤報(bào)道，查看了下allel模塊得滤，確實(shí)沒(méi)有read_vcf函數(shù)：
Python 2.7.18 |Anaconda, Inc.| (default, Nov 25 2022, 06:27:37)
[GCC 11.2.0] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> import allel
>>> allel
<module 'allel' from '/home/miniconda3/envs/xpclr/lib/python2.7/site-packages/scikit_allel-0.20.3-py2.7-linux-x86_64.egg/allel/__init__.py'>
>>> dir(allel)
['AlleleCountsArray', 'AlleleCountsCArray', 'AlleleCountsCTable', 'AlleleCountsChunkedArray', 'AlleleCountsChunkedTable', 'AlleleCountsDaskArray', 'FeatureCTable', 'FeatureChunkedTable', 'FeatureTable', 'GenotypeArray', 'GenotypeCArray', 'GenotypeChunkedArray', 'GenotypeDaskArray', 'HaplotypeArray', 'HaplotypeCArray', 'HaplotypeChunkedArray', 'HaplotypeDaskArray', 'SortedIndex', 'SortedMultiIndex', 'UniqueIndex', 'VariantCTable', 'VariantChunkedTable', 'VariantTable', '__builtins__', '__doc__', '__file__', '__loader__', '__name__', '__package__', '__path__', '__version__', '_bcolz', '_da', 'chunked', 'compat', 'constants', 'io', 'model', 'plot', 'stats', 'util']
但實(shí)際上官網(wǎng)上是有這個(gè)函數(shù)的：
image.png
難道是版本問(wèn)題嗎陨献？裝了下python3版本，發(fā)現(xiàn)確實(shí)是有的：
image.png
試圖重新安裝懂更，舊版本刪除不了：
$ pip uninstall scikit-allel
DEPRECATION: Python 2.7 will reach the end of its life on January 1st, 2020. Please upgrade your Python as Python 2.7 won't be maintained after that date. A future version of pip will dropsupport for Python 2.7. More details about Python 2 support in pip, can be found at https://pip.pypa.io/en/latest/development/release-process/#python-2-support
ERROR: Cannot remove entries from nonexistent file /home/pengjx/miniconda3/envs/xpclr/lib/python2.7/site-packages/easy-install.pth
升級(jí)也升級(jí)不了：
>pip install --upgrade --ignore-installed scikit-allel -i https://pypi.tuna.tsinghua.edu.cn/simple
.......
ERROR: Command errored out with exit status 1:
沒(méi)辦法放棄conda安裝眨业，直接從GitHub安裝：https://github.com/hardingnj/xpclr。
git clone https://github.com/hardingnj/xpclr.git
cd xpclr
python setup.py install
進(jìn)入bin后沮协，可直接運(yùn)行xpclr：
$ xpclr -h
usage: xpclr [-h] --out OUT [--format FORMAT] [--input INPUT] [--gdistkey GDISTKEY] [--samplesA SAMPLESA] [--samplesB SAMPLESB] [--rrate RRATE] [--map MAP] [--popA POPA] [--popB POPB]
--chr CHROM [--ld LDCUTOFF] [--phased] [--verbose VERBOSE] [--maxsnps MAXSNPS] [--minsnps MINSNPS] [--size SIZE] [--start START] [--stop STOP] [--step STEP]
Tool to calculate XP-CLR as per Chen, Patterson, Reich 2010
需要注意的是龄捡，你如果使用xpclr的全路徑是運(yùn)行不了的：
$ /project/xpclr/bin/xpclr
Traceback (most recent call last):
File "/project/xpclr/bin/xpclr", line 4, in <module>
import numpy as np
ImportError: No module named numpy
安裝依賴(lài)包：pip install -r requirements.txt -i https://pypi.tuna.tsinghua.edu.cn/simple。
這時(shí)慷暂，也可以用python（實(shí)際上我當(dāng)前的版本已經(jīng)是Python 3.9.1聘殖，所以xpclr也是支持的。）來(lái)調(diào)用全路徑：
python /project/xpclr/bin/xpclr --format vcf --input /project/04.sweep/sample750_miss0.6_impute/meanDP3.miss0.6.maf0.01.impute.rename.vcf --samplesA /project/04.sweep/sample750_miss0.6_impute/List/w-l-c/Cultivar.list --samplesB /project/04.sweep/sample750_miss0.6_impute/List/w-l-c/Wild.list --chr 1 --maxsnps 600 --size 1000 --step 1000 --out test_out
建議還是將安裝路徑/project/xpclr/bin加入環(huán)境變量（測(cè)試了下行瑞，貌似不用加入環(huán)境變量也可直接調(diào)用奸腺，可能是軟件安裝環(huán)節(jié)已經(jīng)加入），直接用xpclr血久。
原因分析
conda安裝雖然便捷突照，但作者沒(méi)有及時(shí)更新，最早的版本也是在3年前（Python2）氧吐，導(dǎo)致一些包不兼容讹蘑。
image.png
而GitHub版本作者還是在維護(hù)的，所以對(duì)于后來(lái)者推薦之筑舅。
image.png
后續(xù)分析推薦
原版XP-CLR用法參考：https://zhuanlan.zhihu.com/p/145387269
Python版XP-CLR用法參考：http://www.reibang.com/p/9c827a0be66d
