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.

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.

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.

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

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.

MULTIVARIATE EXTREME VALUE DISTRIBUTION AND ITS FISHER INFORMATION MATRIX

The paper is concerned with the basic properties of multivariate extreme value distribution (in the Logistic model). We obtain the characteristic function and recurrence formula of the density function. The explicit algebraic formula for Fisher information matrix is indicated. A simple and accurate procedure for generating random vector from multivariate extreme value distribution is presented.

General Regular Variation of the n-th Order and 2nd Order Edgeworth Expansions of the Extreme Value Distribution （Ⅱ）

In this part II the fundamental inequality of the third order general regular variation is proved and the second order Edgeworth expansion of the distribution of the extreme values is discussed.

Stationary distributions for two-dimensional sticky Brownian motions:Exact tail asymptotics and extreme value distributions

Sticky Brownian motions can be viewed as time-changed semimartingale reflecting Brownian motions,which find applications in many areas including queueing theory and mathematical finance.In this paper,we focus on stationary distributions for sticky Brownian motions.Main results obtained here include tail asymptotic properties in the marginal distributions and joint distributions.The kernel method,copula concept and extreme value theory are the main tools used in our analysis.

Theory of multivariate compound extreme value distribution and its application to extreme sea state prediction

In this paper, a new type of distribution, multivariate compound extreme value distribution (MCEVD), is introduced by compounding a discrete distribution with a multivariate continuous distribution of extreme sea events. In its engineering application the number over certain threshold level per year is fitting to Poisson distribution and the corresponding extreme sea events are fitting to Nested Logistic distribution, then the Poisson-Nested logistic trivariate compound extreme value distribution (PNLTCED) is proposed to predict extreme wave heights, periods and wind speeds in Yellow Sea. The new model gives more stable and reasonable predicted results.

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.

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.

The Pickands' Estimator of the Negative Extreme-value Index

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

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.

Quantile Trends in Temperature Extremes in China

A number of recent studies have examined trends in extreme temperature indices using a linear regression model based on ordinary least-squares. In this study, quantile regression was, for the first time, applied to examine the trends not only in the mean but also in all parts of the distribution of several extreme temperature indices in China for the period 1960–2008. For China as a whole, the slopes in almost all the quantiles of the distribution showed a notable increase in the numbers of warm days and warm nights, and a significant decrease in the number of cool nights. These changes became much faster as the quantile increased. However, although the number of cool days exhibited a significant decrease in the mean trend estimated by classical linear regression, there was no obvious trend in the upper and lower quantiles. This finding suggests that examining the trends in different parts of the distribution of the time-series is of great importance. The spatial distribution of the trend in the 90 th quantile indicated that there was a pronounced increase in the numbers of warm days and warm nights, and a decrease in the number of cool nights for most of China, but especially in the northern and western parts of China, while there was no significant change for the number of cool days at almost all the stations.

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.