On A method for selecting a distributional model

A method for selecting a distributional model for a random variable, given a random sample of observations of it, is studied for various cases. The problems considered include those of choosing between the Weibull and lognormal distributions, between the lognormal and gamma distributions, and between the gamma and Weibull distributions, as well as choosing one of the three. Simulation studies were performed to estimate probabilities of correct selection for the method when it is applied to these problems     

On compounded geometric distributions and their applications

Here, we introduce two-parameter compounded geometric distributions with monotone failure rates. These distributions are derived by compounding geometric distribution and zero-truncated Poisson distribution. Some statistical and reliability properties of the distributions are investigated. Parameters of the proposed distributions are estimated by the maximum likelihood method as well as through the minimum distance method of estimation. Performance of the estimates by both the methods of estimation is compared based on Monte Carlo simulations. An illustration with Air Crash casualties demonstrates that the distributions can be considered as a suitable model under several real situations.     

Reliability study of proportional odds family of discrete distributions

The proportional odds model gives a method of generating new family of distributions by adding a parameter, called tilt parameter, to expand an existing family of distributions. The new family of distributions so obtained is known as Marshall–Olkin family of distributions or Marshall–Olkin extended distributions. In this paper, we consider Marshall–Olkin family of distributions in discrete case with fixed tilt parameter. We study different ageing properties, as well as different stochastic orderings of this family of distributions. All the results of this paper are supported by several examples.     

On the Estimation of Parameters of Kumaraswamy-G Distributions

Cordeiro and de Castro proposed a new family of generalized distributions based on the Kumaraswamy distribution (denoted as Kw-G). Nadarajah et al. showed that the density function of the new family of distributions can be expressed as a linear combination of the density of exponentiated family of distributions. They derived some properties of Kw-G distributions and discussed estimation of parameters using the maximum likelihood (ML) method. Cheng and Amin and Ranneby introduced a new method of estimating parameters based on Kullback–Leibler divergence (the maximum spacing (MSP) method). In this article, the estimates of parameters of Kw-G distributions are obtained using the MSP method. For some special Kw-G distributions, the new estimators are compared with ML estimators. It is shown by simulations and a real data application that MSP estimators have better properties than ML estimators.     

Partially-adaptive robust estimators of location via exponential embedding

A general method is presented for constructing a location estimator which is asymptotically efficient at any two different location-scale families of symmetric distributions as well as at an appropriately defined class of distributions lying in between. The method works by embedding the two families in a comprehensive parametric model and identifying the estimator with the MLE. The case when the families are Normal and Double exponential is examined in detail.     

Lifetime distributions motivated by series and parallel structures

According to Ross, any system can be represented either as a series arrangement of parallel structures or as a parallel arrangement of series structures. Motivated by this, we propose new three-parameter lifetime distributions by compounding geometric, power series, and exponential distributions. The distributions can allow for decreasing, increasing, bathtub-shaped, and upside down bathtub-shaped hazard rates. A mathematical treatment of the new distributions is provided including expressions for their density functions, Shannon and Rényi entropies, mean residual life functions, hazard rate functions, quantiles, and moments. The method of maximum likelihood is used for estimating parameters. Five of the new distributions are studied in detail. Finally, two illustrative data examples and a sensitivity analysis are presented.     

Properties and Inference for a New Class of Skew-t Distributions

In this article we introduce a new generalization of skew-t distributions, which contains the standard skew-t distribution, as a special case. This new class of distributions is an adequate model for modeling some dataset rather than the standard skew-t distributions. This kind of distributions can be represented as a scale-shape mixture of the extended skew-normal distributions. The main properties of this family of distributions are studied and a recurrence relation for the cumulative distribution functions (cdf) of them is presented. We derive the distribution of the order statistics from the trivariate exchangeable t-distribution in terms of our distribution and then an exact expression for the cdf of order statistics is derived. Likelihood inference for this distribution is also examined. The method is illustrated with a numerical example via a simulation study.     

On families of beta- and generalized gamma-generated distributions and associated inference

