code: def get_tukeyQcrit2(k, df, alpha=0.05):
      """
      return critical values for Tukey's HSD (Q)
  
      Parameters
      ----------
      k : int in {2, ..., 10}
          number of tests
      df : int
          degrees of freedom of error term
      alpha : {0.05, 0.01}
          type 1 error, 1-confidence level
  
  
  
      not enough error checking for limitations
      """
      return studentized_range.ppf(1 - alpha, 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 2 for calling a global function with example usage: get_tukeyQcrit2(k, df, alpha) and returns: studentized_range