Getting error when running modkit dmr pair

[agoel13]

Hi,

I am trying to run modkit (v0.3.0) dmr pair using replicates and getting following error

Error msg :
running with replicates, but not matched samples
thread 'main' panicked at src/hmm.rs:330:9:
e_s is not finite -inf
note: run with RUST_BACKTRACE=1 environment variable to display a backtrace

For your reference here is my script:

modkit dmr pair
-a Norm1.bed.gz
-a Norm2.bed.gz
-b Tum1.bed.gz
-b Tum2.bed.gz
-b Tum3.bed.gz
-b Tum4.bed.gz
-o dmr_replicates_pair_segment.bed
--header
--segment dmr_segements.bed
--ref hg38.fa
--base C
--threads 16
--log-filepath dmr_pair_segment.log

While when I run the same script with only one pair of samples and no replicates, it completes without any errors.

PS: Without turning on the segment argument, the above script runs fine.

[ArtRand]

Hello @agoel13,

Sorry about that. Could you try the attached build and tell me if you continue to encounter the crash?

modkit_dev88c2107_centos7_x86_64.tar.gz

[agoel13]

Many thanks. Completed successfully without any error msgs.
When would this version be released? I would then update my pipeline with the released version.

[ArtRand]

@agoel13 Great. Yes this fix and a few others will be released this week. Thanks for testing it!

[agoel13]

Thanks. Looking forward to it.

[ArtRand]

This fix has been released in v0.3.1.