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

chicViewpoint #447

Closed
LeilyR opened this issue Oct 24, 2019 · 4 comments
Closed

chicViewpoint #447

LeilyR opened this issue Oct 24, 2019 · 4 comments

Comments

@LeilyR
Copy link
Collaborator

LeilyR commented Oct 24, 2019

chicViewpoint fails with the following error message:
miniconda3/envs/hicexplorer/lib/python3.6/site-packages/hicexplorer/lib/viewpoint.py", line 571, in computeSumOfDensities
background_nbinom[fixateRange] = nbinom(pBackgroundModel[fixateRange][0], pBackgroundModel[fixateRange][1])
KeyError: 500000

@joachimwolff
Copy link
Collaborator

joachimwolff commented Oct 24, 2019 via email

@LeilyR
Copy link
Collaborator Author

LeilyR commented Oct 24, 2019

The error comes from background model computation. The fixateRange in there is 500,000 and the resolution is 2,000 , so the last computed distance is 500,000 - 2,000 = 498,000 . The solution would be either adding the last distance to the bg model or calling fixateRange -1 here.

@LeilyR
Copy link
Collaborator Author

LeilyR commented Oct 24, 2019

background model:
chicViewpointBackgroundModel -m mergedSamples_merged_bs2kb.corrected.h5
--referencePoints refPoint.bed -t 30
--outFileName bgModel.bed

viewPoint:

chicViewpoint -m mergedSamples_merged_bs2kb.corrected.h5
--referencePoints refPoint.bed -t 30
--backgroundModelFile bgModel.bed
--range 100000 100000 -o interactionFolder

@LeilyR LeilyR mentioned this issue May 18, 2020
19 tasks
@joachimwolff
Copy link
Collaborator

I changed so much in chic-modules in the last half year. Could be fixed by now the current release, or maybe this line of code is not even more there.

Please reopen in case it is still crashing.

Thanks,

Joachim

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

2 participants