code: def kernel_pdf_invgamma(x, sample, bw):
      # Inverse gamma kernel for density, pdf, estimation
      return stats.invgamma.pdf(sample, 1 / bw + 1, scale=x / bw)
  
imports: import numpy
  from scipy import special, stats
label_desc: Use image node_id 13 for calling a global function with example usage: kernel_pdf_invgamma(x, sample, bw) and returns: stats