Next Article in Journal
Design and Analysis of Automated Film Roll Cutter
Previous Article in Journal
Patient Health Monitoring System Using IOT And AI
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Assessing Asymptotic Tail Independence: A Simulation Study †

Centro de Matemática, Universidade do Minho, 4710-057 Braga, Portugal
Presented at the 10th International Conference on Time Series and Forecasting, Gran Canaria, Spain, 15–17 July 2024.
Eng. Proc. 2024, 68(1), 52; https://doi.org/10.3390/engproc2024068052
Published: 18 July 2024
(This article belongs to the Proceedings of The 10th International Conference on Time Series and Forecasting)

Abstract

:
The occurrence of extreme values in one variable can trigger the same in other variables, making it necessary to assess the risk of contagion. The usual dependence measures based on the central part of the data typically fail to assess extreme dependence. Within the scope of EVT, tail dependence measures were developed, such as the Ledford and Tawn coefficient that we discuss here. This is a measure of residual dependence that is particularly important when it comes to analyzing at the tail level where data are scarce. We will consider different estimation methodologies and compare them based on a simulation study. We finish with an application to real data.

1. Introduction

Applications of EVT in the analysis of risky events are diverse, such as finance, insurance, material resistance, quality control, telecommunications, sports, environment, hydrology, biology, and seismology, among others.
The occurrence of extreme observations can have serious consequences, e.g., high water levels in floods, a vast area burned in wildfires, and losses in stock indexes in financial market crashes, among others. Sometimes, the occurrence of extreme values for one variable can contaminate other variables. Typical examples are the financial markets where the fall of one index leads to falls in others as well. Correlation is a measure of dependence between two variables widely used in applications. Although it works well to measure the association between two variables in the central part of the data, it fails to assess dependence within the tails. For instance, in a bivariate Normal with correlation ρ < 1 , if we are far enough in the tail, extreme values seem to occur independently at each margin, regardless of how high a correlation we choose [1].
Extremal dependence in a random pair ( X , Y ) is typically based on the chance for obtaining large values of both variables. It is helpful to remove the influence of marginal aspects by transforming the marginals to a common distribution function (df), e.g., F X ( X ) and F Y ( Y ) , both of which are standard uniform where F X and F Y are the marginal df of X and Y, respectively, and considered continuous. Therefore, differences in distributions are solely due to dependence aspects. As F X ( X ) and F Y ( Y ) are on a common scale, events of the form { F X ( X ) > u } and { F Y ( Y ) > u } , for large values of u, are equally extreme events for each variable, both with their probability approaching zero as u 1 . The tail dependence coefficient (TDC) measures the conditional probability of one variable being extreme, given that the other is extreme too:
λ = lim u 1 P ( F X ( X ) > u | F Y ( Y ) > u )
Observe that the exact independence corresponds to λ = 0 , since P ( F X ( X ) > u | F Y ( Y ) > u ) = P ( F X ( X ) > u ) 0 , as u 1 , whilst perfect dependence means P ( F X ( X ) > u | F Y ( Y ) > u ) = 1 . Indeed, we say that X and Y are tail independent if λ = 0 , and tail dependent if λ > 0 . A positive λ means a strong dependence persisting to the limit. A null λ comprehends exact independence, but also a weak dependence that gradually vanishes as the limit is approached. Such residual dependency can be captured through the rate of convergence of P ( F X ( X ) > u | F Y ( Y ) > u ) towards zero; more precisely, through the asymptotic tail independent coefficient of Ledford and Tawn η ( 0 , 1 ] [2], given by
P ( F X ( X ) > u | F Y ( Y ) > u ) ( 1 u ) 1 / η 1 L ( 1 / ( 1 u ) ) , u 1
where L ( x ) is a slowly varying function, i.e., L ( t x ) / L ( x ) 1 , as x , t > 0 [3]. If η = 1 and L c , 0 < c 1 , the random variables (rvs) are asymptotically dependent. On the other hand, if η < 1 or η =1 and L 0 , then λ = 0 and the rv are asymptotically independent with degree η . The case η = 1 / 2 means near-independence; in particular, exact independence if L 1 . We can also state that 1 / 2 < η < 1 corresponds to a positive association, in the sense that the probability of both rvs exceeding u occurs more frequently than under exact independence, whereas in 0 < η < 1 / 2 , the probability of both rvs exceeding u occur less frequently than under exact independence and, thus, the rvs are negatively associated. Inference on η allows us to distinguish the type of tail (in)dependence: η = 1 means tail dependence, whilst η < 1 means the presence of a residual tail dependence. According to [2], assuming independence and ignoring η can lead to misspecified joint extremes estimations, while considering the dependence between rvs, when a residual dependence in fact occurs, may result in an overestimation.
Observe that
P ( F X ( X ) > u , F Y ( Y ) > u ) = P ( U > ( 1 u ) 1 , V > ( 1 u ) 1 )
where U = ( 1 F X ( X ) ) 1 and V = ( 1 F Y ( Y ) ) 1 are standard Pareto. Considering rv T = min ( U , V ) and t u = ( 1 u ) 1 , by (1) we can state that
P ( T > t u ) t u 1 / η L ( t u ) , t u
and we say that T has a regularly varying tail with index 1 / η . In EVT, this means that η corresponds to the tail index of T, the primary measure in EVT inference. There are many estimation methods in the literature for the tail index. For a survey, see, e.g., [4]. A well-known tail index estimator is the Hill [5]. Given a sample T 1 , , T n , consider
η ^ u ( H ) : = 1 n u j = 1 n u log ( T ( j ) / t u ) ,
where T ( j ) , j = 1 , , n u , are the rvs exceeding the threshold t u . In the Hill estimator, t u corresponds to the n u + 1 up order statistics. In EVT, we can also estimate the tail index through a modeling approach. A common technique is to consider the exceedance values above a high threshold of data, and apply a Generalized Pareto Distribution (GPD) according to the Pickands–Balkema–De Haan theorem [6,7], known as the Peaks Over Threshold (POT) approach (see, e.g., [4] for more details). Here, we apply the POT method on rv T in order to estimate η , and denote this estimator as η ^ u ( G P D ) .
Based on (1), we can rewrite
η η ( u ) log ( 1 u ) log P ( F X ( X ) > u , F Y ( Y ) > u )
and, thus, we can estimate η by taking the empirical counterpart of (3), which we denote as η ^ ( e m p ) . More precisely,
η ^ u ( e m p ) log ( 1 u ) log ( n 1 i = 1 n 1 { F X ( X i ) > u , F Y ( Y i ) > u } ) .
Note that, in practice, the margins F X and F Y are usually unknown, and we can estimate by the empirical df. This is the approach followed for estimators η ^ u ( H ) , η ^ u ( G P D ) and η ^ u ( e m p ) . Observe also that the standard uniform rvs, F X ( X i ) and F Y ( Y i ) , i = 1 , , n , have r-order statistics Beta( r , n r + 1 ) distributed. Thus, replacing the marginal empirical df in (4) with Beta( r , n r + 1 ) df, we have the fourth estimator, which we denote as η ^ u ( b e t a ) . Further details on Beta marginals estimation can be seen in [8].
In Section 2, we analyze the referred estimators through simulation and compare their performance. An application to real data illustrates the methods in Section 3.

