def probit(p): return -sqrt(2.) * erfcinv(2. * p)
def inv_phi(y): return -tt.sqrt(2) * tt.erfcinv(2 * y)
def probit(p): return -sqrt(2) * erfcinv(2 * p)
def probit(p): return -sqrt(2.0) * erfcinv(2.0 * p)
def probit(p): return -sqrt(2.) * erfcinv(2. * p)
def probit(p): return -sqrt(2) * erfcinv(2 * p)