Empirical coverage calculation for P(D) analysis