code: def davies_pvalue(q, w):
      r
  
      q = asarray(atleast_1d(q), float)
      w = asarray(w, float)
  
      re = _pvalue_lambda(_lambda(w), q)
      param = dict()
      param["liu_pval"] = re["p_val_liu"][0]
      param["Is_Converged"] = re["is_converge"][0]
  
      return re["p_value"][0]
docs: r"""Joint significance of statistics derived from chi2-squared distributions.
  
      Parameters
      ----------
      q : float
          Test statistics.
      w : array_like
          Weights of the linear combination.
  
      Returns
      -------
      float
          Estimated p-value.