Normal approximations for the multivariate inverse Gaussian distribution and asymmetric kernel smoothing on $d$-dimensional half-spaces

TL;DR

This paper develops a new asymmetric kernel density estimator on $d$-dimensional half-spaces using the multivariate inverse Gaussian distribution, analyzing its properties, providing a new sampling algorithm, and demonstrating its effectiveness through simulations and an application.

Contribution

It introduces the first asymmetric kernel density estimator supported on half-spaces, based on the MIG distribution, with theoretical analysis, a new sampling algorithm, and practical evaluation.

Findings

The estimator has desirable boundary properties.

The new sampling algorithm outperforms previous methods.

Simulation studies show effective bandwidth selection.

Abstract

This paper introduces a novel density estimator supported on $d$-dimensional half-spaces. It stands out as the first asymmetric kernel density estimator for half-spaces in the literature. Using the multivariate inverse Gaussian (MIG) density from Minami (2003) as the kernel and incorporating locally adaptive parameters, the estimator achieves desirable boundary properties. To analyze its mean integrated squared error (MISE) and asymptotic normality, a local limit theorem and probability metric bounds are established between the MIG and the corresponding multivariate Gaussian distribution with the same mean vector and covariance matrix, which may also be of independent interest. Additionally, a new algorithm for generating MIG random vectors is developed, proving to be faster and more accurate than Minami's algorithm based on a Brownian first-hitting location representation. This algorithm is then used to discuss and compare optimal MISE and likelihood cross-validation bandwidths for the estimator in a simulation study under various target distributions. As an application, the MIG asymmetric kernel is used to smooth the posterior distribution of a generalized Pareto model fitted to large electromagnetic storms.