Next Article in Journal
A Numerical Simulation of Electrical Resistivity of Fiber-Reinforced Composites, Part 1: Brittle Cementitious Concrete
Previous Article in Journal
On the Design of Composite Patch Repair for Strengthening of Marine Plates Subjected to Compressive Loads
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Study of Some Angle Parameter Trigonometric Copulas

by
Christophe Chesneau
Department of Mathematics, University of Caen-Normandie, 14032 Caen, France
Modelling 2022, 3(1), 140-163; https://doi.org/10.3390/modelling3010010
Submission received: 8 February 2022 / Revised: 5 March 2022 / Accepted: 7 March 2022 / Published: 10 March 2022

Abstract

:
Copulas are important probabilistic tools to model and interpret the correlations of measures involved in real or experimental phenomena. The versatility of these phenomena implies the need for diverse copulas. In this article, we describe and investigate theoretically new two-dimensional copulas based on trigonometric functions modulated by a tuning angle parameter. The independence copula is, thus, extended in an original manner. Conceptually, the proposed trigonometric copulas are ideal for modeling correlations into periodic, circular, or seasonal phenomena. We examine their qualities, such as various symmetry properties, quadrant dependence properties, possible Archimedean nature, copula ordering, tail dependences, diverse correlations (medial, Spearman, and Kendall), and two-dimensional distribution generation. The proposed copulas are fleshed out in terms of data generation and inference. The theoretical findings are supplemented by some graphical and numerical work. The main results are proved using two-dimensional inequality techniques that can be used for other copula purposes.

1. Introduction

Multidimensional functions called copulas are important in modeling multivariate random variables and understanding their dependence structures. They find their origin and prime property in the Sklar theorem (see [1]). Mathematically, for any integer n, a n-dimensional copula can be defined as a cumulative distribution function defined on [ 0 , 1 ] n with standard uniform marginal distributions. For the purposes of this article, the definition of a two-dimensional copula in the absolutely continuous case is given below.
Definition 1.
In the absolutely two-dimensional continuous case, the function C : [ 0 , 1 ] 2 [ 0 , 1 ] is a two-dimensional copula if and only if
  • C ( x , 0 ) = C ( 0 , y ) = 0 for any ( x , y ) [ 0 , 1 ] 2 ;
  • C ( x , 1 ) = x and C ( 1 , y ) = y for any ( x , y ) [ 0 , 1 ] 2 ;
  • the two-increasing property holds: 2 C ( x , y ) / ( x y ) 0 for any ( x , y ) [ 0 , 1 ] 2 .
This definition has a certain flexibility, but the main critical point remains the validity of the two-increasing property. On the practical side, copulas are crucial for estimating and interpreting the correlations of measures involved in real or experimental phenomena. Due to the plurality of possible phenomena of interest, a large number of copulas with various attributes are necessary. This has motivated researchers over the years to introduce valuable copulas, and apply them in concrete, real-life scenarios. The Frank, Gumbel–Hougaard, Ali–Mikhail–Haq, Joe, Farlie–Gumbel–Morgenstern, Clayton, Plackett, Raftery, elliptical, Fréchet, Galambos, and Marshall–Olkin are some of the classical copulas. Their definitions are based on motivated transformations of power-polynomial–exponential–logarithmic functions. They have specific qualities that make them of great interest in the context of random dependence modeling. Further details on these copulas can be found in the books of [1,2,3,4].
Recently, a lot of attention has been paid to copulas defined (partially or not) with trigonometric functions. The inclusion of trigonometric functions in this context confers on the copula some oscillating features that are appropriate to model the correlations into phenomena of periodic, circular, or seasonal nature. In particular, they are ideal for analyzing correlations involved in movement data, circular data, and environmental data. The theory of the classical trigonometric copulas can be found in [5,6,7,8,9,10,11]. For practice, we refer to [12,13,14,15,16]. Furthermore, the R package named Cylcop, recently developed by [17], gives the trigonometric (and circular) copulas a new dimension of applicability.
In this article, we present and study trigonometric copulas depending on a tuning angle parameter. They are derived from the very general family of copulas (not especially trigonometric) introduced in [18], but with the importance of the angle parameter in mind. To be more specific, the proposed copulas are defined by the following form:
C ( x , y ) = x y f ( θ ( 1 x ) ( 1 y ) ) ,
where f ( x ) is a simple one-dimensional function involving a trigonometric function, and θ is the angle parameter that only modulates this trigonometric function. Despite the potential for modeling correlations into periodic, circular, or seasonal phenomena, this area of research appears to have received little attention; none of the references [5,6,7,8,9,10,11] consider such an angle parameter approach. Thus, we describe the most intuitive of such angle parameter copulas, and study them on the theoretical side with a maximum of details. In particular, we emphasize the optimal set of values for θ such that the corresponding copula remains valid. We determine the expressions of the corresponding copula density, survival copula, and survival copula density. Some mixed trigonometric copulas are also given. Then, we examine a maximum of theoretical properties of the proposed angle parameter copulas, including symmetry properties, quadrant dependence properties, copula ordering, various expansions, tail dependences, medial correlation, Spearman correlation, Kendall correlation, and two-dimensional distribution generation. Data generation and inference from the proposed copulas are sketched. The theory is illustrated by means of graphics. The proofs of the main results are based on some two-dimensional inequality techniques that can be of independent interest.
The following is the outline for the rest of the article: A cosine angle parameter copula is presented in Section 2, along with its main theoretical properties. In Section 3, an analogue sine copula is introduced and studied. Section 4 contains the conclusions and perspectives.

2. Cosine Angle Parameter Copula

This section is devoted to a simple cosine copula with an angle parameter.

2.1. Definition and Graphics

The following proposition presents the considered cosine copula.
Proposition 1. 
The function C : [ 0 , 1 ] 2 [ 0 , 1 ] defined by
C ( x , y ) = x y cos θ ( 1 x ) ( 1 y ) ,
with θ [ 0 , π / 2 ] is a valid copula.
Proof. 
Let us prove the main points defining an absolutely continuous two-dimensional copula, as recalled in Definition 1.
  • For any x [ 0 , 1 ] , we have C ( x , 0 ) = x × 0 × cos θ ( 1 x ) ( 1 0 ) = 0 , and, for any y [ 0 , 1 ] , C ( 0 , y ) = 0 × y cos θ ( 1 0 ) ( 1 y ) = 0 .
  • For any x [ 0 , 1 ] , we have C ( x , 1 ) = x × 1 × cos θ ( 1 x ) × 0 = x and, similarly, for any y [ 0 , 1 ] , C ( 1 , y ) = 1 × y × cos θ × 0 × ( 1 y ) = y .
  • For any ( x , y ) [ 0 , 1 ] 2 , using standard derivation techniques, simplifications and factorizations, we have
    2 x y C ( x , y ) = 1 θ 2 x y ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y ) 3 θ x y sin θ ( 1 x ) ( 1 y ) + θ ( x + y ) sin θ ( 1 x ) ( 1 y ) .
