code: def _kernel_pdf_gamma(x, sample, bw):
      """Gamma kernel for pdf, without boundary corrected part.
  
      drops `+ 1` in shape parameter
  
      It should be possible to use this if probability in
      neighborhood of zero boundary is small.
  
      """
      return stats.gamma.pdf(sample, x / bw, scale=bw)
  
imports: import numpy
  from scipy import special, stats
label_desc: Use image node_id 9 for calling a global function with example usage: _kernel_pdf_gamma(x, sample, bw) and returns: stats