Extreme value theory applied to the auroral electrojet indices

The study of extreme weather and space events has gained paramount importance in modern society owing to rapid advances in high technology.Understanding and describing exceptional occurrences plays a crucial role in making decisive assessments of their potential impact on technical,economic,and social aspects in various fields.This research focuses on analyzing the hourly values of the auroral electrojet(AE)geomagnetic index from 1957 to 2019 by using the peak over threshold method in extreme value theory.By fitting the generalized Pareto distribution to extreme AE values,shape parameter indices were derived,revealing negative values that establish an upper bound for this time series.Consequently,it became evident that the AE values had reached a plateau,suggesting that extreme events exceeding the established upper limit are rare.As a result,although the need for diligent precautions to mitigate the consequences of such extreme events persists,surpassing the upper limit of AE values becomes increasingly challenging.It is also possible to observe an aurora in the middle-and low-latitude regions during the maximum period of the AE index.

Bivariate Analysis of Pollutants Monthly Maxima in Mexico City Using Extreme Value Distributions and Copula

In the present work, we are interested in studying the joint distributions of pairs of the monthly maxima of the pollutants used by the environmental authorities in Mexico City to classify the air quality in the metropolitan area. In order to obtain the joint distributions a copula will be considered. Since we are analyzing the monthly maxima, the extreme value distributions of Weibull and Fréchet are taken into account. Using these two distributions as marginal distributions in the copula a Bayesian inference was made in order to estimate the parameters of both distributions and also the association parameters appearing in the copula model. The pollutants taken into account are ozone, nitrogen dioxide, sulphur dioxide, carbon monoxide, and particulate matter with diameters smaller than 10 and 2.5 microns obtained from the Mexico City monitoring network. The estimation was performed by taking samples of the parameters generated through a Markov chain Monte Carlo algorithm implemented using the software OpenBugs. Once the algorithm is implemented it is applied to the pairs of pollutants where one of the coordinates of the pair is ozone and the other varies on the set of the remaining pollutants. Depending on the pollutant and the region where they were collected, different results were obtained. Hence, in some cases we have that the best model is that where we have a Fréchet distribution as the marginal distribution for the measurements of both pollutants and in others the most suitable model is the one assuming a Fréchet for ozone and a Weibull for the other pollutant. Results show that, in the present case, the estimated association parameter is a good representation to the correlation parameters between the pair of pollutants analyzed. Additionally, it is a straightforward task to obtain these correlation parameters from the corresponding association parameters.

Analyzing of the ENSO Index Using Extreme Value Theory

We predicted the extreme values of the ENSO index, the Niño3.4 index, and the Southern Oscillation Index (SOI) using extreme value theory. Various diagnostic plots for assessing the accuracy of the Generalized Pareto (GP) model fitted to the Niño3.4 index and SOI are shown, and all four diagnostic plots support the fitted GP model. Because the shape parameter of the Niño3.4 was negative, the Niño3.4 index had a finite upper limit. In contrast, that of the SOI was zero, therefore the SOI did not have a finite upper limit, and there is a possibility that a significant risk will occur. We predicted the maximum return level for the return periods of 10, 20, 50, 100, 350, and 500 years and their respective 95% confidence intervals, CI. The 10-year, and 100-year return levels for Niño3.4 were estimated to be 2.41, and 2.62, with 95% CI [2.22, 2.59], and [2.58, 2.66], respectively. The Niño3.4 index was 2.65 in the 2015/16 super El Niño, which is a phenomenon that occurs once every 500 years. The Niño3.4 index was 2.51 in the 1982/83, and 1997/98 super El Niño, which is a phenomenon that occurs once every 20 years. Recently, a large super El Niño event with a small probability of occurrence has occurred. In response to global warming, the super El Niño events are becoming more likely to occur.

The Pickands' Estimator of the Negative Extreme-value Index

提出一类极值指数为负时的相似于Pickand’s型的新的极值指数估计量 。

Analysis of the USD/JPY and EUR/JPY Exchange Rates Using Multifractal Analysis and Extreme Value Theory

