Nonparametric Copula Density Estimation Using a Petrov–Galerkin Projection

Abstract

Nonparametrical copula density estimation is a meaningful tool for analyzing the dependence structure of a random vector from given samples. Usually kernel estimators or penalized maximum likelihood estimators are considered. We propose solving the Volterra integral equationof the given copula C. In the statistical framework, the copula C is not available and we replace it by the empirical copula of the pseudo samples, which converges to the unobservable copula C for large samples. Hence, we can treat the copula density estimation from given samples as an inverse problem and consider the instability of the inverse operator, which has an important impact if the input data of the operator equation are noisy. The wellknown curse of high dimensions usually results in huge nonsparse linear equations after discretizing the operator equation. We present a Petrov-Galerkin projection for the numerical computation of the linear integral equation. A special choice of test and ansatz functions leads to a very special structure of the linear equations, such that we are able to estimate the copula density also in higher dimensions.

On ill-posedness concepts, stable solvability and saturation

“…Solving an operator equation (1) with A := J d mapping in L 2 ([0, 1] d ) occurs in statistics when y represents a d-dimensional copula, from which an associated d-dimensional copula density x is to be reconstructed, and we refer for details to [22] and [25,Chap. 2.4], but with respect to numerical challenges in the context of this inverse problem to the report [29].…”