2. Estimation of η : A Simulation Study

In our study, we consider different models with different values of η , as follows, where we use copula notation, C ( u , v ) = P ( F X ( X ) < u , F Y ( Y ) < v ) :
  • Bivariate Ali–Mikhail–Haq distribution, where C ( u , v ) = u v ( 1 ξ ( 1 u ) ( 1 v ) ) 1 , with ξ = 1 ( η = 1 / 3 ) (see, e.g., [9]), denoted AMH;
  • Bivariate Frank distribution, where C ( u , v ) = ξ 1 log ( 1 ( 1 e ξ u ) ( 1 e ξ v ) / ( 1 e ξ ) ) , with ξ = 0.5 ( η = 1 / 2 ) (see, e.g., [9]), denoted Frank;
  • Bivariate Normal distribution with ρ = 0.5 ( η = 3 / 4 ) (see, e.g., [10]), denoted BNormal;
  • Bivariate extreme value distribution with a Logistic dependence function, C ( u , v ) = exp ( ( log u ) 1 / α + ( log v ) 1 / α ) α , 0 < α 1 , with α = 0.75 ( η = 1 ) (see, e.g., [2]), denoted Logistic.
We simulate 1000 replicas of each model, with size n = 100 and n = 1000 . We compute estimators η ^ u ( H ) , η ^ u ( G P D ) , η ^ u ( e m p ) and η ^ u ( b e t a ) for thresholds u corresponding to percentiles 0.5 ( 0.01 ) 0.99 . The estimation is conducted using package mev [11] within statistical software R 2020 [12]. We derive the sample mean of estimates in each threshold u along with the Wald 95% confidence interval, which are plotted in Figure 1, Figure 2, Figure 3 and Figure 4, in the case of samples with size n = 1000 , and Figure 5, Figure 6, Figure 7 and Figure 8 for n = 100 . We also calculate the root mean squared error (RMSE) in each threshold u, available in Figure 9 and Figure 10, for sample sizes n = 1000 and n = 100 , respectively. The η ^ u ( G P D ) estimator is based on modeling threshold exceedances through the GPD model and, therefore, fails to be applied at very high thresholds, where the number of observations is small.
In the case n = 1000 , and looking at the estimated means, estimators η ^ u ( H ) and η ^ u ( G P D ) have the best performance. However, concerning the RMSE, η ^ u ( G P D ) presents the worst performance, except in the Logistic model, where η ^ u ( H ) and η ^ u ( G P D ) are the best choices and η ^ u ( e m p ) and η ^ u ( b e t a ) are not recommended. In smaller sample size n = 100 , the confidence intervals become larger and the results are less precise. The RMSE grows slightly for all estimators, except for η ^ u ( G P D ) , where it clearly increases. The low performance of all estimators is in the Logistic model, where η = 1 corresponds to a boundary value.

