Some theory and practical uses of trimmed L-moments

Trimmed L-moments, defined by Elamir and Seheult [2003. Trimmed L-moments. Comput. Statist. Data Anal. 43, 299–314], summarize the shape of probability distributions or data samples in a way that remains viable for heavy-tailed distributions, even those for which the mean may not exist. We derive some further theoretical results concerning trimmed L-moments: a relation with the expansion of the quantile function as a weighted sum of Jacobi polynomials; the bounds that must be satisfied by trimmed L-moments; recurrences between trimmed L-moments with different degrees of trimming; and the asymptotic distributions of sample estimators of trimmed L-moments. We also give examples of how trimmed L-moments can be used, analogously to L-moments, in the analysis of heavy-tailed data. Examples include identification of distributions using a trimmed L-moment ratio diagram, shape parameter estimation for the generalized Pareto distribution, and fitting generalized Pareto distributions to a heavy-tailed data sample of computer network traffic.     

Distributions with maximum entropy subject to constraints on their L-moments or expected order statistics

We find the distribution that has maximum entropy conditional on having specified values of its first r  L$L$-moments. This condition is equivalent to specifying the expected values of the order statistics of a sample of size r. The maximum-entropy distribution has a density-quantile function, the reciprocal of the derivative of the quantile function, that is a polynomial of degree r; the quantile function of the distribution can then be found by integration. This class of maximum-entropy distributions includes the uniform, exponential and logistic, and two new generalizations of the logistic distribution. It provides a new method of nonparametric fitting of a distribution to a data sample. We also derive maximum-entropy distributions subject to constraints on expected values of linear combinations of order statistics.     

The minimum L2 distance estimator for Poisson mixture models

A robust estimator is developed for Poisson mixture models with a known number of components. The proposed estimator minimizes the L₂ distance between a sample of data and the model. When the component distributions are completely known, the estimators for the mixing proportions are in closed form. When the parameters for the component Poisson distributions are unknown, numerical methods are needed to calculate the estimators. Compared to the minimum Hellinger distance estimator, the minimum L₂ estimator can be less robust to extreme outliers, and often more robust to moderate outliers.     

Asymptotic normality of two symmetry test statistics based on the L1-error

In this paper, we propose two new tests to test the symmetry of a distribution. These tests are built up on the asymptotic normality of the L₁-distance to the symmetry of the Kernel and histogram density estimates. A simulation study is carried out to evaluate performances of the kernel based test.     

A note on bias reduction of maximum likelihood estimates for the scalar skew t distribution

The skew t distribution is a flexible parametric family to fit data, because it includes parameters that let us regulate skewness and kurtosis. A problem with this distribution is that, for moderate sample sizes, the maximum likelihood estimator of the shape parameter is infinite with positive probability. In order to try to solve this problem, Sartori (2006) has proposed using a modified score function as an estimating equation for the shape parameter. In this note we prove that the resulting modified maximum likelihood estimator is always finite, considering the degrees of freedom as known and greater than or equal to 2.     

A Nonparametric Shewhart-Type Signed-Rank Control Chart Based on Runs

Measures of distributional symmetry based on quantiles, L-moments, and trimmed L-moments are briefly reviewed, and (asymptotic) sampling properties of commonly used estimators considered. Standard errors are estimated using both analytical and computer-intensive methods. Simulation is used to assess results when sampling from some known distributions; bootstrapping is used on sample data to estimate standard errors, construct confidence intervals, and test a hypothesis of distributional symmetry. Symmetry measures based on 2- or 3-trimmed L-moments have some advantages over other measures in terms of their existence. Their estimators are generally well behaved, even in relatively small samples.     

A New Measure of Kurtosis Adjusted for Skewness

Studies of kurtosis often concentrate on only symmetric distributions. This paper identifies a process through which the standardized measure of kurtosis based on the fourth moment about the mean can be written in terms of two parts: (i) an irreducible component, about L₄, which can be seen to occur naturally in the analysis of fourth moments; (ii) terms that depend only on moments of lower order, in particular including the effects of asymmetry attached to the third moment about the mean. This separation of the effect of skewness allows definition of an improved measure of kurtosis. This paper calculates and discusses examples of the new measure of kurtosis for a range of standard distributions.     

Copula representation of bivariate L-moments: a new estimation method for multiparameter two-dimensional copula models

