code: def tukey_pvalues(std_range, nm, df):
      # corrected but very slow with warnings about integration
      # nm = len(std_range)
      contr = contrast_allpairs(nm)
      corr = np.dot(contr, contr.T) / 2.0
      tstat = (
          std_range / np.sqrt(2) * np.ones(corr.shape[0])
      )  # need len of all pairs
      return multicontrast_pvalues(tstat, corr, df=df)
  
imports: from collections import namedtuple
  from statsmodels.compat.python import lzip, lrange
  import copy
  import math
  import numpy
  from numpy.testing import assert_almost_equal, assert_equal
  from scipy import stats, interpolate
  from statsmodels.iolib.table import SimpleTable
  from statsmodels.stats.multitest import multipletests, _ecdf, fdrcorrection, fdrcorrection_twostage
  from statsmodels.graphics import utils
  from statsmodels.tools.sm_exceptions import ValueWarning
label_desc: Use image node_id 26 for calling a global function with example usage: tukey_pvalues(std_range, nm, df) and returns: multicontrast_pvalues