Open lindenb opened 2 years ago
as far as I can see the error is created by a NAN in MixtureChiSquare.cpp:getLiuPvalue
cdfchn(&which, &p, &q, &x, &l, &ncp, &status, &bound);
with X is NAN
the NAN value is created by a sqrt with a negative value sqrt(VarQ)
in computeIntegrandDavies
modifying if (kappa > lambda.sum() * 10000) {
with if (kappa > lambda.sum() * 10000 || VarQ <=0.0) {
seems to fix the error but I don't know if it is the correct solution. Nevertheles, I'll submit a PR.
Hi, find in attachement where rvtest seems stalled/frozen (skato test ?) while there is a small number of variant.
the test.sh contains:
output:
would it be a bug ? how can I prevent rvtest from freezing ? thanks !
test.zip