We performed a multifractal analysis using wavelet transform to detect the changes in the fractality of the USD/JPY and EUR/JPY exchange rates, and predicted their extreme values using extreme value theory. After the 1997 Asian financial crisis, the USD/JPY and EUR/JPY became multifractal, then the USD/JPY became monofractal and stable, and yen depreciation was observed. However, the EUR/JPY became multifractal and unstable, and a strong yen depreciation was observed. The coherence between the USD/JPY and EUR/JPY was strong between 1995 and 2000. After the 2007-2008 financial crisis, the USD/JPY became monofractal and stable, and yen appreciation was observed. However, the EUR/JPY became multifractal and unstable, and strong yen appreciation was observed. Various diagnostic plots for assessing the accuracy of the GP model fitted to USD/JPY and EUR/JPY are shown, and all the diagnostic plots support the fitted GP model. The shape parameters of USD/JPY and EUR/JPY were close to zero, therefore the USD/JPY and EUR/JPY did not have finite upper limits. We predicted the maximum return level for the return periods of 10, 20, 50, 100, 350, and 500 years and their respective 95% confidence intervals (CI). As a result, the 10-year and 100-year return levels for USD/JPY were estimated to be 149.6 and 164.8, with 95% CI [143.2, 156.0] and [149.4, 180.1], respectively.

The Maximum and Minimum Value of Exponential RandićIndices of Quasi-Tree Graph

The exponential Randić index has important applications in the fields of biology and chemistry. The exponential Randić index of a graph G is defined as the sum of the weights e 1 d( u )d( v ) of all edges uv of G, where d( u ) denotes the degree of a vertex u in G. The paper mainly provides the upper and lower bounds of the exponential Randić index in quasi-tree graphs, and characterizes the extremal graphs when the bounds are achieved.

Projected Changes in Temperature and Precipitation Extremes over China as Measured by 50-yr Return Values and Periods Based on a CMIP5 Ensemble

Future changes in the 50-yr return level for temperature and precipitation extremes over China's Mainland are investigated based on a CMIP5 multi-model ensemble for RCP2.6, RCP4.5 and RCP8.5 scenarios. The following indices are analyzed： TXx and TNn （the annual maximum and minimum of daily annual maximum consecutive 5-day precipitation） and CDD maximum and minimum surface temperature）, RX5day （the （maximum annual number of consecutive dry days）. After first validating the model performance, future changes in the 50-yr return values and return periods for these indices are investigated along with the inter-model spread. Multi-model median changes show an increase in the 50-yr return values of TXx and a decrease for TNn, more specifically, by the end of the 21st century under RCP8.5, the present day 50-yr return period of warm events is reduced to 1.2 yr, while extreme cold events over the country are projected to essentially disappear. A general increase in RX5day 50-yr return values is found in the future. By the end of the 21st century under RCP8.5, events of the present RX5day 50-yr return period are projected to reduce to 〈 10 yr over most of China. Changes in CDD-50 show a dipole pattern over China, with a decrease in the values and longer return periods in the north, and vice versa in the south. Our study also highlights the need for further improvements in the representation of extreme events in climate models to assess the future risks and engineering design related to large-scale infrastructure in China.

Extreme value distribution and reliability of nonlinear stochastic structures

A new approach to evaluate the extreme value distribution （EVD） of the response and reliability of general multi-DOF nonlinear stochastic structures is proposed. The approach is based on the recently developed probability density evolution method, which enables the instantaneous probability density functions of the stochastic responses to be captured. In the proposed method, a virtual stochastic process is first constructed to satisfy the condition that the extreme value of the response equals the value of the constructed process at a certain instant of time. The probability density evolution method is then applied to evaluate the instantaneous probability density function of the response, yielding the EVD. The reliability is therefore available through a simple integration over the safe domain. A numerical algorithm is developed using the Number Theoretical Method to select the discretized representative points. Further, a hyper-ball is imposed to sieve the points from the preceding point set in the hypercube. In the numerical examples, the EVD of random variables is evaluated and compared with the analytical solution. A frame structure is analyzed to capture the EVD of the response and the dynamic reliability. The investigations indicate that the proposed approach provides reasonable accuracy and efficiency.

