Tests for outliers in the inverse Gaussian distribution,with application to first hitting time models

The inverse Gaussian (IG) distribution is often applied in statistical modelling, especially with lifetime data. We present tests for outlying values of the parameters (μ, λ) of this distribution when data are available from a sample of independent units and possibly with more than one event per unit. Outlier tests are constructed from likelihood ratio tests for equality of parameters. The test for an outlying value of λ is based on an F-distributed statistic that is transformed to an approximate normal statistic when there are unequal numbers of events per unit. Simulation studies are used to confirm that Bonferroni tests have accurate size and to examine the powers of the tests. The application to first hitting time models, where the IG distribution is derived from an underlying Wiener process, is described. The tests are illustrated on data concerning the strength of different lots of insulating material.     

The Likelihood Ratio Tests for Homogeneity of Inverse Gaussian Means Under Simple Order and Simple Tree Order

Abstract

The inverse Gaussian (IG) family is now widely used for modeling non negative skewed measurements. In this article, we construct the likelihood-ratio tests (LRTs) for homogeneity of the order constrained IG means and study the null distributions for simple order and simple tree order cases. Interestingly, it is seen that the null distribution results for the normal case are applicable without modification to the IG case. This supplements the numerous well known and striking analogies between Gaussian and inverse Gaussian families     

IG-symmetry and R-symmetry: Interrelations and applications to the inverse Gaussian theory

The two parameter inverse Gaussian (IG) distribution is often more appropriate and convenient for modelling and analysis of nonnegative right skewed data than the better known and now ubiquitous Gaussian distribution. Its convenience stems from its analytic simplicity and the striking similarities of its methodologies with those employed with the normal theory models. These, known as the G–IG analogies, include the concepts and measures of IG-symmetry, IG-skewness and IG-kurtosis, the IG-analogues of the corresponding classical notions and measures. The new IG-associated entities, although well defined and mathematically transparent, are intuitively and conceptually opaque. In this paper, we first elaborate the importance of the IG distribution and of the G–IG analogies. Then we consider the IG-related root-reciprocal IG (RRIG) distribution and introduce a physically transparent, conceptually clear notion of reciprocal symmetry (R-symmetry) and use it to explain the IG-symmetry. We study the moments and mixture properties of the R-symmetric distributions and the relationship of R-symmetry with IG-symmetry and note that RRIG distribution provides a link, in addition to Tweedie's Laplace transform link, between the Gaussian and inverse Gaussian distributions. We also give a structural characterization of the unimodal R-symmetric distributions. This work further expands the long list of G–IG analogies. Several applications including product convolution, monotonicity of power functions, peakedness and monotone limit theorems of R-symmetry are outlined.     

A modified Kolmogorov-Smirnov test for the inverse gaussian density with unknown parameters

The Kolmogorov-Smirnov (KS) test is an empirical distribution function (EDF) based goodness-of-fit test that requires the underlying hypothesized density to be continuous and completely specified. When the parameters are unknown and must be estimated from the data, standard tables of the KS test statistic are not valid. Approximate upper tail percentage points of the KS statistic for the inverse Gaussian (IG) distribution with unknown parameters are tabled in this paper.

A study of the power of the KS test for the IG distribution indicates that the test is able todiscriminate between the IG distribution and distributions such as the uniform and exponentialdistributions that are very different in shape, but is relatively unable to discriminate between the IG distribution and distributions that are similar in shape such as the lognormal and Weibull distributions. In modeling settings the former distinction is typically more important to make than the latter distinction.     

On the bias and mean-square error of order-restricted maximum likelihood estimators

The inverse Gaussian family (IG) (μ,λ) is a versatile family for modelling nonnegative right-skewed data. In this paper, we propose robust methods for testing homogeneity of the scale-like parameters λ_i from k independent IG populations subject to order restrictions. Robustness of the procedures is examined for a variety of IG-symmetric alternatives including lognormal and the recently introduced contaminated inverse Gaussian populations. Our study shows that these inference procedures for the inverse Gaussian scale-like parameters and their properties exhibit striking similarities to those of the scale parameters of the normal distribution.     

Preliminary Goodness-of-Fit Tests for Normality do not Validate the One-Sample Student t

