bandwidth_pi function

hr_kde_pi wraps KernSmooth::dpik to select bandwidth for kernel density estimation the plug-in-the-equation method in two dimensions.