Let us prove the two-increasing property: 2 C ( x , y ) / ( x y ) 0 .
Since θ [ 0 , π / 2 ] , we have θ ( 1 x ) ( 1 y ) [ 0 , π / 2 ] . It follows from the inequality sin ( a ) a cos ( a ) , a [ 0 , π / 2 ] applied with a = θ ( 1 x ) ( 1 y ) that
2 x y C ( x , y ) 1 θ 2 x y ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y ) 3 θ x y sin θ ( 1 x ) ( 1 y ) + θ 2 ( x + y ) ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y ) = 1 + θ 2 ( x + y x y ) ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y ) 3 θ x y sin θ ( 1 x ) ( 1 y ) A 1 + A 2 ,
where
A 1 = θ 2 ( x + y x y ) ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y )
and
A 2 = cos θ ( 1 x ) ( 1 y ) 3 θ x y sin θ ( 1 x ) ( 1 y ) .
Let us prove that A 1 0 and A 2 0 .
For A 1 , we can remark that x + y x y = 1 ( 1 x ) ( 1 y ) [ 0 , 1 ] . So A 1 0 as a product of positive terms.
For A 2 , it follows from the inequality sin ( a ) a , a [ 0 , π / 2 ] applied with a = θ ( 1 x ) ( 1 y ) that
A 2 ψ ( x , y ; θ ) ,
where
ψ ( x , y ; θ ) = cos θ ( 1 x ) ( 1 y ) 3 θ 2 x y ( 1 x ) ( 1 y ) .
We have
θ ψ ( x , y ; θ ) = ( 1 x ) ( 1 y ) 6 θ x y + sin θ ( 1 x ) ( 1 y ) 0 ,
implying that ψ ( x , y ; θ ) is a decreasing function with respect to θ . Therefore, for any θ [ 0 , π / 2 ] , we have
ψ ( x , y ; θ ) ψ ( x , y ; π / 2 ) ,
with
ψ ( x , y ; π / 2 ) = cos π 2 ( 1 x ) ( 1 y ) 3 π 2 4 x y ( 1 x ) ( 1 y ) .
By applying the Kober inequality: cos ( a ) 1 ( 2 / π ) a for a [ 0 , π / 2 ] (see [19,20]), with a = ( π / 2 ) ( 1 x ) ( 1 y ) [ 0 , π / 2 ] , we get
ψ ( x , y ; π / 2 ) ϕ ( x , y ) ,
where
ϕ ( x , y ) = x + y x y 3 π 2 4 x y ( 1 x ) ( 1 y ) .
Let us prove that ϕ ( x , y ) 0 for any ( x , y ) [ 0 , 1 ] 2 . First, ϕ ( x , y ) is a continuous function for ( x , y ) [ 0 , 1 ] 2 , and [ 0 , 1 ] 2 is a compact set, so the maximum and minimum of ϕ ( x , y ) on this set are attained. Let us now perform a critical points analysis. We have
x ϕ ( x , y ) = 1 4 ( 1 y ) ( 3 π 2 ( 2 x 1 ) y + 4 ) y ϕ ( x , y ) = 1 4 ( 1 x ) ( 3 π 2 ( 2 y 1 ) x + 4 ) .
Therefore, we have ϕ ( x , y ) / ( x ) = 0 and ϕ ( x , y ) / ( y ) = 0 if and only if
  • y = 1 and x = 1 ; or
  • y = 1 and 3 π 2 ( 2 y 1 ) x + 4 = 0 which implies that x = 4 / ( 3 π 2 ) [ 0 , 1 ] , so this case is excluded; or
  • x = 1 and 3 π 2 ( 2 x 1 ) y + 4 = 0 which implies that y = 4 / ( 3 π 2 ) [ 0 , 1 ] , so this case is excluded; or
  • 3 π 2 ( 2 x 1 ) y + 4 = 0 and 3 π 2 ( 2 x 1 ) y + 4 = 0 , which implies that x = y and 3 π 2 ( 2 x 1 ) x + 4 = 0 , which is a polynomial of degree 2 with the discriminant equal to Δ = 1 32 / ( 3 π 2 ) < 0 . Therefore, there is no (real) solution.
As a result, there is only one critical point for ϕ ( x , y ) , and it is ( 1 , 1 ) . This point clearly gives a maximum value for ϕ ( x , y ) since ϕ ( 1 , 1 ) = 1 x + y x y ϕ ( x , y ) . Since there is no other critical point, ϕ ( x , y ) has a “two-dimensional decreasing pattern”, and the minimum value of ϕ ( x , y ) is attained on point(s) on the borders of [ 0 , 1 ] 2 . Therefore, for any ( x , y ) [ 0 , 1 ] 2 , we have
ϕ ( x , y ) inf ( x , y ) [ 0 , 1 ] 2 ϕ ( x , y ) = min [ inf x [ 0 , 1 ] ϕ ( x , 0 ) , inf y [ 0 , 1 ] ϕ ( 0 , y ) , inf x [ 0 , 1 ] ϕ ( x , 1 ) , inf y [ 0 , 1 ] ϕ ( 1 , y ) ] = min ( inf x [ 0 , 1 ] x , inf y [ 0 , 1 ] y , inf x [ 0 , 1 ] 1 , inf y [ 0 , 1 ] 1 ) = 0 .
It follows from Equations (2)–(5) that A 2 0 . The two-increasing property is proved.
As a result, C ( x , y ) is a valid two-dimensional copula. This ends the proof of Proposition 1. □
Remark 1.
For the more restrictive case θ [ 0 , π / 4 ] , we can prove the two-increasing property in a more simple and direct manner. We have
2 x y C ( x , y ) = 1 θ 2 x y ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y ) θ ( 3 x y x y ) sin θ ( 1 x ) ( 1 y ) .
Since θ [ 0 , π / 4 ] [ 0 , 1 ] , we have θ 2 x y ( 1 x ) ( 1 y ) 1 , and by applying the inequality cos ( a ) sin ( a ) for any a [ 0 , π / 4 ] with a = θ ( 1 x ) ( 1 y ) [ 0 , π / 4 ] , we get
2 x y C ( x , y ) 1 θ 2 x y ( 1 x ) ( 1 y ) sin θ ( 1 x ) ( 1 y ) θ ( 3 x y x y ) sin θ ( 1 x ) ( 1 y ) .
Since θ [ 0 , 1 ] , we have θ 2 x y ( 1 x ) ( 1 y ) ( 1 x ) ( 1 y ) and θ ( 3 x y x y ) = θ [ x y + x ( y 1 ) + y ( x 1 ) ] x y . Thus
2 x y C ( x , y ) 1 ( 1 x ) ( 1 y ) sin θ ( 1 x ) ( 1 y ) x y sin θ ( 1 x ) ( 1 y ) = [ x ( 1 y ) + y ( 1 x ) ] sin θ ( 1 x ) ( 1 y ) 0 .
The two-increasing property is proved. This ends this alternative proof for θ [ 0 , π / 4 ] .
For the purposes of this study, we call the copula C ( x , y ) in Equation (1) as the cos-copula. It has the feature of having a tuning angle parameter θ that modulates its correlation features. For θ = 0 , we clearly have C ( x , y ) = x y ; the cos-copula is reduced to the independence copula. Furthermore, we can remark that the cos-copula is still valid for θ [ π / 2 , 0 ] because of the evenness of the cosine function. We can, thus, define the cos-copula with θ [ π / 2 , π / 2 ] without loss of generality.
To end the presentation, Figure 1 represents the two-dimensional plot of the cos-copula for selected values of θ .
From Figure 1, we see that the parameter θ skews the triangular shape of the cos-copula.
The next parts are devoted to the functions related to the cos-copula, as well as its main properties.

2.2. Related Functions