Serfling and Xiao [A contribution to multivariate L-moments, L-comoment matrices. J Multivariate Anal. 2007;98:1765–1781] extended the L-moment theory to the multivariate setting. In the present paper, we focus on the two-dimensional random vectors to establish a link between the bivariate L-moments (BLM) and the underlying bivariate copula functions. This connection provides a new estimate of dependence parameters of bivariate statistical data. Extensive simulation study is carried out to compare estimators based on the BLM, the maximum likelihood, the minimum distance and a rank approximate Z-estimation. The obtained results show that, when the sample size increases, BLM-based estimation performs better as far as the bias and computation time are concerned. Moreover, the root-mean-squared error is quite reasonable and less sensitive in general to outliers than those of the above cited methods. Further, the proposed BLM method is an easy-to-use tool for the estimation of multiparameter copula models. A generalization of the BLM estimation method to the multivariate case is discussed.     

Kurtosis of the logistic-exponential survival distribution

ABSTRACT

In this article, the kurtosis of the logistic-exponential distribution is analyzed. All the moments of this survival distribution are finite, but do not possess closed-form expressions. The standardized fourth central moment, known as Pearson’s coefficient of kurtosis and often used to describe the kurtosis of a distribution, can thus also not be expressed in closed form for the logistic-exponential distribution. Alternative kurtosis measures are therefore considered, specifically quantile-based measures and the L-kurtosis ratio. It is shown that these kurtosis measures of the logistic-exponential distribution are invariant to the values of the distribution’s single shape parameter and hence skewness invariant.     

Tables to facilitate fitting lu distributions

This paper is concerned with estimating the parameters of Tadikamalla-Johnson's L_Udistributions based on the method of moments. Tables of the parameters of the L_U distribution are given for selected values of skewness (0.0(0.05) 1.0(0.1)2.0) and for twenty values of kurtosis at intervals of 0.2. The construction and use of these tables is explained with a numerical example.     

On the efficiency of using the sample kurtosis in selecting optimal lpestimators

This paper examines the efficiency of thesample kurtosisin obtaining L_P estimates as an estimates of central tendency for symmetric distributions. Moreover, guidelines are established for determining an optimal value of P based on the kurtosis of the error distribution.     

A lagrange expansion for the parameters of johnson's Su

A series, based on Lagrange's formula, is given for the kurtosis in terms of the skewness and the ω-parameter of S_u. Also a Maclaunn expansion for a function of ω in terms of the skewness and kurtosis is described.     

上证50ETF隐含高阶矩风险对股票收益的预测研究

高阶矩是刻画资产收益涨跌非对称和“尖峰厚尾”现象中不可忽略的系统性风险。本文基于我国上证50ETF期权数据采用无模型方法估计隐含波动率、隐含偏度和隐含峰度,通过自回归滑动平均模型提取期权隐含高阶矩新息(Innovations),将它们作为高阶矩风险的度量,探讨其对股票收益的预测作用。研究表明:①在控制换手率和股息率等变量后,隐含波动率对于上证50指数和市场未来4周的超额收益有显著负向的预测作用;②隐含偏度新息越低,上证50指数和市场的超额收益越高,这种预测能力在未来1周和未来4周均显著,但随着时间的推移,隐含偏度新息的预测能力逐渐下降;③隐含偏度风险对于我国股市横截面收益也有显著的解释能力,投资组合在隐含偏度风险因子上的风险暴露越大即因子载荷值越大,则未来的收益会越低;④隐含峰度新息总体上与股票收益负相关。     

Maximum penalized likelihood estimation for skew-normal and skew-t distributions

The skew-normal and the skew-t distributions are parametric families which are currently under intense investigation since they provide a more flexible formulation compared to the classical normal and t distributions by introducing a parameter which regulates their skewness. While these families enjoy attractive formal properties from the probability viewpoint, a practical problem with their usage in applications is the possibility that the maximum likelihood estimate of the parameter which regulates skewness diverges. This situation has vanishing probability for increasing sample size, but for finite samples it occurs with non-negligible probability, and its occurrence has unpleasant effects on the inferential process. Methods for overcoming this problem have been put forward both in the classical and in the Bayesian formulation, but their applicability is restricted to simple situations. We formulate a proposal based on the idea of penalized likelihood, which has connections with some of the existing methods, but it applies more generally, including the multivariate case.     

Inference Based on Order Statistics from the Generalized Pareto Distribution and Application

ABSTRACT

