error in logic of marginal confidence intervals

[jeff788]

• SOSAT version:
• Python version:
• Operating System:

Description

Currently the confidence intervals are computed by finding the bin where the marginal CDF exceeds (1-confidence range) and the confidence range. This works correctly for confidence ranges > 50% but breaks down for confidence intervals less than or equal to 50%.

What I Did

For confidence interval of 50% the current approach returns the median for both upper and lower limits, which is not correct. For values less than 50% it will return a lower limit that exceeds the upper limit.