3. Application

We are going to analyze the tail dependence between two stock market indexes pairs, CAC40–DAX and CAC40–PSI20, from January 2020 to February 2024, available from https://finance.yahoo.com/world-indices/ (accessed on 1 March 2024). We consider the daily close log-returns of each index, plotted in Figure 11. In order to remove heteroskedasticity, we apply an asymmetric GARCH filter. The scatter-plot in Figure 12 represents the 1037 pairs of residuals, where the association between the variables is evident, both in the central part and in the tail. The estimates of η are plotted in Figure 13 and Figure 14. The value η = 1 seems plausible for the CAC40–DAX pair, particularly in estimators η ^ u ( H ) and η ^ u ( G P D ) , leading to a diagnostic of tail dependence. Regarding CAC40–PSI20, the η estimate is pointed out to be 0.7 , meaning asymptotic tail independence. Therefore, the risk of contagion for the occurrence of extreme values is higher between the CAC40 and DAX indices, when compared with the CAC40–PSI20 pair.

4. Conclusions

In this work, we addressed the inference in tail dependence between two variables, through the Ledford and Tawn coefficient η [2]. This is an important parameter in EVT, as it captures the presence of a possible residual dependence that, once ignored, leads to misspecifications. In the context of EVT, where data are scarce and one needs to extrapolate beyond observations, accuracy in estimations is a particularly sensitive topic. There are different methodologies in the literature for estimating η . In this work, we consider four widely used estimators, and compare their performance in different models. The results obtained constitute an initial guide on the best options for choosing estimators in practical applications. In future work, we intend to continue this study, with more models and variations in η values, and thus improve guidance for users in choosing the best estimation method.

Funding

