Skip to main content

BRIEF RESEARCH REPORT article

Front. Appl. Math. Stat., 29 September 2022
Sec. Dynamical Systems
This article is part of the Research Topic New Trends in Control and Applications of Complex System View all 6 articles

Correlation analysis of financial assets based on asymmetric copula

\nXia Li
Xia Li1*Bing HouBing Hou2
  • 1School of Mathematics and Statistics, Henan University, Kaifeng, China
  • 2School of Mathematics and Statistics, North China University of Water Resources and Electric Power, Zhengzhou, China

Based on the asymmetric copula function, this paper analyzes the static and dynamic correlation between Shanghai Composite Index and Shenzhen Composite Index. Through the static analysis, it is found that the asymmetric copula function is better than Gumbel Copula in describing the distribution characteristics of the top tail dependence between the Shanghai Composite Index and the Shenzhen Composite Index, and the copula correlation coefficient definition based on the asymmetric copula function can well describe the asymmetric dependence between variables. In the time-varying analysis, the paper improves the traditional dynamic evolution equation of the tail-dependence coefficient. Through empirical analysis, the result shows that the improved dynamic evolution equation can better reflect the dynamic evolution process of the tail-dependence coefficient.

Introduction

In investment risk management, diversified investment strategies are generally adopted to avoid risks, and a large number of studies show that selecting assets with low correlation for a portfolio can reduce investment risks. Therefore, it is very important to analyze the correlation between assets in the risk analysis of financial assets.

In origin correlation analysis, Pearson's linear correlation coefficient is generally used to describe the correlation between variables. However, since the calculation of Pearson's correlation coefficient requires the existence of finite variance, financial time series usually follow a thick-tailed distribution, and the variance sometimes does not exist. This greatly limits the application of Pearson's correlation coefficient in the correlation analysis of financial assets.

As a function that connects the joint distribution function and the edge distribution function of random variables, the copula almost contains all the dependent information of random variables. Because it does not limit the type of edge distribution, the copula has been widely used in describing the correlation between financial assets. Common copula-based correlation measures include Spearman's ρ and Kendall's τ.

In recent years, some scholars have done a lot of studies on the application of the copula model and obtained some important results. For example, Romano [1] applied copula theory to study the problem of portfolio risk in the Italian capital market. Rosenberg and Schuerman [2] used the copula theory to analyze the aggregation of market risk, credit risk, and operational risk. Patton [3] studied the impact of asymmetric correlation structure on the performance of two asset portfolios. In addition, Patton [4] proposed the definition of conditional copula and used it to characterize the asymmetric dependence structure between exchange rates. Palaro and Hotta [5] used the conditional copula to dynamically analyze the dependence relationship among financial assets. Nikolai et al. [6] summarized in detail the development and application of copula function in recent years and made a reasonable prospect for its future development. Many Chinese scholars have also conducted numerous studies on the copula and obtained many conclusions [718]. Zhang [7] analyzed the feasibility of the copula function application in financial data analysis. Zhang [8] studied the calculation method of multi-asset VaR based on the copula. Wu et al. [9] applied the copula theory to the risk analysis of the financial asset portfolio. Besides, Wu et al. [11] combined the T GARCH model with multivariate normal copula and multivariate T copula to analyze the portfolio problem among multiple stocks. Wei and Zhang [12] combined the copula function with the GARCH model to analyze the correlation structure between Shanghai and Shenzhen stock markets, and Wei and Zhang [13] used the copula function and the GARCH model to conduct a dynamic analysis of the interdependence between the various sectors of the Shanghai Stock market. Zhang and Li [14] studied the integration risk of the asset portfolio based on copula from the goodness-of-fit. Li et al. [16] measured the financial risk based on Claytoncopula, and Du and Zhang [17] used mixed cane copula to compare the calculation accuracy of VAR of the asset portfolio. Ren and Zhang [18] studied the value at risk of the financial asset portfolio by using multiple Archimedean Copula. Among copula functions, Archimedean Copula has been widely used in the financial field due to its simple construction method and thick-tailed distribution. For example, Clayton Copula [19] with the characteristics of lower tail dependence distribution, Frank Copula [20] with asymptotically independent tail distribution, Gumbel Copula [21] with upper tail dependence, and Joe copula [22] with both upper and lower tails dependence all have been widely used in correlation analysis of financial assets. However, the correlation measures spearman's ρ, Kendall's τ, Blest's correlation coefficient, Kochar, and Gupta's based on the Archimedean copula are all symmetric correlation measures. This is not characteristic of the relationship between financial assets. Generally speaking, for two variables X and Y, the dependence relationship between X and Y is not the same as that between Y and X. In order to describe the dependence structure between variables more accurately, scholars have carried out a large number of studies on the asymmetric copula model between variables. Alfonsi and Brigo [23] proposed a new construction method for asymmetric copula functions based on periodic functions. Khoudragi [24] proposed a method to construct copula based on the product copula function. Liebscher [25] proposed a method for constructing asymmetric copula functions based on copula and Archimedes Copula, respectively.

In this paper, based on the asymmetric copula function that constructed by Eckhard Liebscher, this paper analyzes the static and time-varying dynamic relationship between the Shanghai Composite Index and the Shenzhen Composite Index. Through analysis, it is found that the asymmetric copula can better depict the tail distribution characteristics between variables than the traditional Archimedes copula Gumbel function, and, in the static analysis, the concept of the copula correlation coefficient is introduced. For the asymmetric copula function, copulas correlation coefficients ρ1|2 and ρ2|1, in general, are not equal; the paper engraves the asymmetric relationship between variables through the calculation of the asymmetric copulas correlation coefficient. In addition, in the time-varying analysis, this paper improves the evolution process of the traditional tail-dependence coefficient and makes a comparative analysis between the new evolution process and the traditional evolution process through empirical analysis. The result shows that the improved dynamic evolution equation can better reflect the dynamic evolution process of the tail-dependence coefficient.

The paper is arranged as follows: Part one consists of the Introduction and Part two consists of the construction of the asymmetric copula function. This section introduces the construction method of the asymmetric copula function proposed by Eckhard Liebscher and discusses its tail distribution characteristics. The third part is a correlation analysis of asymmetric copula function. First, this paper gives the parameter estimation method of the asymmetric copula function by using the empirical tail copula function; in this section, the definition of the copula correlation coefficient is introduced, and the estimation value of the copula correlation coefficient is calculated by Monte Carlo simulation so as to explain the asymmetric dependence relationship between variables. Part four is a static empirical analysis based on the asymmetric copula function. The correlation between the Shanghai Composite Index returns and the Shenzhen Composite Index returns is analyzed statically by using the asymmetric copula function. Part five is a dynamic empirical analysis based on asymmetric copula function. The traditional dynamic evolution equation is improved, and the empirical analysis shows that the improved dynamic evolution equation can more accurately describe the dynamic evolution process of the upper- and tail-dependence coefficients. Part six consists of the Conclusion.

Construction method of asymmetric copula function

Archimedes copula function

Let C (u1,,ud)=φ-1[φ(u1)++φ(ud)],; here, ui ∈ [0, 1], φ is a strictly decreasing function for [0, 1] to [0, ∞). According to a document [10], if φ−1 is d- monotonous, namely, (-1)kdkdtkφ-1(t)0, the function C(u1, ⋯, ud) is a copula function.

The Archimedean copula is a symmetric copula function, which can also be expressed in another form:

C(u1,,ud)=φ0-1(φ0(u1)××φ0(ud)),ui[0,1],where φ0:[0,1][0,1],φ0(t)=exp(-φ(t)).