The density associated with the cos-copula is the function c : [ 0 , 1 ] 2 [ 0 , + ) given by
c ( x , y ) = 2 x y C ( x , y ) = 1 θ 2 x y ( 1 x ) ( 1 y ) cos θ ( 1 x ) ( 1 y ) θ ( 3 x y x y ) sin θ ( 1 x ) ( 1 y ) .
It is worth noting that c ( 0 , 0 ) = cos ( θ ) which can be negative for some values of θ . In particular, C ( x , y ) is not a copula for θ ( π / 2 , π ] . In this sense, if we restrict our attention to the interval [ 0 , π ] , [ 0 , π / 2 ] is the optimal set of values for θ .
As usual, the copula density can be involved in various moment-type measures, and estimation methods (see, for instance, Refs. [13,17]).
Figure 2 represents the two-dimensional plot of the cos-copula density for selected values of θ .
Thus, we see how the parameter θ affects the overall shape of the cos-copula density.
As a last important function, the survival cos-copula is the function C ^ : [ 0 , 1 ] 2 [ 0 , 1 ] defined by
C ^ ( x , y ) = x + y 1 + C ( 1 x , 1 y ) = x + y 1 + ( 1 x ) ( 1 y ) cos ( θ x y ) = x y ( 1 x ) ( 1 y ) 1 cos ( θ x y ) .
It defines a valid copula, which is also a new trigonometric copula with an angle parameter represented by θ .
Figure 3 represents the two-dimensional plot of the survival cos-copula for selected values of θ .
The density associated with the survival cos-copula is given by
c ^ ( x , y ) = 2 x y C ^ ( x , y ) = 1 θ 2 x y ( 1 x ) ( 1 y ) cos ( θ x y ) θ ( 3 x y 2 x 2 y + 1 ) sin ( θ x y ) .
To have an idea of its shapes, Figure 4 shows the two-dimensional plot of the survival cos-copula density for selected values of θ .
Because any convex linear combination of copulas is a copula (see [1]), the cos-copula can be used in a variety of new mixed copulas. For instance, we can consider the following mixed cos-copulas:
  • Mixed copula 1: For any angle parameters θ 1 [ 0 , π / 2 ] and θ 2 [ 0 , π / 2 ] and λ [ 0 , 1 ] , by setting C ( x , y ) = C ( x , y ; θ ) , a possible mixed copula is given as
    C ( x , y ) = λ C ( x , y ; θ 1 ) + ( 1 λ ) C ( x , y ; θ 2 ) = λ x y cos θ 1 ( 1 x ) ( 1 y ) + ( 1 λ ) x y cos θ 2 ( 1 x ) ( 1 y ) .
  • Mixed copula 2: A second example is
    C ( x , y ) = 1 2 C ( x , y ) + y C ( 1 x , y ) = 1 2 x y cos θ ( 1 x ) ( 1 y ) + y ( 1 x ) y cos θ x ( 1 y ) .
They are also new angle parameter trigonometric copulas by construction.

2.3. Properties

We now list the important properties of the cos-copula C ( x , y ) as specified in Equation (1).
  • As already mentioned before:
    For θ = 0 , it is clear that C ( x , y ) = x y . Therefore, the cos-copula is reduced to the independence copula.
    If we restrict our attention to the interval [ 0 , π ] , the set [ 0 , π / 2 ] is the optimal set of values for θ for validating C ( x , y ) as a copula.
  • For any θ [ 0 , π / 2 ] , we have C ( x , y ) x y . Hence, the cos-copula satisfies the negative quadrant dependence property (see [21]).
  • The cos-copula is symmetric since C ( x , y ) = C ( y , x ) for any ( x , y ) [ 0 , 1 ] 2 .
  • The cos-copula can be expressed under various analytical forms. Two of them are given below:
    In terms of simple cosine-sine functions, we can write
    C ( x , y ) = x y [ sin ( θ ) sin ( θ x ) sin ( θ y ) sin ( θ x y ) + cos ( θ ) cos ( θ x ) cos ( θ y ) cos ( θ x y ) + sin ( θ ) sin ( θ x ) cos ( θ y ) cos ( θ x y ) + sin ( θ ) cos ( θ x ) sin ( θ y ) cos ( θ x y ) cos ( θ ) sin ( θ x ) sin ( θ y ) cos ( θ x y ) sin ( θ ) cos ( θ x ) cos ( θ y ) sin ( θ x y ) + cos ( θ ) sin ( θ x ) c o s ( θ y ) sin ( θ x y ) + cos ( θ ) cos ( θ x ) sin ( θ y ) sin ( θ x y ) ]
    We thus see the intrinsic analytical complexity into the cos-copula.
    In terms of power series, by using the cosine series expansion and binomial formula, we get
    C ( x , y ) = x y k = 0 + ( 1 ) k ( 2 k ) ! θ 2 k ( 1 x ) 2 k ( 1 y ) 2 k = k = 0 + = 0 2 k m = 0 2 k ( 1 ) k + + m ( 2 k ) ! θ 2 k 2 k 2 k m x + 1 y m + 1 .
    In particular, upon differentiation with respect to x and y on the interior of the domain of convergence, one has
    c ( x , y ) = k = 0 + = 0 2 k m = 0 2 k ( 1 ) k + + m ( 2 k ) ! θ 2 k 2 k 2 k m ( + 1 ) ( m + 1 ) x y m .
    This expansion can be used in a variety of mathematical applications, such as determining various moment-type measurements.
  • By arbitrary taking θ = π / 2 , we notice that
    C 1 4 , C 1 2 , 1 3 = 0.01925131 0.02047885 = C C 1 4 , 1 2 , 1 3 .
    As a result, the cos-copula is not Archimedean (see [1]). In other words, there is no generator function ψ : [ 0 , 1 ] [ 0 , + ) such that C ( x , y ) = ψ [ 1 ] [ ψ ( x ) + ψ ( y ) ] , where ψ [ 1 ] ( x ) denotes the pseudo-inverse of ψ ( x ) .
  • The cos-copula is not radially symmetric since there clearly exists ( x , y ) such that C ^ ( x , y ) C ( x , y ) .
  • As any copula, the Fréchet–Hoeffding bounds can be expressed as follows: For any ( x , y ) [ 0 , 1 ] 2 , we have max ( x + y 1 , 0 ) C ( x , y ) min ( x , y ) .
  • Thanks to the Kober inequality, the following inequality holds:
    C ( x , y ) C x , y ; 2 π θ ,
    where C ( x , y ; λ ) = x y + λ x y ( 1 x ) ( 1 y ) with λ [ 1 , 1 ] refers to the Farlie–Gumbel–Morgenstern (FGM) copula (see [1]). Since ( 2 / π ) θ [ 1 , 1 ] , the cos-copula and FGM copula are involved in a complete copula ordering.
  • For any θ [ 0 , π / 2 ] , the two following results are obtained:
    λ L = lim x 0 C ( x , x ) x = lim x 0 x cos θ ( 1 x ) 2 = 0
    and
    λ U = lim x 1 1 2 x + C ( x , x ) 1 x = lim x 1 1 2 x + x 2 cos θ ( 1 x ) 2 1 x = lim x 1 1 2 x + x 2 1 x = lim x 1 ( 1 x ) = 0 .
    Hence, the cos-copula has no tail dependence (see [1]).
  • The medial correlation of the cos-copula is defined by
    M = 4 C 1 2 , 1 2 1 = cos θ 4 1 .
    It is clearly a decreasing and negative function with respect to θ for θ [ 0 , π / 2 ] , with M = 0 for θ = 0 and M = 0.07612047 for θ = π / 2 . Figure 5 represents the medial correlation for θ [ 0 , 1.5 ] .
    Thus, the cos-copula has a weak medial correlation with M [ 0.076 , 0 ] .
  • A useful dependence measure based on copula is the Spearman rho (see [1]). The Spearman rho of the cos-copula, as an example of copula, is defined by
    ρ = 12 0 1 0 1 C ( x , y ) d x d y 3 .
