Combining local and global smoothing in multivariate density estimation