I was wondering if scipy/numpy has the inverse cumulative normal function, ie the function f in this expression f(scipy.stats.norm.cdf(1.2)) = 1.2 or more generally, a function f which fits the criteria f(scipy.stats.norm.cdf(x)) = x There is a distribution called invnorm, but I am not sure of how to use it.