Based on well-known mathematical techniques, the following proposition provides a mathematical expression for this measure.
Proposition 2.
The Spearman rho of the cos-copula can be expressed as
ρ = 0 , for θ = 0 12 θ 2 Ci ( θ ) + θ Si ( θ ) log ( θ ) + cos ( θ ) γ 1 3 , for θ ( 0 , π / 2 ] ,
where Ci ( x ) denotes the cosine integral defined by Ci ( x ) = x [ cos ( t ) / t ] d t , Si ( x ) denotes the sine integral referred by Si ( x ) = 0 x [ sin ( t ) / t ] d t , and γ 0.5772 is the Euler–Mascheroni constant.
Proof. 
  • For θ = 0 , we have C ( x , y ) = x y and, by Equation (7), it is immediate that ρ = 0 .
  • For θ ( 0 , π / 2 ] , based on Equation (7), only the integral term requires attention. By using a step-by-step integration and some properties of the special functions Ci ( x ) and Si ( x ) , we obtain
    0 1 0 1 C ( x , y ) d x d y = 0 1 0 1 x y cos θ ( 1 x ) ( 1 y ) d x d y = 1 θ 2 0 1 y θ x ( y 1 ) sin θ ( 1 x ) ( 1 y ) + cos θ ( 1 x ) ( 1 y ) ( 1 y ) 2 x = 0 x = 1 d y = 1 θ 2 0 1 y 1 cos θ ( 1 y ) ( 1 y ) 2 d y = 1 θ 2 log ( 1 y ) Ci [ θ ( 1 y ) ] θ Si [ θ ( 1 y ) ] + 1 cos [ θ ( 1 y ) ] 1 y y = 0 y = 1 = 1 θ 2 Ci ( θ ) + θ Si ( θ ) log ( θ ) + cos ( θ ) γ 1 .
The desired result follows immediately.
The proof of Proposition 2 ends. □
Thus, ρ is a decreasing function with respect to θ for θ [ 0 , π / 2 ] , with ρ = 0 for θ = 0 and ρ = 0.1 for θ = π / 2 . Figure 6 represents the Spearman rho for θ [ 0 , 1.5 ] .
In light of the above results, the cos-copula is adapted to model weak negative correlations.
Remark 2.
Based on Equations (6) and (8), upon integration over [ 0 , 1 ] 2 and some arrangements, the following series formula holds:
k = 0 + = 0 2 k m = 0 2 k ( 1 ) k + + m ( 2 k ) ! θ 2 ( k + 1 ) 2 k 2 k m 1 ( + 2 ) ( m + 2 ) = Ci ( θ ) + θ Si ( θ ) log ( θ ) + cos ( θ ) γ 1 .
Of course, this result is of independent interest to our copula study. It is, however, fascinating to see the versatility in the nature of the functions involved, i.e., logarithmic, special integral function, cosine function, and Euler–Mascheroni constant.
  • In complement of the Spearman rho, we can present the Kendall tau of the cos-copula. It is defined by
    τ = 4 0 1 0 1 C ( x , y ) c ( x , y ) d x d y 1 .
The closed form expression for τ is unmanageable due to the complexity of the product function C ( x , y ) c ( x , y ) . We can, however, state that it is a decreasing function with respect to θ for θ [ 0 , π / 2 ] , with τ = 0 for θ = 0 and τ = 0.06641903 for θ = π / 2 . Figure 7 represents the Kendall tau for θ [ 0 , 1.5 ] .
The small values of τ confirm the fact that the cos-copula is ideal to model weak negative correlations.
  • The cos-copula opens some interesting perspectives in distribution theory and modeling. The most immediate of these perspectives is the creation of simple and new two-dimensional distributions with cumulative distribution functions of the following form: H ( x , y ) = C ( F ( x ) , G ( y ) ) , so
    H ( x , y ) = F ( x ) G ( y ) cos θ ( 1 F ( x ) ) ( 1 G ( y ) ) ,
    where F ( x ) and G ( x ) denote two cumulative distribution functions. This gives two-dimensional trigonometric distributions, which seem slightly underexplored in the literature. By considering the exponential distribution as the parent, we may define the cos-two-exponential distribution by the following cumulative distribution:
    H ( x , y ; α , β ) = ( 1 e α x ) ( 1 e β y ) cos θ e ( α x + β y ) , ( x , y ) [ 0 , + ) 2 ,
    and H ( x , y ; α , β ) = 0 for ( x , y ) [ 0 , + ) 2 . To understand the importance of the trigonometric distributions in theory and practice for the uni-dimensional case, we may refer to [22,23].

2.4. Data Generation and Inference

The most straightforward method of generating random data (or values) from a distribution defined by a copula is what might be termed the inverse conditional method. For any positive integer n, with the aim of generating n data from the cos-copula, this method may be described as follows:
  • Generate n data ( s 1 , t 1 ) , , ( s n , t n ) from a random vector ( S , T ) , where S and T are independent random variables with the uniform distribution over ( 0 , 1 ) .
  • Choose a value of θ [ 0 , π / 2 ] .
  • Consider the following “conditional function”:
    C ( x , y ) = x C ( x , y ) = y cos θ ( 1 x ) ( 1 y ) + θ x ( 1 y ) sin θ ( 1 x ) ( 1 y ) .
  • For any i = 1 , , n , compute u i such that C ( s i , u i ) = t i .
  • Then ( s 1 , u 1 ) , , ( s n , u n ) are n data generated from the cos-copula defined with the chosen θ .
Other methods exist (see [24]). Such simulated data can be used for computational tests, or estimation purposes.
On the other hand, in a data analysis scenario, the angle parameter θ is generally unknown. Its evaluation is, thus, of interest for precise data fitting, or at least to know if θ is close to 0 corresponding to the independent case, or close to π / 2 , corresponding to the most highly correlated case. For this evaluation, it can be estimated from n data ( x 1 , y 1 ) , , ( x n , y n ) , that are susceptible to coming from the cos-copula distribution, by the maximum likelihood method; θ is thus estimated by
θ ˜ = argmax θ [ 0 , π / 2 ] i = 1 n c ( x i , y i ) = argmax θ [ 0 , π / 2 ] i = 1 n log [ c ( x i , y i ) ] ,
provided that it is unique. This estimation method ensures satisfying qualities to θ ˜ , such as underlying consistency and asymptotic normality, which are the basis for the construction of confidence intervals and statistical tests (see, for instance, Refs. [13,17]). At this point, concrete applications of the aforementioned techniques to real-world datasets remain a possibility.

3. Sine Angle Parameter Copula

This section completes the findings of the previous section; it is devoted to a simple sine copula with an angle parameter. This copula can be viewed as a parametric generalization of one copula introduced in [18]. Indeed, in ([18], Example 9), it is proved that the function C : [ 0 , 1 ] 2 [ 0 , 1 ] defined by
C ( x , y ) = x y 1 + sin π 4 ( 1 x ) ( 1 y ) ,
is a valid copula. Among the questions that arise are:
  • Can we replace the angle-value π / 4 with a tuning parameter and, if so, what is its “optimal values set”?
  • What are the related functions of such a copula?
  • What are its theoretical properties?
The answers to these questions are given below.

3.1. Definition and Graphics

The following proposition presents the considered angle parameter sine copula.
Proposition 3.
The function C : [ 0 , 1 ] 2 [ 0 , 1 ] defined by
C ( x , y ) = x y 1 + sin θ ( 1 x ) ( 1 y ) ,
with θ [ 1 , 1 ] is a valid copula.
Proof. 
Let us prove the main points defining an absolutely continuous two-dimensional copula, as recalled in Definition 1.
  • For any x [ 0 , 1 ] , we have C ( x , 0 ) = x × 0 × 1 + sin θ ( 1 x ) ( 1 0 ) = 0 , and, for any y [ 0 , 1 ] , C ( 0 , y ) = 0 × y 1 + sin θ ( 1 0 ) ( 1 y ) = 0 .
  • For any x [ 0 , 1 ] , we have C ( x , 1 ) = x × 1 × 1 + sin θ ( 1 x ) × 0 = x , similarly, for any y [ 0 , 1 ] , C ( 1 , y ) = 1 × y 1 + sin θ × 0 × ( 1 y ) = y .
  • For any ( x , y ) [ 0 , 1 ] 2 , using standard derivation techniques, we have
    2 x y C ( x , y ) = 1 + 1 θ 2 x y ( 1 x ) ( 1 y ) sin θ ( 1 x ) ( 1 y ) + θ ( 3 x y x y ) cos θ ( 1 x ) ( 1 y ) .
    Let us now study the sign of the above function by distinguishing the case θ [ 0 , 1 ] and the case θ [ 1 , 0 ) .
    Case θ [ 0 , 1 ] : We can write
    2 x y C ( x , y ) = A 1 + A 2 ,
    where
    A 1 = 1 + θ ( 3 x y x y ) cos θ ( 1 x ) ( 1 y )
    and
    A 2 = 1 θ 2 x y ( 1 x ) ( 1 y ) sin θ ( 1 x ) ( 1 y ) .
    Let us prove that A 1 0 and A 2 0 .
    For A 1 , let us first remark that
    3 x y x y = ( 1 x ) ( 1 y ) + 2 x y 1 1 .
    Therefore, since θ [ 0 , 1 ] and cos ( a ) [ 0 , 1 ] for any a [ 0 , 1 ] , we have
    A 1 1 θ cos θ ( 1 x ) ( 1 y ) 1 θ 0 .
    For A 2 , since θ [ 0 , 1 ] , we have θ ( 1 x ) ( 1 y ) [ 0 , 1 ] [ 0 , π / 2 ] , implying that sin θ ( 1 x ) ( 1 y ) 0 , and 1 θ 2 x y ( 1 x ) ( 1 y ) 0 . It follows that A 2 0 . Hence
    2 x y C ( x , y ) = A 1 + A 2 0 .
    The two-increasing property is proved.
