Introduction
Copula can express the nonlinear dependency between two or more variables, in comparison with correlation that can capture linear dependency between variables. It is a technique used to implement the joint distribution of the many variables (two or more). This has many wide applications in modern era of science like hydrology and geophysics (Valle & Kaplan, 2019) (Liu et al., 2018) (Salvadori et al., 2007), transportation research(Ma et al., 2017) (Bhat & Eluru, 2009) , medicine (Kuss et al., 2014) (Lapuyade-Lahorgue et al., 2017) (Gomes et al., 2019) , engineering (Kilgore & Thompson, 2011), biology (Konigorski et al., 2014) (Dokuzoğlu & Purutçuoğlu, 2017), climate research (Schölzel & Friederichs, 2008) (Oppenheimer et al., 2016), economics (Oh & Patton, 2018) (Kole et al., 2007) (De Lira Salvatierra & Patton, 2015). Sklar (Abe Sklar, 1973) showed that a copula exists for any multivariate distribution, such that the joibt distribution equals the copula applied to the marginal.
Definition 1 (
Sklar’s theorem 1): For any bivariate distribution function F with marginal
and
, there exists a copula such that
If F is continuous, then the copula C is unique. Otherwise, it is uniquely determined on the . The converse is also true, for any copula and univariate distribution functions and , the function is a bivariate distribution function with margins and .
A bivariate distribution function
with marginal
and
is said to be generated by an Archimedean copula, if it can be written in the form of equation (1)
The generator should satisfy the following conditions (sufficient conditions):
- (1)
is continuous, strictly decreasing and convex function mapping onto , in other words , ( decreasing) and ( convex)
- (2)
- (3)
Definition 2 : a bivariate copula is a non-decreasing and a right-continuous bivariate function mapping
which satisfies the following conditions
- (1)
( Grounded)
- (2)
and ( uniform marginal)
- (3)
(2- increasing) . This last property is equal to the joint PDF of the copula.
So the copula links the marginal distributions through the dependency parameter.
The singular part of the copula is defined by (Genest & Mackay, 1986) in the following Theorem 2: the ratio between the generator and the first derivative of the generator evaluated at zero does not equal zero.
Theorem: the distribution generated by has a singular component iff in that case, with probability =
This paper is structured in 4 sections. Section 1 discusses the first copula, its derivation, properties and figures for its PDF and CDF. Section 2 explains the second copula with the same items. Section 3 explore the third copula in the same manner. Section 4 illustrate the discussion and conclusion with future work
Section 1: (First Copula)
Let’s say the inverse generator is
The generator should fulfill the sufficient conditions:
- (1)
- (2)
- (3)
-
This ensures that the generator is a decreasing function.
- (4)
This ensures that the generator is convex for .
for bivariate distribution:
For a copula to be valid copula it should fulfill the boundary condition, the marginal uniformity and 2- increasing condition. And this is equivalent for the inverse generator the necessary conditions to be fulfilled are the following:
Proposition 1 : The boundary conditions and marginal uniformity:
, boundary condition.
, the marginal uniformity.
Proposition 2 : for copula to be valid it should be 2-increasing, in other words,
Although the second derivative is positive at the following intervals: so alpha can be defined on the interval .
Proposition 3 : this copula is absolutely continuous copula and it has no singular part.
Proof: to test for singularity:
As long as this limit is zero at u=0 so it has no singular part and it is absolutely continuous copula.
Proposition 4 : Kendall tau for this copula is τ = 1 − α
If indicating independence.
Although the second derivative is positive at the following intervals: considering that tau is defined within the interval so alpha can be defined on the interval . In this interval the Kendall tau can show positive dependency and independency.
Preposition 5 : This copula is product copula at alpha =1, in other word, if alpha=1, , this is the product copula indicating independency.
The following figures illustrates the Copula surface with different alpha values
Figure 1.
shows the generator function (decreasing and convex) at alpha = 0.6.
Figure 1.
shows the generator function (decreasing and convex) at alpha = 0.6.
Figure 2.
shows the joint PDF copula (copula density) at alpha=0.1.
Figure 2.
shows the joint PDF copula (copula density) at alpha=0.1.
Figure 3.
shows the joint CDF copula (Copula) at alpha =0.1.
Figure 3.
shows the joint CDF copula (Copula) at alpha =0.1.
Figure 4.
shows the joint pdf copula ( copula density) at alpha=0.6.
Figure 4.
shows the joint pdf copula ( copula density) at alpha=0.6.
Figure 5.
shows the joint cdf copula ( Copula) at alpha=0.6.
Figure 5.
shows the joint cdf copula ( Copula) at alpha=0.6.
Figure 6.
shows the joint pdf copula ( copula density) at alpha =1.
Figure 6.
shows the joint pdf copula ( copula density) at alpha =1.
Figure 7.
shows joint cdf Copula ( Copula) at alpha=1.
Figure 7.
shows joint cdf Copula ( Copula) at alpha=1.
The limitation for this copula is that it does not cover the negative dependency between variables.
Section 2: (Second Copula)
Let’s say the inverse generator is
The generator should fulfill the sufficient conditions:
- (1)
- (2)
- (3)
-
This ensures that the generator is a decreasing function.
- (4)
This ensures that the generator is convex at
for bivariate distribution:
For a copula to be a valid copula it should fulfill the boundary condition, the marginal uniformity and 2- increasing condition. And this is equivalent for necessary conditions of the inverse generator to be fulfilled which are the following:
Proposition 6 : The boundary conditions and marginal uniformity:
, boundary condition.
, the marginal uniformity.
Proposition 7 : for copula to be valid it should be 2-increasing, in other words
The second derivative is positive at
Preposition 8 : this copula is absolutely continuous copula and it has no singular part.
Proof: to test for singularity:
As long as this limit is zero at u=0 so the copula has no singular part and it is absolutely continuous copula.
Preposition 9 : Kendall tau for this copula is τ = 1 − 2α2
If indicating positive dependency.
If indicating negative dependency.
To be product copula indicating independency, the alpha parameter should be square root of 2 at which the copula is invalid. So this copula cannot represent independency.
The following figures illustrates the Copula surface with different alpha values
Figure 8.
shows the generator function (decreasing and convex) at alpha = 1.
Figure 8.
shows the generator function (decreasing and convex) at alpha = 1.
Figure 9.
shows the joint CDF copula (copula) at alpha=0.1.
Figure 9.
shows the joint CDF copula (copula) at alpha=0.1.
Figure 10.
shows the joint PDF copula (copula density) at alpha =0.1.
Figure 10.
shows the joint PDF copula (copula density) at alpha =0.1.
Figure 11.
shows the joint CDF copula (copula) at alpha =0.4.
Figure 11.
shows the joint CDF copula (copula) at alpha =0.4.
Figure 12.
shows the joint PDF copula (copula density) at alpha =0.4.
Figure 12.
shows the joint PDF copula (copula density) at alpha =0.4.
Figure 13.
shows the joint CDF copula (copula) at alpha =0.6.
Figure 13.
shows the joint CDF copula (copula) at alpha =0.6.
Figure 14.
shows the joint PDF copula (copula density) at alpha =0.6.
Figure 14.
shows the joint PDF copula (copula density) at alpha =0.6.
Figure 15.
shows the joint CDF copula (copula) at alpha =0.888.
Figure 15.
shows the joint CDF copula (copula) at alpha =0.888.
Figure 16.
shows the joint PDF copula (copula density) at alpha =0.888.
Figure 16.
shows the joint PDF copula (copula density) at alpha =0.888.
Figure 17.
shows the joint CDF copula (copula) at alpha =0.988.
Figure 17.
shows the joint CDF copula (copula) at alpha =0.988.
Figure 18.
shows the joint PDF copula (copula density) at alpha =0.988.
Figure 18.
shows the joint PDF copula (copula density) at alpha =0.988.
Figure 19.
shows the joint CDF copula (copula) at alpha =1.
Figure 19.
shows the joint CDF copula (copula) at alpha =1.
Figure 20.
shows the joint PDF copula (copula density) at alpha =1.
Figure 20.
shows the joint PDF copula (copula density) at alpha =1.
The limitation of this copula it does not cover independency.
Section 3: (Third Copula)
This copula depends on the frailty method for generating the copula. Assuming two individuals have survival time distributed as exponential
and
.They have exponential baseline hazard function with
. The survival function is shown in equation (2) and the hazard function in equation ()
The Cox proportional hazard model uses the hazard function as shown in equation (3)
where the Z are the explanatory variables in survival analysis and b(t) is the baseline hazard function. And B is the vector of regression coefficients. It is proportional because all information is contained in the multiplicative factor
, this is called the frailty model when some explanatory variables (Z) and hence the factor
are unobserved. And the factor
is called the frailty parameter. Integrating and exponentiating the negative hazard
for 2 explanatory variables for example
, to sum up
The marginal distribution for a single life time T is obtained by taking expectation over the potential values of
that is
, assuming the
is the PDF of the frailty variable so
Let y be a random variable distributed as Median Based Unit Rayleigh (BMUR)in equation (4) previously introduced by Iman Attia (Iman M. Attia, 2024), transform this variable to be defined on the interval from zero to infinity as shown below:
Let
Now this w is the frailty variable
which is distributed as unbounded MBUR defined on the interval
.
We can think of this expectation as a function of the frailty variable, it is averaging the probability to survive beyond specific time given that frailty. It is the expectation of a function of the frailty not the variable frailty itself so this function which is the baseline survival function should be multiplied by the PDF of the frailty variable. This way the frailty variable is integrated out. So this expectation is considered as the joint CDF of the two random variables T1 and T2. How to get this time? By inverting this expectation. Let us call this expectation u then solving second degree polynomial to get its root.
This (t) is the generator and (u) is the inverse generator. So the time (t) is a function of (u)
The generator should fulfill the sufficient conditions:
- (1)
- (2)
- (3)
-
This ensures that the generator is a decreasing function.
- (4)
This ensures that the generator is convex at
for bivariate distribution:
For a copula to be a valid copula it should fulfill the boundary condition, the marginal uniformity and 2- increasing conditions. And this is equivalent for necessary conditions of the inverse generator to be fulfilled which are the following:
Proposition 10 : The boundary conditions and marginal uniformity:
, boundary condition.
, the marginal uniformity.
The same is true for v , if v=1 so
Proposition 11 : for copula to be valid it should be 2-increasing in other words
Proposition 12 : this copula is absolutely continuous copula and it has no singular part.
Proof: to test for singularity:
As long as this limit is zero at u=0 so the copula has no singular part and it is absolutely continuous copula.
Proposition 13 : Kendall tau for this copula is τ = 1 − 2α2
The following figures show the generator and the CDF and PDF of this third copula
Figure 21.
shows the generator function (decreasing and convex) at alpha = 0.1.
Figure 21.
shows the generator function (decreasing and convex) at alpha = 0.1.
Figure 22.
shows the joint CDF copula (copula density) at alpha =0.1.
Figure 22.
shows the joint CDF copula (copula density) at alpha =0.1.
Figure 23.
shows the joint PDF copula (copula density) at alpha =0.1.
Figure 23.
shows the joint PDF copula (copula density) at alpha =0.1.
Figure 23.
shows the joint CDF copula (copula) at alpha =0.6.
Figure 23.
shows the joint CDF copula (copula) at alpha =0.6.
Figure 24.
shows the joint PDF copula (copula density) at alpha =0.6.
Figure 24.
shows the joint PDF copula (copula density) at alpha =0.6.
Figure 25.
shows the joint CDF copula (copula) at alpha =1.
Figure 25.
shows the joint CDF copula (copula) at alpha =1.
Figure 26.
shows the joint PDF copula (copula density) at alpha =1.
Figure 26.
shows the joint PDF copula (copula density) at alpha =1.
Figure 27.
shows the joint CDF copula (copula) at alpha =10.
Figure 27.
shows the joint CDF copula (copula) at alpha =10.
Figure 28.
shows the joint PDF copula (copula density) at alpha =1.
Figure 28.
shows the joint PDF copula (copula density) at alpha =1.
To sum up for this copula, the generator and inverse generator fulfill the conditions. The copula density fulfills the criteria but the values of the joint CDF and DPF are the same for all values of the dependency parameter. Kendal tau is fixed at 0.32 and it does not depend on the dependency parameter.
Section 4: Conclusions and Future Work
Although the generator and inverse generator fulfilled the criteria as well as the copula density but for each copula discussed previously, there were limitations. These limitation may be solved by mixture of copulas or better re-parameterization of these copula once again.
Ethics Approval and Consent to Participate
Not applicable.
Consent for publication
Not applicable
Availability of Data and Material
Not applicable. Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.
Competing Interests
The author declares no competing interests of any type.
Funding
No funding resource. No funding roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript are declared.
Authors’ Contribution
AI carried the conceptualization by formulating the goals, aims of the research article, formal analysis by applying the statistical, mathematical and computational techniques to synthesize and analyze the hypothetical data, carried the methodology by creating the model, software programming and implementation, supervision, writing, drafting, editing, preparation, and creation of the presenting work.
References
- Abe Sklar. Random Variables, Joint distribution functions, and copulas. Kybernetika 1973, 9, 449–460. [Google Scholar]
- Bhat, C. R. , & Eluru, N. A copula-based approach to accommodate residential self-selection effects in travel behavior modeling. Transportation Research Part B: Methodological 2009, 43, 749–765. [Google Scholar] [CrossRef]
- De Lira Salvatierra, I. , & Patton, A. J. Dynamic copula models and high frequency data. Journal of Empirical Finance 2015, 30, 120–135. [Google Scholar] [CrossRef]
- Dokuzoğlu, D. , & Purutçuoğlu, V. Comprehensive Analyses of Gaussian Graphical Model under Different Biological Networks. Acta Physica Polonica A 2017, 132, 1106–1111. [Google Scholar] [CrossRef]
- Genest, C. , & Mackay, J. The Joy of Copulas: Bivariate Distributions with Uniform Marginals. The American Statistician 1986, 40, 280–283. [Google Scholar] [CrossRef]
- Gomes, M. , Radice, R., Camarena Brenes, J., & Marra, G. Copula selection models for non-Gaussian outcomes that are missing not at random. Statistics in Medicine 2019, 38, 480–496. [Google Scholar] [CrossRef] [PubMed]
- Iman, M. Attia. A Novel Unit Distribution Named asMedian Based Unit Rayleigh (MBUR):Properties and Estimations. Preprints.Org, 7 October. [CrossRef]
- Kilgore, R. T. , & Thompson, D. B. Estimating Joint Flow Probabilities at Stream Confluences by Using Copulas. Transportation Research Record: Journal of the Transportation Research Board 2011, 2262, 200–206. [Google Scholar] [CrossRef]
- Kole, E. , Koedijk, K., & Verbeek, M. Selecting copulas for risk management. Journal of Banking & Finance 2007, 31, 2405–2423. [Google Scholar] [CrossRef]
- Konigorski, S. , Yilmaz, Y. E., & Bull, S. B. Bivariate genetic association analysis of systolic and diastolic blood pressure by copula models. BMC Proceedings 2014, 8, S72. [Google Scholar] [CrossRef]
- Kuss, O. , Hoyer, A., & Solms, A. Meta-analysis for diagnostic accuracy studies: A new statistical model using beta-binomial distributions and bivariate copulas. Statistics in Medicine 2014, 33, 17–30. [Google Scholar] [CrossRef] [PubMed]
- Lapuyade-Lahorgue, J. , Xue, J.-H., & Ruan, S. Segmenting Multi-Source Images Using Hidden Markov Fields With Copula-Based Multivariate Statistical Distributions. IEEE Transactions on Image Processing 2017, 26, 3187–3195. [Google Scholar] [CrossRef] [PubMed]
- Liu, Z. , Guo, S., Xiong, L., & Xu, C.-Y. Hydrological uncertainty processor based on a copula function. Hydrological Sciences Journal 2018, 63, 74–86. [Google Scholar] [CrossRef]
- Ma, X. , Luan, S., Du, B., & Yu, B. Spatial Copula Model for Imputing Traffic Flow Data from Remote Microwave Sensors. Sensors 2017, 17, 2160. [Google Scholar] [CrossRef]
- Oh, D. H. , & Patton, A. J. Time-Varying Systemic Risk: Evidence From a Dynamic Copula Model of CDS Spreads. Journal of Business & Economic Statistics 2018, 36, 181–195. [Google Scholar] [CrossRef]
- Oppenheimer, M. , Little, C. M., & Cooke, R. M. Expert judgement and uncertainty quantification for climate change. Nature Climate Change 2016, 6, 445–451. [Google Scholar] [CrossRef]
- Salvadori, G., Michele, C. D., Kottegoda, N. T., & Rosso, R. (2007). Extremes in Nature: An Approach Using Copulas (Vol. 56). Springer Netherlands. [CrossRef]
- Schölzel, C. , & Friederichs, P. Multivariate non-normally distributed random variables in climate research – introduction to the copula approach. Nonlinear Processes in Geophysics 2008, 15, 761–772. [Google Scholar] [CrossRef]
- Valle, D. , & Kaplan, D. Quantifying the impacts of dams on riverine hydrology under non-stationary conditions using incomplete data and Gaussian copula models. Science of The Total Environment 2019, 677, 599–611. [Google Scholar] [CrossRef] [PubMed]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).