rseg-diff: ../common/TwoStateScaleSplitResolveMixture.cpp:57: Assertion `val_weight > 0' failed

[MIRIAM GAGLIARDI]

Hi,

I ran RSEG with the following script:

#!/bin/bash

#PBS -o /home/userA/RSEG_protein1/RSEG_MG_30-6-15.out

#PBS -e /home/userA/RSEG_protein1/RSEG_MG_30-6-15.err

/home/userA/RSEG_protein1/rseg-0.4.8/bin/rseg-diff -p 3 -c /home/userA/RSEG_protein1/Hg19_chr_length.sorted.bed -o -out /home/userA/RSEG_protein1/RSEG_patient1_protein1 -score Posterior_scores_file_protein1_patien1_RSEG -readcount readcounts_file_protein1_patient1_RSEG -boundary domain_boundary_file_protein1_patien1_RSEG -boundary-score boundary_transition_scores_file_protein1_patient1_RSEG -i 20 -v -mode 2 -d ./deadzones-k50-hg19.bed new_ChIP_protein1_patient1_UNIQUE.sorted.toRSEG.bed new_ChIP_INPUT_patien1_UNIQUE.sorted.toRSEG.bed

but I got the following error:

rseg-diff: ../common/TwoStateScaleSplitResolveMixture.cpp:57: void initialize_values(const std::vector<double, std::allocator<double> >&, const std::vector<double, std::allocator<double> >&, const std::vector<double, std::allocator<double> >&, const std::vector<double, std::allocator<double> >&, SplitDistro&, SplitDistro&): Assertion `val_weight > 0' failed.

/var/spool/torque/mom_priv/jobs/18268.lilligridbio.na.iac.cnr.it.SC: line 7: 28990 Aborted                 /home/userA/RSEG_protein1/rseg-0.4.8/bin/rseg-diff -p 3 -c /home/userA/RSEG_protein1/Hg19_chr_length.sorted.bed -o -out /home/userA/RSEG_protein1/RSEG_patient1_protein1 -score Posterior_scores_file_protein1_patien1_RSEG -readcount readcounts_file_protein1_patient1_RSEG -boundary domain_boundary_file_protein1_patient1_RSEG -boundary-score boundary_transition_scores_file_protein1_patient1_RSEG -i 20 -v -mode 2 -d ./deadzones-k50-hg19.bed new_ChIP_protein1_patient1_UNIQUE.sorted.toRSEG.bed new_ChIP_INPUT_patient1_UNIQUE.sorted.toRSEG.bed

and in the .out file there was write:

[LOADING_DATA] chromosomes

[LOADING_DATA] reads

[LOADING_DATA] deadzones

[SELECTING BIN SIZE] bin size =  1550

[ESTIMATING PARAMETERS]

I searched some explanation about this error on the web and in the similar topics in this group but I didn't find any solution...

Do you have any suggestion?

Thanks a lot,

Miriam