Further approximate pearson percentage points and Cornish-Fisher

A set of 19-point rational fraction approximants for Pearson distributions at levels a = 0.025, 0.25, 0.50, 0.75, and 0.975 is added to a previous study (Biometrika, 1979). The relation to a Cornish-Fisher expansion is discussed and advice proffered on interpolation for a new deviate or new level.     

A higher order approximation to a percentage point of the non–central t-distribution

A higher order approximation formula for a percentage point of the noncentral t–distribution with v degrees of freedom is given up to the order o(v^-3), using the Cornish-Fisher expansion for the statistic based on a lin-ear combination of a normal random variable and a chi-random variable. The upper confidence limit and the confidence interval for the non–centrality parameter are given. Numerical results are also obtained.     

Approximate percentiles for ratios of f-variates

The need for percentiles of ratios of independent F-variates arises in many types of applications. Some tabulated values are available mostly for equal or nearly equal degrees-of-freedom cases, but it is not possible to provide comprehensive tables for the numerous combinations of possible degrees-of-freedom and probability levels. A simple approximation for the percentiles based on the Cornish-Fisher expansion technique is provided which is adequate for larger degrees-of-freedom and a more complicated approximation is provided which is suitable for all cases. The accuracy of the approximation is shown over a wide range of typical cases     

The cornish-fisher expansion for estimating percentage points: a numerical perspective

The expansion, in standard form, consists of some 66 terms involving polyno - mials in a normal deviate, and cumulants and cumulant products to order ten. An assumed order of magnitude reduces these terms to eight groups. Sign patterns in the terms are not obvious. We take a number of Pearson densities and assess from the expansions a set of standard percentiles (1%, 5%, 95%, 99%). Validity of the as-sessments is pivoted on two alternative models:(i) the Bowman-Shenton algorithm for percentage points of Pearson densities, (ii) the 4-moment Johnson translation model. This approach has wide application since the models have proved to be remarkably reliable when compared, and also when compared with simulation as-sessments. A brief account is given of acceleration of convergence for the series, but there seems to be no analogue of the Padè or Levin algorithms.

The Cornish-Fisher application to the Fisher z-statistic is studied and the cumu- lants defined in general. Irwin's expression for the density of means from Pearson Type II is recalled. There is an error in the Cornish-Fisher treatment of the z-statistic but this is one which has its source in the write-up. Again the Irwin density in the general case has a factor missing.     

Approximations to the chi-square distribution

The chi-square distribution arises frequently in applied statistics.Associated with the chi-square random variable with v degrees of freedom are two interdependent variables: the probability integral and the percentage point.Given one of these variables,the other can be obtained from chi-square tables for selected values.In order to overcome the inconvenience of statistical tables and interpolation,many approximations have been suggested.The computational difficulty and accuracy of various approximations is compared.     

Approximating the difference of two t-variables for all degrees of freedom using truncated variables

A scaled t‐distribution is used to approximate the distribution of a linear combination of two independent t‐variables for any number of degrees of freedom, and in particular for low degrees of freedom where moments do not exist. The approximation is the method‐of‐moments solution to the analogous problem with truncated t‐variables. The approximation exists for all degrees of freedom, is very accurate for more than two degrees of freedom, and performs as well as other approximations of this form when they exist.     

A modification of the Fisher-Cornish approximation for the student t percentiles

An asymptotic expansion of the Student t distribution is derived by expanding the standardized Student t distribution in terms of the normal distribution. This expansion is inverted to obtain corresponding asymptotic expansions for the Student t percentiles as functions of the standard normal percentiles0 Using the first two, three or four terms of these expansions, we get approximations of the Student t percentiles which are generally more accurate than the approximations given by Fisher and Cornish(1960) and Koehler (1983).An approximation of the distribution function obtained from this expansion is compared with the approximations discussed by Ling (1978) andfound to be more accurate for moderate degrees of freedom.     

Exact and resampling probability values for the Piccarreta nominal–ordinal index of association

Exact, resampling, and Pearson type III permutation methods are provided to compute probability values for Piccarreta's nominal–ordinal index of association. The resampling permutation method provides good approximate probability values based on the proportion of resampled test statistic values equal to or greater than the observed test statistic value.     

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.     

从统计学史解析“自由度”概念

自由度是统计学中一个十分重要而又长期没有被圆满解释的概念。对此,从统计学史角度,对皮尔逊、费歇尔有关自由度问题争论原始文献细致考察,彻底澄清了自由度概念的内涵及与其相关的统计思想,弥补了Fienberg、Stigler与陈希孺已有解释的缺陷。研究表明：皮尔逊关于卡方检验中无论总体分布已知还是其来自于样本推断统计量都具有同一分布的错误判断,导致卡方检验的准确性出现偏差,这种偏差虽被同时代少数几个统计学家察觉但他们却无法解释其根源。费歇尔提出自由度概念并结合n维几何、假设检验与最大似然方法的论证不仅修正了皮尔逊的错误,也完善了从样本统计量估计总体参数的数理逻辑。     

Higher-order Bartlett-type adjustment

This paper deals with Bartlett-type adjustment which makes all the terms up to order n^−k in the asymptotic expansion vanish, where k is an integer k ⩾ 1 and n depends on the sample size. Extending Cordeiro and Ferrari (1991, Biometrika, 78, 573–582) for the case of k = 1, we derive a general formula of the kth-order Bartlett-type adjustment for the test statistic whose kth-order asymptotic expansion of the distribution is given by a finite linear combination of chi-squared distribution with suitable degrees of freedom. Two examples of the second-order Bartlett-type adjustment are given. We also elucidate the connection between Bartlett-type adjustment and Cornish-Fisher expansion.     