One of the most basic topics in many introductory statistical methods texts is inference for a population mean, μ. The primary tool for confidence intervals and tests is the Student t sampling distribution. Although the derivation requires independent identically distributed normal random variables with constant variance, σ², most authors reassure the readers about some robustness to the normality and constant variance assumptions. Some point out that if one is concerned about assumptions, one may statistically test these prior to reliance on the Student t. Most software packages provide optional test results for both (a) the Gaussian assumption and (b) homogeneity of variance. Many textbooks advise only informal graphical assessments, such as certain scatterplots for independence, others for constant variance, and normal quantile–quantile plots for the adequacy of the Gaussian model. We concur with this recommendation. As convincing evidence against formal tests of (a), such as the Shapiro–Wilk, we offer a simulation study of the tails of the resulting conditional sampling distributions of the Studentized mean. We analyze the results of systematically screening all samples from normal, uniform, exponential, and Cauchy populations. This pretest does not correct the erroneous significance levels and makes matters worse for the exponential. In practice, we conclude that graphical diagnostics are better than a formal pretest. Furthermore, rank or permutation methods are recommended for exact validity in the symmetric case.     

Moments of nonparametric probability density functions of entropy estimators applied to testing the inverse Gaussian distribution

The inverse Gaussian (IG) distribution is widely used to model data and then it is important to develop efficient goodness of fit tests for this distribution. In this article, we introduce some new test statistics for examining the IG goodness of fit based on correcting moments of nonparametric probability density functions of entropy estimators. These tests are consistent against all alternatives. Critical points and power of the tests are explored by simulation. We show that the proposed tests are more powerful than competitor tests. Finally, the proposed tests are illustrated by real data examples.     

A Bayesian analysis for the Wilcoxon signed-rank statistic

A Bayesian analysis is provided for the Wilcoxon signed-rank statistic (T⁺). The Bayesian analysis is based on a sign-bias parameter φ on the (0, 1) interval. For the case of a uniform prior probability distribution for φ and for small sample sizes (i.e., 6 ? n ? 25), values for the statistic T⁺ are computed that enable probabilistic statements about φ. For larger sample sizes, approximations are provided for the asymptotic likelihood function P(T⁺|φ) as well as for the posterior distribution P(φ|T⁺). Power analyses are examined both for properly specified Gaussian sampling and for misspecified non Gaussian models. The new Bayesian metric has high power efficiency in the range of 0.9–1 relative to a standard t test when there is Gaussian sampling. But if the sampling is from an unknown and misspecified distribution, then the new statistic still has high power; in some cases, the power can be higher than the t test (especially for probability mixtures and heavy-tailed distributions). The new Bayesian analysis is thus a useful and robust method for applications where the usual parametric assumptions are questionable. These properties further enable a way to do a generic Bayesian analysis for many non Gaussian distributions that currently lack a formal Bayesian model.     

Survival analysis for the inverse Gaussian distribution with the Gibbs sampler

This paper describes a comprehensive survival analysis for the inverse Gaussian distribution employing Bayesian and Fiducial approaches. It focuses on making inferences on the inverse Gaussian (IG) parameters μ and λ and the average remaining time of censored units. A flexible Gibbs sampling approach applicable in the presence of censoring is discussed and illustrations with Type II, progressive Type II, and random rightly censored observations are included. The analyses are performed using both simulated IG data and empirical data examples. Further, the bootstrap comparisons are made between the Bayesian and Fiducial estimates. It is concluded that the shape parameter ( ϕ=λ/μ) of the inverse Gaussian distribution has the most impact on the two analyses, Bayesian vs. Fiducial, and so does the size of censoring in data to a lesser extent. Overall, both these approaches are effective in estimating IG parameters and the average remaining lifetime. The suggested Gibbs sampler allowed a great deal of flexibility in implementation for all types of censoring considered.     

An efficient correction to the density-based empirical likelihood ratio goodness-of-fit test for the inverse Gaussian distribution

The inverse Gaussian (IG) distribution is widely used to model positively skewed data. An important issue is to develop a powerful goodness-of-fit test for the IG distribution. We propose and examine novel test statistics for testing the IG goodness of fit based on the density-based empirical likelihood (EL) ratio concept. To construct the test statistics, we use a new approach that employs a method of the minimization of the discrimination information loss estimator to minimize Kullback–Leibler type information. The proposed tests are shown to be consistent against wide classes of alternatives. We show that the density-based EL ratio tests are more powerful than the corresponding classical goodness-of-fit tests. The practical efficiency of the tests is illustrated by using real data examples.     

A new look at the inverse Gaussian distribution with applications to insurance and economic data

