Model-based clustering with non-elliptically contoured distributions

The majority of the existing literature on model-based clustering deals with symmetric components. In some cases, especially when dealing with skewed subpopulations, the estimate of the number of groups can be misleading; if symmetric components are assumed we need more than one component to describe an asymmetric group. Existing mixture models, based on multivariate normal distributions and multivariate t distributions, try to fit symmetric distributions, i.e. they fit symmetric clusters. In the present paper, we propose the use of finite mixtures of the normal inverse Gaussian distribution (and its multivariate extensions). Such finite mixture models start from a density that allows for skewness and fat tails, generalize the existing models, are tractable and have desirable properties. We examine both the univariate case, to gain insight, and the multivariate case, which is more useful in real applications. EM type algorithms are described for fitting the models. Real data examples are used to demonstrate the potential of the new model in comparison with existing ones.     

Penalized likelihood-ratio test for finite mixture models with multinomial observations

Due to the irregularity of finite mixture models, the commonly used likelihood-ratio statistics often have complicated limiting distributions. We propose to add a particular type of penalty function to the log-likelihood function. The resulting penalized likelihood-ratio statistics have simple limiting distributions when applied to finite mixture models with multinomial observations. The method is especially effective in addressing the problems discussed by Chernoff and Lander (1995). The theory developed and simulations conducted show that the penalized likelihood method can give very good results, better than the well-known C(α) procedure, for example. The paper does not, however, fully explore the choice of penalty function and weight. The full potential of the new procedure is to be explored in the future.     

Modeling lifetime of sequential r-out-of-n systems with independent and heterogeneous components

This article deals with an extension of sequential order statistics which is useful for describing system lifetimes with independent but heterogeneous components. Explicit expressions for marginal distributions as well means of system lifetimes are derived. Some special cases and illustrative examples are also investigated.     

A fiducial-based approach to the one-way ANOVA in the presence of nonnormality and heterogeneous error variances

In this study, we propose a new test for testing the equality of the treatment means in one-way ANOVA when the usual normality and the homogeneity of variances assumptions are not met. In developing the proposed test, we benefit from the Fisher's fiducial inference [1–3]. Distribution of the error terms is assumed to be long-tailed symmetric (LTS) which includes the normal distribution as a limiting case. Modified maximum likelihood (MML) estimators are used in the test statistics rather than the traditional least squares (LS) estimators, since LS estimators have very low efficiencies under nonnormal distributions, see Tiku [4] for the details of MML methodology. An extensive Monte Carlo simulation study is done to compare the efficiency of the proposed test with the corresponding test based on normal theory, see Li et al. [5]. Finally, we give a real life example to show the applicability of the proposed methodology.     

Testing goodness-of-fit for some lifetime distributions with conventional Type-I censoring

A goodness-of-fit test procedure is proposed for some lifetime distributions when the available data are subject to Type-I censoring. The proposed method extends the test procedure of Pakyari and Balakrishnan to other lifetime distributions. The extension to Weibull and log-normal models is studied in details. The new test recovers the nominal level of significance and exhibits more power in comparison to the existing tests for several alternative distributions by means of Monte Carlo simulations. Finally, a real dataset is considered for illustrative purposes.     

A unified view on skewed distributions arising from selections

Parametric families of multivariate nonnormal distributions have received considerable attention in the past few decades. The authors propose a new definition of a selection distribution that encompasses many existing families of multivariate skewed distributions. Their work is motivated by examples that involve various forms of selection mechanisms and lead to skewed distributions. They give the main properties of selection distributions and show how various families of multivariate skewed distributions, such as the skew‐normal and skew‐elliptical distributions, arise as special cases. The authors further introduce several methods of constructing selection distributions based on linear and nonlinear selection mechanisms.     

Gompertz-power series distributions

ABSTRACT

In this article, we introduce the Gompertz power series (GPS) class of distributions which is obtained by compounding Gompertz and power series distributions. This distribution contains several lifetime models such as Gompertz-geometric (GG), Gompertz-Poisson (GP), Gompertz-binomial (GB), and Gompertz-logarithmic (GL) distributions as special cases. Sub-models of the GPS distribution are studied in details. The hazard rate function of the GPS distribution can be increasing, decreasing, and bathtub-shaped. We obtain several properties of the GPS distribution such as its probability density function, and failure rate function, Shannon entropy, mean residual life function, quantiles, and moments. The maximum likelihood estimation procedure via a EM-algorithm is presented, and simulation studies are performed for evaluation of this estimation for complete data, and the MLE of parameters for censored data. At the end, a real example is given.     

