code: def mu_adj_score(y, loc, scale, shape, skewness):
      try:
          if (y - loc) >= 0:
              return (
                  ((shape + 1.0) * power(y - loc, 2))
                  / float(
                      power(skewness, 2) * shape * exp(scale)
                      + power(y - loc, 2)
                  )
              ) - 1.0
          else:
              return (
                  ((shape + 1.0) * power(y - loc, 2))
                  / float(
                      power(skewness, -2) * shape * exp(scale)
                      + power(y - loc, 2)
                  )
              ) - 1.0
      except:
          return -1.0
  
imports: from numpy import abs, exp, power, array, sqrt, pi
  from scipy.special import gamma
label_desc: Use image node_id 4 for calling the SkewBetatScore obj's underlying member method code with example usage: obj.mu_adj_score(y, loc, scale, shape, skewness) and returns: unknown, unknown, unknown