A general family of univariate distributions generated by beta random variables, proposed by Jones, has been discussed recently in the literature. This family of distributions possesses great flexibility while fitting symmetric as well as skewed models with varying tail weights. In a similar vein, we define here a family of univariate distributions generated by Stacy’s generalized gamma variables. For these two families of univariate distributions, we discuss maximum entropy characterizations under suitable constraints. Based on these characterizations, an expected ratio of quantile densities is proposed for the discrimination of members of these two broad families of distributions. Several special cases of these results are then highlighted. An alternative to the usual method of moments is also proposed for the estimation of the parameters, and the form of these estimators is particularly amenable to these two families of distributions.     

Bayesiam simultaneous estimation for several multinomial distributions

A Bayesian method is proposed for estimating the cell probabilities of several multinomial distributions. Parameters of different distributions are taken to be a priori exchangeable. The prior specification is based upon mixtures of a hierarchical distribution, referred to as the multivariate “Dirichlet-Dirichlet” distribution. The analysis is facilitated by a multinomial approximation relating to the multinomial-Dirichlet distribution. The posterior estimates depend upon measures of entropy for the various distributions and shrink the individual observed proportions towards values obtained by pooling the data across the distributions. As well as incorporating prior information they are particularly useful when some of the cell frequencies are zero. We use them to investigate a numerical classification of males of various vocations, according to cause of death.     

On the Harris extended family of distributions

A new method for generating new classes of distributions based on the probability-generating function is presented in Aly and Benkherouf [A new family of distributions based on probability generating functions. Sankhya B. 2011;73:70–80]. In particular, they focused their interest to the so-called Harris extended family of distributions. In this paper, we provide several general results regarding the Harris extended models such as the general behaviour of the failure rate function. We also derive a very useful representation for the Harris extended density function as an absolutely convergent power series of the survival function of the baseline distribution. Additionally, some stochastic order relations are established and limiting distributions of sample extremes are also considered for this model. These general results are illustrated in several special Harris extended models. Finally, we discuss estimation of the model parameters by the method of maximum likelihood and provide an application to real data for illustrative purposes.     

On the discrete analogues of continuous distributions

In this paper, a new method is proposed for generating discrete distributions. A special class of the distributions, namely, the T-geometric family contains the discrete analogues of continuous distributions. Some general properties of the T-geometric family of distributions are obtained. A member of the T-geometric family, namely, the exponentiated-exponential–geometric distribution is defined and studied. Various properties of the exponentiated-exponential–geometric distribution such as the unimodality, the moments and the probability generating function are discussed. The method of maximum likelihood estimation is proposed for estimating the model parameters. Three real data sets are used to illustrate the applications of the exponentiated-exponential–geometric distribution.     

A use of mixtures of two normal distributions in a classification problem

An assumption made in the classification problem is that the distribution of the data being classified has the same parameters as the data used to obtain the discriminant functions. A method based on mixtures of two normal distributions is proposed as method of checking this assumption and modifying the discriminant functions accordingly. As a first step, the case considered in this paper, is that of a shift in the mean of one or two univariate normal distributions with all other parameters remaining fixed and known. Calculations based on the asymptotic the proposed method works well even for small shifts.     

On the discrete analogues of continuous distributions

In this paper, a new method is proposed for generating discrete distributions. A special class of the distributions, namely, the T-geometric family contains the discrete analogues of continuous distributions. Some general properties of the T-geometric family of distributions are obtained. A member of the T-geometric family, namely, the exponentiated-exponential–geometric distribution is defined and studied. Various properties of the exponentiated-exponential–geometric distribution such as the unimodality, the moments and the probability generating function are discussed. The method of maximum likelihood estimation is proposed for estimating the model parameters. Three real data sets are used to illustrate the applications of the exponentiated-exponential–geometric distribution.     

Reliability computation under dynamic stress–strength modeling with cumulative stress and strength degradation

Reliability function is defined under suitable assumptions for dynamic stress–strength scenarios where strength degrades and stress accumulates over time. Methods for numerical evaluation of reliability are suggested under deterministic strength degradation and cumulative damage due to shocks arriving according to a point process, in particular a Poisson process, using simulation method and inversion theorem. These methods are specifically useful in the scenarios where damage distributions do not possess closure property under convolution. The method is also extended for non-identical, dependent damage distributions as well as for random strength degradation. Results from inversion method is compared with known approximate methods and also verified by simulation. As it turns out, the simulation method seems to have an edge in terms of computational burden and has much wider domain of applicability.     

Approximate MLEs of the parameters of location-scale models under type II censoring

