code: def pairwise_tukeyhsd(endog, groups, alpha=0.05):
      """
      Calculate all pairwise comparisons with TukeyHSD confidence intervals
  
      Parameters
      ----------
      endog : ndarray, float, 1d
          response variable
      groups : ndarray, 1d
          array with groups, can be string or integers
      alpha : float
          significance level for the test
  
      Returns
      -------
      results : TukeyHSDResults instance
          A results class containing relevant data and some post-hoc
          calculations, including adjusted p-value
  
      Notes
      -----
      This is just a wrapper around tukeyhsd method of MultiComparison
  
      See Also
      --------
      MultiComparison
      tukeyhsd
      statsmodels.sandbox.stats.multicomp.TukeyHSDResults
      """
  
      return MultiComparison(endog, groups).tukeyhsd(alpha=alpha)
  
label_desc: Use image node_id 1 for calling a global function with example usage: pairwise_tukeyhsd(endog, groups, alpha) and returns: MultiComparison