code: def get_tukey_pvalue(k, df, q):
      """
      return adjusted p-values for Tukey's HSD
  
      Parameters
      ----------
      k : int in {2, ..., 10}
          number of tests
      df : int
          degrees of freedom of error term
      q : scalar, array_like; q >= 0
          quantile value of Studentized Range
  
      """
      return studentized_range.sf(q, k, 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 3 for calling a global function with example usage: get_tukey_pvalue(k, df, q) and returns: studentized_range