A Novel Asymmetric Distribution with Power Tails

In this article, we propose a four-parameter asymmetric doubly Pareto-uniform (DPU) distribution with support (?∞, ∞) whose density and cumulative distribution functions are constructed by seamlessly concatenating the left and right Pareto tails with a uniform central part. Properties of the distribution are described and a maximum likelihood estimation (MLE) procedure for its parameters is obtained. Two illustrative examples of the MLE procedure are provided. The first example utilizes an i.i.d. sample of standardized log-differences of bi-monthly 30-year U.S. conventional mortgage interest rates (1971–2004). The second example deals with the height of 100 female Australian athletes.     

A Comparison of Bayes Factors for Separated Models: Some Simulation Results

Some alternative Bayes Factors: Intrinsic, Posterior, and Fractional have been proposed to overcome the difficulties presented when prior information is weak and improper prior are used. Additional difficulties also appear when the models are separated or non nested. This article presents both simulation results and some illustrative examples analysis comparing these alternative Bayes factors to discriminate among the Lognormal, the Weibull, the Gamma, and the Exponential distributions. Simulation results are obtained for different sample sizes generated from the distributions. Results from simulations indicates that these alternative Bayes factors are useful for comparing non nested models. The simulations also show some similar behavior and that when both models are true they choose the simplest model. Some illustrative example are also presented.     

Tests of goodness of fit based on Phi-divergence

In this paper, we introduce a general goodness of fit test based on Phi-divergence. Consistency of the proposed test is established. We then study some special cases of tests for normal, exponential, uniform and Laplace distributions. Through Monte Carlo simulations, the power values of the proposed tests are compared with some known competing tests under various alternatives. Finally, some numerical examples are presented to illustrate the proposed procedure.     

A lifetime distribution motivated by parallel and series structures

A new three-parameter distribution with decreasing, increasing, and bathtub-shaped hazard rates obtained by compounding geometric, power series, and exponential distributions is introduced. It includes some well-known distributions as particular cases. Various mathematical properties of the new distribution as well as details of the maximum likelihood estimation and a sensitivity analysis for its parameters are presented. Finally, two real data applications are presented.     

Tail index varying coefficient model

This paper deals with a new class of tail index varying coefficient models with the random covariate under Pareto-type distributions. To estimate the unknown coefficient functions, we develop an estimation procedure via a local polynomial maximum likelihood techniques. The asymptotic normality of the estimated coefficient functions under some mild regularity conditions are established. Two numerical examples and one application are used to illustrate the performance of the proposed procedure.     

A unified approach to estimating parameters in some generalized poission distributions

A family of distributions labelled as Poisson v Katz is formulated, which includes, as particular or limiting cases, the Negative Binomial, Neyman Type A, Poisson v Pascal, and Poisson v Binomial. Thus, while analyzing data, estimating the parameters in the Poisson v Katz family obviates the necessity of having to choose from among the particular or limiting cases. In this article minimum chi-square estimators are presented and their asymptotic relative efficiency obtained. An example is presented to illustrate the procedure     

Testing linear inequality constraints in the standard linear model

In this paper we are concerned with the problem of testing whether the â-parameters of the standard linear model satisfy the linear equality constraints R = r when they are known to satisfy the corresponding linear inequality constraints Râ ? r. In particular we will show that the exact finite sample null distributions of the Likelihood Ratio, Wald and Kuhn-Tucker

statistics are known when R is of full row rank but not known when R has less than full row rank. The less than full row rank problem has not been discussed previously but it is of considerable potential importance.

This paper contains several simple numerical examples which illustrate the computational details of the tests     

Nonparametric test for the homogeneity of the overall variability