Log-location-scale distributions are widely used parametric models that have fundamental importance in both parametric and semiparametric frameworks. The likelihood equations based on a Type II censored sample from location-scale distributions do not provide explicit solutions for the para-meters. Statistical software is widely available and is based on iterative methods (such as, Newton Raphson Algorithm, EM algorithm etc.), which require starting values near the global maximum. There are also many situations that the specialized software does not handle. This paper provides a method for determining explicit estimators for the location and scale parameters by approximating the likelihood function, where the method does not require any starting values. The performance of the proposed approximate method for the Weibull distribution and Log-Logistic distributions is compared with those based on iterative methods through the use of simulation studies for a wide range of sample size and Type II censoring schemes. Here we also examine the probability coverages of the pivotal quantities based on asymptotic normality. In addition, two examples are given.     

A comparative study of the K-means algorithm and the normal mixture model for clustering: Bivariate homoscedastic case

The K-means algorithm and the normal mixture model method are two common clustering methods. The K-means algorithm is a popular heuristic approach which gives reasonable clustering results if the component clusters are ball-shaped. Currently, there are no analytical results for this algorithm if the component distributions deviate from the ball-shape. This paper analytically studies how the K-means algorithm changes its classification rule as the normal component distributions become more elongated under the homoscedastic assumption and compares this rule with that of the Bayes rule from the mixture model method. We show that the classification rules of both methods are linear, but the slopes of the two classification lines change in the opposite direction as the component distributions become more elongated. The classification performance of the K-means algorithm is then compared to that of the mixture model method via simulation. The comparison, which is limited to two clusters, shows that the K-means algorithm provides poor classification performances consistently as the component distributions become more elongated while the mixture model method can potentially, but not necessarily, take advantage of this change and provide a much better classification performance.     

Implied distributions in multiple change point problems

A method for efficiently calculating exact marginal, conditional and joint distributions for change points defined by general finite state Hidden Markov Models is proposed. The distributions are not subject to any approximation or sampling error once parameters of the model have been estimated. It is shown that, in contrast to sampling methods, very little computation is needed. The method provides probabilities associated with change points within an interval, as well as at specific points.     

Calculation of the Prokhorov distance by optimal quantization and maximum flow

Calculating exact values of the Prokhorov metric for the set of probability distributions on a metric space is a challenging problem. In this paper probability distributions are approximated by finite-support distributions through optimal or quasi-optimal quantization, in such a way that exact calculation of the Prokhorov distance between a distribution and a quantizer can be performed. The exact value of the Prokhorov distance between two quantizers is obtained by solving an optimization problem through the Simplex method. This last value is used to approximate the Prokhorov distance between the two initial distributions, and the accuracy of the approximation is measured. We illustrate the method on various univariate and bivariate probability distributions. Approximation of bivariate standard normal distributions by quasi-optimal quantizers is also considered.     

Bayesian Estimation of the Parameters of Bivariate Exponential Distributions

In this article, we develop an empirical Bayesian approach for the Bayesian estimation of parameters in four bivariate exponential (BVE) distributions. We have opted for gamma distribution as a prior for the parameters of the model in which the hyper parameters have been estimated based on the method of moments and maximum likelihood estimates (MLEs). A simulation study was conducted to compute empirical Bayesian estimates of the parameters and their standard errors. We use moment estimators or MLEs to estimate the hyper parameters of the prior distributions. Furthermore, we compare the posterior mode of parameters obtained by different prior distributions and the Bayesian estimates based on gamma priors are very close to the true values as compared to improper priors. We use MCMC method to obtain the posterior mean and compared the same using the improper priors and the classical estimates, MLEs.     

A Bayesian approach on the two-piece scale mixtures of normal homoscedastic nonlinear regression models

In this application note paper, we propose and examine the performance of a Bayesian approach for a homoscedastic nonlinear regression (NLR) model assuming errors with two-piece scale mixtures of normal (TP-SMN) distributions. The TP-SMN is a large family of distributions, covering both symmetrical/ asymmetrical distributions as well as light/heavy tailed distributions, and provides an alternative to another well-known family of distributions, called scale mixtures of skew-normal distributions. The proposed family and Bayesian approach provides considerable flexibility and advantages for NLR modelling in different practical settings. We examine the performance of the approach using simulated and real data.KEYWORDS: Gibbs sampling, MCMC method, nonlinear regression model, scale mixtures of normal family, two-piece distributions