Practical Representations of Copula and Joint Density Estimates

- A moment-based approximation methodology for estimating a copula density from bivariate observations is introduced. The resulting simple representation of the copula density is suitable for reporting purpose or carrying out further algebraic manipulation. Empirical copula density functions will also be determined from kernel density estimates. A technique for obtaining a joint density from marginal density estimates and a copula density is proposed as well. The Bernstein copula density approximants will be utilized for comparison purposes. The results are applied to two stocks’ closing prices and a stock’s price and its running maximum. In the latter case, the model is related to a Brownian motion process.


Introduction
Copulas are principally utilized for modelling dependencies in multivariate distributions.Their properties have been increasingly exploited in numerous types of scientific investigations; for instance, the reader may refer to Chao et al. [1], Carrera et al. [2] and Chen et al. [3] for recent advances in the area of artificial intelligence.The key idea behind copulas is that the joint distribution of two or more variables can be represented in terms of their marginal distributions and a specific correlation structure.As a measure of dependence, they have for instance found applications in reliability theory, signal processing, geodesy, hydrology and medicine.Results involving empirical bivariate copula densities are discussed in this paper.
Let ( 1 ,  2 ) be the joint cumulative distribution function of random variables  1 and  2 having continuous marginal distribution functions  1 ( 1 ) and  2 ( 2 ) .According to Sklar [4], there exists a unique bivariate copula :  2 ↦  (the unit interval) such that where (⋅,⋅) is a joint cumulative distribution function having uniform marginals.Conversely, for any continuous cumulative distribution function  1 ( 1 ) and  2 ( 2 ) and any copula  , the function  defined in Equation ( 1) is a joint distribution function with marginal distributions  1 and  2 .Sklar's theorem provides a technique for constructing copulas.Indeed, the function is a bivariate copula, where the quasi-inverse   −1 for  = 1, 2 is defined by We shall denote the probability density function corresponding to the copula ( 1 ,  2 ) by The following relationship between the joint density function (⋅,⋅) and the copula density function (⋅,⋅) can readily be obtained from Equation (1): where  1 ( 1 ) and  2 ( 2 ) respectively denote the marginal density functions of  1 and  2 .Accordingly, the copula density function can be expressed as follows: The proposed approaches to estimating copula density functions are summarized in Sections 2 and 3.The estimation of a joint density function in terms of marginal distributions and a copula density estimate is discussed in Section 4. Three illustrative examples are provided in the last section.

Copula Density Based on Kernel Density Estimation
On applying the kernel density estimation (kde) method to a given bivariate dataset, one can obtain an estimate of the bivariate probability density of .In light of Equation ( 6), the copula density can be represented as follows: where   can be estimated by a bivariate kde and    (⋅) denotes the quantile function.
The marginal density function of each variable can be obtained by determining their respective kernel density estimates.The inverse cdfs   1 and   2 can be determined in polynomial form by making use of a moment-based method or the method of least squares.

A Moment-based Bivariate Polynomial Approximation of the Copula Density
Once a copula density is determined from Equation (7), it can be approximated by the product of a base density and a bivariate polynomial whose coefficients are obtained from the joint moments associated with the copula density.
The proposed procedure for achieving this is described in the next result.
The degree  used in the polynomial adjustment should be selected so that   provides an accurate approximation to the estimate of the copula density, which can be determined for instance by evaluating their integrated squared differences.
Suppose that observations are available on the random variables  1 and  2 .The marginal densities  1 ( 1 ) and  2 ( 2 ) are estimated and a Bernstein's copula density of high order is determined.A joint density estimate of ( 1 ,  2 ) can then be obtained via Equation (5).

Applications
The results are applied to two stocks' closing prices as well as a stock's price and its running maximum.

Copula Density Estimation Methodologies Applied to Two Stocks
The two stocks selected are GOOG (Alphabet Inc.) and AAPL (Apple Inc.).The bivariate data are the daily closing prices of (GOOG, AAPL) in 2019.Each component of the data has been standardized.The joint KDE of , the marginal densities of  1 and  2 and the corresponding cdf inverses are displayed in Figures 1-5.The resulting KDE based (Section 2) and momentbased (Section 3) copula densities plotted in Figures 6 and 7 are seen to be similar.

Determination of a Joint Density Estimate from the Marginal Density Estimates and the Copula Density
Using the two stocks' standardized data, a Bernstein's copula density of degree 250, which is shown in Figure 8, is utilized as an estimate of the underlying copula density.The density estimate (Figure 9) resulting from the application of the methodology described in Section 4 and a KDE (Figure 10) exhibit very similar features.

Copula Associated with a Brownian Motion Process and Its Running Maximum
The data consists of the daily closing prices of AC.TO (Air Canada) during 2019.To relate the data to a standard Wiener process, the first data point should be zero, the differences between successive observations should ideally often change signs and have a variance of one, and there should be one unit of time between successive observations.Hence the following transformation is used.
The   's and   's are plotted in Figures 11 and 12.The resulting copula density as obtained by applying the methodology described in Section 2, is shown in Figure 13.

Conclusion
As copula density estimates are usually expressed in complicated forms, the bivariate polynomial approximation that is proposed in this paper ought to prove more suitable for reporting purposes.Approximations by means of Bernstein polynomials and the kernel density estimation approach are discussed as well.Additionally, a flexible technique for estimating joint density functions is introduced.The proposed methodologies were successfully applied to two stocks' closing prices as well as a set of observations and its running maximum.

Figure 7 .
Figure 7. Moment-based bivariate polynomial estimate of the copula density.

Figure 11 .
Figure 11.List Plot of the   's.

Figure 12 .
Figure 12.List Plot of the   's.