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

segment-cnv and plot-cnv issue #65

Open
BERICP01 opened this issue Apr 26, 2024 · 0 comments
Open

segment-cnv and plot-cnv issue #65

BERICP01 opened this issue Apr 26, 2024 · 0 comments

Comments

@BERICP01
Copy link

BERICP01 commented Apr 26, 2024

Hello!
I am facing two problems while running neoloopfinder

  1. when I run segment-cnv I continue to get this error with most of my cool matrix at different resolutions:

ARGUMENT LIST:
CNV Profile = /gpfs/data/HernandoLab/home/bericp01/HiC/results/CNV/501mel_sgCHD7-2_NeoLoopFinder/501mel_sgCHD7-2_25k.CNV-profile.bedGraph
Number of States = None
CBS P-value Cutoff = 1e-05
Maximum HMM-CBS distance = 4
Minimum Copy Number difference = 0.4
Minimum Segment Size = 3
Ploidy = 2
Bin Size = 25000
Output Path = /gpfs/data/HernandoLab/home/bericp01/HiC/results/CNV/501mel_sgCHD7-2_NeoLoopFinder/501mel_sgCHD7-2_25k.CNV-seg.bedGraph
Number of Processes = 10
Log file name = cnv-seg.log
root INFO @ 04/26/24 18:51:53: Loading CNV profile ...
neoloop.cnv.segcnv INFO @ 04/26/24 18:51:54: Training HMM with 3 states...
Traceback (most recent call last):
File "/gpfs/data/HernandoLab/home/bericp01/miniconda_bericp01/envs/NeoLoopFinder/bin/segment-cnv", line 111, in run
work.segment(
File "/gpfs/data/HernandoLab/home/bericp01/miniconda_bericp01/envs/NeoLoopFinder/lib/python3.10/site-packages/neoloop/cnv/segcnv.py", line 167, in segment
model = HiddenMarkovModel.from_samples(
File "pomegranate/hmm.pyx", line 3705, in pomegranate.hmm.HiddenMarkovModel.from_samples
ValueError: need at least one array to concatenate

  1. all my plots are not centered to the 0 but rather always below zero irrespective if are cancer or normal cells. Any clue?
    501mel_sgCHD7-2_25k CNV genome-wide

Any help will be really appreciate!!!
Thanks
Pietro

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

1 participant