powered by
kernel function
k(z, h = 1, type = "gaussian")
evaluate kernel at
bandwidth
either "gaussian" or "epanechnikov"
k(z/h)