code: def quantiles(x, qlist=(2.5, 25, 50, 75, 97.5)):
      
  
      # Make a copy of trace
      x = x.copy()
  
      # For multivariate node
      if x.ndim > 1:
          # Transpose first, then sort, then transpose back
          sx = sort(x.T).T
      else:
          # Sort univariate node
          sx = sort(x)
  
      try:
          # Generate specified quantiles
          quants = [sx[int(len(sx) * q / 100.0)] for q in qlist]
  
          return dict(zip(qlist, quants))
  
      except IndexError:
          print_("Too few elements for quantile calculation")
docs: Returns a dictionary of requested quantiles from array
  
      :Arguments:
        x : Numpy array
            An array containing MCMC samples
        qlist : tuple or list
            A list of desired quantiles (defaults to (2.5, 25, 50, 75, 97.5))