Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Running FitHiChIP error for .hic data #91

Open
cauyrd opened this issue Feb 17, 2023 · 2 comments
Open

Running FitHiChIP error for .hic data #91

cauyrd opened this issue Feb 17, 2023 · 2 comments

Comments

@cauyrd
Copy link

cauyrd commented Feb 17, 2023

Hello, I have a .hic (h3k27ac HiChIP) data and a .bed.gz (H3k27ac peaks) downloaded from GEO and run FitHiChIP below:

./FitHiChIP_HiCPro.sh -C configfile_BiasCorrection_CoverageBias

Then I got errors below:

================ Processing input .hic file : /home/FitHiChIP/TestData/GSM5229037_VCaP_DMSO_4h_H3K27Ac.hic =================

Input - InpHiCFile : /home/FitHiChIP/TestData/GSM5229037_VCaP_DMSO_4h_H3K27Ac.hic
Input - OutFile : /home/FitHiChIP/TestData/results/Matrix_BinSize5000/VCaP_DMSO.interactions.initial_temp.bed
Input - TARGETRES : 5000
genomeID : hg38
[2500000, 1000000, 500000, 250000, 100000, 50000, 25000, 10000, 5000, 1000]
processing chromosome : ALL
Error finding block data
../FitHiChIP_HiCPro.sh: line 1185: 31124 Floating point exception(core dumped) $PythonExec ./src/Process_contactMat_HiC.py --hicfile ${InpHiCFile} --outfile ${Interaction_Initial_File_temp} --chromsizefile ${ChrSizeFile} --resolution ${BIN_SIZE}

Any suggestions why I got this error? You can find the GEO accession number based on the file name: GSM5229037 (.hic) and GSM5228990 (.bed.gz)

@ay-lab
Copy link
Owner

ay-lab commented Mar 4, 2023

Hi @cauyrd
Could you please edit the line 52 of the script "src/Process_contactMat_HiC.py" by the following line?
if chrom.name != "All" or chrom.name != "ALL":

Let me know if it works. I'll then update the release.

@danidigi
Copy link

Hi all,
I have the same error when running FitHiChIP using my .hic file as input. Did your change work?
Thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants