A NEW APPROXIMATION FOR FISHER'S Z

As the sample size increases, the coefficient of skewness of the Fisher's transformation z= tanh^-1r, of the correlation coefficient decreases much more rapidly than the excess of its kurtosis. Hence, the distribution of standardized z can be approximated more accurately in terms of the t distribution with matching kurtosis than by the unit normal distribution. This t distribution can, in turn be subjected to Wallace's approximation resulting in a new normal approximation for the Fisher's z transform. This approximation, which can be used to estimate the probabilities, as well as the percentiles, compares favorably in both accuracy and simplicity, with the two best earlier approximations, namely, those due to Ruben (1966) and Kraemer (1974). Fisher (1921) suggested approximating distribution of the variance stabilizing transform z=(1/2) log ((1 +r)/(1r)) of the correlation coefficient r by the normal distribution with mean = (1/2) log ((1 + p)/(lp)) and variance =l/(n3). This approximation is generally recognized as being remarkably accurate when ||Gr| is moderate but not so accurate when ||Gr| is large, even when n is not small (David (1938)). Among various alternatives to Fisher's approximation, the normalizing transformation due to Ruben (1966) and a t approximation due to Kraemer (1973), are interesting on the grounds of novelty, accuracy and/or aesthetics. If r?= r/√ (1r²) and r?|Gr = |Gr/√(1|Gr²), then Ruben (1966) showed that (1) g_n (r,|Gr) ={(2n5)/2}^1/2r?r{(2n3)/2}^1/2r?|GR, {1 + (1/2)(r?r²+r?|Gr²)}^1/2 is approximately unit normal. Kraemer (1973) suggests approximating (2) t_n (r, |Gr) = (r|GR¹) √ (n2), √(11r²) √(1|Gr²) by a Student's t variable with (n2) degrees of freedom, where after considering various valid choices for |Gr¹ she recommends taking |Gr¹= |Gr*, the median of r given n and |Gr.     

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.     

Inference for the Generalized Normal Laplace Distribution

The generalized normal Laplace distribution has been used in financial modeling because of its skewness and excess kurtosis. To estimate its parameters, we use a method based on minimizing the quadratic distance between the real and imaginary parts of the empirical and theoretical characteristic functions. The quadratic distance estimator (QDE) obtained is shown to be robust, consistent, and with an asymptotically normal distribution. The goodness-of-fit test statistics presented follow an asymptotic chi-square distribution. The performance of the QDE is illustrated by simulation results and an application to financial data.     

Bimodal skew-symmetric normal distribution

ABSTRACT

We introduce a new parsimonious bimodal distribution, referred to as the bimodal skew-symmetric Normal (BSSN) distribution, which is potentially effective in capturing bimodality, excess kurtosis, and skewness. Explicit expressions for the moment-generating function, mean, variance, skewness, and excess kurtosis were derived. The shape properties of the proposed distribution were investigated in regard to skewness, kurtosis, and bimodality. Maximum likelihood estimation was considered and an expression for the observed information matrix was provided. Illustrative examples using medical and financial data as well as simulated data from a mixture of normal distributions were worked.     

On distribution of AIC in linear regression models

This paper investigates an asymptotic distribution of the Akaike information criterion (AIC) and presents its characteristics in normal linear regression models. The bias correction of the AIC has been studied. It may be noted that the bias is only the mean, i.e., the first moment. Higher moments are important for investigating the behavior of the AIC. The variance increases as the number of explanatory variables increases. The skewness and kurtosis imply a favorable accuracy of the normal approximation. An asymptotic expansion of the distribution function of a standardized AIC is also derived.     

The Mean Squared Error of the Instrumental Variables Estimator When the Disturbance Has an Elliptical Distribution

This paper generalizes Nagar's (1959) approximation to the finite sample mean squared error (MSE) of the instrumental variables (IV) estimator to the case in which the errors possess an elliptical distribution whose moments exist up to infinite order. This allows for types of excess kurtosis exhibited by some financial data series. This approximation is compared numerically to Knight's (1985) formulae for the exact moments of the IV estimator under nonnormality. We use the results to explore two questions on instrument selection. First, we complement Buse's (1992) analysis by considering the impact of additional instruments on both bias and MSE. Second, we evaluate the properties of Andrews's (1999) selection method in terms of the bias and MSE of the resulting IV estimator.     

Characterizing parameters of multivariate elliptical distributions*

This paper defines new parameters characterizing multivariate elliptical distributions. Mardia's coefficient of multivariate kurtosis is shown to be essentially one of these parameters. A simple relation is established between centered multivariate product moments and second moments of the variables. The general results are verified on the contaminated normal distribution as an example.     

The null distribution of multivariate kurtosis

The paper introduces a x²-approximation to multivariate kurtosis b_2,punder normality. It requires calculating the third moment of b_2,pwhich is obtained. We compare the approximation with simulated percentage points and the normal approximation, and find it to be adequate for p=l and 2. For p=3, the simple average of this estimate and the normal approximation is found to be generally superior to either approximation on its own. For p=4, the normal approximation is best for non-extreme values of ∝     

Estimation of variance of mean using known coefficient of variation

An optimum unbiased estimator of the variance of mean is given It is defined as a function of the mean and itscustomary unbiased variance estimator, utilizing known coefficient of variation, skewness and kurtosis of the underlying distributions. Exact results are obtained. Normal and large sample cases receive particular treatment. The proposed variance estimator is generally more efficient than the customary variance estimator; its relative efficiency becomes appreciably higher for smaller coefficient of variation, smaller sample (in the normal case at least), higher negative skewness, or higher positive skewness with sufficiently large kurtosis. The empirical findings are reassuring and supportive.     

The two-sided power distribution for the treatment of the uncertainty in PERT

The aim of this paper is to include the Two-Sided Power (TSP) distribution in the PERT methodology making use of the advantages that this four-parameter distribution offers. In order to be completely determined, a distribution of this type needs, the same as the beta distribution, a new datum apart from the three usual values a (pessimistic), m (most likely) and b (optimistic). To solve this question, when using the beta distribution in the PERT context, we are looking for the maximum similarity with the normal and so it is required that the distribution has the same variance as the normal or its same kurtosis, giving rise to the constant variance and mesokurtic families, respectively. Nevertheless, while this approach can be only applied to the beta distribution for some values in the range of the standardized mode, in the case of the TSP distribution this methodology leads always to a solution. A detailed analysis comparing the beta and TSP distribution based on their PERT means and variances is presented indicating better results for the second. We are very grateful for the comments and suggestions of two anonymous referees.     

A powerful and interpretable alternative to the Jarque–Bera test of normality based on 2nd-power skewness and kurtosis,using the Rao's score test on the APD family

We introduce the 2nd-power skewness and kurtosis, which are interesting alternatives to the classical Pearson's skewness and kurtosis, called 3rd-power skewness and 4th-power kurtosis in our terminology. We use the sample 2nd-power skewness and kurtosis to build a powerful test of normality. This test can also be derived as Rao's score test on the asymmetric power distribution, which combines the large range of exponential tail behavior provided by the exponential power distribution family with various levels of asymmetry. We find that our test statistic is asymptotically chi-squared distributed. We also propose a modified test statistic, for which we show numerically that the distribution can be approximated for finite sample sizes with very high precision by a chi-square. Similarly, we propose a directional test based on sample 2nd-power kurtosis only, for the situations where the true distribution is known to be symmetric. Our tests are very similar in spirit to the famous Jarque–Bera test, and as such are also locally optimal. They offer the same nice interpretation, with in addition the gold standard power of the regression and correlation tests. An extensive empirical power analysis is performed, which shows that our tests are among the most powerful normality tests. Our test is implemented in an R package called PoweR.     

Brownian–Laplace Motion and Its Use in Financial Modelling

Brownian-Laplace motion is a Lévy process which has both continuous (Brownian) and discontinuous (Laplace motion) components. The increments of the process follow a generalized normal Laplace (GNL) distribution which exhibits positive kurtosis and can be either symmetrical or exhibit skewness. The degree of kurtosis in the increments increases as the time between observations decreases. This and other properties render Brownian-Laplace motion a good candidate model for the motion of logarithmic stock prices. An option pricing formula for European call options is derived and it is used to calculate numerically the value of such an option both using nominal parameter values (to explore its dependence upon them) and those obtained as estimates from real stock price data.     

Estimating kurtosis and confidence intervals for the variance under nonnormality

Exact confidence intervals for variances rely on normal distribution assumptions. Alternatively, large-sample confidence intervals for the variance can be attained if one estimates the kurtosis of the underlying distribution. The method used to estimate the kurtosis has a direct impact on the performance of the interval and thus the quality of statistical inferences. In this paper the author considers a number of kurtosis estimators combined with large-sample theory to construct approximate confidence intervals for the variance. In addition, a nonparametric bootstrap resampling procedure is used to build bootstrap confidence intervals for the variance. Simulated coverage probabilities using different confidence interval methods are computed for a variety of sample sizes and distributions. A modification to a conventional estimator of the kurtosis, in conjunction with adjustments to the mean and variance of the asymptotic distribution of a function of the sample variance, improves the resulting coverage values for leptokurtically distributed populations.     

Effects of some design factors on the distribution of similarity indices in cluster analysis

