Error M0 > 0 are not all TRUE

[kdkorthauer]

Hi,

I am running into the error M0 > 0 are not all TRUE using nulltype = "empirical" with the additional warning Warning message: In sqrt(-1/{ : NaNs produced. Looks like this is thrown from the following line of the efron() function:

sig = as.numeric(sqrt(-1/{ 2 * b0[3] }))

I’m not sure what this line is doing, or why it expects that b0[3] is negative, so I'm not sure how to fix.

Here’s a toy reproducible example where it occurs - when the test statistics look a bit bimodal (but I’m not sure if this is the only type of scenario):

set.seed(57)
res <- FDRreg::FDRreg(c(rnorm(500,1,0.8), rnorm(500,-1,0.8)), 
               features = model.matrix( ~  splines::bs(rnorm(1000), df = 3) - 1),
               nulltype = 'empirical’)

Please let me know if this can be resolved. Thanks!

Best,
Keegan