Extended generalised Pareto models for tail estimation

The most popular approach in extreme value statistics is the modelling of threshold exceedances using the asymptotically motivated generalised Pareto distribution. This approach involves the selection of a high threshold above which the model fits the data well. Sometimes, few observations of a measurement process might be recorded in applications and so selecting a high quantile of the sample as the threshold leads to almost no exceedances. In this paper we propose extensions of the generalised Pareto distribution that incorporate an additional shape parameter while keeping the tail behaviour unaffected. The inclusion of this parameter offers additional structure for the main body of the distribution, improves the stability of the modified scale, tail index and return level estimates to threshold choice and allows a lower threshold to be selected. We illustrate the benefits of the proposed models with a simulation study and two case studies.     

中国居民收入分布专题实证研究——居民收入分布变迁测度及其影响因素分解

基于核密度估计方法,给出居民收入分布变迁测度及其影响因素分解的相对收入分布方法。使用中国健康和营养调查的微观居民收入数据,实证分析中国居民收入分布及其变迁的特征、过程及其影响因素。实证结果显示:中国居民收入分布变迁的速度与经济体制改革进程相对应;经济增长因素是引起居民收入分布变迁的主要因素,在提高总体居民收入水平的同时加剧了居民收入分布的极化程度;收入分配因素对居民收入分布变迁的作用相对较弱,但有一定的减贫或缩小居民收入差距作用。     

公平分配与居民福利——收入差距对中国居民恩格尔系数的影响分析

居民收入水平不同,恩格尔系数也不同,因此收入分配对社会总体恩格尔系数的影响是不确定的。基于马斯洛层次需求理论证明边际食品消费倾向递减,并且当消费函数为线性或近似线性时,收入差距扩大会降低恩格尔系数。基于省际面板数据的实证分析表明,中国居民边际食品消费倾向递减且消费函数为(近似)线性,居民收入差距扩大降低了居民恩格尔系数(基尼系数升高对恩格尔系数下降的贡献率为22%),从而出现"分配越不平等,居民总体生活水平越高"的矛盾现象。因此在使用总体指标评价居民生活水平时,收入分配指标必须得到同等的重视;在经济发展中,公平和效率必须兼顾才能实现居民福利最大化。     

指数族分布中参数的经验贝叶斯估计

指数族分布是一类应用广泛的分布类,包括了泊松分布、Gamma分布、Beta分布、二项分布等常见分布.在非寿险中,索赔额或索赔次数过程常常被假定服从指数族分布,由于风险的非齐次性,指数族分布中的参数θ也为随机变量,假定服从指数族共轭先验分布.此时风险参数的估计落入了Bayes框架,风险参数θ的Bayes估计被表达“信度”形式.然而,在实际运用中,由于先验分布与样本分布中仍然含有结构参数,根据样本的边际分布的似然函数估计结构参数,从而获得风险参数的经验Bayes估计,最后证明了该经验Bayes估计是渐近最优的.     

Generation of multivariate distributions by vertical density representation

Troutt (1991,1993) proposed the idea of the vertical density representation (VDR) based on Box-Millar method. Kotz, Fang and Liang (1997) provided a systematic study on the multivariate vertical density representation (MVDR). Suppose that we want to generate a random vector X[d]Rⁿthat has a density function ?(x). The key point of using the MVDR is to generate the uniform distribution on [D]_?(v) = {x :?(x) = v} for any v > 0 which is the surface in RⁿIn this paper we use the conditional distribution method to generate the uniform distribution on a domain or on some surface and based on it we proposed an alternative version of the MVDR(type 2 MVDR), by which one can transfer the problem of generating a random vector X with given density f to one of generating (X, X_n+i) that follows the uniform distribution on a region in Rⁿ⁺¹defined by ?. Several examples indicate that the proposed method is quite practical.     

Progressive type II censored order statistics from exponential distributions

In the model of progressive type II censoring, point and interval estimation as well as relations for single and product moments are considered. Based on two-parameter exponential distributions, maximum likelihood estimators (MLEs), uniformly minimum variance unbiased estimators (UMVUEs) and best linear unbiased estimators (BLUEs) are derived for both location and scale parameters. Some properties of these estimators are shown. Moreover, results for single and product moments of progressive type II censored order statistics are presented to obtain recurrence relations from exponential and truncated exponential distributions. These relations may then be used to compute all the means, variances and covariances of progressive type II censored order statistics based on exponential distributions for arbitrary censoring schemes. The presented recurrence relations simplify those given by Aggarwala and Balakrishnan (1996)     

On optimising the estimation of high quantiles of a probability distribution

One of the major aims of one-dimensional extreme-value theory is to estimate quantiles outside the sample or at the boundary of the sample. The underlying idea of any method to do this is to estimate a quantile well inside the sample but near the boundary and then to shift it somehow to the right place. The choice of this “anchor quantile” plays a major role in the accuracy of the method. We present a bootstrap method to achieve the optimal choice of sample fraction in the estimation of either high quantile or endpoint estimation which extends earlier results by Hall and Weissman (1997) in the case of high quantile estimation. We give detailed results for the estimators used by Dekkers et al. (1989). An alternative way of attacking problems like this one is given in a paper by Drees and Kaufmann (1998).     

Robust estimation in the generalized poisson model

A computationally simple method of robust estimation in the generalized Poisson model is presented. Estimators are proved to be optimal in the sense of local minimax testing, conditionally on the explanatory variable. Results of a Monte Carlo experiment are supplemented where robust and efficient estimators are compared.     

The asymptotic distribution of the partial attributable risk in cross-sectional studies

The partial attributable risk (PAR) has been introduced as a tool for partitioning the responsibility for causing an adverse event between various risk factors. It has arisen from epidemiology, but it is also a valid general risk allocation concept, which can, for example, be applied to data from customer satisfaction surveys. So far, a variance formula for the PAR has been missing so that the confidence intervals were not directly available. This paper provides the asymptotic normal distribution for the PAR determined from a cross-sectional study.     

Kernel density estimation for heavy-tailed distributions using the champernowne transformation

When estimating loss distributions in insurance, large and small losses are usually split because it is difficult to find a simple parametric model that fits all claim sizes. This approach involves determining the threshold level between large and small losses. In this article, a unified approach to the estimation of loss distributions is presented. We propose an estimator obtained by transforming the data set with a modification of the Champernowne cdf and then estimating the density of the transformed data by use of the classical kernel density estimator. We investigate the asymptotic bias and variance of the proposed estimator. In a simulation study, the proposed method shows a good performance. We also present two applications dealing with claims costs in insurance.