Robust Estimation of the Generalized Pareto Distribution

One approach used for analyzing extremes is to fit the excesses over a high threshold by a generalized Pareto distribution. For the estimation of the shape and scale parameters in the generalized Pareto distribution, under some restrictions on the value of the scale parameter, maximum likelihood, method of moments and probability weighted moments' estimators are available. However, these are not robust estimators. In this paper we implement a robust estimation procedure known as the method of medians (He and Fung, 1999) to estimate the parameters in the generalized Pareto distribution. The asymptotic distribution of our estimator is normal for any value of the shape parameter except –1.     

�׶�II����������ֵķ��յ�����ļ��ο���ͼ

In recent years, statistical process control (SPC) has been widely used to monitor the performance of clinical practitioners, such as surgeons and general practitioners. In this paper, two risk-adjusted geometric control charts namely cumulative sum (CUSUM) and weighted likelihood ratio test (WLRT) are proposed to monitor surgery performance in phase II. The performance of the proposed control charts is evaluated and compared by simulation experiments for different shift values in the parameters of a risk-adjusted logistic regression model in terms of the average run length (ARL) criterion. The results show that all methods work well in the sense that they can effectively detect shifts in the process parameters.     

Economic design of an attribute control chart for monitoring mean life based on multiple deferred state sampling

In this paper, we design an attribute np control chart using multiple deferred state (MDS) sampling under Weibull distribution based on time truncated life test. This chart is constructed for monitoring the variation of mean life of the product in a manufacturing process. The optimal parameters of MDS sampling and the control limit coefficients are determined so that the in‐control average run length (ARL) is as close as to the target ARL. The optimal parameters of MDS sampling are sample size and number of successive subgroups required for declaring the current state of process. Out‐of‐control ARL is considered as a measure of the performance of proposed chart and reported with determined optimal parameters for various shift constants. The out‐of‐control ARL of the proposed chart obtained under various distributions is compared with each other. The performance of proposed control chart is compared with the performance of the existing control chart designed under single sampling. In addition, the economic design of proposed chart using variable sampling interval scheme is discussed, and sensitivity analysis on expected costs is also investigated.     

Estimation of Parameters of the Generalized Pareto\Distribution by the Least Squares

Traditional estimations of parameters of the generalized Pareto distribution (GPD) are generally constrained by the shape parameter of GPD. Such as: the method-of-moments (MOM), the probability-weighted moments (PWM), L-moments (LM), the maximum likelihood estimation (MLE) and so on. In this paper we use the fact that GPD can be transformed into the exponential distribution and use the results of parameters estimation for the exponential distribution, than we propose parameters estimators of the two-parameter or three-parameter GPD by the least squares method. Some asymptotic results are provided and the proposed method not constrained by the shape parameter of GPD. A simulation study is carried out to evaluate the performance of the proposed method and to compare them with other methods suggested in this paper. The simulation results indicate that the proposed method performs better than others in some common situation.     

Improved threshold diagnostic plots for extreme value analyses

A crucial aspect of threshold-based extreme value analyses is the level at which the threshold is set. For a suitably high threshold asymptotic theory suggests that threshold excesses may be modelled by a generalized Pareto distribution. A common threshold diagnostic is a plot of estimates of the generalized Pareto shape parameter over a range of thresholds. The aim is to select the lowest threshold above which the estimates are judged to be approximately constant, taking into account sampling variability summarized by pointwise confidence intervals. This approach doesn’t test directly the hypothesis that the underlying shape parameter is constant above a given threshold, but requires the user subjectively to combine information from many dependent estimates and confidence intervals. We develop tests of this hypothesis based on a multiple-threshold penultimate model that generalizes a two-threshold model proposed recently. One variant uses only the model fits from the traditional parameter stability plot. This is particularly beneficial when many datasets are analysed and enables assessment of the properties of the test on simulated data. We assess and illustrate these tests on river flow rate data and 72 series of significant wave heights.     

Generalized fiducial confidence intervals for extremes

The generalized Pareto distribution is relevant to many situations when modeling extremes of random variables. In particular, peaks over threshold data approximately follow the generalized Pareto distribution. We use a fiducial framework to perform inference on the parameters and the extreme quantiles of the generalized Pareto. This inference technique is demonstrated both when the threshold is a known and unknown parameter. Assuming the threshold is a known parameter resulted in fiducial intervals with good empirical properties and asymptotically correct coverage. Likewise, our simulation results suggest that the fiducial intervals and point estimates compare favorably to the competing methods seen in the literature. The proposed intervals for the extreme quantiles when the threshold is unknown also have good empirical properties regardless of the underlying distribution of the data. Comparisons to a similar Bayesian method suggest that the fiducial intervals have better coverage and are similar in length with fewer assumptions. In addition to simulation results, the proposed method is applied to a data set from the NASDAQ 100. The data set is analyzed using the fiducial approach and its competitors for both cases when the threshold is known and unknown. R code for our procedure can be downloaded at .     

