A new family of life distributions based on the elliptically contoured distributions

In this article we propose a new family of life distributions, generated from an elliptically contoured distribution, in which the density and some of its properties are obtained. Explicit expressions for the density are found for a large number of specific elliptical distributions, such as Pearson type VII, t, Cauchy, Kotz type, normal, Bessel, Laplace and logistic.     

A test for bivariate two sample location problem

A consistent test for difference in locations between two bivariate populations is proposed, The test is similar as the Mann-Whitney test and depends on the exceedances of slopes of the two samples where slope for each sample observation is computed by taking the ratios of the observed values. In terms of the slopes, it reduces to a univariate problem, The power of the test has been compared with those of various existing tests by simulation. The proposed test statistic is compared with Mardia's(1967) test statistics, Peters-Randies(1991) test statistic, Wilcoxon's rank sum test. statistic and Hotelling' T² test statistic using Monte Carlo technique. It performs better than other statistics compared for small differences in locations between two populations when underlying population is population 7(light tailed population) and sample size 15 and 18 respectively. When underlying population is population 6(heavy tailed population) and sample sizes are 15 and 18 it performas better than other statistic compared except Wilcoxon's rank sum test statistics for small differences in location between two populations. It performs better than Mardia's(1967) test statistic for large differences in location between two population when underlying population is bivariate normal mixture with probability p=0.5, population 6, Pearson type II population and Pearson type VII population for sample size 15 and 18 .Under bivariate normal population it performs as good as Mardia' (1967) test statistic for small differences in locations between two populations and sample sizes 15 and 18. For sample sizes 25 and 28 respectively it performs better than Mardia's (1967) test statistic when underlying population is population 6, Pearson type II population and Pearson type VII population     

On the efficiency property of the linear estimators in the bivariate pearson type vii distribution

In this paper the efficiency property of the estimators of the parameters of the bivariate Pearson type VII distribution is studied inside the family of linear estimators, assuming that the sample is constituted by dependent random vectors. It is proven that, although there are not efficient linear estimators, the sample mean and the sample covariance matrix (affected by an unbiasedness weighting) are unbiased linear estimators of minimum distance to the Cramér-Rao lower bound. Finally, a numerical simulation example shows that the proposed estimators are computationally feasible.     

Moments of truncated normal/independent distributions

In this work we have considered the problem of finding the moments of a doubly truncated member of the class of normal/independent distributions. We obtained a general result and then use it to derive the moments in the case of doubly truncated versions of Pearson type VII distribution, slash distribution, contaminated normal distribution, double exponential distribution and variance gamma distribution. We also give an application of some actuarial data.     

A new family of multivariate heavy-tailed distributions with variable marginal amounts of tailweight: application to robust clustering

We propose a family of multivariate heavy-tailed distributions that allow variable marginal amounts of tailweight. The originality comes from introducing multidimensional instead of univariate scale variables for the mixture of scaled Gaussian family of distributions. In contrast to most existing approaches, the derived distributions can account for a variety of shapes and have a simple tractable form with a closed-form probability density function whatever the dimension. We examine a number of properties of these distributions and illustrate them in the particular case of Pearson type VII and t tails. For these latter cases, we provide maximum likelihood estimation of the parameters and illustrate their modelling flexibility on simulated and real data clustering examples.     

On influence diagnostic in univariate elliptical linear regression models

We discuss in this paper the assessment of local influence in univariate elliptical linear regression models. This class includes all symmetric continuous distributions, such as normal, Student-t, Pearson VII, exponential power and logistic, among others. We derive the appropriate matrices for assessing the local influence on the parameter estimates and on predictions by considering as influence measures the likelihood displacement and a distance based on the Pearson residual. Two examples with real data are given for illustration.     

SINGULAR ELLIPTICAL DISTRIBUTION: DENSITY AND APPLICATIONS

ABSTRACT

The paper present an explicit expression for the density of a n-dimensional random vector with a singular Elliptical distribution. Based on this, the densities of the generalized Chi-squared and generalized t distributions are derived, examining the Pearson Type VII distribution and Kotz Type distribution (as specific Elliptical distributions). Finally, the results are applied to the study of the distribution of the residuals of an Elliptical linear model and the distribution of the t-statistic, based on a sample from an Elliptical population.     