Case θ [ 1 , 0 ) : For this case, we develop a strategy different to the previous case. We can write
2 x y C ( x , y ) = B 1 + B 2 ,
where
B 1 = θ 2 x y ( 1 x ) ( 1 y ) sin θ ( 1 x ) ( 1 y )
and
B 2 = 1 + sin θ ( 1 x ) ( 1 y ) + θ ( 3 x y x y ) cos θ ( 1 x ) ( 1 y ) .
Let us prove that B 1 0 and B 2 0 .
For B 1 , since θ [ 1 , 0 ) , we have sin θ ( 1 x ) ( 1 y ) 0 , which implies that B 1 0 .
For B 2 , since θ [ 1 , 0 ) , the following inequality holds: sin ( a ) a for any a > 0 , which implies that sin ( θ a ) = sin [ ( θ ) a ] θ a for any a > 0 . It follows from this inequality applied with a = θ ( 1 x ) ( 1 y ) , the inequality
θ ( 3 x y x y ) = θ [ x y + x ( y 1 ) + y ( x 1 ) ] θ x y ,
and the inequality: cos ( a ) 1 for any a R , applied with a = θ ( 1 x ) ( 1 y ) , that
B 2 1 + θ ( 1 x ) ( 1 y ) + x y cos θ ( 1 x ) ( 1 y ) 1 ( 1 x ) ( 1 y ) + x y = x + y 2 x y = x ( 1 y ) + y ( 1 x ) 0 .
Hence
2 x y C ( x , y ) = B 1 + B 2 0 .
The two-increasing property is proved.
As a result, C ( x , y ) is a valid two-dimensional copula. This ends the proof of Proposition 3.
For the purposes of this paper, we call the copula C ( x , y ) in Equation (9) as the sin-copula. Like the cos-copula, it has the feature to have a tuning angle parameter θ . Clearly, for θ = 0 , we have C ( x , y ) = x y ; the sin-copula is reduced to the independence copula, and for θ = π / 4 , it is reduced to the copula in ([18], Example 9).
To end the presentation, Figure 8 represents the two-dimensional plot of the sin-copula for selected values of θ .
From Figure 8, we see that the parameter θ skews the triangular shape of the sin-copula.
The functions of the sin-copula, as well as its key features, are discussed in the following parts.

3.2. Related Functions

The density associated with the sin-copula is the function c : [ 0 , 1 ] 2 [ 0 , + ) given by
c ( x , y ) = 2 x y C ( x , y ) = 1 + 1 θ 2 x y ( 1 x ) ( 1 y ) sin θ ( 1 x ) ( 1 y ) + θ ( 3 x y x y ) cos θ ( 1 x ) ( 1 y ) .
It is worth noting that c ( 1 , 0 ) = 1 θ , which is negative for θ > 1 , and c ( 1 , 1 ) = 1 + θ , which is negative for θ < 1 ; C ( x , y ) is not a copula for θ R [ 1 , 1 ] . In this sense, [ 1 , 1 ] is the optimal set of values for θ .
Figure 9 represents the two-dimensional plot of the sin-copula density for selected values of θ .
As a result of Figure 9, we can see how θ influences the overall form of the sin-copula density.
As a last important function, the survival sin-copula is the function C ^ : [ 0 , 1 ] 2 [ 0 , 1 ] defined by
C ^ ( x , y ) = x + y 1 + C ( 1 x , 1 y ) = x y + ( 1 x ) ( 1 y ) sin ( θ x y ) .
It establishes a valid copula, which is also a new trigonometric copula with θ as an angle parameter.
Figure 10 represents the two-dimensional plot of the survival sin-copula for selected values of θ .
The copula density associated with the survival sin-copula is given by
c ^ ( x , y ) = 2 x y C ^ ( x , y ) = 1 + 1 θ 2 x y ( 1 x ) ( 1 y ) sin ( θ x y ) + θ ( 3 x y 2 x 2 y + 1 ) cos ( θ x y ) .
To have an idea of its shapes, Figure 11 shows its two-dimensional plot for selected values of θ .
On the other hand, the sin-copula can be used in a variety of new mixed copulas. For instance, we can consider the following mixed sin-copulas:
  • Mixed copula 1: For any angle parameters θ 1 [ 1 , 1 ] and θ 2 [ 1 , 1 ] , and λ [ 0 , 1 ] , by setting C ( x , y ) = C ( x , y ; θ ) , we can consider
    C ( x , y ) = λ C ( x , y ; θ 1 ) + ( 1 λ ) C ( x , y ; θ 2 ) = x y + x y λ sin θ 1 ( 1 x ) ( 1 y ) + ( 1 λ ) x y sin θ 2 ( 1 x ) ( 1 y ) .
  • Mixed copula 2: Similarly, for any angle parameters θ 1 [ 0 , π / 2 ] and θ 2 [ 1 , 1 ] , and λ [ 0 , 1 ] , by setting C ( x , y ) = C ( x , y ; θ ) and C ( x , y ) = C ( x , y ; θ ) , we can set
    C ( x , y ) = λ C ( x , y ; θ 1 ) + ( 1 λ ) C ( x , y ; θ 2 ) = λ x y cos θ 1 ( 1 x ) ( 1 y ) + ( 1 λ ) x y 1 + sin θ 2 ( 1 x ) ( 1 y ) .
  • Mixed copula 3: Another example is
    C § ( x , y ) = 1 2 C ( x , y ) + y C ( 1 x , y ) = x y + 1 2 x y sin θ ( 1 x ) ( 1 y ) + y ( 1 x ) y sin θ x ( 1 y ) .
They are also new angle parameter trigonometric copulas by construction.

3.3. Properties