联合检测过程位置参数和尺度参数的非参数Cusum控制图

用于检测生产过程的多数传统控制图都假定过程的受控分布是已知的,并假定数据服从正态分布。然而在很多情况下,由于没有足够的数据来估计过程的分布,这种假定就变得不现实,而非参数控制图却不需要任何关于分布的特殊形式的假定。另外,多数的已有控制图都是使用两个单独的均值与方差控制图来同时检测生产过程.本文中,我们提出一个新的基于Cramer-von-Mises(CvM)检验的非参数累积和控制图(称为CvM图)来同时检测过程位置参数和尺度参数。文中给出了基于不同受控平均运行长度(ARL)下的CvM图的控制限,通过步长的均值、方差及分位数来研究控制图的性能表现。最后用一个实例来说明CvM图的实际应用。     

监测Lévy稳定过程均值变点的平均运行时间估计

平均运行长度(时间)(ARL)是判断一控制图监测变点效果好坏的一个重要工具.本文主要研究Lévy稳定过程的均值变点监测问题.我们给出了三个控制图,即EWMA,GEWMA和GLR的ARL估计,并通过数值模拟比较了4个控制图监测均值变点的效果和差异.     

监测L\'{e}vy稳定过程均值变点的平均运行时间估计

平均运行长度(时间) (ARL)是判断一控制图监测变点效果好坏的一个重要工具\bd 本文主要研究L\'{e}vy 稳定过程的均值变点监测问题\bd 我们给出了三个控制图, 即EWMA, GEWMA和GLR的ARL 估计, 并通过数值模拟比较了4个控制图监测均值变点的效果和差异.     

ARL计算方法综述

质量控制图的研究在统计过程控制(Statistical Process Control,简称SPC)中占有很重要的位置,且在实际中得到了很好的应用,也取得了很好的经济效益和社会效益。而ARL(AverageRun Length)和ATS(Average Time to Signal)分别是评价各种静态和动态控制图性质好坏的一个常用的关键指标。到现在为止,关于各种控制图ARL或ATS的计算大致有三种方法:马氏链法、积分方程法和随机模拟,本文将针对各种不同的控制图,对上述前两种方法做一个综述性的总结,以方便研究者和实际工作参考使用。     

Generalizing the Pareto to the log-Pareto model and statistical inference

In this article we introduce a full-fledged statistical model of log-Pareto distribution functions (dfs) parametrized by two shape parameters and a scale parameter. Pareto dfs can be regained in the limit by varying parameters of log-Pareto dfs, whence the log-Pareto model can be regarded as an extension of the Pareto model. Log-Pareto dfs are first of all obtained by means of exponential transformations of Pareto dfs. We also indicate an iterated application of such a procedure. A class of generalized log-Pareto dfs is considered as well. In addition, power-pot (p-pot) stable dfs – related to p-max stable dfs – are introduced and log-Pareto dfs are identified as special cases. A modification of a quick (systematic) estimator is proposed as an initial estimator for the numerical computation of the maximum likelihood estimator (MLE) in the 3-parameter model.      

高质量过程控制技术的比较研究与分析

随着质量改进活动的不断开展,现代制造过程中的不合格项在逐渐降低。在这种情况下,常规的休哈特型计数控制图往往是失效的。为了监控高质量的过程运行,一种方法是采用累计合格品计数(CCC)图;另一种方法是采用几何Q图,本文首先分析了这两种控制图的基本原理;进而以平均链长(ARL)和探测过程发生漂移的概率为准则,系统分析和比较了这两种控制图的性能,仿真结果表明,在大多数情况下,这两种控制图具有相似的性能;最后,通过实例说明了这两种控制图的应用,并给出了若干建议。     

An evaluation of the multivariate dispersion charts with estimated parameters under non‐normality

Various charts such as |S|, W, and G are used for monitoring process dispersion. Most of these charts are based on the normality assumption, while exact distribution of the control statistic is unknown, and thus limiting distribution of control statistic is employed which is applicable for large sample sizes. In practice, the normality assumption of distribution might be violated, while it is not always possible to collect large sample size. Furthermore, to use control charts in practice, the in‐control state usually has to be estimated. Such estimation has a negative effect on the performance of control chart. Non‐parametric bootstrap control charts can be considered as an alternative when the distribution is unknown or a collection of large sample size is not possible or the process parameters are estimated from a Phase I data set. In this paper, non‐parametric bootstrap multivariate control charts |S|, W, and G are introduced, and their performances are compared against Shewhart‐type control charts. The proposed method is based on bootstrapping the data used for estimating the in‐control state. Simulation results show satisfactory performance for the bootstrap control charts. Ultimately, the proposed control charts are applied to a real case study.     