In this article, we derive exact explicit expressions for the single, double, triple, and quadruple moments of order statistics from the generalized Pareto distribution (GPD). Also, we obtain the best linear unbiased estimates of the location and scale parameters (BLUE's) of the GPD. We then use these results to determine the mean, variance, and coefficients of skewness and kurtosis of certain linear functions of order statistics. These are then utilized to develop approximate confidence intervals for the generalized Pareto parameters using Edgeworth approximation and compare them with those based on Monte Carlo simulations. To show the usefulness of our results, we also present a numerical example. Finally, we give an application to real data.     

Behaviour of skewness, kurtosis and normality tests in long memory data

We establish the limiting distributions for empirical estimators of the coefficient of skewness, kurtosis, and the Jarque–Bera normality test statistic for long memory linear processes. We show that these estimators, contrary to the case of short memory, are neither ${\sqrt{n}}We establish the limiting distributions for empirical estimators of the coefficient of skewness, kurtosis, and the Jarque–Bera normality test statistic for long memory linear processes. We show that these estimators, contrary to the case of short memory, are neither ?n{\sqrt{n}}-consistent nor asymptotically normal. The normalizations needed to obtain the limiting distributions depend on the long memory parameter d. A direct consequence is that if data are long memory then testing normality with the Jarque–Bera test by using the chi-squared critical values is not valid. Therefore, statistical inference based on skewness, kurtosis, and the Jarque–Bera normality test, needs a rescaling of the corresponding statistics and computing new critical values of their nonstandard limiting distributions.     

Direct density estimation of L-estimates via characteristic functions with applications

In this note we propose a new and novel kernel density estimator for directly estimating the probability and cumulative distribution function of an L-estimate from a single population based on utilizing the theory in Knight (1985) in conjunction with classic inversion theory. This idea is further developed for a kernel density estimator for the difference of L-estimates from two independent populations. The methodology is developed via a “plug-in” approach, but it is distinct from the classic bootstrap methodology in that it is analytically and computationally feasible to provide an exact estimate of the distribution function and thus eliminates the resampling related error. The asymptotic and finite sample properties of our estimators are examined. The procedure is illustrated via generating the kernel density estimate for the Tukey's trimean from a small data set.     

THE RELATIONSHIPS BETWEEN SKEWNESS AND KURTOSIS

Theoretical considerations of kurtosis, whether of partial orderings of distributions with respect to kurtosis or of measures of kurtosis, have tended to focus only on symmetric distributions. With reference to historical points and recent work on skewness and kurtosis, this paper defines anti-skewness and uses it as a tool to discuss the concept of kurtosis in asymmetric univariate distributions. The discussion indicates that while kurtosis is best considered as a property of symmetrised versions of distributions, symmetrisation does not simply remove skewness. Skewness, anti-skewness and kurtosis are all inter-related aspects of shape. The Tukey g and h family and the Johnson Su family are considered as examples.     

Higher order moments of order statistics from the Lindley distribution and associated inference

ABSTRACT

The Lindley distribution is an important distribution for analysing the stress–strength reliability models and lifetime data. In many ways, the Lindley distribution is a better model than that based on the exponential distribution. Order statistics arise naturally in many of such applications. In this paper, we derive the exact explicit expressions for the single, double (product), triple and quadruple moments of order statistics from the Lindley distribution. Then, we use these moments to obtain the best linear unbiased estimates (BLUEs) of the location and scale parameters based on Type-II right-censored samples. Next, we use these results to determine the mean, variance, and coefficients of skewness and kurtosis of some certain linear functions of order statistics to develop Edgeworth approximate confidence intervals of the location and scale Lindley parameters. In addition, we carry out some numerical illustrations through Monte Carlo simulations to show the usefulness of the findings. Finally, we apply the findings of the paper to some real data set.     

Testing for normality in linear regression models

The importance of the normal distribution for fitting continuous data is well known. However, in many practical situations data distribution departs from normality. For example, the sample skewness and the sample kurtosis are far away from 0 and 3, respectively, which are nice properties of normal distributions. So, it is important to have formal tests of normality against any alternative. D'Agostino et al. [A suggestion for using powerful and informative tests of normality, Am. Statist. 44 (1990), pp. 316–321] review four procedures Z ²(g ₁), Z ²(g ₂), D and K ² for testing departure from normality. The first two of these procedures are tests of normality against departure due to skewness and kurtosis, respectively. The other two tests are omnibus tests. An alternative to the normal distribution is a class of skew-normal distributions (see [A. Azzalini, A class of distributions which includes the normal ones, Scand. J. Statist. 12 (1985), pp. 171–178]). In this paper, we obtain a score test (W) and a likelihood ratio test (LR) of goodness of fit of the normal regression model against the skew-normal family of regression models. It turns out that the score test is based on the sample skewness and is of very simple form. The performance of these six procedures, in terms of size and power, are compared using simulations. The level properties of the three statistics LR, W and Z ²(g ₁) are similar and close to the nominal level for moderate to large sample sizes. Also, their power properties are similar for small departure from normality due to skewness (γ₁≤0.4). Of these, the score test statistic has a very simple form and computationally much simpler than the other two statistics. The LR statistic, in general, has highest power, although it is computationally much complex as it requires estimates of the parameters under the normal model as well as those under the skew-normal model. So, the score test may be used to test for normality against small departure from normality due to skewness. Otherwise, the likelihood ratio statistic LR should be used as it detects general departure from normality (due to both skewness and kurtosis) with, in general, largest power.