code: def kernel_pdf_invgauss(x, sample, bw):
      # Inverse gaussian kernel for density, pdf, estimation
      m = x
      lam = 1 / bw
      return stats.invgauss.pdf(sample, m / lam, scale=lam)
  
imports: import numpy
  from scipy import special, stats
label_desc: Use image node_id 15 for calling a global function with example usage: kernel_pdf_invgauss(x, sample, bw) and returns: stats