Bayesian approaches for analyzing earthquake catastrophic risk

Extreme value theory has been widely used in analyzing catastrophic risk. The theory mentioned that the generalized Pareto distribution (GPD) could be used to estimate the limiting distribution of the excess value over a certain threshold; thus the tail behaviors are analyzed. However, the central behavior is important because it may affect the estimation of model parameters in GPD, and the evaluation of catastrophic insurance premiums also depends on the central behavior. This paper proposes four mixture models to model earthquake catastrophic loss and proposes Bayesian approaches to estimate the unknown parameters and the threshold in these mixture models. MCMC methods are used to calculate the Bayesian estimates of model parameters, and deviance information criterion values are obtained for model comparison. The earthquake loss of Yunnan province is analyzed to illustrate the proposed methods. Results show that the estimation of the threshold and the shape and scale of GPD are quite different. Value-at-risk and expected shortfall for the proposed mixture models are calculated under different confidence levels.     

CUSUM,GLR,GEWMA以及RFCuscore控制图监测稳定过程均值漂移的效果比较(英)

本文给出了累积和控制图(CUSUM)监测稳定过程均值漂移的平均运行长度(ARL)的区间估计,并采用数字模拟的方法对CUSUM,GLR,GEWMA以及RFCuscore四种控制图监测稳定过程均值漂移的效果进行比较,结果显示CUSUM效果最好.     

具有可变样本容量的非正态EWMA控制图

为了提高控制图的监控效率,本文研究非正态分布下,EWMA控制图的可变样本容量设计问题。首先利用Burr分布近似各种非正态分布,构造可变样本容量的非正态EWMA控制图;其次运用马尓科夫链法计算可变样本容量非正态EWMA控制图的平均运行长度;然后与传统的非正态EWMA控制图进行比较得出:当过程中出现小波动时,可变样本容量的非正态EWMA控制图能够更快地发现过程中的异常波动,具有较小的平均运行长度,其监控效率明显优于传统的非正态EWMA控制图。     

Comparison of the cuscore, GLRT and cusum control charts for detecting a dynamic mean change

Although statistical process control (SPC) techniques have been focused mostly on detecting constant mean shifts, dynamic and time-varying process changes frequently occur in the monitoring of feedback-controlled and autocorrelated processes. In this research, the performances of cumulative score (Cuscore), generalized likelihood ratio test (GLRT), and cumulative sum (CUSUM) charts in detecting a dynamic mean change that finally approaches a steady-state value are compared. Theoretical results in average run length (ARL) comparison are provided. From the theretical study we find that, when the steady-state value is greater or less than a critical value,Rδ/2+δ/2, the Cuscore and CUSUM charts have a different performance in detecting the mean change. We prove also that the GLRT has the best performance among the three charts in detecting any mean change for which the steady-state value is not equal to δ or δR, when the in-control ARL is large.     

Testing the Gumbel Hypothesis Via the Pot-Method

We consider an i.i.d. sample, generated by some distribution function, which belongs to the domain of attraction of an extreme value distribution with unknown shape and scale parameters. We treat the scale parameter as a nuisance parameter and establish for the hypothesis of Gumbel domain of attraction an asymptotically optimal test based on those observations among the sample, which exceed a given threshold sequence. Asymptotic optimality is achieved along certain contiguous extreme value alternatives within the concept of local asymptotic normality (LAN). Adaptive test procedures exist under restrictive assumptions. The finite sample size behavior of the proposed test is studied by simulations and it is compared to that of a test based on the sample coefficient of variation.     

Extreme value modeling under power normalization

In this paper the mathematical modeling of extremes under power normalization is developed. An estimate of the shape parameter within the generalized extreme value distribution under power normalization is suggested. The statistical inference about the upper tail of a distribution function by using the power normalization is studied. Two models for generalized Pareto distribution under power normalization (GPDP) are given. Estimates for the shape and scale parameters within these GPDP’s are obtained. Finally, a simulation study illustrates and corroborates theoretical results.     

监测均值变动的多重控制图方法

本文证明了当受控平均运行长度充分大时, 多重控制图有两个优点: 一是相比较GLR (广义似然比) 和GEWMA (广义指数权重移动平均)控制图它可以大大降低运算的复杂性; 二是能够较快地监测均值变化的大小. 数值模拟也表明: 多重控制图不仅优于其构成的单个控制图, 而且在监测未知的均值变动方面也优于单个的CUSUM, EWMA, 多重EWMA和GLR控制图.