In the Archimedean copulas, Clayton copula, Frank copula, and Gumbel copula are widely used in the financial field due to their thick-tailed distribution.

Tail copula function and tail-dependence coefficient

Definition 2.1. Set H is a distribution function, the corresponding copula is C, if for each point in ¯+2=[0,]2{(,)}, the following limits ΛL(x,y)=limtt·C(xt,yt) exists, then, the function ΛL(x, y): ¯2 is called the lower-tail copula function associated with H.

Likewise, if the ΛU(x,y)=limttC¯(xtlimt,yt)  limit exists, we say ΛU(x, y) the upper tail copula function.

Tail-dependence estimation is not an easy task, especially for non-standard distributions. It is for this reason that we are prompted to consider tail-dependence coefficients, or, rather, tail-dependence coefficients are a special case of tail copula functions.

Definition 2.2. (The tail-dependence coefficient) A random vector (X, Y) is called upper tail dependence; if ΛU(1, 1) exists and λU=ΛU(1,1)=tC¯(1t,1t) > 0, λU is called the upper tail-dependence coefficient; if λU = 0, (X, Y) is said to be upper tail independent. Similarly, the lower tail-dependence coefficient is defined by λL = ΛL(1, 1). If λL > 0(= 0), the lower tail-dependence coefficient exists, and λL = 0 says (X, Y) is the lower tail asymptotically independent.

In addition, for λU=limtt·C¯(1t,1t)=limtC¯(1t,1t)1t=limu0C¯(u,u)u=limu1C¯(1-u,1-u)1-u;λL=limtt·C(xt,yt)=limtC(xt,yt)1t=limu0C(u,u )u.

Therefore, this definition is essentially consistent with the definition of the tail-dependence coefficient by Joe [22].

The upper- and lower-tail-dependence coefficients are the values of the upper- and lower-tail-dependence functions at a point. By comparison, the tail-dependence function can more intuitively reflect the dependence structure between variables.

As the dependence relationship between financial variables is not fixed, the dependence coefficient and even the dependence structure between variables may change with the change of time. Therefore, many scholars have carried out a large number of dynamic analyses on the dependence structure between financial variables. In the past dynamic correlation analysis, we mostly used the symmetric copula function model to describe the dependence relationship, so the dependence relationship is also mostly symmetric, but, in practice, the dependence relationship between different variables is not symmetric. Specifically, for X and Y variables, the dependence relationship between X and Y is not equal to the dependence relationship between Y and X. Based on this, this paper uses the asymmetric copula function to dynamically analyze the correlation between financial variables. First, the asymmetric copula function model that was selected for analysis is presented.

Asymmetric copula function model

According to the expression of Archimedean copulas C(u1,,ud)=φ0-1(φ0(u1)××φ0(ud)),ui[0,1],φ0:[0,1][0,1], we can see, since the functions acting on each ui are the same, the Archimedean copula is characterized by symmetry. In order to construct the asymmetric copula function, we can compound different functions φ0(ui) for different ui. Based on this idea, Eckhard Liebscher generalized the Archimedes copula function in literature [25] and proposed a construction method of the asymmetric copula function.

Theorem 2.1 Let C(u1,,ud)=φ(1mj=1mhj1(φ-1(u1))hjd(φ-1(ud))), here, φ(d)exists,φ′(u) > 0, φ(i)(u) ≥ 0, i = 2, ⋯, d, u ∈ [0, 1], let φ(0) = 0, φ(1) = 1. Assume that, for each j ∈ {1, ⋯, m} and k ∈ {1, ⋯, d}, hjk is derivable and strictly increasing function on [0, 1] to [0, 1], and hjk(0)=0,hjk(1)=1,1mj=1mhjk(x)=x,k=1,,d,x[0,1], so C(u1, ⋯, ud) can be concluded to be an absolutely continuous copula function.

Proof: Fist,

C(u1,1,,1)=φ(1mj=1mhj1(φ-1(u1)) )= u1,C(0,u2,,ud)=φ(0 )=0.

Then, we show that du1u2udC(u)0 for all u; i.e., the copula has a density. Let φ~(·)=φ(·m). By induction, we can prove that

ku1u2ukC(u)=v=1kφ˜(v)(j=1mhj1(φ1(u1))hjd(φ1(ud)))×M1Mμ is a decomposition of {1,,k}avk(M1,,Mμ)×i=1μuMi(j=1mhj1(φ1(u1))hjd(φ1(ud)))()

For k = 1, ⋯, d. The notation uM means lui1uil, where M = {i1, ⋯, il}. In formula (*), avk(M1, ⋯, Mμ) denotes an integer depending on k, v, M1, ⋯, Mμ. Obviously, (*) is true for k = 1. Differentiating both sides of (11) w·r·t·uk+1, the right-hand side gets the form of the right-hand side of (11) with k replaced by k + 1. Therefore, by induction, we obtain that (11) is valid for all 1 ≤ kd. Since mui1uim(j=1mhj1(φ-1(u1))hjd(φ-1(ud)))0 for different i1, ⋯, im by assumption, we have uM(j=1mhj1(φ-1(u1))hjd(φ-1(ud)))0 for M ⊂ {1, ⋯, k}. Remember that φ~(v)(u)0, by assumption. Thus, the proof of this theorem is complete.