Analysis of Japan and World Records in the 100 m Dash Using Extreme Value Theory

Extreme value theory provides methods to analyze the most extreme parts of data. We predicted the ultimate 100 m dash records for men and women for specific periods using the generalized extreme value (GEV) distribution. The various diagnostic plots, which assessed the accuracy of the GEV model, were well fitted to the 100 m records in the world and Japan, validating the model. The men’s world record had a shape parameter of -0.250 with a 95% confidence interval of [-0.391, -0.109]. The 100 m record had a finite limit and a calculated upper limit was 9.46 s. The return level estimates for the men’s world record were 9.74, 9.62, and 9.58 s with a 95% confidence interval of [9.69, 9.79], [9.54, 9.69], and [9.48, 9.67] for 10-, 100- and 350-year return periods, respectively. In one year, the probability of occurrence for a new world record of men, 9.58 s (Usain Bolt), was 1/350, while that for women, 10.49 s (Florence Griffith-Joyner), was about 1/100, confirming it was more difficult for men to break records than women.

Importance of Generalized Logistic Distribution in Extreme Value Modeling

We consider a problem from stock market modeling, precisely, choice of adequate distribution of modeling extremal behavior of stock market data. Generalized extreme value (GEV) distribution and generalized Pareto (GP) distribution are the classical distributions for this problem. However, from 2004, [1] and many other researchers have been empirically showing that generalized logistic (GL) distribution is a better model than GEV and GP distributions in modeling extreme movement of stock market data. In this paper, we show that these results are not accidental. We prove the theoretical importance of GL distribution in extreme value modeling. For proving this, we introduce a general multivariate limit theorem and deduce some important multivariate theorems in probability as special cases. By using the theorem, we derive a limit theorem in extreme value theory, where GL distribution plays central role instead of GEV distribution. The proof of this result is parallel to the proof of classical extremal types theorem, in the sense that, it possess important characteristic in classical extreme value theory, for e.g. distributional property, stability, convergence and multivariate extension etc.

Analyzing the Annual Maximum Magnitude of Earthquakes in Japan by Extreme Value Theory

One of the most important and interesting issues associated with the earthquakes is the long-term trend of the extreme events. Extreme value theory provides methods for analysis of the most extreme parts of data. We estimated the annual maximum magnitude of earthquakes in Japan by extreme value theory using earthquake data between 1900 and 2019. Generalized extreme value (GEV) distribution was applied to fit the extreme indices. The distribution was used to estimate the probability of extreme values in specified time periods. The various diagnostic plots for assessing the accuracy of the GEV model fitted to the magnitude of maximum earthquakes data in Japan gave the validity of the GEV model. The extreme value index, ξ was evaluated as −0.163, with a 95% confidence interval of [−0.260, −0.0174] by the use of profile likelihood. Hence, the annual maximum magnitude of earthquakes has a finite upper limit. We obtained the maximum return level for the return periods of 10, 20, 50, 100 and 500 years along with their respective 95% confidence interval. Further, to get a more accurate confidence interval, we estimated the profile log-likelihood. The return level estimate was obtained as 7.83, 8.60 and 8.99, with a 95% confidence interval of [7.67, 8.06], [8.32, 9.21] and [8.61, 10.0] for the 10-, 100- and 500-year return periods, respectively. Hence, the 2011 off the Pacific coast of Tohoku Earthquake, which was the largest in the observation history of Japan, had a magnitude of 9.0, and it was a phenomenon that occurs once every 500 year.

Extreme value distributions of mixing two sequences with different MDA＇s

Suppose {Xi, i≥l } and {Yi, i≥l } are two independent sequences with distribution functions Fx(x) and Fr(x), respectively. Zi,n is the combination of Xi and Y1 with a probability pn for each i with 1≤i≤<_n. The extreme value distribution Gz(x) of this particular triangular array of the i.i.d, random variables Z1,n, Z2.n ,..., Zn.n is discussed. We found a new form of the extreme value distribution A^A(px)A(x)(0＜p＜1), which is not max-stable. It occurs if Fx(x) and Fr(x) belong to the same MDA(A). Gz(x) does not exist as mixture forms of the different types of extreme value distributions.

