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.

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.

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.

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

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

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.

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(?.

More general results on mixed extreme value distributions

The sequences {Zi,n, 1≤i≤n}, n≥1 are multi-nomial distribution among i.i.d, random variables {X1,i, i≥1}, {X2,i, i≥1 } {Xm,i, i≥1 }. The extreme value distribution Gz（x） of this particular triangular array of i.i,d, random variables Z1,n, Z2 n,...,Zn,n is discussed. A new type of not max-stable extreme value distributions which are Fréchet mixture, Gumbel mixture and Weibull mixture has been found if Fj,…… Fm belong to the same MDA. Whether mixtures of different types of extreme value distributions exist or not and the more general case are discussed in this paper. We found that Gz（x） does not exist as mixture forms of the different types of extreme value distributions after we investigated all cases.

A class of not max-stable extreme value distributions

The sequences {Zi , 1≤i≤n}, n≥1 have multi-nomial distribution among i.i.d. random variables {X1, , i≥1}, {X2, , ,n i i i≥1}, …, {Xm , i≥1}. The extreme value distribution GZ(x) of this particular triangular array of i.i.d. random variables Z1, , Z2, , …, ,i n n r ?1 Zn is discussed in this paper. We found a new type of not max-stable extreme value distributions, i) GZ (x) = ,n ∏Φα Ai(x)×Φαr (x); i i=1 r ?1 r?1 ii) GZ (x) = ∏Ψα Ai(x)×Ψαr (x); iii) GZ (x) = ∏Λ Ai(λix)×Λ(x), r≥2, 0<α1≤α2≤…≤αr and λi∈(0,1] for i, 1≤i≤r?1 which occur if i i=1 i=1 Fj, …, Fm belong to the same MDA.

Limiting Distribution of Extreme Values for FGM Random Sequences

This paper mainly study extreme values of FGM random sequences.We prove a technique theorem by the dependence structure of FGM sequences,and further obtain the limiting distributions of maxima and k-th largest for stationary FGM random sequences.

Applications of Bootstrap in Analyzing General Extreme Value Distributions

The bootstrap method is one of the new ways of studying statistical math which this article uses but is a major tool for studying and evaluating the values of parameters in probability distribution.Our research is concerned overview of the theory of infinite distribution functions.The tool to deal with the problems raised in the paper is the mathematical methods of random analysis(theory of random process and multivariate statistics).In this article,we introduce the new function to find out the bias and standard error with jackknife method for Generalized Extreme Value distributions.

The Pickands' Estimator of the Negative Extreme-value Index

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

The Convergence Rate of Fréchet Distribution under the Second-Order Regular Variation Condition

In this article we consider the asymptotic behavior of extreme distribution with the extreme value index γ>0 . The rates of uniform convergence for Fréchet distribution are constructed under the second-order regular variation condition.

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.

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.

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.

Red Tide Frequency Analysis Using the Extreme Value Theory

In recent years, the red tide erupted frequently, and caused a great economic loss. At present, most literatures emphasize the academic research on the growth mechanism of red tide alga. In order to find out the characters of red tide in detail and improve the precision of forecast, this paper gives some new approaches to dealing with the red tide. By the extreme values, we deal with the red tide frequency analysis and get the estimation of T-times red tide level U (T), which is the level once the consistence of red tide alga exceeds on the average in a period of T times.

Analysis of Network Traffic with Extreme Value Theory

It is very im portant to analyze network traffic in the network control and management. In thi s paper, extreme value theory is first introduced and a model with threshold met hods is proposed to analyze the characteristics of network traffic. In this mode l, only some traffic data that is greater than threshold value is considered. Th en the proposed model with the trace is simulated by using S Plus software. The modeling results show the network traffic model constructed from the extreme va lue theory fits well with that of empirical distribution. Finally, the extreme v alue model with the FARIMA(p,d,q) modeling is compared. The anal ytical results illustrate that extreme value theory has a good application foreg round in the statistic analysis of network traffic. In addition, since only some traffic data which is greater than the threshold is processed, the computation overhead is reduced greatly.

General Regular Variation of n-th Order and the 2nd Order Edgeworth Expansion of the Extreme Value Distribution (Ⅰ)

In Part Ⅰ the concept of the general regular variation of n-th order is proposed and its construction is discussed. The uniqueness of the standard expression and the higher order regularity of the auxiliary functions are proved.

A Newly-Discovered GPD-GEV Relationship Together with Comparing Their Models of Extreme Precipitation in Summer

It has been theoretically proven that at a high threshold an approximate expression for a quantile of GEV （Generalized Extreme Values） distribution can be derived from GPD （Generalized Pareto Distribution）. Afterwards, a quantile of extreme rainfall events in a certain return period is found using L-moment estimation and extreme rainfall events simulated by GPD and GEV, with all aspects of their results compared. Numerical simulations show that POT （Peaks Over Threshold）-based GPD is advantageous in its simple operation and subjected to practically no effect of the sample size of the primitive series, producing steady high-precision fittings in the whole field of values （including the high-end heavy tailed）. In comparison, BM （Block Maximum）-based GEV is limited, to some extent, to the probability and quantile simulation, thereby showing that GPD is an extension of GEV, the former being of greater utility and higher significance to climate research compared to the latter.

A New Method for Determining Threshold in Using PGCEVD to Calculate Return Values of Typhoon Wave Height

In using the PGCEVD （Poisson-Gumbel Compound Extreme Value Distribution） model to calculate return values of typhoon wave height, the quantitative selection of the threshold has blocked its application. By analyzing the principle of the threshold selection of PGCEVD model and in combination of the change point statistical methods, this paper proposes a new method for quantitative calculation of the threshold in PGCEVD model. Eleven samples from five engineering points in several coastal waters of Guangdong and Hainan, China, are calculated and analyzed by using PGCEVD model and the traditional Pearson type III distribution （P-III） model, respectively. By comparing the results of the two models, it is shown that the new method of selecting the optimal threshold is feasible. PGCEVD model has more stable results than that of P-III model and can be used for the return wave height in every direction.