Variable selection in elliptical linear mixed model

Variable selection in elliptical Linear Mixed Models (LMMs) with a shrinkage penalty function (SPF) is the main scope of this study. SPFs are applied for parameter estimation and variable selection simultaneously. The smoothly clipped absolute deviation penalty (SCAD) is one of the SPFs and it is adapted into the elliptical LMM in this study. The proposed idea is highly applicable to a variety of models which are set up with different distributions such as normal, student-t, Pearson VII, power exponential and so on. Simulation studies and real data example with one of the elliptical distributions show that if the variable selection is also a concern, it is worthwhile to carry on the variable selection and the parameter estimation simultaneously in the elliptical LMM.     

SCALE MIXTURES DISTRIBUTIONS IN STATISTICAL MODELLING

This paper presents two types of symmetric scale mixture probability distributions which include the normal, Student t, Pearson Type VII, variance gamma, exponential power, uniform power and generalized t (GT) distributions. Expressing a symmetric distribution into a scale mixture form enables efficient Bayesian Markov chain Monte Carlo (MCMC) algorithms in the implementation of complicated statistical models. Moreover, the mixing parameters, a by-product of the scale mixture representation, can be used to identify possible outliers. This paper also proposes a uniform scale mixture representation for the GT density, and demonstrates how this density representation alleviates the computational burden of the Gibbs sampler.     

Bayesian analysis of censored linear regression models with scale mixtures of normal distributions

As is the case of many studies, the data collected are limited and an exact value is recorded only if it falls within an interval range. Hence, the responses can be either left, interval or right censored. Linear (and nonlinear) regression models are routinely used to analyze these types of data and are based on normality assumptions for the errors terms. However, those analyzes might not provide robust inference when the normality assumptions are questionable. In this article, we develop a Bayesian framework for censored linear regression models by replacing the Gaussian assumptions for the random errors with scale mixtures of normal (SMN) distributions. The SMN is an attractive class of symmetric heavy-tailed densities that includes the normal, Student-t, Pearson type VII, slash and the contaminated normal distributions, as special cases. Using a Bayesian paradigm, an efficient Markov chain Monte Carlo algorithm is introduced to carry out posterior inference. A new hierarchical prior distribution is suggested for the degrees of freedom parameter in the Student-t distribution. The likelihood function is utilized to compute not only some Bayesian model selection measures but also to develop Bayesian case-deletion influence diagnostics based on the q-divergence measure. The proposed Bayesian methods are implemented in the R package BayesCR. The newly developed procedures are illustrated with applications using real and simulated data.     

Coverage‐adjusted Confidence Intervals for a Binomial Proportion

We consider the classic problem of interval estimation of a proportion p based on binomial sampling. The ‘exact’ Clopper–Pearson confidence interval for p is known to be unnecessarily conservative. We propose coverage adjustments of the Clopper–Pearson interval that incorporate prior or posterior beliefs into the interval. Using heatmap‐type plots for comparing confidence intervals, we show that the coverage‐adjusted intervals have satisfying coverage and shorter expected lengths than competing intervals found in the literature.     

The X2-goodness-of-fit tests with moment type estimators

In the x²-goodness-of-fit test the underlying null hypothesis usually involves unknown parameters. In this article we study the asymptotic distribution of the Pearson statistic when the unknown parameters are estimated by a moment type estimator based on the ungrouped data. As is expected the usual Pearson statistic is no longer asymptotically x²-distributed in this situation. We propose a statistic [Qcirc] which under certain regularity conditions is asymptotically x²-distributed. We also propose a statistic Q? for the goodness-of-fit test when the class boundaries are random. The asymptotic powers of [Qcirc] and [Qcirc]? tests are discussed.     

Statistical characteristics of the surface wind speeds over Kenya