Insurance and economic data are often positive, and we need to take into account this peculiarity in choosing a statistical model for their distribution. An example is the inverse Gaussian (IG), which is one of the most famous and considered distributions with positive support. With the aim of increasing the use of the IG distribution on insurance and economic data, we propose a convenient mode-based parameterization yielding the reparametrized IG (rIG) distribution; it allows/simplifies the use of the IG distribution in various branches of statistics, and we give some examples. In nonparametric statistics, we define a smoother based on rIG kernels. By construction, the estimator is well-defined and does not allocate probability mass to unrealistic negative values. We adopt likelihood cross-validation to select the smoothing parameter. In robust statistics, we propose the contaminated IG distribution, a heavy-tailed generalization of the rIG distribution to accommodate mild outliers. Finally, for model-based clustering and semiparametric density estimation, we present finite mixtures of rIG distributions. We use the EM algorithm to obtain maximum likelihood estimates of the parameters of the mixture and contaminated models. We use insurance data about bodily injury claims, and economic data about incomes of Italian households, to illustrate the models.     

A comparative simulation study of the small sample powers of several goodness fit test

Eight goodness of fit tests are compared with respect to their simulated small sample power of detecting an inbreeding alternative to the Hardy-Weinberg null hypothesis. The Pearson's x ² test is found to be most powerful, and the small rample levels of this test are close to the nominal (x ²) significance levels. The use of conditional expectations, rather than expected frequencies based on ML estimates, increases the power and improves thc x ² fit to the true significance level. The small sample powers are also compared to the asymptotic (Pitman) pourer, based on the noncenlral x ² distribution.     

A Test for Slope Heterogeneity in Fixed Effects Models

Typical panel data models make use of the assumption that the regression parameters are the same for each individual cross-sectional unit. We propose tests for slope heterogeneity in panel data models. Our tests are based on the conditional Gaussian likelihood function in order to avoid the incidental parameters problem induced by the inclusion of individual fixed effects for each cross-sectional unit. We derive the Conditional Lagrange Multiplier test that is valid in cases where N → ∞ and T is fixed. The test applies to both balanced and unbalanced panels. We expand the test to account for general heteroskedasticity where each cross-sectional unit has its own form of heteroskedasticity. The modification is possible if T is large enough to estimate regression coefficients for each cross-sectional unit by using the MINQUE unbiased estimator for regression variances under heteroskedasticity. All versions of the test have a standard Normal distribution under general assumptions on the error distribution as N → ∞. A Monte Carlo experiment shows that the test has very good size properties under all specifications considered, including heteroskedastic errors. In addition, power of our test is very good relative to existing tests, particularly when T is not large.     

On the use of priors in goodness‐of‐fit tests

Priors are introduced into goodness‐of‐fit tests, both for unknown parameters in the tested distribution and on the alternative density. Neyman–Pearson theory leads to the test with the highest expected power. To make the test practical, we seek priors that make it likely a priori that the power will be larger than the level of the test but not too close to one. As a result, priors are sample size dependent. We explore this procedure in particular for priors that are defined via a Gaussian process approximation for the logarithm of the alternative density. In the case of testing for the uniform distribution, we show that the optimal test is of the U‐statistic type and establish limiting distributions for the optimal test statistic, both under the null hypothesis and averaged over the alternative hypotheses. The optimal test statistic is shown to be of the Cramér–von Mises type for specific choices of the Gaussian process involved. The methodology when parameters in the tested distribution are unknown is discussed and illustrated in the case of testing for the von Mises distribution. The Canadian Journal of Statistics 47: 560–579; 2019 © 2019 Statistical Society of Canada     

A hierarchical Bayesian regression model for the uncertain functional constraint using screened scale mixtures of Gaussian distributions

This paper considers a hierarchical Bayesian analysis of regression models using a class of Gaussian scale mixtures. This class provides a robust alternative to the common use of the Gaussian distribution as a prior distribution in particular for estimating the regression function subject to uncertainty about the constraint. For this purpose, we use a family of rectangular screened multivariate scale mixtures of Gaussian distribution as a prior for the regression function, which is flexible enough to reflect the degrees of uncertainty about the functional constraint. Specifically, we propose a hierarchical Bayesian regression model for the constrained regression function with uncertainty on the basis of three stages of a prior hierarchy with Gaussian scale mixtures, referred to as a hierarchical screened scale mixture of Gaussian regression models (HSMGRM). We describe distributional properties of HSMGRM and an efficient Markov chain Monte Carlo algorithm for posterior inference, and apply the proposed model to real applications with constrained regression models subject to uncertainty.     