Extreme value distributions of mixing two sequences with the same MDA

Suppose {Xi, i 1} and {Yi, i 1} are two independent sequences with distribution functions ()XFx and ()YFx, respectively. Zi,n is the combination of Xi and Yi with a probability np for each i with 1 in. The extreme value distribution GZ(x) of this particular triangular array of the i.i.d. random variables Z1,n, Z2,n, ,L Zn,n is discussed. We found a new form of the extreme value distributions i) 12()()AxxaaFF and ii) 12()()AxxaaYY (a1<a2), which are not max-stable. It occurs if FX and FY belong to the same MDA(? or MDA(?.

Empirical Analysis of Value-at-Risk Estimation Methods Using Extreme Value Theory

This paper investigates methods of value-at-risk (VaR) estimation using extreme value theory (EVT). It compares two different estimation methods, 'two-step subsample bootstrap' based on moment estimation and maximum likelihood estimation (MLE), according to their theoretical bases and computation procedures. Then, the estimation results are analyzed together with those of normal method and empirical method. The empirical research of foreign exchange data shows that the EVT methods have good characters in estimating VaR under extreme conditions and 'two-step subsample bootstrap' method is preferable to MLE.

Network Traffic Based on GARCH-M Model and Extreme Value Theory

GARCH-M ( generalized autoregressive conditional heteroskedasticity in the mean) model is used to analyse the volatility clustering phenomenon in mobile communication network traffic. Normal distribution, t distribution and generalized Pareto distribution assumptions are adopted re- spectively to simulate the random component in the model. The demonstration of the quantile of network traffic series indicates that common GARCH-M model can partially deal with the "fat tail" problem. However, the "fat tail" characteristic of the random component directly affects the accura- cy of the calculation. Even t distribution is based on the assumption for all the data. On the other hand, extreme value theory, which only concentrates on the tail distribution, can provide more ac- curate result for high quantiles. The best result is obtained based on the generalized Pareto distribu- tion assumption for the random component in the GARCH-M model.

Estimation of Poisson-Generalized Pareto Compound Extreme Value Distribution by Probability-Weighted Moments and Empirical Analysis

This paper puts forward a Poisson-generalized Pareto (Poisson-GP) distribution. This new form of compound extreme value distribution expands the existing application of compound extreme value distribution, and can be applied to predicting financial risk, large insurance settlement and high-grade earthquake, etc. Compared with the maximum likelihood estimation (MLE) and compound moment estimation (CME), probability-weighted moment estimation (PWME) is used to estimate the parameters of the distribution function. The specific formulas are presented. Through Monte Carlo simulation with sample sizes 10, 20, 50, 100, 1 000, it is concluded that PWME is an efficient method and it behaves steadily. The mean square errors (MSE) of estimators by PWME are much smaller than those of estimators by CME, and there is no significant difference between PWME and MLE. Finally, an example of foreign exchange rate is given. For Dollar/Pound exchange rates from 1990-01-02 to 2006-12-29, this paper formulates the distribution function of the largest loss among the investment losses exceeding a certain threshold by Poisson-GP compound extreme value distribution, and obtains predictive values at different confidence levels.

Estimation of Dynamic VaR in Chinese Stock Markets Based on Time Scale and Extreme Value Theory

The accuracy and time scale invariance of value-at-risk （VaR） measurement methods for different stock indices and at different confidence levels are tested. Extreme value theory （EVT） is applied to model the extreme tail of standardized residual series of daily/weekly indices losses, and parametric and nonparametric methods are used to estimate parameters of the general Pareto distribution （GPD）, and dynamic VaR for indices of three stock markets in China. The accuracy and time scale invariance of risk measurement methods through back-testing approach are also examined. Results show that not all the indices accept time scale invariance; there are some differences in accuracy between different indices at various confidence levels. The most powerful dynamic VaR estimation methods are EVT-GJR-Hill at 97.5% level for weekly loss to Shanghai stock market, and EVT-GARCH-MLE （Hill） at 99.0% level for weekly loss to Taiwan and Hong Kong stock markets, respectively.

Evaluation of the prognostic nutritional index for the prognosis of Chinese patients with high/extremely high-risk prostate cancer after radical prostatectomy

BACKGROUND The incidence of prostate cancer(PCa)is on the rise in China.The risk level of patients with PCa is associated with disease-free survival rate at 10 years after radical prostatectomy.Predicting prognosis in advance according to the degree of risk can provide a reference for patients,especially treatment options and postoperative adjuvant treatment measures for high-risk/extremely high-risk patients.AIM To explore the predictive value of the prognostic nutritional index(PNI)for biological recurrence in Chinese patients with high/extremely high-risk PCa after radical prostatectomy.METHODS The biochemical test results and clinical data of 193 patients who underwent radical prostatectomy for the first time from January 2015 to December 2020 were retrospectively collected.The PNI value of peripheral blood within 1 wk before surgery was calculated,and during the follow-up period,prostate-specific antigen≥0.2 ng/mL was considered to have biological recurrence.The receiver operating characteristic(ROC)curve was used to calculate the optimal critical value and area under the curve(AUC)of the patients.According to the critical value,the progression-free survival of the high PNI group and low PNI group was compared.The independent influencing factors of the patients’prognosis were obtained by the Cox proportional hazards regression model.RESULTS The non-biological recurrence rates at 1,3,and 5 years were 92.02%,84.05%,and 74.85%,respectively.The optimal critical value for PNI to predict biological recurrence was 46.23,and the AUC was 0.789(95%confidence interval:0.651-0.860;P<0.001).The sensitivity and specificity were 82.93%and 62.30%,respectively.In accordance with the optimal critical value of the ROC curve(46.23),193 patients were further divided into a high PNI group(PNI≤46.23,n=108)and low PNI group(PNI>46.23,n=85).The incidence of postoperative complications in the high PNI group was lower than that in the low PNI group(21.18%vs 38.96%).Kaplan-Meier survival analysis showed that the overall survival rate at 5 years in the low PNI group was 87.96%(13/108),which was lower than that in the high PNI group(61.18%,33/85;P<0.05).Low PNI[hazard ratio(HR)=1.74;P=0.003]and positive incisal margin status(HR=2.14;P=0.001)were independent predictors of biological recurrence in patients with high/extremely high-risk PCa.CONCLUSION The PNI has predictive value for the prognosis of patients with high/extremely high-risk PCa,and is an independent prognostic factor.Patients with low PNI value have a shorter time of nonbiological recurrence after prostatectomy.It is expected that the combined prediction of other clinicopathological data will further improve the accuracy and guide postoperative adjuvant therapy to improve the quality of prognosis.

Extreme Values Approach in Food Risk Modeling

In this paper we estimate the incubation period of a possible pathology following exposure to dioxins during a poor diet. The tools developed for this purpose include the probabilistic extremal model and the stochastic behavior of the distribution tails of the contamination. We propose a cumulative distribution function for a random variable that follows both a Gaussian distribution and a GPD. A global optimization method is also explored for the efficient estimation of parameters of GPD.

Analysis of Design Rainstorm Hydrograph Based on Asymmetrical Extreme Value Copula

In this paper,the maximum 1-hour rainfall( rain peak),the maximum 6-hour rainfall and the maximum 24-hour rainfall in the Caojiang River basin from 1967 to 2013 were taken as samples. The typical typhoon rainstorm hydrograph of joint distribution of rainfall in three periods was constructed based on the asymmetric Archimedean Gumbel-Hougaard extreme value Copula. The main conclusions were as follows:( 1) the design rainstorm value in the Caojiang River basin calculated by using the joint distribution of rainfall in three periods was larger than the design rainstorm value of the joint distribution in two periods and that of a single period. The design rainstorm process hydrograph amplified at the same frequency had the optimal overall effect,which provided a new idea and method for studying the design rainfall patterns.( 2) According to the maximum 24-hour rainfall,the risk rate of the multi-peak rainstorm process that the main peak was in the back was the highest,and the constructed typical design rainstorm process hydrograph was the most representative.( 3) " OR" joint return period of rainfall combination in three periods as the design criteria of a watershed was applicable to responding to the risk of rainfall and flood in this watershed.