def random_normal(): """ returns a random draw from normally distributed data""" return inverse_normal_cdf(random.random())
def normal_lower_bound(probability, mu=0, sigma=1): """ returns the z for which P(Z > z) = probability """ return inverse_normal_cdf(1-probability, mu, sigma)