Generalized Gumbel distribution

A generalization of the Gumbel distribution is presented to deal with general situations in modeling univariate data with broad range of skewness in the density function. This generalization is derived by considering a logarithmic transformation of an odd Weibull random variable. As a result, the generalized Gumbel distribution is not only useful for testing goodness-of-fit of Gumbel and reverse-Gumbel distributions as submodels, but it is also convenient for modeling and fitting a wide variety of data sets that are not possible to be modeled by well-known distributions. Skewness and kurtosis shapes of the generalized Gumbel distribution are illustrated by constructing the Galton’s skewness and Moor’s kurtosis plane. Parameters are estimated by using maximum likelihood method in two different ways due to the fact that the reverse transformation of the proposed distribution does not change its density function. In order to illustrate the flexibility of this generalization, wave and surge height data set is analyzed, and the fitness is compared with Gumbel and generalized extreme value distributions.     

Information Loss in Grade Point Conversion

Some information gets lost when numerical scores evaluating performances are converted into letter grades. We propose to measure this information loss through the proportion of variance lost due to grouping. We study various properties of this measure, including its invariance in location and scale equivariant families. The information loss typically decreases with an increase in the number of levels of letter grades. However, it is not appropriate to have too many levels. The optimum number of levels may be determined, either by visual inspection when the information loss becomes marginal/stable, or by minimizing the sum of the information loss and a penalty term, the latter being taken as linear in the number of levels. We also address the problem of determining the groups, or equivalently, the boundaries so that the information loss is minimized, given a fixed number of groups. Finding these optimal boundaries is a computationally intensive exercise even for moderate size data, unless the number of groups is very small. We recommend an alternative way by fitting an appropriate probability distribution. When the probabilistic nature of the data is known, the boundary points turn out to be the solutions to a system of equations; however these solutions may not necessarily have any closed form. We derive the exact or approximate solutions of these equations when the composite scores follow a probability distribution belonging to the Uniform, Triangular, and Gaussian family.     

Gamma mixture of generalized error distribution

Abstract

A new symmetric heavy-tailed distribution, namely gamma mixture of generalized error distribution is defined by scaling generalized error distribution with gamma distribution, its probability density function, k-moment, skewness and kurtosis are derived. After tedious calculation, we also give the Fisher information matrix, moment estimators and maximum likelihood estimators for the parameters of gamma mixture of generalized error distribution. In order to evaluate the effectiveness of the point estimators and the stability of Fisher information matrix, extensive simulation experiments are carried out in three groups of parameters. Additionally, the new distribution is applied to Apple Inc. stock (AAPL) data and compared with normal distribution, F-S skewed standardized t distribution and generalized error distribution. It is found that the new distribution has better fitting effect on the data under the Akaike information criterion (AIC). To a certain extent, our results enrich the probability distribution theory and develop the scale mixture distribution, which will provide help and reference for financial data analysis.     

A new weighted distribution as an extension of the generalized half-normal distribution with applications

In this study, a new extension of generalized half-normal (GHN) distribution is introduced. Since this new distribution can be viewed as weighted version of GHN distribution, it is called as weighted generalized half-normal (WGHN) distribution. It is shown that WGHN distribution can be observed as a single constrained and hidden truncation model. Therefore, the new distribution is more flexible than the GHN distribution. Some statistical properties of the WGHN distribution are studied, i.e. moments, cumulative distribution function, hazard rate function are derived. Furthermore, maximum likelihood estimation of the parameters is considered. Some real-life data sets taken from the literature are modelled using the WGHN distribution. It is seen that for these data sets the WGHN distribution provides better fitting than the GHN and slashed generalized half-normal (SGHN) distributions.     

Consistency of Bernstein polynomial posteriors

A Bernstein prior is a probability measure on the space of all the distribution functions on [0, 1]. Under very general assumptions, it selects absolutely continuous distribution functions, whose densities are mixtures of known beta densities. The Bernstein prior is of interest in Bayesian nonparametric inference with continuous data. We study the consistency of the posterior from a Bernstein prior. We first show that, under mild assumptions, the posterior is weakly consistent for any distribution function P ₀ on [0, 1] with continuous and bounded Lebesgue density. With slightly stronger assumptions on the prior, the posterior is also Hellinger consistent. This implies that the predictive density from a Bernstein prior, which is a Bayesian density estimate, converges in the Hellinger sense to the true density (assuming that it is continuous and bounded). We also study a sieve maximum likelihood version of the density estimator and show that it is also Hellinger consistent under weak assumptions. When the order of the Bernstein polynomial, i.e. the number of components in the beta distribution mixture, is truncated, we show that under mild restrictions the posterior concentrates on the set of pseudotrue densities. Finally, we study the behaviour of the predictive density numerically and we also study a hybrid Bayes–maximum likelihood density estimator.     

Bivariate inverse Weibull distribution

ABSTRACT

