Error with rseg diff : "gsl/gamma.c ERROR" or "Assertion val_weight>0 failed"

[julie]

Hi,
I've tried to use rseg-diff with this command line :

/opt/rseg/bin/rseg-diff -c /usr/share/rseg/genomes/mouse.mm9.genome -d /usr/share/rseg/deadzones/deadzones-k36-mm9.bed -i 2 -mode 3 fileA.bed fileB.bed -o out_domains.bed -readcount out_counts.bed -boundary out_boundary_domains.bed -boundary-score out_boundary_scores.bed -bin-step 50 -s 20000 -fragment-length 0 -v

and this command failed with two different errors if I run several time :

This is the first type of error:

gsl: ../gsl/gamma.c:1185: ERROR: error
Default GSL error handler invoked.
Abandon (core dumped)

This is the second type of error:
rseg-diff: ../common/ThreeStateScaleSplitResolveMixture.cpp:65: void initialize_values(const std::vector<double>&, const std::vector<double>&, const std::vector<double>&, const std::vector<double>&, SplitDistro&, SplitDistro&, SplitDistro&): Assertion `val_weight > 0' failed.
Abandon (core dumped)

Is it because my files are not enough differentialy enriched ? or is it another problem?
Thanks

Julie

[Li Shen]

I'm having the same kind of problem. Has this been resolved? Is it because of large input file sizes?

- L

[Song, Qiang]

Hi Shen,

I have not figure out the exact cause of this issue since I don't have a test dataset to reproduce the error.
Is your dataset accessible by us?

Regards,
Song Qiang

--
You received this message because you are subscribed to the Google Groups "RSEG Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to rseg-support...@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.