heterogen error in function boost::math::cdf

[Welmoed van Zuiden]

Hello,

I am trying to run heterogen analysis comparing two groups, and run into the following error:

2025-11-03 19:44:46,245 (PID:4179269) - INFO - Calculating statistics for LSVs 24732 to 27480 out of 95656 (progress: 25.9%-28.7%)
2025-11-03 19:44:54,549 (PID:4179269) - INFO - Calculating statistics for LSVs 27480 to 30228 out of 95656 (progress: 28.7%-31.6%)
terminate called after throwing an instance of 'boost::wrapexcept<std::domain_error>'
  what():  Error in function boost::math::cdf(const normal_distribution<double>&, double): Random variate x is -nan, but must be finite!

This has not happened before, and when I compare two different groups within the same data, it works fine. I'm using majiq v2.5.7.dev1 and the data is from human subjects, so in the majiq build, each subject was listed as a different experiment.

Would love to understand what is going on!

[San Jewell]

Hello <unknown name>, 

This is a new error. Can you describe what you mean by "compare two different groups in the same data", vs what you are trying to do, where it fails? A redacted outline/shell script of steps to reproduce the issue would be helpful. 

I also notice that this function call is only involved in the calculation of the 'mannwhitneyu' stat output. While we debug this issue, if you would like to continue other work, and you do not depend heavily on this stat, you can call majiq heterogen with other stats using the --stats argument, which should sidestep the issue in the short term. 

Thanks, 

-San