A note on the power of the f-test for testing linear hypothesis with elliptic t errors

The power of the classical .F-test for testing the regression coefficient of a general linear model with elliptic t error variable depends on the degrees of freedom of the t- distribution. In this note it is shown that the power of the F-test based on t-distribution is greater than the normal based test at smaller level of significance.     

Higher Order Asymptotic Cumulants of Studentized Estimators in Covariance Structures

In covariance structure analysis, the Studentized pivotal statistic of a parameter estimator is often used since the statistic is asymptotically normally distributed with mean zero and unit variance. For more accurate asymptotic distribution, the first and third asymptotic cumulants can be used to have the single-term Edgeworth, Cornish-Fisher, and Hall type asymptotic expansions. In this paper, the higher order asymptotic variance and the fourth asymptotic cumulant of the statistic are obtained under nonnormality when the partial derivatives of a parameter estimator with respect to sample variances and covariances up to the third order and the moments of the associated observed variables up to the eighth order are available. The result can be used to have the two-term Edgeworth expansion. Simulations are performed to see the accuracy of the asymptotic results in finite samples.     

A system for computing joint probabilities from jensen's bivariate f distribution

A system of subroutines is presented for efficient computation of joint probabilities from Jensen's bivariate F distribution. Any valid set of parameters is permitted, whereas previous work was limited to the special case of equal numerator degrees of freedom and equal canonical correlations in the underlying multinormal distribution. The use of joint Probabilities from Jensen's bivariate F distribution is demonstrated via an application to two-way ANOVA without interaction.     

The Cornish-Fisher expansion for a class of statistics in first order autoregression

ABSTRACT

This article is concerned with the derivation and study of the Cornish-Fisher expansion for a wide class of estimators of the parameter in the first order autoregressive process. Second and third order Cornish-Fisher approximations to the quantile of the distribution of the corresponding asymptotically normal standardized statistic are stated explicitly and their accuracy is examined, both theoretically and numerically, by comparing them with the exact value of the quantile obtained by Monte Carlo simulation.     

Simple approximations to the behrens—fisher distribution

Dayal and Dickey (1977) have published in this journal a rather efficient numerical integration procedure for the product of k Student t-densities, and point out the evaluation of Behrens-Fisher (BF) densities as an important special case. The present note adds to this three simple normal approximations to Behrens-Fisher tail probabilities, that will save computer time for someone using the Dayal-Dickey results, and even allow evaluation on a desk calculator for moderately large degrees of freedom.

A direct normal approximation (method U) will be too coarse unless both degrees of freedom are large. A combination of the Peizer-Pratt (1968) approximation to the t-distribution and the Patil (1965) t-approximation to the BF distribution turns out to be very accurate. For very small degrees of freedom it may still be refined by an adhoc correction presented below. Other approximations and expansions turn out to be less satisfactory than the present trio. It facilitates a quick evaluation of BF probabilities and quantiles on a small computer or even a pocket calculator.     

New monte carlo results on the robustness of the anova f,w and f statistics

Because the usual F test for equal means is not robust to unequal variances, Brown and Forsythe (1974a) suggest replacing F with the statistics F or W which are based on the Satterthwaite and Welch adjusted degrees of freedom procedures. This paper reports practical situations where both F and W give * unsatisfactory results. In particular, both F and W may not provide adequate control over Type I errors. Moreover, for equal variances, but unequal sample sizes, W should be avoided in favor of F (or F ), but for equal sample sizes, and possibly unequal variances, W was the only satisfactory statistic. New results on power are included as well. The paper also considers the effect of using F or W only after a significant test for equal variances has been obtained, and new results on the robustness of the F test are described. It is found that even for equal sample sizes as large as 50 per treatment group, there are practical situations where the F test does not provide adequately control over the probability of a Type I error.     

Estimation of Weibull parameters from common percentiles

Estimation of Weibull distribution shape and scale parameters is accomplished through use of symmetrically located percentiles from a sample. The process requires algebraic solution of two equations derived from the cumulative distribution function. Three alternatives examined are compared for precision and variability with maximum likelihood (MLE) and least squares (LS) estimators. The best percentile estimator (using the 10th and 90th) is inferior to MLE in variability and to one least squares estimator in accuracy and variability to a small degree. However, application of a correction factor related to sample size improves the percentile estimator substantially, making it more accurate than LS.     

Additional bowman-shenton approximate percentage points for pearson distributions based on pearson type vi

The Bowman and Shenton approximate percentage points for Pearson distributions are extended to include some cases for which the skewness β¹ exceeds four.Use is made of the reciprocity property existing between the Pearson Type VI distribution and the Type I distribution     

speed of covergence to the extreme value distributions on their probability ploting parers

The speed of convergence of the distribution of the normalized maximum, of a sample of independent and identically distributed random variables, to its asymptotic distribution is considered in this article. Assuming that the cumulative distribution function of the random variables is known, the error committed replacing the actual distribution of the normalized maximum by its asympotic distribution is studied. Instead of using the arithmetical scale of probabilities, we measure the difference between the actual and asympotic distribution in terms of the double-log scale used for building the probability plotting paper for the the latter. We demonstrate that the difference between the double-log values corresponding to two probabilities in the upper tail is almost exactly equal to the logarithm of the distribution may not be uniform in this double-log scale and that the difference between the actual and asymptotic distributions, on the probebility plotting paper, may be a logarithmic, power, or even exponential function in the upper tail when the latter distribution is of Fisher-Tippett type I, but that difference is at most logarithmic in the upper tail for type II and III distributions. This fact is exploited to obtain transformed variables that converge tothe asymptotic distribution faster than the original variable on the probabilites plotting paper