The main features of the sin-copula C ( x , y ) as described in Equation (9) are now listed.
  • As already mentioned before:
    For θ = 0 , it is clear that C ( x , y ) = x y . Therefore, the sin-copula is reduced to the independence copula.
    The set [ 1 , 1 ] is the optimal set of values for θ for validating C ( x , y ) as a copula.
  • For any θ [ 1 , 0 ] , we have C ( x , y ) x y , so the negative quadrant dependence property is satisfied. Similarly, for any θ [ 0 , 1 ] , we have C ( x , y ) x y , so the positive quadrant dependence property is satisfied (see [21]).
  • The sin-copula is symmetric since C ( x , y ) = C ( y , x ) for any ( x , y ) [ 0 , 1 ] 2 .
  • The sin-copula can be expressed under various analytical forms. Two of them are given below:
    In terms of simple cosine-sine functions, we can write
    C ( x , y ) = x y [ 1 + sin ( θ ) cos ( θ x ) cos ( θ y ) cos ( θ x y ) sin ( θ ) sin ( θ x ) sin ( θ y ) cos ( θ x y ) + sin ( θ ) sin ( θ x ) cos ( θ y ) sin ( θ x y ) + sin ( θ ) cos ( θ x ) sin ( θ y ) sin ( θ x y ) cos ( θ ) sin ( θ x ) cos ( θ y ) cos ( θ x y ) cos ( θ ) cos ( θ x ) sin ( θ y ) cos ( θ x y ) + cos ( θ ) cos ( θ x ) cos ( θ y ) sin ( θ x y ) cos ( θ ) sin ( θ x ) sin ( θ y ) sin ( θ x y ) ] .
    As a result, we can observe that the sin-copula has inherent analytical complexity.
    In terms of power series, by using the cosine series expansion and binomial formula, we get
    C ( x , y ) = x y + x y k = 0 + ( 1 ) k ( 2 k + 1 ) ! θ 2 k + 1 ( 1 x ) 2 k + 1 ( 1 y ) 2 k + 1 = x y + k = 0 + = 0 2 k + 1 m = 0 2 k + 1 ( 1 ) k + + m ( 2 k + 1 ) ! θ 2 k + 1 2 k + 1 2 k + 1 m x + 1 y m + 1 .
    In particular, upon differentiation with respect to x and y on the interior of the domain of convergence, one has
    c ( x , y ) = 1 + k = 0 + = 0 2 k + 1 m = 0 2 k + 1 ( 1 ) k + + m ( 2 k + 1 ) ! θ 2 k + 1 2 k + 1 2 k + 1 m ( + 1 ) ( m + 1 ) x y m .
    This expansion can be used to determine various moment-type measurements in a range of mathematical applications.
  • By arbitrary taking θ = 1 , we notice that
    C 1 4 , C 1 2 , 1 3 = 0.01071536 0.01119522 = C C 1 4 , 1 2 , 1 3 .
    As a result, the sin-copula is not Archimedean (see [1]).
  • The sin-copula is not radially symmetric since there exists ( x , y ) such that C ^ ( x , y ) C ( x , y ) .
  • As any copula, the Fréchet-Hoeffding bounds can be expressed as follows: For any ( x , y ) [ 0 , 1 ] 2 , we have max ( x + y 1 , 0 ) C ( x , y ) min ( x , y ) .
  • Thanks to the inequality: sin ( a ) a and the Jordan inequality: sin ( a ) ( 2 / π ) a for a [ 0 , π / 2 ] (see [20]), we have a copula ordering between the sin-copula and FGM copula:
    For θ [ 0 , 1 ] , we have
    C x , y ; 2 π θ C ( x , y ) C ( x , y ; θ ) .
    For θ [ 1 , 0 ] , the contrary holds:
    C x , y ; 2 π θ C ( x , y ) C ( x , y ; θ ) .
  • The following relationship between the cos-copula and sin-copula holds:
    C ( x , y ) C ( x , y ) = 2 2 x y sin θ ( 1 x ) ( 1 y ) 2 sin θ ( 1 x ) ( 1 y ) 2 + π 4 .
    Therefore, the following ordering results are established:
    For θ [ 1 , 0 ] , since the first sine term in Equation (12) is negative, and the second one is positive, we have
    C ( x , y ) C ( x , y ) .
    For θ [ 0 , 1 ] , since the first sine term in Equation (12) is positive, and the second one too, we have
    C ( x , y ) C ( x , y ) .
  • For any θ [ 0 , π / 2 ] , the two following results are obtained:
    λ L = lim x 0 C ( x , x ) x = lim x 0 x 1 + sin θ ( 1 x ) 2 = 0
    and
    λ U = lim x 1 1 2 x + C ( x , x ) 1 x = lim x 1 1 2 x + x 2 1 + sin θ ( 1 x ) 2 1 x = lim x 1 ( 1 x ) ( 1 + θ x 2 ) = 0 .
    Hence, the sin-copula has no tail dependence.
  • The medial correlation of the sin-copula is defined by
    M = 4 C 1 2 , 1 2 1 = sin θ 4 .
    It is clearly an increasing function with respect to θ for θ [ 1 , 1 ] , with M = 0.247404 for θ = 1 and M = 0.247404 for θ = 1 . Figure 12 represents the medial correlation for θ [ 1 , 1 ] .
    The possible values of this medial correlation are not negligible; we have M [ 0.25 , 0.25 ] . Hence, the sin-copula has a certain flexibility in this regard.
  • The Spearman rho of the sin-copula, as an example of copula, is defined by
    ρ = 12 0 1 0 1 C ( x , y ) d x d y 3 .
    Using well-known mathematical methods, the following assertion provides a mathematical expression for this measure.
Proposition 4.
The Spearman rho of the sin-copula can be expressed as
ρ = 0 , for θ = 0 12 θ 2 θ Ci ( θ ) + Si ( θ ) + θ ( log ( θ ) + γ 2 ) + sin ( θ ) , for θ ( 0 , 1 ] 12 θ 2 θ Ci ( θ ) + Si ( θ ) θ ( log ( θ ) + γ 2 ) sin ( θ ) , for θ [ 1 , 0 ) .
Proof. 
  • For θ = 0 , we have C ( x , y ) = x y and, by Equation (13), it is immediate that ρ = 0 .
  • For θ ( 0 , 1 ] , still based on the definition of ρ in Equation (13), we have
    ρ = 12 0 1 0 1 C ( x , y ) d x d y 3 = 12 0 1 0 1 x y sin θ ( 1 x ) ( 1 y ) d x d y .
    By using a step-by-step integration, we obtain
    0 1 0 1 x y sin θ ( 1 x ) ( 1 y ) d x d y = 0 1 0 1 x y sin θ ( 1 x ) ( 1 y ) d x d y = 1 θ 2 0 1 y sin θ ( 1 x ) ( 1 y ) θ x ( y 1 ) cos θ ( 1 x ) ( 1 y ) ( 1 y ) 2 x = 0 x = 1 d y = 1 θ 2 0 1 y θ ( 1 y ) sin θ ( 1 y ) ( 1 y ) 2 d y = 1 θ 2 0 1 θ Ci [ θ ( 1 y ) ] Si [ θ ( 1 y ) ] θ y θ log ( 1 y ) sin θ ( 1 y ) 1 y y = 0 y = 1 d y = 1 θ 2 θ Ci ( θ ) + Si ( θ ) + θ ( log ( θ ) + γ 2 ) + sin ( θ ) .
    Immediately, the intended result occurs.
  • For θ [ 1 , 0 ) , thanks to the oddity of the sine function, we can write
    ρ = 12 0 1 0 1 x y sin θ ( 1 x ) ( 1 y ) d x d y = 12 0 1 0 1 x y sin ( θ ) ( 1 x ) ( 1 y ) d x d y .
    Since θ ( 0 , 1 ] , the expression of ρ can be transposed with θ instead of θ , with the minus in factor of the overall expression.