Transformations for Testing the Fit of the Inverse-Gaussian Distribution

Abstract

Use of the MVUE for the inverse-Gaussian distribution has been recently proposed by Nguyen and Dinh [Nguyen, T. T., Dinh, K. T. (2003). Exact EDF goodnes-of-fit tests for inverse Gaussian distributions. Comm. Statist. (Simulation and Computation) 32(2):505–516] where a sequential application based on Rosenblatt's transformation [Rosenblatt, M. (1952). Remarks on a multivariate transformation. Ann. Math. Statist. 23:470–472] led the authors to solve the composite goodness-of-fit problem by solving the surrogate simple goodness-of-fit problem, of testing uniformity of the independent transformed variables. In this note, we observe first that the proposal is not new since it was proposed in a rather general setting in O'Reilly and Quesenberry [O'Reilly, F., Quesenberry, C. P. (1973). The conditional probability integral transformation and applications to obtain composite chi-square goodness-of-fit tests. Ann. Statist. I:74–83]. It is shown on the other hand that the results in the paper of Nguyen and Dinh (2003) are incorrect in their Sec. 4, specially the Monte Carlo figures reported. Power simulations are provided here comparing these corrected results with two previously reported goodness-of-fit tests for the inverse-Gaussian; the modified Kolmogorov–Smirnov test in Edgeman et al. [Edgeman, R. L., Scott, R. C., Pavur, R. J. (1988). A modified Kolmogorov-Smirnov test for inverse Gaussian distribution with unknown parameters. Comm. Statist. 17(B): 1203–1212] and the A ² based method in O'Reilly and Rueda [O'Reilly, F., Rueda, R. (1992). Goodness of fit for the inverse Gaussian distribution. T Can. J. Statist. 20(4):387–397]. The results show clearly that there is a large loss of power in the method explored in Nguyen and Dinh (2003) due to an implicit exogenous randomization.     

Power and Sample Size Determination When the Alternative Hypotheses are Given in Quantiles

The power of normal-theory tests about means depends on a noncentrality parameter which is a function of the unknown parameter σ. In order to calculate power and to solve sample-size problems based on power, differences between hypothesized and alternative values of the means are frequently selected as a multiple of σ, a choice which eliminates σ from the noncentrality parameter and permits a solution. Perhaps a more natural (but equivalent) way to express alternatives is to give one or more means as the quantile of order p (say Q_p ) of a distribution with another mean. As we will demonstrate, this kind of alternative also eliminates σ from the problem.     

News about Members

In this article a natural extension of the beta-binomial distribution is developed. Forced binary choice situations are modeled such that each individual has a probability p of knowing the correct answer. (This probability is distributed f(p) across the population.) Hence each individual will guess at the correct answer with probability 1 – p. The observable random variable R, the total number of correct answers (both by knowing and guessing) out of k trials has a rather complicated distribution. However, when f(p) is distributed beta with parameters m and n, the distribution P(r; k, m, n) can be expressed in terms of the well-known Gaussian hypergeometric function. This distribution has implications for true-false tests, taste tests, and virtually every other forced binary choice situation.     

On bivariate transformation of scale distributions

ABSTRACT

Elsewhere, I have promoted (univariate continuous) “transformation of scale” (ToS) distributions having densities of the form 2g(Π^?1(x)) where g is a symmetric distribution and Π is a transformation function with a special property. Here, I develop bivariate (readily multivariate) ToS distributions. Univariate ToS distributions have a transformation of random variable relationship with Azzalini-type skew-symmetric distributions; the bivariate ToS distribution here arises from marginal variable transformation of a particular form of bivariate skew-symmetric distribution. Examples are given, as are basic properties—unimodality, a covariance property, random variate generation—and connections with a bivariate inverse Gaussian distribution are pointed out.     

Testing Equality of Means Under Order Restrictions Using Multiple Contrasts

The likelihood ratio test for equality of ordered means is known to have power characteristics that are generally superior to those of competing procedures. Difficulties in implementing this test have led to the development of alternative approaches, most of which are based on contrasts. While orthogonal contrasts can be chosen to simplify the distribution theory, we propose a class of tests that is easy to implement even if the contrasts used are not orthogonal. An overall measure of significance may be obtained by using Fisher's combination statistic to combine the dependent p-values arising from these contrasts. This method can be easily implemented for testing problems involving unequal sample sizes and any partial order, and has power properties that compare well with those of the likelihood ratio test and other contrast-based tests.