Recently it is observed that the inverse Weibull (IW) distribution can be used quite effectively to analyse lifetime data in one dimension. The main aim of this paper is to define a bivariate inverse Weibull (BIW) distribution so that the marginals have IW distributions. It is observed that the joint probability density function and the joint cumulative distribution function can be expressed in compact forms. Several properties of this distribution such as marginals, conditional distributions and product moments have been discussed. We obtained the maximum likelihood estimates for the unknown parameters of this distribution and their approximate variance– covariance matrix. We perform some simulations to see the performances of the maximum likelihood estimators. One data set has been re-analysed and it is observed that the bivariate IW distribution provides a better fit than the bivariate exponential distribution.     

Fitting with Matrix-Exponential Distributions

Abstract

It is well known that general phase-type distributions are considerably overparameterized, that is, their representations often require many more parameters than is necessary to define the distributions. In addition, phase-type distributions, even those defined by a small number of parameters, may have representations of high order. These two problems have serious implications when using phase-type distributions to fit data. To address this issue we consider fitting data with the wider class of matrix-exponential distributions. Representations for matrix-exponential distributions do not need to have a simple probabilistic interpretation, and it is this relaxation which ensures that the problems of overparameterization and high order do not present themselves. However, when using matrix-exponential distributions to fit data, a problem arises because it is unknown, in general, when their representations actually correspond to a distribution. In this paper we develop a characterization for matrix-exponential distributions and use it in a method to fit data using maximum likelihood estimation. The fitting algorithm uses convex semi-infinite programming combined with a nonlinear search.     

The Marshall-Olkin logistic-exponential distribution

We introduce a three-parameter extension of the exponential distribution which contains as sub-models the exponential, logistic-exponential and Marshall-Olkin exponential distributions. The new model is very flexible and its associated density function can be decreasing or unimodal. Further, it can produce all of the four major shapes of the hazard rate, that is, increasing, decreasing, bathtub and upside-down bathtub. Given that closed-form expressions are available for the survival and hazard rate functions, the new distribution is quite tractable. It can be used to analyze various types of observations including censored data. Computable representations of the quantile function, ordinary and incomplete moments, generating function and probability density function of order statistics are obtained. The maximum likelihood method is utilized to estimate the model parameters. A simulation study is carried out to assess the performance of the maximum likelihood estimators. Two actual data sets are used to illustrate the applicability of the proposed model.     

The beta generalized exponential distribution

A new distribution called the beta generalized exponential distribution is proposed. It includes the beta exponential and generalized exponential (GE) distributions as special cases. We provide a comprehensive mathematical treatment of this distribution. The density function can be expressed as a mixture of generalized exponential densities. This is important to obtain some mathematical properties of the new distribution in terms of the corresponding properties of the GE distribution. We derive the moment generating function (mgf) and the moments, thus generalizing some results in the literature. Expressions for the density, mgf and moments of the order statistics are also obtained. We discuss estimation of the parameters by maximum likelihood and obtain the information matrix that is easily numerically determined. We observe in one application to a real skewed data set that this model is quite flexible and can be used effectively in analyzing positive data in place of the beta exponential and GE distributions.     

The McDonald arcsine distribution: a new model to proportional data

We propose a new three-parameter continuous model called the McDonald arcsine distribution, which is a very competitive model to the beta, beta type I and Kumaraswamy distributions for modelling rates and proportions. We provide a mathematical treatment of the new distribution including explicit expressions for the density function, moments, generating and quantile functions, mean deviations, two probability measures based on the Bonferroni and Lorenz curves, Shannon entropy, Rényi entropy and cumulative residual entropy. Maximum likelihood is used to estimate the model parameters and the expected information matrix is determined. An application of the proposed model to real data shows that it can give consistently a better fit than other important statistical models.     

The Touchard distribution

We present a novel model, which is a two-parameter extension of the Poisson distribution. Its normalizing constant is related to the Touchard polynomials, hence the name of this model. It is a flexible distribution that can account for both under- or overdispersion and concentration of zeros that are frequently found in non-Poisson count data. In contrast to some other generalizations, the Hessian matrix for maximum likelihood estimation of the Touchard parameters has a simple form. We exemplify with three data sets, showing that our suggested model is a competitive candidate for fitting non-Poisson counts.     

Pareto Poisson–Lindley distribution with applications

A new lifetime distribution is introduced based on compounding Pareto and Poisson–Lindley distributions. Several statistical properties of the distribution are established, including behavior of the probability density function and the failure rate function, heavy- and long-right tailedness, moments, the Laplace transform, quantiles, order statistics, moments of residual lifetime, conditional moments, conditional moment generating function, stress–strength parameter, Rényi entropy and Song's measure. We get maximum-likelihood estimators of the distribution parameters and investigate the asymptotic distribution of the estimators via Fisher's information matrix. Applications of the distribution using three real data sets are presented and it is shown that the distribution fits better than other related distributions in practical uses.     

Segmentation of circular data