This article investigates the effects of number of clusters, cluster size, and correction for chance agreement on the distribution of two similarity indices, namely, Jaccard and Rand indices. Skewness and kurtosis are calculated for the two indices and their corrected forms then compared with those of the normal distribution. Three clustering algorithms are implemented: complete linkage, Ward, and K-means. Data were randomly generated from bivariate normal distributions with specified means and variance covariance matrices. Three-way ANOVA is performed to assess the significance of the design factors using skewness and kurtosis of the indices as responses. Test statistics for testing skewness and kurtosis and observed power are calculated. Simulation results showed that independent of the clustering algorithms or the similarity indices used, the interaction effect cluster size x number of clusters and the main effects of cluster size and number of clusters were found always significant for skewness and kurtosis. The three way interaction of cluster size x correction x number of clusters was significant for skewness of Rand and Jaccard indices using all clustering algorithms, but was not significant using Ward's method for both Rand and Jaccard indices, while significant for Jaccard only using complete linkage and K-means algorithms. The correction for chance agreement was significant for skewness and kurtosis using Rand and Jaccard indices when complete linkage method is used. Hence, such design factors must be taken into consideration when studying distribution of such indices.     

Generalized Tukey-type distributions with application to financial and teletraffic data

Constructing skew and heavy-tailed distributions by transforming a standard normal variable goes back to Tukey (Exploratory data analysis. Addison-Wesley, Reading, 1977) and was extended and formalized by Hoaglin (In: Data analysis for tables, trends, and shapes. Wiley, New York, 1983) and Martinez and Iglewicz (Commun Statist Theory Methods 13(3):353–369, 1984). Applications of Tukey’s GH distribution family—which are composed by a skewness transformation G and a kurtosis transformation H—can be found, for instance, in financial, environmental or medical statistics. Recently, alternative transformations emerged in the literature. Rayner and MacGillivray (Statist Comput 12:57–75, 2002b) discuss the GK distributions, where Tukey’s H-transformation is replaced by another kurtosis transformation K. Similarly, Fischer and Klein (All Stat Arch, 88(1):35–50, 2004) advocate the J-transformation which also produces heavy tails but—in contrast to Tukey’s H-transformation—still guarantees the existence of all moments. Within this work we present a very general kurtosis transformation which nests H-, K-and an approximation to the J-transformation and, hence, permits to discriminate between them. Applications to financial and teletraffic data are given.     

Correlated Errors in the Parameters Estimation of the ARFIMA Model: A Simulated Study

Kurtosis, usually as measured by the standardised fourth central moment, has been examined on a number of occasions by observing the effect of contaminating the distribution, that is, mixing in another distribution. However, superficial treatment can lead, and indeed has led, to misunderstandings. This paper considers, firstly for a symmetric distribution contaminated at two points symmetrically placed around its centre and then for a mixture of two continuous symmetric distributions, the behaviour of three measures of kurtosis. This is done in general and not just as the mixing proportion tends to zero as in the influence function approach. It is seen that when both scale and kurtosis change, the latter is not necessarily intuitive. It is also illustrated that parameter interpretation in terms of distributional properties such as shape can be misleading without the use of the appropriate distributional partial ordering     

Generalized modified slash distribution with applications

Abstract

In this article a generalization of the modified slash distribution is introduced. This model is based on the quotient of two independent random variables, whose distributions are a normal and a one-parameter gamma, respectively. The resulting distribution is a new model whose kurtosis is greater than other slash distributions. The probability density function, its properties, moments, and kurtosis coefficient are obtained. Inference based on moment and maximum likelihood methods is carried out. The multivariate version is also introduced. Two real data sets are considered in which it is shown that the new model fits better to symmetric data with heavy tails than other slash extensions previously introduced in literature.     

A comparison of some approximations to the Wilcoxon-Mann-Whitney distribution

The following approximations to the exact distribution of the Wilcoxon rank sum test (Mann-Whitney U-test) are compared at or near significance levels .05, .025, .01 and .005: the normal approximation with continuity correction, the simple and complete version of the Edgeworth series approximation proposed by Fix and Hodges (1955), Buckle, Kraft and van Eeden’s (1969a) uniform approximation and Iman’s (1976) recently proposed approximation. The comparison takes into account simplicity of application as well as closeness, and some preference rules are suggested for different user objectives.     

A note on the accuracy of fisher's approximation to the large sample variance of an intraclass correlationa

The adequacy of Fisher's approximation to the large sample variance of an intraclass correlation is investigated in the context of family studies. It is found that the approximation is highly accurate in samples of moderately large size (≧ 30 families), and can also be used for significance-testing under a broad range of circumstances. The exact sampling of distribution of the intraclass correlation coefficient is also derived.