The stated proposition is proved. □
The measure ρ is an increasing function with respect to θ for θ [ 1 , 1 ] , with ρ = 0 for θ = 0 , ρ = 0.3283896 for θ = 1 , and ρ = 0.3283896 for θ = 1 . Figure 13 represents the Spearman rho for θ [ 1 , 1 ] .
In light of the above results, the sin-copula is adapted to model moderate correlations, which may be negative or positive.
Remark 3.
Similarly to Remark 2, a fascinating series formula, not repertoried to the best of our knowledge, can be proved. Based on Equations (6) and (8), upon integration over [ 0 , 1 ] 2 , the following formula holds:
k = 0 + = 0 2 k + 1 m = 0 2 k + 1 ( 1 ) k + + m ( 2 k + 1 ) ! θ 2 k + 3 2 k + 1 2 k + 1 m 1 ( + 2 ) ( m + 2 ) = θ Ci ( θ ) + Si ( θ ) + θ ( log ( θ ) + γ 2 ) + sin ( θ ) .
  • In complement of the Spearman rho, we can present the Kendall tau of the sin-copula. It is defined by
    τ = 4 0 1 0 1 C ( x , y ) c ( x , y ) d x d y 1 .
    The complexity of the product function C ( x , y ) c ( x , y ) makes the closed form expression for τ unmanageable. We can, however, show that it is an increasing function with respect to θ for θ [ 1 , 1 ] , with τ = 0 for θ = 0 , τ = 0.2185653 for θ = 1 and τ = 0.2185653 for θ = 1 . Figure 14 represents the Kendall tau for θ [ 1 , 1 ] .
    The wide range of values of τ confirm the fact that the sin-copula is ideal to model moderate correlations.
  • Similarly to the cos-copula, the sin-copula opens up several fascinating possibilities, such as the development of simple and new two-dimensional distributions with cumulative distribution functions of the form: H ( x , y ) = C ( F ( x ) , G ( y ) ) , so
    H ( x , y ) = F ( x ) G ( y ) 1 + sin θ ( 1 F ( x ) ) ( 1 G ( y ) ) ,
    where F ( x ) and G ( x ) denote two cumulative distribution functions. As a result, probabilistic or statistical modeling in this context becomes more feasible.

3.4. Data Generation and Inference

The data generation method described in Section 2.4 can be configured for the sin-copula. For any positive integer n, we can generate n data from the sin-copula by proceeding as follows:
  • Generate n data ( s 1 , t 1 ) , , ( s n , t n ) from a random vector ( S , T ) , where S and T are independent random variables with the uniform distribution over ( 0 , 1 ) .
  • Choose a value for θ [ 1 , 1 ] .
  • Consider the following “conditional function”:
    C ( x , y ) = x C ( x , y ) = y 1 + sin θ ( 1 x ) ( 1 y ) θ x ( 1 y ) cos θ ( 1 x ) ( 1 y ) .
  • For any i = 1 , , n , compute u i such that C ( s i , u i ) = t i .
  • Then ( s 1 , u 1 ) , , ( s n , u n ) are n data generated from the sin-copula defined with the chosen value of θ .
Other techniques are available (see [24]).
In a data analysis scenario, the angle parameter θ is usually unknown. Its estimation is thus useful for exact data fitting, or at the very least for determining if θ is close to 0, corresponding to the independent case, or close to 1 or 1, corresponding to the most highly negatively or positively correlated situations, respectively. For this evaluation, it can be estimated from n data ( x 1 , y 1 ) , , ( x n , y n ) , that are susceptible to coming from the sin-copula distribution, by the maximum likelihood method; θ is thus estimated by
θ ˜ = argmax θ [ 1 , 1 ] i = 1 n c ( x i , y i ) = argmax θ [ 1 , 1 ] i = 1 n log [ c ( x i , y i ) ] .
Concrete applications of the aforementioned methodologies to real-world datasets are thus possible.

4. Conclusions and Perspectives

4.1. Conclusions

We have introduced and studied several trigonometric copulas that have the features to depend on a tuning angle parameter. We have shown that they are quite simple from the mathematical point of view, and possess interesting properties.
For the first copula, called cos-copula, we have demonstrated that: (i) it extends the independence copula; (ii) it has the negative quadrant property; (iii) it is symmetric; (iv) it is not Archimedean; (v) it is not radially symmetric; (vi) a special ordering exists between it and the FGM copula; (vii) it has no tail dependence; (viii) the medial correlation belongs to the interval [ 0.0762 , 0 ] ; (ix) the Spearman rho belongs to the interval [ 0.1 , 0 ] ; (x) the Kendall tau belongs to the interval [ 0.067 , 0 ] ; (xi) it can serve to create a plethora of trigonometric two-dimensional distributions. Thus, the cos-copula is adapted to model weak negative correlations. It is thus not adapted to model moderate or large correlations. The corresponding copula density, survival copula, and survival copula density have been expressed, as well as two mixed copula versions.
For the second copula, called the sin-copula, its properties can be listed as follows: (i) it extends the independence copula; (ii) it has the negative and positive quadrant properties; (iii) it is symmetric; (iv) it is not Archimedean; (v) it is not radially symmetric; (vi) comprehensive ordering exists between it and both the FGM copula and cos-copula; (vii) it has no tail dependence; (viii) the medial correlation belongs to the interval [ 0.25 , 0.25 ] ; (ix) the Spearman rho belongs to the interval [ 0.33 , 0.33 ] ; (x) the Kendall tau belongs to the interval [ 0.22 , 0.22 ] ; (xi) it can also serve to create a plethora of trigonometric two-dimensional distributions. Thus, the sin-copula is adapted to model moderate negative or positive correlations. The corresponding copula density, survival copula, and survival copula density have been expressed, as well as three mixed copula versions.
The shapes of the principal functions and their related features have been visually observed using graphics.

4.2. Perspectives

Thus, the first elements for the promotion of the proposed copulas are in this article. Perspectives of further research includes the following points:
  • Following the spirit of some power-extended FGM copulas (see [25]), one can think of considering some extensions of the cos-copula and sin-copula of the forms:
    C ( x , y ) = x y cos θ ( 1 x α ) c ( 1 y β ) c
    and
    C ( x , y ) = x y 1 + sin θ ( 1 x α ) c ( 1 y β ) c ,
    respectively, where α , β and c are newly introduced shape parameters. However, the possible values of α , β and c such that C ( x , y ) and C ( x , y ) are valid copulas remain to discover.
  • The n-dimensional versions of the cos-copula and sin-copula, which can be defined as C : [ 0 , 1 ] n [ 0 , 1 ] and C : [ 0 , 1 ] n [ 0 , 1 ] , respectively, where
    C ( x 1 , , x n ) = i = 1 n x i cos θ i = 1 n ( 1 x i )
    and
    C ( x 1 , , x n ) = i = 1 n x i 1 + sin θ i = 1 n ( 1 x i ) ,
    respectively, deserve to be investigated for n 3 . Especially, the possible values for θ in this case need to be determined.
  • The use of the sin-copula and cos-copula for data analysis, following the methodology in [13,15], is an important perspective. Furthermore, the development of an R package is envisageable, inspired by the research design of Cylcop developed by [17].
  • Last but not least, other simple angle parameter copulas can be created on the basis of this study. One could think of C : [ 0 , 1 ] 2 [ 0 , 1 ] defined by
    C ( x , y ) = x y 1 + tan θ ( 1 x ) ( 1 y ) ,
    but the optimal values of θ such that C ( x , y ) is a valid copula remain unknown. It is proven that θ [ 0 , π / 4 ] satisfies the required conditions (in an unpublished work), but this set is thought to be less than optimal.

Funding

This research received no external funding.

Acknowledgments

