def infid_vs_th(x): cdf = ro_CDF(x=[x, x], **bvn) return (1 - np.abs(cdf[0] - cdf[1])) / 2
def CDF(x): return ro_CDF(x=x, **bvn)