The research of the author was partially financed by Portuguese Funds through FCT (Fundação para a Ciência e a Tecnologia) within the Projects UIDB/00013/2020 and UIDP/00013/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original data presented in the study are openly available at https://finance.yahoo.com/world-indices/.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Embrechts, P.; McNeil, A.; Straumann, D. Correlation and dependency in risk management: Properties and pitfalls. In Risk Management: Value at Risk and Beyond; Dempster, M.A.H., Ed.; Cambridge University Press: Cambridge, UK, 2002; pp. 176–223. [Google Scholar]
  2. Ledford, A.W.; Tawn, J.A. Statistics for near independence in multivariate extreme values. Biometrika 1996, 83, 169–187. [Google Scholar] [CrossRef]
  3. Gnedenko, B.V. Sur la distribution limite du terme maximum d’une série aléatoire. Ann. Math. 1943, 44, 423–453. [Google Scholar] [CrossRef]
  4. Beirlant, J.; Goegebeur, Y.; Segers, J.; Teugels, J. Statistics of Extremes: Theory and Applications; Wiley: Chichester, UK, 2004. [Google Scholar]
  5. Hill, B. A simple general approach to inference about the tail of a distribution. Ann. Stat. 1975, 3, 1163–1174. [Google Scholar] [CrossRef]
  6. Balkema, A.A.; de Haan, L. Residual life time at great age. Ann. Probab. 1974, 2, 792–804. [Google Scholar] [CrossRef]
  7. Pickands, J., III. Statistical inference using extreme order statistics. Ann. Stat. 1975, 3, 119–131. [Google Scholar]
  8. Segers, J.; Sibuya, M.; Tsukahara, H. The empirical beta copula. J. Multivar. Anal. 2017, 155, 35–51. [Google Scholar] [CrossRef]
  9. Goegebeur, Y.; Guillou, A. Asymptotically unbiased estimation of the coefficient of tail dependence. Scand. J. Stat. 2013, 40, 174–189. [Google Scholar] [CrossRef]
  10. Draisma, G.; Drees, H.; Ferreira, A.; Haan, L.D. Bivariate tail estimation: Dependence in asymptotic independence. Bernoulli 2004, 10, 251–280. [Google Scholar] [CrossRef]
  11. Belzile, L.; Wadsworth, J.L.; Northrop, P.J.; Grimshaw, S.D.; Huser, R. mev: Multivariate Extreme Value Distributions. R Package Version 1.13.1. 2020. Available online: https://CRAN.R-project.org/package=mev (accessed on 1 March 2024).
  12. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