Circular data – data whose values lie in the interval [0,2π) – are important in a number of application areas. In some, there is a suspicion that a sequence of circular readings may contain two or more segments following different models. An analysis may then seek to decide whether there are multiple segments, and if so, to estimate the changepoints separating them. This paper presents an optimal method for segmenting sequences of data following the von Mises distribution. It is shown by example that the method is also successful in data following a distribution with much heavier tails.     

The Poisson Birnbaum–Saunders model with long-term survivors

In this paper, we propose a cure rate survival model by assuming that the number of competing causes of the event of interest follows the Poisson distribution and the time to event has the Birnbaum–Saunders (BS) distribution. We define the Poisson BS distribution and provide two useful representations for its density function which facilitate to obtain some mathematical properties. Two closed-form expressions for the moments of the new distribution are given. We estimate the parameters of the model with cure rate using maximum likelihood. For different parameter settings, sample sizes and censoring percentages, several simulations are performed. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and present some ways to perform a global influence study. We analyse a real data set from the medical area.     

Gamma-Maxwell distribution

A new two-parameter distribution, the gamma-Maxwell distribution, isproposed and studied. We generate the new distribution using the gamma-G generator of distributions. The proposal distribution can be seen as an extension of the Maxwell distribution with more flexibility in terms of the distribution asymmetry and kurtosis. We study some probability properties, discuss maximum-likelihood estimation and present a real data application indicating that the new distribution can improve the ordinary Maxwell distribution in fitting real data.     

Bayes estimation and prediction of the two-parameter gamma distribution

In this article, the Bayes estimates of two-parameter gamma distribution are considered. It is well known that the Bayes estimators of the two-parameter gamma distribution do not have compact form. In this paper, it is assumed that the scale parameter has a gamma prior and the shape parameter has any log-concave prior, and they are independently distributed. Under the above priors, we use Gibbs sampling technique to generate samples from the posterior density function. Based on the generated samples, we can compute the Bayes estimates of the unknown parameters and can also construct HPD credible intervals. We also compute the approximate Bayes estimates using Lindley's approximation under the assumption of gamma priors of the shape parameter. Monte Carlo simulations are performed to compare the performances of the Bayes estimators with the classical estimators. One data analysis is performed for illustrative purposes. We further discuss the Bayesian prediction of future observation based on the observed sample and it is seen that the Gibbs sampling technique can be used quite effectively for estimating the posterior predictive density and also for constructing predictive intervals of the order statistics from the future sample.     

A generalization of the beta–binomial distribution

Summary.  The paper describes a distribution generated by the Gaussian hypergeometric function that may be seen as a generalization of the beta–binomial distribution. It is expressed as a generalized beta mixture of a binomial distribution. This new mixing distribution allows the existence of a mode and an antimode, which is very useful for fitting some data sets. Two examples illustrate the greater versatility of the new distribution compared with the beta–binomial distribution.     

Linear models with an inverse gaussian poisson error distribution

The inverse Gaussian-Poisson (two-parameter Sichel) distribution is useful in fitting overdispersed count data. We consider linear models on the mean of a response variable, where the response is in the form of counts exhibiting extra-Poisson variation, and assume an IGP error distribution. We show how maximum likelihood estimation may be carried out using iterative Newton-Raphson IRLS fitting, where GLIM is used for the IRLS part of the maximization. Approximate likelihood ratio tests are given.     

A new skew-bimodal distribution with applications

The modeling and analysis of experiments is an important aspect of statistical work in a wide variety of scientific and technological fields. We introduce and study the odd log-logistic skew-normal model, which can be interpreted as a generalization of the skew-normal distribution. The new distribution can be used effectively in the analysis of experiments data since it accommodates unimodal, bimodal, symmetric, bimodal and right-skewed, and bimodal and left-skewed density function depending on the parameter values. We illustrate the importance of the new model by means of three real data sets in analysis of experiments.     

The beta exponentiated Weibull distribution

The Weibull distribution is one of the most important distributions in reliability. For the first time, we introduce the beta exponentiated Weibull distribution which extends recent models by Lee et al. [Beta-Weibull distribution: some properties and applications to censored data, J. Mod. Appl. Statist. Meth. 6 (2007), pp. 173–186] and Barreto-Souza et al. [The beta generalized exponential distribution, J. Statist. Comput. Simul. 80 (2010), pp. 159–172]. The new distribution is an important competitive model to the Weibull, exponentiated exponential, exponentiated Weibull, beta exponential and beta Weibull distributions since it contains all these models as special cases. We demonstrate that the density of the new distribution can be expressed as a linear combination of Weibull densities. We provide the moments and two closed-form expressions for the moment-generating function. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The density of the order statistics can also be expressed as a linear combination of Weibull densities. We obtain the moments of the order statistics. The expected information matrix is derived. We define a log-beta exponentiated Weibull regression model to analyse censored data. The estimation of the parameters is approached by the method of maximum likelihood. The usefulness of the new distribution to analyse positive data is illustrated in two real data sets.