code: def gaussian_fdr_threshold(x, alpha=0.05):
      """Return FDR threshold given normal variates
  
      Given an array x of normal variates, this function returns the
      critical p-value associated with alpha.
      x is explicitly assumed to be normal distributed under H_0
  
      Parameters
      -----------
      x: ndarray
          input data
      alpha: float, optional
          desired significance
  
      Returns
      -------
      threshold : float
          threshold, given as a Gaussian critical value
      """
      pvals = st.norm.sf(x)
      pth = fdr_threshold(pvals, alpha)
      return st.norm.isf(pth)
  
file_path: nipy/nipy/algorithms/statistics/empirical_pvalue.py
imports: import numpy
  import scipy.stats
  from numpy.linalg import pinv
label_desc: Use image node_id 3 for calling a global function with example usage: gaussian_fdr_threshold(x, alpha) and returns: st