Example 2.1. Let φ(t)=exp(-α(1-t)1r)-e-α1-e-α, so φ1(t)=1(ln(eα+tteαα)r, then set h1k=xζk,ζk[1,2],h2k=2xxζk, where r ∈ (1, +∞), to be simple, let d = 2, m = 2. According to the theorem, φ(t) and hij(x) can be composited to construct the asymmetric copula function.

C*(u1,u2)=φ(12j=12hj1(φ-1(u1))hj2(φ-1(u2)))

In the following, the upper tail copulaΛU(x, y) and the upper- and lower tail-dependence coefficients of the copula C*(u1,u2) are calculated.

First, the upper tail copula function is evaluated as follows:

ΛU(x,y)=x+y-limtC(1-xt,1-yt)(1-xt)x                  -limtC(1-xt,1-yt)(1-yt)y              =x+y-x·21r-1-y·21r-1λU(r)=ΛU(1,1)=1 + 1-21r-1-21r-1=2-21r

The lower-tail-dependence coefficient λL = 0 can be calculated by using the same method, that is, the copula function constructed in the example is upper tail dependence and lower tail asymptotically independent.

In order to compare the tail dependence structure of asymmetric copula and Archimedean copula, this paper analyzes the tail-dependence function of Gumbel copula, which also has the characteristics of upper-tail-dependence distribution.

Gumbel copula is an Archimedean copula that is generated by φ(t)=(-lnt)1α,α(0,1], C(u,v)=exp(-((-lnu)1α+(-lnv)1α)α), the upper-tail-dependence function of C(u, v) is calculated by the following:

ΛU(x,y)=limttC¯(xt,yt)                  =limtt(xt+yt-1+C(1-xt,1-yt))                  =x+y+limtt·C(1-xt,1-yt)                  =x+y+limtC(1-xt,1-yt)-11t                  =x+y-2α-1·x-2α-1·y

So, λU(r)=ΛU(1,1)=1 + 1-2α-1-2α-1=2-2α,; the upper-tail-dependence function and the upper-tail-dependence coefficient of Gumbel copula have the same form as the asymmetric copula function  C*(u1,u2 ).

Correlation analysis of the asymmetric copula function

Parameter estimation method of the asymmetric copula function

Due to the complex structure, the commonly used parameter estimation methods, such as maximum likelihood estimation, marginal distribution inference method, and maximum likelihood estimation, are not easy to calculate for the asymmetric copula function. In this paper, we present a parameter estimation method of asymmetric copula function based on empirical tail copula function.

First, the concept of rank is given.

Definition 3.1. Let(x1, y1), (x2, y2), ⋯, (xN, yN) be N samples of a random variable (X, Y); arrange x1, x2, ⋯, xN in order, from small to large; the position of xn in the queue, namely, rank rn, is called its rank. Similarly, the rank of yn iny1, y2, ⋯, yNis sn, where n = 1, 2, ⋯, N.

Suppose (X, Y), (X1, Y1), ⋯, (Xm, Ym) are random vectors with distribution functions F, edge distributions G, H, and joint distributions C, in the following non-parametric estimates of upper- and lower-tail-dependent functions are given.

Let Cm represent the empirical copula function, defined as follows:

Cm(u,v)=Fm(Gm-1(u),Hm-1(v)),(u,v)[0,1]2 

Here, Fm, Gm, Hm represent the empirical distribution functions corresponding to F, G, H. Similarly, experience survival can be defined as copula Cm¯(u,v).

Cm¯(u,v)=H¯(Fm¯-1(u),Gm¯-1(v)),(u,v)[0,1]2

Hm¯(x,y)=1mj=1mI{xj>x,yj>y},Fm¯(u)=1Fm(u),Gm¯(v)=1Gm(v), so Fm¯1(u)=sup{x|F¯(x)u}, and sup {Φ} = −∞

Let Rm1j and Rm2j represent the rank of Xj and Yj, respectively, j = 1, 2, ⋯, m, and then the estimation of tail-dependent function as the following:

Λ^L,m(x,y)=mkCM(kxm,kym)1kj=1mI{Rm1jkx,Rm2jky}Λ^U,m(x,y)=mkC¯M(kxm,kym)                    1kj=1mI{Rm1jm-kx,Rm2jm-ky}

They are called the experience tail copulas function; λ^U=Λ^U(1,1),λ^L=Λ^L(1,1)are called the experience tail correlation coefficient.

Since the empirical tail copula function is calculated according to the empirical copula function, the empirical copula function is only related to the sample data. For the same group of sample data, when the asymmetric copula function or Archimedes copula function is used to describe the interdependent structure between variables, an empirical tail copula can be used as an estimate of their tail copulas. Through previous analysis, we know the upper-tail-dependence function of asymmetric copulas as  ΛUC*(x,y)=x+y-x·21r-1,and the upper-tail-dependence coefficient λUC*(x,y)=ΛU(1,1)=2-21r, the upper-tail-dependence function of Gumbel copulas as ΛUG(x,y)=x+y-x·2α-1-x·2α-1, the upper-tail-dependence coefficient λUG(x,y)=ΛU(1,1)=2-2α, because the empirical tail-dependence coefficient λ^U,λ^Lcan also be estimated as either λUC or λUG. Since  λUC*(r)=2-21r, we can get an estimation of r that is calculated from λUC*. In addition to using λU G to estimate λUC*, we can also directly estimate λUC* with empirical tail-dependence coefficients, and also get the parameter estimates of the asymmetric copula function.

Copula correlation coefficient

Nelsen [26] gives Spearman's ρ and Kendall's τ expressions.

Definition 3.2. Let the random variable (U, V) with copula C(U, V), expressions of Spearman's ρ and Kendall's τ as follows:

ρ(C)=120101C(u,v)dudv-3,τ(C)             =40101C(u,v)dC(u,v)-1

In correlation analysis, Pearson's correlation coefficient describes the linear correlation between variables, while Spearman's ρ and Kendall's τ describe the non-linear correlation between variables, but cannot describe the asymmetric relationship between variables. Therefore, in order to describe the asymmetric correlation between variables, many scholars have carried out numerous studies. Sungur [27] proposed a correlation coefficient based on the copula function. Dabrawska [28] analyzed a correlation measure based on regression. Dette et al. [29] proposed a regression-dependent non-parametric measure based on copula function on the basis of regression correlation measure. Shih and Emura [30] improved the copula correlation coefficient proposed by Sungur and gave an expression of the copula correlation coefficient based on convolution operation. Based on the copula correlation coefficient proposed by Jia-Han Shih, this paper analyzes the asymmetric dependence relation of the asymmetric copula function that is presented in this paper.

Sungur [27] pointed out, ρ2|1(C) and ρ1|2(C) respectively express directional dependence of V on Uand directional dependence of U on V, where ρ2|1(C)=Var{E(V|U)}Var(V),ρ1|2(C)=Var{E(U|V)}Var(U),ρ2|1(C),ρ1|2(C)[0,1]. As can be seen from the definition that ρ2|1(C) measure the part of the variation in V that is explained by U and ρ1|2(C) measure the part of the variation in U that is explained by V, if ρ2|1(C) > ρ1|2(C) , then it means that U explains more change in V than V explains the change in U.

Based on Sungur's research, Jia - Han Shih improved the definition of the type of ρ2|1(C) and ρ1|2(C), and put forward the following two methods of definition.

Definition 3.3. Let the random variable (U, V) have copula C(U, V), according to Spearman's ρ of the convolution copula CT*C or C*CT, the copula correlation coefficient can be defined as follows:  ρ2|1(C)=ρ(CT*C) and ρ1|2(C)=ρ(C*CT).

where convolution is defined as follows:

(CT*C)(u,v)=012CT(u,t)1(t,v)dt, u,v[ 0,1]

Since Parsow et al. [31] proved that if C1, C2 are copula, then C1*C2(u, v) is also copula, that is, the copula function is closed under convolution operation. Therefore, the above definition indicates that the copula correlation coefficient ρ2|1(C1|2(C) is Spearman's ρ of copula function obtained after convolution operation of CT and C. Thus, the calculation problem of the copula correlation coefficient is transformed into Spearman's ρ of convolution copula function. According to the definition, if C is asymmetric, then ρ2|1(C) ≠ ρ1|2(C). If C is symmetrical, then ρ2|1(C) = ρ1|2(C).

According to the definition of Spearman's ρ, we can obtain:

ρ2|1(C)=ρ(CT*C)=120101CT*Cdudv-3ρ1|2(C)=ρ(C*CT)=120101C*CTdudv-3

If the C is symmetrical, CT = C, then ρ2|1(C)=120101CT*Cdudv-3 = 120101C*CTdudv-3=ρ1|2(C),. There is no conclusion that the Archimedean copula is closed under convolution; in other words, for an Archimedean copula C (u, v), convolving with itself is not necessarily the same type of Archimedean copula. The definition of Spearman's ρ is ρ(C)=120101C(u,v)dudv-3,; therefore, ρ2|1(C) and ρ(C) are two dependent indicators from different angles. Even for the FGM copula function closed under convolution operation, that is, for the FGM copula, CT*C and C*CT are still FGM copula, but the obtained copula function and the original copula function have different parameters, and the dependence relationship is also different.

The copula correlation coefficient gives us a way to define the asymmetric correlation between variables, but, for the general asymmetric copula function, expressions of Cu(v) and Cv(u) are very complex and difficult to calculate. Therefore, the literature [31] presents a method to calculate estimates of copula correlation coefficients using the quasi-Monte Carlo simulation.

Estimation of the copula correlation coefficient

Sungur [27, 33] mentioned that, if (U, V) is a random vector with uniform edge distribution of [0, 1] and copula C, then in the case of given U = u, the conditional distribution function of V is represented by CU(V), and then CU(V)=P(Vv|U=u)=C(U,V)u. The copula regression function of V on U can be expressed with rV|UC(u), and rV|UC(u)=EC(V|U=u)=1-01CU(V)dV. Similarly, CV(U) =P(Uu|V=v)=C(U,V)v,rU|VC(v)=EC(U|V=v)=1-01CV(U)dU, and then,

ρ1|2(C)=Var{E(U|V)}Var(U)=Var{rU|Vc(v)}Var(U)                =E[(rU|VC(v)12)2]112=12E[(rU|VC(v))2]3ρ2|1(C)=Var{E(V|U)}Var(V)=Var{rV|Uc(u)}Var(V)                =E[(rV|UC(u)12)2]112=12E[(rV|UC(u))2]3.

For asymmetric copula functions, CU(V) and CV(U) calculations are complicated.

Kim and Kim [32] proposed the copula correlation coefficient estimation method by using the method of quasi-Monte Carlo simulation, and the specific steps are as follows:

Step 1. First, the pseudo-Monte Carlo simulation method is used to generate two independent random sequences on [0,1], (u1, ⋯, un), (v1, ⋯, vn).

Step 2. Analog data were used to calculate CU(vn) andCV(un ),

CU(vn)=  C(U,V,α)U|V=vn,CV(un)=C(U,V,α)V|U=un.

Step 3. According to CU(vn) and CV(un), the estimation of the copulas regression function  rV|UC(u),rU|VC(v)can be calculated; specific methods are as follows:

rV|UC(u)r~V|UC(u)=1-1ni=1nCU(vi),rU|VC(v)r~U|VC(v)                =1-1ni= 1nCV(ui).

Step 4. The estimation of the copula correlation coefficient was calculated according to the estimation of the copula regression function.

ρ1|2ρ~1|2=12ni=1n(r~U|VC(ui))2-3,ρ2|1ρ~2|1                      =12ni=1n(r~V|UC(vi))2-3

The above is the method to estimate the copula correlation coefficient. For the asymmetric copula function mentioned in this paper, we can use the above method to estimate its copula correlation coefficient so as to explain the asymmetric correlation between variables.

The following is a static analysis of the correlation between financial assets using the asymmetric copula function that is introduced in Example 3.1. The specific steps are as follows:

Step 1. First, the edge distribution of financial time series is modeled. Since the fluctuation of financial time series is characterized by time variation, aggregation, skew, peak, and thick tail, we can choose ARCH, GARCH, and SV models to describe the distribution characteristics of financial time series to describe the edge distribution.

Step 2. According to the GARCH model selected by the sample data, the original time series is transformed by probability integral, and the marginal distribution is transformed into (0, 1) uniform distribution, and the transformed variables are represented by (u, v).

Step 3. Parameter estimation and copula function selection. First, the maximum likelihood estimation method is used to estimate the parameter r in copula function from three common Archimedean copulas, Gumbel, Clayton, and Frank, and the copula with the best fitting degree is selected by AI Ck.

Step 4. According to the tail distribution characteristics of the copula selected in Step 3, the asymmetric copula function with the same distribution characteristics is selected, and the parameter estimates of the asymmetric copula function are given according to the tail dependence coefficients of the copula selected in Step 3.

Step 5. By calculating the copula function selected in Step 3 and the square Euclidean distance between the asymmetric copula function selected in Step 4 and the empirical copula function, the goodness of fit between the two copula functions was judged.

Step 6. Finally, according to the estimation methods of copula correlation coefficients ρ1|2 and ρ2|1 that have been introduced before, copulas correlation coefficient estimations of the asymmetric copulas are calculated, which were selected in Step 4, thus explaining the asymmetric correlation between financial variables.

Static empirical analysis based on asymmetric copula function

Data selection and statistical description

In order to analyze the relationship between the Shanghai stock market and the Shenzhen stock market, the Shanghai Composite Index and the Shenzhen Composite Index are selected as the research objects (represented by X1, X2 respectively). The stock price is defined as the daily closing price Pt. The time span of sample selection is from 31 July 2015 to 27 September 2021, with a total of 1,500 valid data. The data source is the Sina Securities network.

The return rate of assets Xion the t trading day Rjt = 100(ln pntlnPnt−1). The following is a correlation analysis of the return rates R1t, R2t of financial assets X1, X2. First, the statistical characteristics analysis table of return rates is presented, as shown in Table 1. The data show that the return rates of both the Shanghai Composite Index and the Shenzhen Composite index are right-biased, and both show peaks. The ADF test statistics are significant at the significance level of 0.01, so both sequences are stationary.

TABLE 1
www.frontiersin.org

Table 1. Analysis of statistical characteristics of the return rate.

Selection of edge sequence

Since the conditional distribution of financial time series shows the distribution characteristics of wave aggregation, peak, and thick tail, while t-distribution and GED distribution show the distribution characteristics of the thick tail, the GARCH model can well describe the fluctuation law of financial time series. Therefore, GARCH-T or GARCH-GED can be selected to describe the distribution characteristics of the peak, thick tail, and fluctuation aggregation of financial time series. The GARCH-t and GARCH-GED models were used to simulate and analyze the sample data with R software. The simulation result shows that the GARCH (1, 1)-t model can describe the fluctuation rule of each return rate series well. The specific model is as follows:

         Rnt=εnt,n=1,2,t=1,2,,T         εnt=hnt12ξnt         hnt=Wn+αnεn,t-12+βnhn,t-1(ξ1t,ξ2t)~C(Tv1(ξ1t),Tv2(ξ2t))

where Tv1(·) and Tv2(·) express the normalized t-distribution of mean 0, variance 1, the freedom degrees of v1 and v2, respectively, namely, v1v1-2ξ1t~t(v1),v2v2-2ξ2t~t(v2).

Parameter estimation results of the GARCH model with two edge sequences are shown in Table 2.

TABLE 2
www.frontiersin.org

Table 2. Parameter estimation of the GARCH model.

According to the k-s probability value, it can be seen that, by using the estimated conditional marginal distribution, the sequence after the probability integral transformation of the original financial time series can be considered to obey (0, 1) uniform distribution, indicating that the GARCH (1,1)-t model can better simulate the marginal distribution of various financial series.

Selection of copula model

The following is the correlation analysis of the two sequences obtained after the probability integral transformation of the original financial sequence using Clayton, Gumbel, and Frank Archimedes copulas, respectively. The analysis results are as follows:

The Akaike Information Criterion (AIC) is a standard to measure the excellence of statistical model fitting. The smaller the AIC value is, the better the model is. Generally, the model with the smallest AIC is selected. Bayesian Information Criterion (BIC) refers to BIC that is similar to AIC, which is used for model selection. The smaller the value, the better the model is.

According to AIC and BIC criteria in Table 3, the Gumbel copula is selected to describe the interdependence structure between the Shanghai Composite Index and the Shenzhen Composite index. According to the results of the analysis in 2.3, λU(α)=2-2α, and according to 4-3,α = 0.2857, so λU = 0.7810 , λL = 0, we can say that there is a strong upper tail dependence between the Shanghai Composite Index and the Shenzhen Composite Index.

TABLE 3
www.frontiersin.org

Table 3. Parameter estimation of the copula function.

Selection of asymmetric copula functions

As we know from the previous analysis, the return rates of the Shanghai Composite Index and the Shenzhen Composite Index have a strong upper tail dependence relationship. Here, we select the asymmetric copula function with the distribution characteristics of upper tail dependence for analysis. According to the previous analysis, the asymmetric copula C*(u,v)=φ(12j=12hj1(φ-1(u))hj2(φ-1(v))), where φ(t)=exp(-α(1-t)1r)-e-α1-e-α,r(1,+),h11(x)=x32,h21(x)=2x-x32,h12(x)=x2,h22(x)=2x-x2, namely, c*(u,v)=φ(12((φ1(u))32(φ1(v))2+(2(φ1(u))(φ1(u))32)(2(φ1(v))(φ1(v))2))), and Gumbel copula are identical in form with the tail-dependence function and the tail-dependence coefficient, and, according to the previous discussion, the upper-tail-dependence coefficient of C*(u, v) is λUC*=2-21r,  and it is identical in form with Gumbel copula. According to the analysis results in 4.3, the upper- and lower-tail- dependence coefficients of Gumbel copula are, respectively, λU = 0.7810 , λL = 0, so λUC*and λUG can be estimated by empirical tail-dependence functions and empirical tail-dependence coefficients. Therefore, the value λUG can be used as an estimate of λUC*, namely, 2-21r=0.78,so r^= 3.5.

In order to compare the advantages and disadvantages of the asymmetric copula C*‵, Gumbel copula in describing dependent structures, the square Euclidean distance between asymmetric copula C*, Gumbel copula, and empirical copula is calculated respectively; the result is as follows: (and Ci* is the asymmetric copula function with α = i. For example, C1* indicates α = 1).

This is only the case where α = 1, 2... 10. It can be seen from the Euclidean distance in Table 4 that when α = 5, 6, 7, 8, 9, 10, the Euclidean distance between the asymmetric copula C* and the empirical distribution function is smaller than that between Gumbel and the empirical distribution function, and that, when α = 6 and α = 7, d2 minimizes. Both equal 0.0640. The smaller the Euclidean distance is, the higher the fitting degree of the original data is. Therefore, when describing the interdependent structure of financial variables, asymmetric copula shows higher advantages than Archimedean copula. In addition, the asymmetric copula function can also describe the asymmetric interdependent structure between variables, and more truly reflect the interdependent relationship between variables.

TABLE 4
www.frontiersin.org

Table 4. d2 of C* (u,v) and Gumbel.

According to the above analysis, the upper-tail-dependence coefficient of the asymmetric copula is only related to parameter r. The larger the value of r is, the larger the value of the tail-dependence coefficient will be. Therefore, we can say that parameter R is a parameter representing the strength of tail dependence. In addition, the Euclidean distance between the copula function and the empirical copula function can be changed by changing the parameter α, and the empirical copula can reflect the overall dependence pattern of the data. Therefore, it can be said that the parameter α is the parameter representing the overall dependence structure. The different parameters of asymmetric copula provide convenience for a more accurate description of the correlation between variables.

Copula correlation coefficient of the asymmetric copula

The biggest advantage of using the asymmetric copula function to describe the interdependent structure of financial variables is that the asymmetric interdependent relationship between financial variables can be explained by the copula correlation coefficient.

According to the discussion in 3.3, at first, the stochastic simulation method is used to generate two independent random sequences (u1, ⋯, un), (v1, ⋯, vn),ui, viϵ [0, 1], for some fixed U, and then calculate C*U(vi)=C*(U,V)U|V=vi. Since

C*(U,V)=φ-1(12(φ(u))32(φ(v))2+12(2φ(u)-(φ(u))32)(2φ(v)-(φ(v))2)),We can getC*(U,V)U|V=vi=(34φ12(u)φ2(vi)+12(2-32φ12(u))(2(φ(vi)-φ2(vi))))(-ln(e-α+u-ue-α))r-1exp(-α(1-s)1r)(α(1-s)1r)r-1(e-α+u-ue-α)rV|UC(u)r~V|UC(u)=1-1ni=1nCU(vi)rU|VC(v)r~U|VC*(v)=1-1ni=1nCV(ui),ρ2|1ρ~2|1=12ni=1n(r~V|UC*(vi))2-3

where r = 3.5; the result can be obtained by Matlab,  ρ~2|1=0.3536. This result shows that 35.36% of the changes in the Shenzhen Composite Index can be explained by the copula regression on the Shanghai Composite Index.

The same method should be used to calculate  ρ~2|1. First, for some fixed V, we calculated  CV*(ui)=C*(U,V)V|U=ui,

C*(U,V)V|U=ui=(φ32(ui)φ(v)+ (2φ(ui)-φ32(ui)) (1-φ(vi)))(-ln(e-α+v-ve-α))r-1exp(-α(1-s)1r)(α(1-s)1r)r-1(e-α+v-ve-α)r~U|VC*(v)=1-1ni=1nCV(ui),ρ~1|2=12ni=1n(r~U|VC*(vi))2-3,where r=3.5,φ-1(t)=exp(-α(1-t)1r)-e-α1-e-α,φ(t)=1-(-ln(e-α+t-te-α )α)r.

Using Matlab to calculate available ρ~1|2=0.8450, the results show that changes in the Shanghai index, 84.5% of the part can be explained by the Shenzhen index on the Shenzhen index Copula regression. The different values of ρ˜1|2N˙ρ˜2|1 illustrate the different relationships between the Shanghai index and the Shenzhen index. It provides conditions for correctly describing the dependence between variables.

Time-varying analysis based on the asymmetric copula function

Because the correlation between variables will change with the change of time, scholars have conducted a lot of research on the dynamic copula model, and some results have been obtained. Deng [34] gave the estimation of Var based on a time-varying copula. Wei and Zhang [35] analyzed the correlation mode and dynamic correlation structure among financial markets, respectively, based on the copula theory and the time series model. Wang et al. [36, 37] analyzed the time-varying correlation model and the variable structure correlation model between financial assets based on Spearman's ρ respectively. However, the symmetric Archimedean copula and elliptic copula were used in most previous analyses. In this paper, we used the asymmetric copula function to analyze the time dependence of variables.

The basic idea of time-varying analysis is to determine the evolution process of the population parameter, according to the dynamic evolution process of the dependent indexes and the correspondence between the dependent indexes and the population parameter. According to the previous analysis, the upper-tail-dependence coefficient λU(r)=2-21r, and the lower-tail-dependence coefficient λL = 0 for the asymmetric copula C*. Since there is a one-to-one correspondence between the upper-tail-dependence coefficient and the population parameter r, this paper uses the dynamic evolution process of the upper-tail-dependence coefficient to analyze the time variation of asymmetric copula C*.

Improvement of the dynamic evolution process of the upper-tail-dependence coefficient

Patton [38] was the first to study the time-dependent copula model and proposed that the dynamic evolution process of the parameter ρ of the binary normal copula function could be described by an equation similar to ARMA (1, 10). For tail-dependence coefficients, the commonly used dynamic evolution equation is as follows:

λmU=Λ(γU+βUλm-1U+αUWm)    (5-1)
λmL=Λ(γL+βLλm-1L+αLWm)    (5-2)

In the above equation, the function Λ(·) is the logistic transformation function, defined as Λ(x)=11+e-x. The reason why Λ(x) is defined is to ensure that the upper-tail correlation coefficient and the lower-tail correlation coefficient are in the interval of (0, 1) at any time. The right side of the equation includes an autoregressive term and an exogenous variable Wm. There are two common forms of the exogenous variable at present: one is Wm=1qi=1q|um-j-vm-j|,; another is Wm=1qi=1qum-jvm-j; q means the lagging term, where Wm is the mean of the absolute value of the difference between ui and vi in the lagging q period, and Wm is the mean of the product of dependent random variables ui and vi in the lagging q period. If X and Y are positively correlated, U and V are probability integral transformations of X and Y, so U and V are also positively correlated; therefore, the variable will be the main diagonal on the unit square, because |utvt|, and the point (ut, vt) is proportional to the main diagonal of the shortest distance. Therefore, Wmselects the mean of the absolute value of the difference between ut and vtin the lag q period as an exogenous variable, and because the tail-dependence coefficient describes the dependence relationship of variables at the extreme value, so Wm_m selects the mean of the product of two dependent variables in the lag q period as an exogenous variable.

The dynamic evolution process of the tail dependence coefficients above describes the local evolution process of the copula function model to a certain extent. The autoregressive term on the right-hand side of the equation can explain the sustainability of the tail distribution, but in the process of exogenous variables calculation, only using the lag q period data. These data are not necessarily tail data, and the evolution process of the upper tail and the lower tail is uniform in form, so it is not very good to reflect the tail data influence on the tail-dependence coefficient based on this idea. This paper attempts to improve the above evolution equation into the following form:

λm*U=Λ(γU*+βU*λm-1*U+αu*Wm*U)    (5-3)
λm*L=Λ(γL*+βL*λm-1L+ αL*Wm*L)    (5-4)

The function Λ (·) is still the logistic transformation function, defined as Λ (x)=11+e-x

Wm*U=1qj=1q|um-j-vm-j|I(um-j>34,vm-j>34 ),Wm*L=1qj=1q|um-j-vm-j|I(um-j<14,vm-j<14)

Equation 5-3 represents the dynamic evolution equation of the upper-tail-dependence coefficient, and Equation 5-4 represents the dynamic evolution equation of the lower-tail-dependence coefficient.

Exogenous variables can also be taken in the following form:

Wm*U=1qj=1qum-jvm-jI(um-j>34,vm-j>34)Wm*L=1qj=1qum-jvm-jI(um-j<14,vm-j<14)

In an exogenous variable expression, an indicator function

I(umj>34,vmj>34)={1,    umj>34,vmj>340,    elseI(umj<14,vmj<14)={1,    umj<14,vmj<140,    else

By multiplying the indicator function, we can better characterize the dynamic distribution of the variable in the tail.

Here, we still take asymmetric copula C* as an example to conduct a dynamic analysis of the correlation between financial assets. The analysis steps are as follows:

Step 1. First, the edge distribution of financial time series is modeled.

Step 2. According to the GARCH model selected by the sample data, the original time series is transformed by probability integral, and the marginal distribution is transformed into (0, 1) uniform distribution, and the transformed variables are expressed by (U, V).

Step 3. The statistics of the upper-tail-dependence coefficient should be calculated according to the sample data, and the upper-tail-dependence coefficient sample statistics are as follows:

λU=i=m950I(ui>34,vi>34)i=m950I(ui>34),m=1,2,,900

Step 4. According to the sample data distribution obtained in Step 3, the parameters of the following four dynamic evolution equations are calculated.

λ1m*U=Λ(γ1u*+β1u*λ1m-1*U+α1u*W1m*U)

where W1m*U=1qj=1q|um-j-vm-j|I(um-j>34,vm-j>34 );

λ2m*U=Λ(γ2u*+β2u*λ2m-1*U+α2u*W2m*U)

where W2m*U=1qj=1qum-jvm-jI(um-j>34,vm-j>34 );

λ1mU=Λ(γ1U+β1Uλ1m-1U+α1UW1mU)

where W1mU=1q0j=1q0|um-j-vm-j |;

λ2mU=Λ(γ2U+β2Uλ2m-1U+α2UW2mU)where W2mU=1q0j=1q0(um-jvm-j)

Improved dynamic evolution equation, the exogenous variables of considering only the tail ui>34,vi>34 ; thus, the amount of data involved in calculating reduces; in order to improve the fitting accuracy, the lag period q of the improved dynamic evolution equation of the lag period is chosen to 20, and the lag periods of the original dynamic evolution equation, in general, are chosen to 10.

Step 5. λmU is estimated according to the obtained dynamic evolution equations, and the square root of the mean square error is calculated for the estimated values and sample values of the four evolution equations.

RMSE =(1N(λ^-λ)2)12. The advantages and disadvantages of four kinds of dynamic evolution equations are judged. In this paper, the dynamic evolution equation, which minimizes the square root of the mean square error, is selected for further analysis.

Step 6. Finally, by using the selected dynamic evolution equation, according to the one-to-one correspondence between the upper-tail-dependence coefficient and the population parameter, the estimated value of the population parameter is obtained so as to dynamically analyze the dependence relationship between the Shanghai Composite Index return rate and the Shenzhen Composite Index return rate.

Empirical analysis

The Shanghai Composite Index and the Shenzhen Component Index are still selected as research objects (represented by X1 and X2, respectively). Stock The stock price is defined as the daily closing price Pt,. The time span of sample selection is from 31 July 2015 to 27 September 2021, with a total of 1,500 valid data. The return rate of financial asset Xi on the t trading day Rjt = 100 (ln PntlnPnt−1). According to the previous analysis, the GARCH (1,1) -t model can be used to describe the probability distribution of each return rate series. The GARCH (1,1) -t model can be used to describe the probability distribution of each return rate series. The sequence (ui, vi) is obtained by the probability integral transformation of the original sequence X1 and X2, i = 1, 2, ⋯, 1, 500.

1. The sample sequence (ui, vi) should e selected, i = 1, 2, ⋯, 900, and the statistical value of the upper- tail- dependence coefficient λmU should be calculated by using the statistics defined in Step 3, as shown below:

λU=i=m900I(ui>34,vi>34)i=m900I(ui>34),m= 1,2,,850

where the numerator and the denominator are sliding window data. In the above equation, the maximum value of m is 850. In the definition of 5.1 λU, the maximum value of m is 900. The reason for this choice in the paper is that it is found in the data analysis process that, if the maximum value of m is selected as 850, the advantages of describing the dynamic dependence relationship between variables are more obvious according to the improved dynamic evolution equation. In order to fully demonstrate the advantages of the improved dynamic evolution equation in theory and practice, the maximum value of m is selected as 850 in the empirical analysis.

According to Figure 1, λUincreases with the increase of m, indicating that, for the sample sequence (ui, vi), i = 1, 2, ⋯, 900, the larger the m is, the further the position of the sliding window is, the further the data involved in the calculation are, and the smaller the amount of data is. At this point, when (ui>34),i=m,m+1,, 900, the number of (vi>34) also increases. As shown in Figure 1, the statistical value of λU reaches the maximum when the sliding window is located at m = 700, indicating that, under the conditions of (ui>34),i=700, 701,, 900. The number of data that simultaneously satisfies (vi>34) is largest. Through calculation, it is found that the statistical values of λmU swing around 0.8. In the previous static analysis of the Shanghai Composite Index and the Shenzhen Component Index, the estimated value of the upper-tail-dependence coefficient is 0.7809863, close to 0.8. Therefore, the estimated value of the upper-tail-dependence coefficient in the static analysis is basically consistent with the results of the dynamic analysis.

FIGURE 1
www.frontiersin.org

Figure 1. A dynamic diagram of upper- tail-dependence coefficients statistics.

2. The estimated values of the four model parameters are obtained using the statistics of λmU, as shown in the table below:

3. The four dynamic evolution equations were used to estimate the upper-tail-dependence coefficients, and the square root of the mean square error (RMSE) of the estimated value and the sample value were calculated. RMSEs were expressed as d1*,d2*,d1,d2, respectively. The calculated results are shown in the following table:

It can be seen from the table that the RMSEs of the two improved dynamic evolution equations are both smaller than that of the original dynamic evolution equation, indicating that the improved dynamic evolution equation can more accurately depict the dynamic evolution process of the upper-tail-dependence coefficient compared with the original dynamic evolution equation.

In fact, according to the values of (ui, vi), it can be seen that some pairs of (ui, vi) do not belong to the upper tail data. It is obviously not appropriate to calculate the upper-tail-dependence coefficients directly based on all values of (ui, vi). It is theoretically more feasible to select the upper-tail data by multiplying the function of expression to calculate the upper-tail-dependence coefficients. According to Table 9, d1* is 3 × 10−5smaller than d1, and d2* is 2 × 10−5 smaller than d2,; the difference is relatively small. The reasons for the small difference are as follows: (1) By multiplying the exogenous variable by an explicit function, the value of the exogenous variable will be smaller than that of the original exogenous variable; however, according to Tables 58 it can be seen that the absolute values of the coefficient α1u* and  α2u* of the exogenous variable in the improved dynamic evolution equation increase. Therefore, after multiplying the improved exogenous variable and its coefficient, the larger and smaller parts will cancel each other. There is not much difference between di* and di(i = 1, 2) on the surface. (2) The difference between di* and di(i = 1, 2) becomes smaller due to a large amount of data. Although the difference between di* and di is not large, theoretically, the improved dynamic evolution equation is more reasonable, and the value of di*  is smaller, and the dynamic evolution process of the tail-dependence coefficient can be described more accurately.

TABLE 5
www.frontiersin.org

Table 5. Parameters estimates of the time-varying model for the upper-tail-dependence coefficients  λ1m*U.

TABLE 6
www.frontiersin.org

Table 6. Parameters estimates of the time-varying model for the upper-tail-dependence coefficients  λ2m*U.

TABLE 7
www.frontiersin.org

Table 7. Parameters estimates of the time-varying model for the upper-tail-dependence coefficients  λ1mU.

TABLE 8
www.frontiersin.org

Table 8. Parameters estimates of the time-varying model for the upper-tail-dependence coefficients  λ2mU.

In addition, among the two improved dynamic evolution equations, the RMSE corresponding to the exogenous variable W1m *U is smaller than that corresponding to the exogenous variable W2m*U, indicating that it is more appropriate to choose the exogenous variable W1m*U. However, according to Table 9, the difference between d1* and d2* is very small. It shows that there is no significant difference between the two types of exogenous variables in the dynamic analysis of the Shanghai Stock index and the Shenzhen Stock index.

TABLE 9
www.frontiersin.org

Table 9. RMSE values of estimated values and sample values.

4. The selected dynamic evolution equation should be used to calculate the estimated value of λuand draw the dynamic evolution diagram of the estimated value and the sample value, as shown in the figure below:

As can be seen from Figure 2, the estimated value and the sample value calculated by using the dynamic evolution equation are highly consistent in both trends and fluctuation, which fully demonstrates that dynamic evolution Equation 5-3 can accurately depict the dynamic evolution process of the upper-tail-dependence coefficient.

FIGURE 2
www.frontiersin.org

Figure 2. A dynamic diagram of statistics and estimates of upper-tail-dependence coefficients.

5. According to the correspondence between the upper-tail-dependence coefficient λuand the population parameter r, λu(r)=2-21r, the dynamic estimate of the population parameter r is obtained, as shown below Figure 3.

FIGURE 3
www.frontiersin.org

Figure 3. A dynamic diagram of the population parameters estimation.

Conclusion

In the correlation analysis of financial assets, symmetric Archimedes copula and elliptic copula functions are mostly used for analysis, and the dependence indicators between variables are also mostly symmetric, which does not conform to the relationship characteristics between financial assets. Based on this, this paper tries to conduct static and dynamic analyses on the correlation of financial assets based on the asymmetric copula. In the static analysis, through the comparative analysis of the asymmetric copula function and the Gumbel copula function, it is found that the asymmetric copula function has more advantages in depicting tail dependence between financial assets. Furthermore, this paper also introduced the concept of the asymmetric copulas correlation coefficient. For the asymmetric copula function, copula correlation coefficients ρ1|2 and ρ2|1 are not equal; and asymmetric dependence between variables can be explained by the copula correlation coefficient, which also shows the great advantages of the asymmetric copula. In the dynamic analysis, the traditional dynamic evolution equation is improved through the empirical analysis, which shows that the improved dynamic evolution equation more accurately portrays the tail-dependence dynamic evolution process, since there is a one-to-one correspondence between the upper-tail-dependence coefficient λu and the population parameter; according to the dynamic evolution equations of λu, the population parameter can be estimated. Thus, the interdependence between financial assets is analyzed asymmetrically and dynamically.

Through empirical analysis, it can be seen that asymmetric copula has more advantages in describing the dependence relationship between financial assets compared with the commonly used Archimedean copula. The asymmetric copula is more suitable for the dependence structure between financial variables. However, for asymmetric copula function, it is difficult to estimate the parameters due to the complex form. In this paper, the Archimedean copula is first used to analyze the correlation between the Shanghai Stock index and the Shenzhen Stock index, and the empirical dependence coefficient of the tail is used to estimate the parameters of the asymmetric copula. Therefore, the asymmetric analysis is developed on the basis of the symmetric analysis. Although the asymmetric copula function is superior to the symmetric copula in application, it is difficult to make further progress in asymmetric copula analysis without some conclusions obtained from symmetric copula analysis. The analysis of asymmetric dependence between variables is one of the advantages of the asymmetric copula, but the calculation of the asymmetric dependence coefficient defined on the basis of the asymmetric copula is more difficult. In this paper, the method of Monte Carlo simulation is used to estimate the dependence coefficient. Therefore, how to estimate the parameters and calculate the dependence coefficient of asymmetric copula more easily is the problem to be studied in the following paper.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fams.2022.1005956/full#supplementary-material

References

1. Romano C. Calibrating and Simulating Copula Functions: An Application to the Italian Stock Market. Working paper, CIDEM. (2002).

Google Scholar

2. Rosenberg JV, Schuermann T. A general approach to integrated risk management with skewed, fat-tailed risks. J Financ Econ. (2006) 79:569–614. doi: 10.1016/j.jfineco.2005.03.001

CrossRef Full Text | Google Scholar

3. Patton AJ. Modelling asymmetric exchange rate dependence. Int Econ Rev. (2006) 47:527–56. doi: 10.1111/j.1468-2354.2006.00387.x

CrossRef Full Text | Google Scholar

4. Patton AJ. Application of copula theory in financial econometrics. Doctor Thesis, San Diego: Department of Economics, University of California, San Diego. (2002).

Google Scholar

5. Palaro HP, Hotta LK. Using conditional copula to estimate value at risk. J Data Sci. (2006) 4:93–115. doi: 10.6339/JDS.2006.04(1).226

CrossRef Full Text | Google Scholar

6. Nikolai K, Anjos U, Mendes BVM. Copulas: a review and recent developments. Stochastic Models. (2006) (4):617–60. doi: 10.1080/15326340600878206

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Zhang YT. Copula technology and financial risk analysis. Stat Res. (2002) 4:48–51. doi: 10.19343/j.cnki.11-1302/c.2002.04.011

CrossRef Full Text

8. Zhang MH. Research on the Copula measurement method of Value at risk of multiple financial assets. J Quant Tech Econ. (2004) 21:67–70. doi: 10.13653/j.cnki.jqte.2004.04.011

CrossRef Full Text

9. Wu ZX, Chen M, Ye WY, Miao BQ. Analysis of foreign exchange portfolio risk based on Copula. Chin J Manage Sci. (2004) 12:1–5. doi: 10.16381/j.cnki.issn:1003-207x.2004.04.001

CrossRef Full Text

10. Wang PJ, Wen GH, Huang TW, Yu WW, Lü YZ. Asymptotical neuro-adaptive consensus of multi-agent systems with a high dimensional leader and directed switching topology. In: IEEE Transactions on Neural Networks and Learning Systems. (2022). doi: 10.1109/TNNLS.2022.3156279

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wu ZX, Chen M, Ye WY, Miao BQ. Portfolio risk analysis based on Copula GARCH. Syst Eng Theor Pract. (2006) 26:45–52. doi: 10.3321/j.issn:1000-6788.2006.03.007

CrossRef Full Text

12. Wei YH, Zhang SY. Correlation analysis of financial market: copula GARCH model and its application. Syst Eng. (2004) 22:7–12. doi: 10.3969/j.issn:1001-4098.2004.04.002

CrossRef Full Text

13. Wei YH, Zhang SY. Research on dynamic correlation structure of financial market. J Syst Eng. (2006) 21:313–7. doi: 10.3969/j.issn:1001-45781.2006.03.015

CrossRef Full Text

14. Zhang JQ, Li X. Integrated risk measurement of asset portfolio and its application: research on VaR based on optimal fitting copula function. Syst Eng Theor Pract. (2008) 6:23–6. doi: 10.3321/j.issn:1000-6788.2008.06.002

CrossRef Full Text

15. Wang PJ, Wen GH, Huang TW Yu WW, Ren Y. Observer-based consensus protocol for directed switching networks with a leader of nonzero inputs. IEEE Trans Cybern. (2022) 52:630–40. doi: 10.1109/TCYB.2020.2981518

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Li M, Liu XB, He H, Yao, P. Measurement of portfolio VaR based on clayton Copula-GARCH model. Pract Understand Mathem. (2013) 6:51–4. doi: 10.3969/j.issn:1000-0984.2013.06.017

CrossRef Full Text

17. Du ZP, Zhang XF. Comparative study of two kinds of mixed rattan copula model: based on VaR of foreign exchange asset portfolio. Technl Econs Manag Res. (2014) 1:27–31.

18. Ren XL, Zhang SY. Analysis of investment risk based on kernel estimation and multiple archimedes copula. Manage Sci. (2007) 20:92–97. doi: 10.3969/j.issn:1672-0334.2007.05.012

CrossRef Full Text

19. Clayton DG. A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence. Biometrika. (1978) 65(1):141–151. doi: 10.1093/biomet/65.1.141

CrossRef Full Text | Google Scholar

20. Frank MJ. On the simultaneous associativity of F(x,y) and x+y-F(x,y). Aequ Math. (1979) 19:194–226. doi: 10.1007/BF02189866

CrossRef Full Text | Google Scholar

21. Gumbel EJ. Distributions des valeurs extrêmes en plusiers dimensions. Publ. Inst. Statist. Univ. Paris. (1960) 9:171–173.

Google Scholar

22. Joe H. Multivariate Models and Dependence Concepts. London: Chapman & Hall. (1997). doi: 10.1201/9780367803896

CrossRef Full Text

23. Alfonse AE, Brigo D. New families of Copulas based on periodic function. Comn Statist Theory Methods. (2005) 34:1437–47. doi: 10.1081/STA-200063351

CrossRef Full Text | Google Scholar

24. Khoudraji A. Contributions à l' étude des Copulas et à la modélisation des valeurs extrêmes bivaviées. Ph.D Thesis, Canada: Université Lavab Quebec. (1995).

Google Scholar

25. Liebscher E. Construction of asymmetric multivariate Copulas. J Multivar Anal. (2008) 99:2234–50. doi: 10.1016/j.jmva.2008.02.025

CrossRef Full Text | Google Scholar

26. Nelsen RB. An Introduction to Copulas. New York: Springer. (2006).

Google Scholar

27. Sungur EA. Some observations on Copula regression functions. Stat Theory Methods. (2005) 34:1967–78. doi: 10.1080/03610920500201244

CrossRef Full Text | Google Scholar

28. Dabrowska D. Regression-based orderings and measures of stochastic dependence. Math Oper Forsch Stat Ser Stat. (1981) 12:317–25. doi: 10.1080/02331888108801592

CrossRef Full Text | Google Scholar

29. Dette H, Siburg KF, Stoimenov PA, A Copula-based non-parametric measure of regression dependence. Scand J Stat. (2013) 40:21–41. doi: 10.1111/j.1467-9469.2011.00767.x

CrossRef Full Text | Google Scholar

30. Shih JH, Emura T. On the Copula Correlation ratio and its generalization. J Multivar Anal. (2021) 182:104708. doi: 10.1016/j.jmva.2020.104708

CrossRef Full Text | Google Scholar

31. Darsow WF, Nguyen B, Olsen ET. Copulas and Markov processes. Illinois J Math. (1992) 36:600–42. doi: 10.1215/ijm/1255987328

CrossRef Full Text | Google Scholar

32. Kim DY, Kim JM. Analysis of directional dependence using asymmetric Copula-based regression models. J Stat Comput Simul. (2014) 84:1990–2010. doi: 10.1080/00949655.2013.779696

CrossRef Full Text | Google Scholar

33. Sungur EA. A note on directional dependence in regression setting. Comm Statist Theory Methods. (2005) 34:1957–65. doi: 10.1080/03610920500201228

CrossRef Full Text | Google Scholar

34. Deng GM. The VaR Estimating on Time-varying Copula. Syst Eng. (2007) 25:28–33. doi: 10.3969/j.issn:1001-4098.2007.08.005

CrossRef Full Text

35. Wei YH, Zhang SY. Copula Theory and its Application in Financial Analysis. Beijing: Qing dynasty Hua University Press. (2008).

36. Wang Q, Wang L, Cheng SJ. Analysis of the Interdependence of Shanghai and Shenzhen stock markets based on time-varying Copula Model. Statist Decis. (2010) 19:139–41. doi: 10.13546/j.cnki.tjyjc.2010.19.002

CrossRef Full Text

37. Wang Q, Wang L, He P. Spearman based time-varying Copula module Simulation and Application of type 2. Mathem Statist Manage. (2011) 1:76–84. doi: 10.13860/j.cnki.sltj.2011.01.007

CrossRef Full Text

38. Patton AJ. Modeling time-varying exchange rate dependence using conditional copula. Working Paper of Department of Economics, University of California, San Diego. (2001). doi: 10.2139/ssrn.275591

CrossRef Full Text | Google Scholar

Keywords: asymmetric copula function, copula correlation coefficient, upper tail dependence coefficient, time-varying analysis, dynamic evolution process

Citation: Li X and Hou B (2022) Correlation analysis of financial assets based on asymmetric copula. Front. Appl. Math. Stat. 8:1005956. doi: 10.3389/fams.2022.1005956

Received: 28 July 2022; Accepted: 25 August 2022;
Published: 29 September 2022.

Edited by:

Peijun Wang, Anhui Normal University, China

Reviewed by:

Liang Wang, Yunnan Normal University, China
Fengyun Zhang, Jining University, China

Copyright © 2022 Li and Hou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xia Li, lxhoub@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.