In this study attempts were made to fit several statistical distributions to the surface maximum, minimum and mean daily wind speed records at 24 Kenyan sites. The statistical distributions fitted included the 2 and 3 parameter Lognormal, Pearson type III, Log Pearson type III, the 2 and 3 parameter Weibull distributions. The various parameters for these distributions were obtained from the methods of moments and maximum likelihood. The goodness of fit of the various distributions were investigated at each of the 24 sites using the Kolmogorov-Smirnov and x2 tests. The corresponding standard errors were also computed.     

Approximate mrpp p-values obtained from four exact moments

A F0RTRAN-77 subroutine for a general version of multi-response permutation procedures (MRPP) is described. The exact four moments are employed in conjunction with the Pearson type I, type III, and type VI distributions to calculate the associated P-values.     

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.     

Bayesian (mean) most powerful tests

A fundamental theorem in hypothesis testing is the Neyman‐Pearson (N‐P) lemma, which creates the most powerful test of simple hypotheses. In this article, we establish Bayesian framework of hypothesis testing, and extend the Neyman‐Pearson lemma to create the Bayesian most powerful test of general hypotheses, thus providing optimality theory to determine thresholds of Bayes factors. Unlike conventional Bayes tests, the proposed Bayesian test is able to control the type I error.     

On the use of curve fitting to model the error of the Cornish-Fisher expansion of the pearson type-VI distribution

The Cornish-Fisher expansion of the Pearson type VI distribution is known to be reasonably accurate when both degrees of freedom are relatively large (say greater than or equal to 5). However, when either or both degrees of freedom are less than 5, the accuracy of the computed percentage point begins to suffer; in some cases severely. To correct for this, the error surface in the degrees of freedom plane is modeled by least squares curve fitting for selected levels of tail probability (.025, .05, and .10) which can be used to adjust the percentage point obtained from the usual Cornish-Fisher expansion. This adjustment procedure produces a computing algorithm that computes percentage points of the Pearson type VI distribution at the above probability levels, accurate to at least + 1 in 3 digits in approximately 11 milliseconds per subroutine call on an IBM 370/145. This adjusted routine is valid for both degrees of freedom greater than or equal to 1.     

A Comparative Study of Some Modified Chi-Squared Tests

Some recent results in the theory and applications of modified chi-squared goodness-of-fit tests are briefly discussed. It seems that for the first time power of modified chi-squared type tests for the logistic and three-parameter Weibull distributions based on moment type estimators is studied. Power of different modified tests against some alternatives for equiprobable fixed or random grouping intervals, and for Neyman–Pearson classes is investigated. It is shown that power of test statistic essentially depends on the quantity of Fisher's sample information this statistic uses. Some recommendations on implementing modified chi-squared type tests are given.     

Asymptotic convergence in distribution for spearman's rank correlation coefficient

The Edgeworth expansion for the distribution function of Spearman's rank correlation coefficient may be used to show that the rates of convergence for the normal and Pearson type II approximations are l/nand l/n² respectively. Using the Edgeworth expansion up to terms involving the sixth moment of the exact distribution allows an approximation with an error of order l/n³.     

The Pearson Diffusions: A Class of Statistically Tractable Diffusion Processes

Abstract.  The Pearson diffusions form a flexible class of diffusions defined by having linear drift and quadratic squared diffusion coefficient. It is demonstrated that for this class explicit statistical inference is feasible. A complete model classification is presented for the ergodic Pearson diffusions. The class of stationary distributions equals the full Pearson system of distributions. Well-known instances are the Ornstein–Uhlenbeck processes and the square root (CIR) processes. Also diffusions with heavy-tailed and skew marginals are included. Explicit formulae for the conditional moments and the polynomial eigenfunctions are derived. Explicit optimal martingale estimating functions are found. The discussion covers GMM, quasi-likelihood, non-linear weighted least squares estimation and likelihood inference too. The analytical tractability is inherited by transformed Pearson diffusions, integrated Pearson diffusions, sums of Pearson diffusions and Pearson stochastic volatility models. For the non-Markov models, explicit optimal prediction-based estimating functions are found. The estimators are shown to be consistent and asymptotically normal.