使用HiC-Pro分析完Hi-C数据之后,打算用HiCExplorer对TAD进行分析绘图。
按照HiCExplorer官网建议,
“The fastet way to obtain Python 3.6 or 3.7 together with numpy and scipy is via the Anaconda Scientific Python Distribution. Just download the version that’s suitable for your operating system and follow the directions for its installation. All of the requirements for HiCExplorer can be installed in Anaconda with:
$ conda install hicexplorer -c bioconda -c conda-forge”
尝试

conda create -n hicexplorer python=3.7 #创建软件安装环境
conda activate hicexplorer 
conda install numpy -c conda-forge
conda install scipy -c conda-forge
conda install hicexplorer -c bioconda -c conda-forge

numpy、scipy安装都很顺利,但到hicexplorer时就卡住了。一直卡在Solving environment,安装失败。
开始以为是网络的问题,设置了.condarc的proxy,但并没有用。
后来删掉环境,重新建立,没有指定python3的具体版本,让conda默认安装了3.10,并且numpy、scipy、hicexplorer整合为一行命令。

conda create -n hicexplorer python=3 #创建软件安装环境
conda activate hicexplorer
conda install numpy scipy hicexplorer -c conda-forge -c bioconda

等了一会,总算是安装成功了。

HiCExplorer要求contact matrix以h5格式存放,因此需要对HiC-Pro的结果进行转换。

hicConvertFormat \
    -m /mnt/data/GSE164_HiC_out/hic_results/matrix/SRR164/iced/2000/SRR164_2000_iced.matrix \
    -o SRR164_2000 \
    --bedFileHicpro /mnt/data/GSE164_HiC_out/hic_results/matrix/SRR164/raw/2000/SRR164_2000_abs.bed \
    --inputFormat hicpro --outputFormat h5

输出SRR164_2000.h5文件,就可以进行HiCExplorer后续的分析和绘图了。

标签: Hi-C, HiC-Pro, HiCExplorer