Figure 1. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model AMH where η = 1 / 3 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 1. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model AMH where η = 1 / 3 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g001
Figure 2. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model Frank where η = 1 / 2 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 2. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model Frank where η = 1 / 2 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g002
Figure 3. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model BNormal where η = 3 / 4 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 3. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model BNormal where η = 3 / 4 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g003
Figure 4. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model Logistic where η = 1 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 4. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model Logistic where η = 1 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g004
Figure 5. Estimated sample means (solid line) for each threshold u, in samples of size n = 100 , along with Wald 95% confidence interval (dashed line), in model AMH where η = 1 / 3 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 5. Estimated sample means (solid line) for each threshold u, in samples of size n = 100 , along with Wald 95% confidence interval (dashed line), in model AMH where η = 1 / 3 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g005
Figure 6. Estimated sample means (solid line) for each threshold u, in samples of size n = 100 , along with Wald 95% confidence interval (dashed line), in model Frank where η = 1 / 2 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 6. Estimated sample means (solid line) for each threshold u, in samples of size n = 100 , along with Wald 95% confidence interval (dashed line), in model Frank where η = 1 / 2 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g006
Figure 7. Estimated sample means (solid line) for each threshold u, in samples of size n = 100 , along with Wald 95% confidence interval (dashed line), in model BNormal where η = 3 / 4 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 7. Estimated sample means (solid line) for each threshold u, in samples of size n = 100 , along with Wald 95% confidence interval (dashed line), in model BNormal where η = 3 / 4 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g007
Figure 8. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model Logistic where η = 1 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Figure 8. Estimated sample means (solid line) for each threshold u, in samples of size n = 1000 , along with Wald 95% confidence interval (dashed line), in model Logistic where η = 1 (horizontal line), of η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right).
Engproc 68 00052 g008
Figure 9. Estimated RMSE for each threshold u, in samples of size n = 1000 , for models η ^ u ( H ) (dotted line), η ^ u ( G P D ) (dot-dash line), η ^ u ( e m p ) (solid line) and η ^ u ( b e t a ) (dashed line), for models AMH (top-left), Frank (top-right), BNormal (bottom-left) and Logistic (bottom-right).
Figure 9. Estimated RMSE for each threshold u, in samples of size n = 1000 , for models η ^ u ( H ) (dotted line), η ^ u ( G P D ) (dot-dash line), η ^ u ( e m p ) (solid line) and η ^ u ( b e t a ) (dashed line), for models AMH (top-left), Frank (top-right), BNormal (bottom-left) and Logistic (bottom-right).
Engproc 68 00052 g009
Figure 10. Estimated RMSE for each threshold u, in samples of size n = 100 , for models η ^ u ( H ) (dotted line), η ^ u ( G P D ) (dot-dash line), η ^ u ( e m p ) (solid line) and η ^ u ( b e t a ) (dashed line), for models AMH (top-left), Frank (top-right), BNormal (bottom-left) and Logistic (bottom-right).
Figure 10. Estimated RMSE for each threshold u, in samples of size n = 100 , for models η ^ u ( H ) (dotted line), η ^ u ( G P D ) (dot-dash line), η ^ u ( e m p ) (solid line) and η ^ u ( b e t a ) (dashed line), for models AMH (top-left), Frank (top-right), BNormal (bottom-left) and Logistic (bottom-right).
Engproc 68 00052 g010
Figure 11. Daily close log-returns of CAC40 (left), DAX (middle) and PSI20 (right) indexes, in the period between January 2020 and February 2024 .
Figure 11. Daily close log-returns of CAC40 (left), DAX (middle) and PSI20 (right) indexes, in the period between January 2020 and February 2024 .
Engproc 68 00052 g011
Figure 12. Scatter-plot of filtered daily close log-returns: CAC40–DAX (left) and CAC40–PSI20 (right).
Figure 12. Scatter-plot of filtered daily close log-returns: CAC40–DAX (left) and CAC40–PSI20 (right).
Engproc 68 00052 g012
Figure 13. CAC40–DAX: Estimates given by η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right) for thresholds u between 0.5 and 0.99 and respective 95% Wald confidence intervals.
Figure 13. CAC40–DAX: Estimates given by η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right) for thresholds u between 0.5 and 0.99 and respective 95% Wald confidence intervals.
Engproc 68 00052 g013
Figure 14. CAC40–PSI20: Estimates given by η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right) for thresholds u between 0.5 and 0.99 and respective 95% Wald confidence intervals.
Figure 14. CAC40–PSI20: Estimates given by η ^ u ( H ) (top-left), η ^ u ( G P D ) (top-right), η ^ u ( e m p ) (bottom-left) and η ^ u ( b e t a ) (bottom-right) for thresholds u between 0.5 and 0.99 and respective 95% Wald confidence intervals.
Engproc 68 00052 g014
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ferreira, M. Assessing Asymptotic Tail Independence: A Simulation Study. Eng. Proc. 2024, 68, 52. https://doi.org/10.3390/engproc2024068052

AMA Style

Ferreira M. Assessing Asymptotic Tail Independence: A Simulation Study. Engineering Proceedings. 2024; 68(1):52. https://doi.org/10.3390/engproc2024068052

Chicago/Turabian Style

Ferreira, Marta. 2024. "Assessing Asymptotic Tail Independence: A Simulation Study" Engineering Proceedings 68, no. 1: 52. https://doi.org/10.3390/engproc2024068052

APA Style

Ferreira, M. (2024). Assessing Asymptotic Tail Independence: A Simulation Study. Engineering Proceedings, 68(1), 52. https://doi.org/10.3390/engproc2024068052

Article Metrics

Back to TopTop