In this paper, we propose a nonparametric test for homogeneity of overall variabilities for two multi-dimensional populations. Comparisons between the proposed nonparametric procedure and the asymptotic parametric procedure and a permutation test based on standardized generalized variances are made when the underlying populations are multivariate normal. We also study the performance of these test procedures when the underlying populations are non-normal. We observe that the nonparametric procedure and the permutation test based on standardized generalized variances are not as powerful as the asymptotic parametric test under normality. However, they are reliable and powerful tests for comparing overall variability under other multivariate distributions such as the multivariate Cauchy, the multivariate Pareto and the multivariate exponential distributions, even with small sample sizes. A Monte Carlo simulation study is used to evaluate the performance of the proposed procedures. An example from an educational study is used to illustrate the proposed nonparametric test.     

Inference of the derivative of nonparametric curve based on confidence distribution

Abstract

This paper focuses on inference based on the confidence distributions of the nonparametric regression function and its derivatives, in which dependent inferences are combined by obtaining information about their dependency structure. We first give a motivating example in production operation system to illustrate the necessity of the problems studied in this paper in practical applications. A goodness-of-fit test for polynomial regression model is proposed on the basis of the idea of combined confidence distribution inference, which is the Fisher’s combination statistic in some cases. On the basis of this testing results, a combined estimator for the p-order derivative of nonparametric regression function is provided as well as its large sample size properties. Consequently, the performances of the proposed test and estimation method are illustrated by three specific examples. Finally, the motivating example is analyzed in detail. The simulated and real data examples illustrate the good performance and practicability of the proposed methods based on confidence distribution.     

Confidence intervals based on L-moments for quantiles of the GP and GEV distributions with application to market-opening asset prices data

In a ground-breaking paper published in 1990 by the Journal of the Royal Statistical Society, J.R.M. Hosking defined the L-moment of a random variable as an expectation of certain linear combinations of order statistics. L-moments are an alternative to conventional moments and recently they have been used often in inferential statistics. L-moments have several advantages over the conventional moments, including robustness to the the presence of outliers, which may lead to more accurate estimates in some cases as the characteristics of distributions. In this contribution, asymptotic theory and L-moments are used to derive confidence intervals of the population parameters and quantiles of the three-parametric generalized Pareto and extreme-value distributions. Computer simulations are performed to determine the performance of confidence intervals for the population quantiles based on L-moments and to compare them to those obtained by traditional estimation techniques. The results obtained show that they perform well in comparison to the moments and maximum likelihood methods when the interest is in higher quantiles, or even best. L-moments are especially recommended when the tail of the distribution is rather heavier and the sample size is small. The derived intervals are applied to real economic data, and specifically to market-opening asset prices.     

On a Mixture Model for Directional Data on the Sphere

We consider mixtures of general angular central Gaussian distributions as models for multimodal directional data. We prove consistency of the maximum‐likelihood estimates of model parameters and convergence of their numerical approximations based on an expectation–maximization algorithm. Then, we focus on mixtures of special angular central Gaussian distributions and discuss the details of a fast numerical algorithm, which allows to fit multimodal distributions to massive data, occurring, for example, in the study of the microstructure of materials. We illustrate the applicability with some data from fibre composites and from ceramic foams.     

Goodness-of-fit-based outlier detection for Phase I monitoring

ABSTRACT

Nonparametric charts are useful in statistical process control when there is a lack of or limited knowledge about the underlying process distribution. Most existing approaches in the literature of Phase I monitoring assume that outliers have the same distributions as the in-control sample but only differ in location or scale parameters, they may not be effective with distributional changes. This article develops a new procedure based on the integration of the classical Anderson–Darling goodness-of-fit test and the stepwise isolation method. Our proposed procedure is efficient in detecting potential shifts in location, scale, or shape, and thus it offers robust protection against variation in various underlying distributions. The finite sample performance of our method is evaluated through simulations and is compared with that of available outlier detection methods for Phase I monitoring.     

On discrete distributions generated from drawing balls with bivariate labels

In this article, we consider urn models under three types of sampling schemes in terms of the probability-generating functions. The tools are developed for the evaluation of the distributions arising from the urn models along with some examples. Furthermore, the distributions are investigated by making use of the Bell polynomials. The results presented here provide a wide framework for developing the theory of urn models. As examples, we propose new class of probability models, which are called multiple-player problems and examine their properties. Finally, we treat the parameter estimation problem in the waiting time distributions with a numerical example.