We thank the two reviewers and the associate editor for their in-depth comments on the first version of the article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Nelsen, R. An Introduction to Copulas, 2nd ed.; Springer Science+Business Media, Inc.: Singapore, 2006. [Google Scholar]
  2. Durante, F.; Sempi, C. Principles of Copula Theory; CRS Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  3. Emura, T.; Chen, Y.-H. Analysis of Survival Data with Dependent Censoring, Copula-Based Approaches; JSS Research Series in Statistics; Springer: Singapore, 2018. [Google Scholar]
  4. Joe, H. Dependence Modeling with Copulas; CRS Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  5. Alfonsi, A.; Brigo, D. New families of copulas based on periodic functions. Commun. Stat.-Theory Methods 2005, 34, 1437–1447. [Google Scholar] [CrossRef] [Green Version]
  6. Amblard, C.; Girard, S. Symmetry and dependence properties within a semiparametric family of bivariate copulas. J. Nonparametr. Stat. 2002, 14, 715–727. [Google Scholar] [CrossRef] [Green Version]
  7. Chesneau, C. A study of the power-cosine copula. Open J. Math. Anal. 2021, 5, 85–97. [Google Scholar] [CrossRef]
  8. Chesneau, C. On new types of multivariate trigonometric copulas. Appl. Math. 2021, 1, 3–17. [Google Scholar] [CrossRef]
  9. Chesneau, C. A note on a simple polynomial-sine copula. Asian J. Math. Appl. 2022, 2, 1–14. [Google Scholar]
  10. Durante, F. A new class of symmetric bivariate copulas. J. Nonparametr. Stat. 2006, 18, 499–510. [Google Scholar] [CrossRef]
  11. Jones, M.C.; Pewsey, A.; Kato, S. On a class of circulas: Copulas for circular distributions. Ann. Inst. Stat. Math. 2015, 67, 843–862. [Google Scholar] [CrossRef]
  12. Davy, M.; Doucet, A. Copulas: A new insight into positive time-frequency distributions. IEEE Signal Process. Lett. 2003, 10, 215–218. [Google Scholar] [CrossRef]
  13. Hodel, F.H.; Fieberg, J.R. Circular-linear copulae for animal movement data. bioRxiv 2021. [Google Scholar] [CrossRef]
  14. Knockaert, L. A class of positive isentropic time-frequency distributions. IEEE Signal Process. Lett. 2002, 9, 22–25. [Google Scholar] [CrossRef] [Green Version]
  15. Susam, S.O. A new family of archimedean copula via trigonometric generator function. Gazi Univ. J. Sci. 2020, 33, 795–802. [Google Scholar] [CrossRef]
  16. Xiao, Q.; Zhou, S.W. Modeling correlated wind speeds by trigonometric Archimedean copulas. In Proceedings of the 11th International Conference on Modelling, Identification and Control (ICMIC2019); Lecture Notes in Electrical Engineering; Wang, R., Chen, Z., Zhang, W., Zhu, Q., Eds.; Springer: Singapore, 2020; Volume 582. [Google Scholar]
  17. Hodel, F.H.; Fieberg, J.R. Cylcop: An R package for circular-linear copulae with angular symmetry. bioRxiv 2021. [Google Scholar] [CrossRef]
  18. Manstavičius, M.; Bagdonas, G. A class of bivariate independence copula transformations. Fuzzy Sets Syst. 2021, 428, 58–79. [Google Scholar] [CrossRef]
  19. Kober, H. Approximation by integral functions in the complex domain. Trans. Am. Math. Soc. 1944, 56, 7–31. [Google Scholar] [CrossRef] [Green Version]
  20. Qi, F.; Hao, Q.-D. Refinements and sharpenings of Jordan’s and Kober’s inequality. Math. Inform. Q. 1998, 8, 116–120. [Google Scholar]
  21. Lehmann, E.L. Some concepts of dependence. Ann. Math. Stat. 1966, 37, 1137–1153. [Google Scholar] [CrossRef]
  22. Chesneau, C.; Artault, A. On a comparative study on some trigonometric classes of distributions by the analysis of practical data sets. J. Nonlinear Model. Anal. 2021, 3, 225–262. [Google Scholar]
  23. Tomy, L.; Satish, G. A review study on trigonometric transformations of statistical distributions. Biomet. Biostat. Int. J. 2021, 10, 130–136. [Google Scholar]
  24. Johnson, M.E. Multivariate Statistical Simulation; J. Wiley & Sons: New York, NY, USA, 1987. [Google Scholar]
  25. Bairamov, I.; Kotz, S. Dependence structure and symmetry of Huang-Kotz FGM distributions and their extensions. Metrika 2002, 56, 55–72. [Google Scholar] [CrossRef]
Figure 1. Representations of the cos-copula for (a) θ = 0.2 and (b) θ = 1.5 .
Figure 1. Representations of the cos-copula for (a) θ = 0.2 and (b) θ = 1.5 .
Modelling 03 00010 g001
Figure 2. Representations of the cos-copula density for (a) θ = 0.2 and (b) θ = 1.5 .
Figure 2. Representations of the cos-copula density for (a) θ = 0.2 and (b) θ = 1.5 .
Modelling 03 00010 g002
Figure 3. Representations of the survival cos-copula for (a) θ = 0.2 and (b) θ = 1.5 .
Figure 3. Representations of the survival cos-copula for (a) θ = 0.2 and (b) θ = 1.5 .
Modelling 03 00010 g003
Figure 4. Representations of the survival cos-copula density for (a) θ = 0.2 and (b) θ = 1.5 .
Figure 4. Representations of the survival cos-copula density for (a) θ = 0.2 and (b) θ = 1.5 .
Modelling 03 00010 g004
Figure 5. Representation of the medial correlation of the cos-copula for θ [ 0 , 1.5 ] .
Figure 5. Representation of the medial correlation of the cos-copula for θ [ 0 , 1.5 ] .
Modelling 03 00010 g005
Figure 6. Representation of the Spearman rho of the cos-copula for θ [ 0 , 1.5 ] .
Figure 6. Representation of the Spearman rho of the cos-copula for θ [ 0 , 1.5 ] .
Modelling 03 00010 g006
Figure 7. Representation of the Kendall tau of the cos-copula for θ [ 0 , 1.5 ] .
Figure 7. Representation of the Kendall tau of the cos-copula for θ [ 0 , 1.5 ] .
Modelling 03 00010 g007
Figure 8. Representations of the sin-copula for (a) θ = 0.2 and (b) θ = 1 .
Figure 8. Representations of the sin-copula for (a) θ = 0.2 and (b) θ = 1 .
Modelling 03 00010 g008
Figure 9. Representations of the sin-copula density for (a) θ = 0.2 and (b) θ = 1 .
Figure 9. Representations of the sin-copula density for (a) θ = 0.2 and (b) θ = 1 .
Modelling 03 00010 g009
Figure 10. Representations of the survival sin-copula for (a) θ = 0.2 and (b) θ = 1 .
Figure 10. Representations of the survival sin-copula for (a) θ = 0.2 and (b) θ = 1 .
Modelling 03 00010 g010
Figure 11. Representations of the survival sin-copula density for (a) θ = 0.2 and (b) θ = 1 .
Figure 11. Representations of the survival sin-copula density for (a) θ = 0.2 and (b) θ = 1 .
Modelling 03 00010 g011
Figure 12. Representation of the medial correlation of the sin-copula for θ [ 1 , 1 ] .
Figure 12. Representation of the medial correlation of the sin-copula for θ [ 1 , 1 ] .
Modelling 03 00010 g012
Figure 13. Representation of the Spearman rho of the sin-copula for θ [ 1 , 1 ] .
Figure 13. Representation of the Spearman rho of the sin-copula for θ [ 1 , 1 ] .
Modelling 03 00010 g013
Figure 14. Representation of the Kendall tau of the sin-copula for θ [ 1 , 1 ] .
Figure 14. Representation of the Kendall tau of the sin-copula for θ [ 1 , 1 ] .
Modelling 03 00010 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chesneau, C. Theoretical Study of Some Angle Parameter Trigonometric Copulas. Modelling 2022, 3, 140-163. https://doi.org/10.3390/modelling3010010

AMA Style

Chesneau C. Theoretical Study of Some Angle Parameter Trigonometric Copulas. Modelling. 2022; 3(1):140-163. https://doi.org/10.3390/modelling3010010

Chicago/Turabian Style

Chesneau, Christophe. 2022. "Theoretical Study of Some Angle Parameter Trigonometric Copulas" Modelling 3, no. 1: 140-163. https://doi.org/10.3390/modelling3010010

APA Style

Chesneau, C. (2022). Theoretical Study of Some Angle Parameter Trigonometric Copulas. Modelling, 3(1), 140-163. https://doi.org/10.3390/modelling3010010

Article Metrics

Back to TopTop