Goodness‐of‐Fit Tests for Bivariate and Multivariate Skew‐Normal Distributions

Abstract. Goodness‐of‐fit tests are proposed for the skew‐normal law in arbitrary dimension. In the bivariate case the proposed tests utilize the fact that the moment‐generating function of the skew‐normal variable is quite simple and satisfies a partial differential equation of the first order. This differential equation is estimated from the sample and the test statistic is constructed as an L ₂ ‐type distance measure incorporating this estimate. Extension of the procedure to dimension greater than two is suggested whereas an effective bootstrap procedure is used to study the behaviour of the new method with real and simulated data.     

On batch queueing systems: A combinatorial approach

We consider batch queueing systems M/M^H/1$M/M^H/1$ and M^H/M/1$M^H/M/1$ with catastrophes. The transient probability functions of these queueing systems are obtained by a Lattice Path Combinatorics approach that utilizes randomization and dual processes. Steady state distributions are also determined. Generalization to systems having batches of different sizes are discussed.     

Data-driven sensitivity analysis to detect missing data mechanism with applications to structural equation modelling

Missing data are a common problem in almost all areas of empirical research. Ignoring the missing data mechanism, especially when data are missing not at random (MNAR), can result in biased and/or inefficient inference. Because MNAR mechanism is not verifiable based on the observed data, sensitivity analysis is often used to assess it. Current sensitivity analysis methods primarily assume a model for the response mechanism in conjunction with a measurement model and examine sensitivity to missing data mechanism via the parameters of the response model. Recently, Jamshidian and Mata (Post-modelling sensitivity analysis to detect the effect of missing data mechanism, Multivariate Behav. Res. 43 (2008), pp. 432–452) introduced a new method of sensitivity analysis that does not require the difficult task of modelling the missing data mechanism. In this method, a single measurement model is fitted to all of the data and to a sub-sample of the data. Discrepancy in the parameter estimates obtained from the the two data sets is used as a measure of sensitivity to missing data mechanism. Jamshidian and Mata describe their method mainly in the context of detecting data that are missing completely at random (MCAR). They used a bootstrap type method, that relies on heuristic input from the researcher, to test for the discrepancy of the parameter estimates. Instead of using bootstrap, the current article obtains confidence interval for parameter differences on two samples based on an asymptotic approximation. Because it does not use bootstrap, the developed procedure avoids likely convergence problems with the bootstrap methods. It does not require heuristic input from the researcher and can be readily implemented in statistical software. The article also discusses methods of obtaining sub-samples that may be used to test missing at random in addition to MCAR. An application of the developed procedure to a real data set, from the first wave of an ongoing longitudinal study on aging, is presented. Simulation studies are performed as well, using two methods of missing data generation, which show promise for the proposed sensitivity method. One method of missing data generation is also new and interesting in its own right.     

Nonparametric Random Variate Generation Using a Piecewise-Linear Cumulative Distribution Function

The standard approach to solving the interpolation problem for a trace-driven simulation involving a continuous random variable is to construct a piecewise-linear cdf that fills in the gaps between the data values. Some probabilistic properties of this estimator are derived, and three extensions to the standard approach (matching moments, weighted values, and right-censored data) are presented, along with associated random variate generation algorithms. The algorithm is a nonparametric blackbox variate generator requiring only observed data from the user.     

Stability of the roots of random algebraic polynomials

In this paper we consider the behavior of the roots of random algebraic polynomials. A code was developed which generates a sample of random algebraic polynomials, calculates the roots of each sample polynomial, and then calculates the averages of the roots. Finally, the roots of the deterministic algebraic polynomial whose coefficients are the averages of the sample coefficients are calculated. These data are then tabulated and graphically displayed. The relationship between the averages of the roots of the sample polynomials and the roots of the average polynomial is discussed.     

Relative efficiency of a parameter for a M/G/1 queueing system based on reduced and full likelihood functions

The present paper derives the relative efficiency of a parameter for the M/G/1 queueing system based on reduced and full likelihood functions. Monte Carlo simulations were carried out to study the finite sample properties for estimating the parameters of a M/G/1 queueing system. The simulation runs were conducted using various traftic intensities with increaseing sample sizes. The simulation results indicate that the loss in efficiency is quite small due to the use of a reduced likelihood function approach for estimating the parameter instead of the full likelihood, even for a moderate sample size of 50     

Multi-response permutation procedures for a priori classifications

An exact permutation test for analyzing and/or dredging multi-response data at the ordinal or higher levels is presented. The associated test statistic is based on the average distance (or any specified norm) between points within a priori disjoint subgroups of a finite population of points in an r-dimensional space (corresponding to r measured responses from each object in a finite population of objects). Alternative approximate tests based on the beta and normal distributions are provided. Two detailed examples utilizing actual social science data are considered, including comparisons of the approximate tests. An additional example describes the behavior of these tests under a variety of conditions, including extreme data configurations     

Testing the Fit of Gamma Distributions Using the Empirical Moment Generating Function

ABSTRACT

This article presents goodness-of-fit tests for two and three-parameter gamma distributions that are based on minimum quadratic forms of standardized logarithmic differences of values of the moment generating function and its empirical counterpart. The test statistics can be computed without reliance to special functions and have asymptotic chi-squared distributions. Monte Carlo simulations are used to compare the proposed test for the two-parameter gamma distribution with goodness-of-fit tests employing empirical distribution function or spacing statistics. Two data sets are used to illustrate the various tests.