Variable bandwidth kernel density estimation