code: def probit(x):
      """Probit function that ensures result is in (0, 1)"""
      eps = np.finfo(float).eps
      result = 0.5 + 0.5 * pt.erf(x / pt.sqrt(2))
      result = pt.switch(pt.eq(result, 0), eps, result)
      result = pt.switch(pt.eq(result, 1), 1 - eps, result)
  
      return result
  
file_path: bambi/bambi/backend/links.py
imports: import numpy
  import pytensor.tensor
label_desc: Use image node_id 1 for calling a global function with example usage: probit(x) and returns: result