ashr
ashr copied to clipboard
doesn't handle zero sebetahats correctly
For example, ash(betahat = c(1, 2, 3, 4), sebetahat = c(0, 0, 1, 1))
gives me posterior means of zero for the first two observations.
We might want to handle infinite SEs too, possibly by throwing an error. (We can get infinite SEs in flash
.)
calc_pm
automatically sets these cases (covered by get_exclusions
) to have a mean of zero and a SD of zero. We should change this: zero SEs should return mean = x, sd = 0; infinite SEs and NAs should return prior mean and sd.
Issue #59 can be addressed as part of the same fix.
Infinite SEs also screw up the returned log likelihood.