function [lower_bound, upper_bound] = ModWald(numerator,denominator)

p_prime = (numerator+2)/(denominator+4);
dev = 1.96*sqrt((p_prime*(1-p_prime))/(denominator+4));
lower_bound = max(p_prime - dev,0);   % Forces the value to be nonzero
upper_bound = p_prime + dev;