Eliciting expert judgements about a set of proportions

Eliciting expert knowledge about several uncertain quantities is a complex task when those quantities exhibit associations. A well-known example of such a problem is eliciting knowledge about a set of uncertain proportions which must sum to 1. The usual approach is to assume that the expert's knowledge can be adequately represented by a Dirichlet distribution, since this is by far the simplest multivariate distribution that is appropriate for such a set of proportions. It is also the most convenient, particularly when the expert's prior knowledge is to be combined with a multinomial sample since then the Dirichlet is the conjugate prior family. Several methods have been described in the literature for eliciting beliefs in the form of a Dirichlet distribution, which typically involve eliciting from the expert enough judgements to identify uniquely the Dirichlet hyperparameters. We describe here a new method which employs the device of over-fitting, i.e. eliciting more than the minimal number of judgements, in order to (a) produce a more carefully considered Dirichlet distribution and (b) ensure that the Dirichlet distribution is indeed a reasonable fit to the expert's knowledge. The method has been implemented in a software extension of the Sheffield elicitation framework (SHELF) to facilitate the multivariate elicitation process.     

Discrete Linnik Weibull distribution

Abstract

We introduce a new family of distributions using truncated discrete Linnik distribution. This family is a rich family of distributions which includes many important families of distributions such as Marshall–Olkin family of distributions, family of distributions generated through truncated negative binomial distribution, family of distributions generated through truncated discrete Mittag–Leffler distribution etc. Some properties of the new family of distributions are derived. A particular case of the family, a five parameter generalization of Weibull distribution, namely discrete Linnik Weibull distribution is given special attention. This distribution is a generalization of many distributions, such as extended exponentiated Weibull, exponentiated Weibull, Weibull truncated negative binomial, generalized exponential truncated negative binomial, Marshall-Olkin extended Weibull, Marshall–Olkin generalized exponential, exponential truncated negative binomial, Marshall–Olkin exponential and generalized exponential. The shape properties, moments, median, distribution of order statistics, stochastic ordering and stress–strength properties of the new generalized Weibull distribution are derived. The unknown parameters of the distribution are estimated using maximum likelihood method. The discrete Linnik Weibull distribution is fitted to a survival time data set and it is shown that the distribution is more appropriate than other competitive models.     

Some properties of the dirichlet-multinomial distribution and its use in prior elicitation

Two results on the unimodality of the Dirichlet-multinomial distribution are proved, and a further result is alos proved on the identifiability of mixtures of multinomial distributions. These properties are used in developing a method for eliciting a Dirchlet prior distribution. The elicitation method is based on the mode, and region around the mode, of the Dirichlet-multinomial predictive distribution.     

Household budget-share distributions and welfare implications: an application of multivariate distributional statistics

In this paper the consequences of considering the household ‘food share’ distribution as a welfare measure, in isolation from the joint distribution of itemized budget shares, is examined through the unconditional and conditional distribution of ‘food share’ both parametrically and nonparametrically. The parametric framework uses Dirichlet and Beta distributions, while the nonparametric framework uses kernel smoothing methods. The analysis, in a three commodity setup (‘food’, ‘durables’, ‘others’), based on household level rural data for West Bengal, India, for the year 2009–2010 shows significant underrepresentation of households by the conventional unconditional ‘food share’ distribution in the higher range of food budget shares that correspond to the lower end of the income profile. This may have serious consequences for welfare measurement.     

Bayesian nonparametric survival analysis using mixture of Burr XII distributions

ABSTRACT

Recently, the Bayesian nonparametric approaches in survival studies attract much more attentions. Because of multimodality in survival data, the mixture models are very common. We introduce a Bayesian nonparametric mixture model with Burr distribution (Burr type XII) as the kernel. Since the Burr distribution shares good properties of common distributions on survival analysis, it has more flexibility than other distributions. By applying this model to simulated and real failure time datasets, we show the preference of this model and compare it with Dirichlet process mixture models with different kernels. The Markov chain Monte Carlo (MCMC) simulation methods to calculate the posterior distribution are used.     

THE WRAPPED t FAMILY OF CIRCULAR DISTRIBUTIONS

This paper considers the three‐parameter family of symmetric unimodal distributions obtained by wrapping the location‐scale extension of Student's t distribution onto the unit circle. The family contains the wrapped normal and wrapped Cauchy distributions as special cases, and can be used to closely approximate the von Mises distribution. In general, the density of the family can only be represented in terms of an infinite summation, but its trigonometric moments are relatively simple expressions involving modified Bessel functions. Point estimation of the parameters is considered, and likelihood‐based methods are used to fit the family of distributions in an illustrative analysis of cross‐bed measurements. The use of the family as a means of approximating the von Mises distribution is investigated in detail, and new efficient algorithms are proposed for the generation of approximate pseudo‐random von Mises variates.     

Nonparametric Mixtures Based on Skew-normal Distributions: An Application to Density Estimation

This article addresses the density estimation problem using nonparametric Bayesian approach. It is considered hierarchical mixture models where the uncertainty about the mixing measure is modeled using the Dirichlet process. The main goal is to build more flexible models for density estimation. Extensions of the normal mixture model via Dirichlet process previously introduced in the literature are twofold. First, Dirichlet mixtures of skew-normal distributions are considered, say, in the first stage of the hierarchical model, the normal distribution is replaced by the skew-normal one. We also assume a skew-normal distribution as the center measure in the Dirichlet mixture of normal distributions. Some important results related to Bayesian inference in the location-scale skew-normal family are introduced. In particular, we obtain the stochastic representations for the full conditional distributions of the location and skewness parameters. The algorithm introduced by MacEachern and Müller in 1998 MacEachern, S.N., Müller, P. (1998). Estimating mixture of Dirichlet Process models. J. Computat. Graph. Statist. 7(2):223–238.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar] is used to sample from the posterior distributions. The models are compared considering simulated data sets. Finally, the well-known Old Faithful Geyser data set is analyzed using the proposed models and the Dirichlet mixture of normal distributions. The model based on Dirichlet mixture of skew-normal distributions captured the data bimodality and skewness shown in the empirical distribution.     

Bayesian design and analysis of composite endpoints in clinical trials with multiple dependent binary outcomes

The author considers studies with multiple dependent primary endpoints. Testing hypotheses with multiple primary endpoints may require unmanageably large populations. Composite endpoints consisting of several binary events may be used to reduce a trial to a manageable size. The primary difficulties with composite endpoints are that different endpoints may have different clinical importance and that higher‐frequency variables may overwhelm effects of smaller, but equally important, primary outcomes. To compensate for these inconsistencies, we weight each type of event, and the total number of weighted events is counted. To reflect the mutual dependency of primary endpoints and to make the weighting method effective in small clinical trials, we use the Bayesian approach. We assume a multinomial distribution of multiple endpoints with Dirichlet priors and apply the Bayesian test of noninferiority to the calculation of weighting parameters. We use composite endpoints to test hypotheses of superiority in single‐arm and two‐arm clinical trials. The composite endpoints have a beta distribution. We illustrate this technique with an example. The results provide a statistical procedure for creating composite endpoints. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

On representing second-order uncertainty in multi-state systems via moments of mixtures of Dirichlet distributions†

Second-order probabilities have been proposed as representations of the uncertainty in the parameters of probabilistic models such as Bayesian belief networks. We investigate conditions under which second-order probabilities can be represented in terms of their marginal moments. We show that certain combinations of marginal means and variances do not correspond to any valid second-order joint distribution. By fitting a Dirichlet mixture to marginal mean and variance information, we derive sufficient conditions for a valid second-order joint distribution to exist.     

On the Dirichlet distributions on symmetric matrices

In this paper, we introduce a generalization of the Dirichlet distribution on symmetric matrices which represents the multivariate version of the Connor and Mosimann generalized real Dirichlet distribution. We establish some properties concerning this generalized distribution. We also extend to the matrix Dirichlet distribution a remarkable characterization established in the real case by Darroch and Ratcliff.     

Zero-modified power series distribution and its Hurdle distribution version

This paper presents a new family of distributions for count data, the so called zero-modified power series (ZMPS), which is an extension of the power series (PS) distribution family, whose support starts at zero. This extension consists in modifying the probability of observing zero of each PS distribution, enabling the new zero-modified distribution to appropriately accommodate data which have any amount of zero observations (for instance, zero-inflated or zero-deflated data). The Hurdle distribution version of the ZMPS distribution is presented. PS distributions included in the proposed ZMPS family are the Poisson, Generalized Poisson, Geometric, Binomial, Negative Binomial and Generalized Negative Binomial distributions. The paper also describes the properties and particularities of the new distribution family for count data. The distribution parameters are estimated via maximum likelihood method and the use of the new family is illustrated in three real data sets. We emphasize that the new distribution family can accommodate sets of count data without any previous knowledge on the characteristic of zero-inflation or zero-deflation present in the data.     

The roles of prior catchability assumptions and sample features on bayes estimates of the number of classes in a population

The object of this paper is to explain the role played by the catchability and sampling in the Bayesian estimation of k, the unknown number of classes in a multinomial population. It is shown that the posterior distribution of k increases as the capture probabilities of the classes become more unequal, and that the posterior distribution of k increases with the number of classes observed in the sample and decreases with the sample size. Moreover, it is shown that the posterior mean of k is consistent.     

A Bivariate Geometric Distribution with Applications to Reliability

This article studies a bivariate geometric distribution (BGD) as a plausible reliability model. Maximum likelihood and Bayes estimators of parameters and various reliability characteristics are obtained. Approximations to the mean, variance, and Bayes risk of these estimators have been derived using Taylor's expansion. A Monte-Carlo simulation study has been performed to compare these estimators. At the end, the theory is illustrated with a real data set example of accidents.     

The Exact Likelihood Ratio Test for Equality of Two Normal Populations

Testing the equality of two independent normal populations is a perfect case of the two-sample problem, yet it is not treated in the main text of any textbook or handbook. In this article, we derive the exact distribution of the likelihood ratio test and implement this test with an R function. This article has supplementary materials online.     

Posterior Distributions for the Gini Coefficient Using Grouped Data

When available data comprise a number of sampled households in each of a number of income classes, the likelihood function is obtained from a multinomial distribution with the income class population proportions as the unknown parameters. Two methods for going from this likelihood function to a posterior distribution on the Gini coefficient are investigated. In the first method, two alternative assumptions about the underlying income distribution are considered, namely a lognormal distribution and the Singh–Maddala (1976) income distribution. In these cases the likelihood function is reparameterized and the Gini coefficient is a nonlinear function of the income distribution parameters. The Metropolis algorithm is used to find the corresponding posterior distributions of the Gini coefficient from a sample of Bangkok households. The second method does not require an assumption about the nature of the income distribution, but uses (a) triangular prior distributions, and (b) beta prior distributions, on the location of mean income within each income class. By sampling from these distributions, and the Dirichlet posterior distribution of the income class proportions, alternative posterior distributions of the Gini coefficient are calculated.     

On integrals of dirichlet distributions and their applications

The purpose of this paper is to extend results obtained by Sahai and Anderson (1973) and Hurlburt and Spiegel (1976) computationally as well as theoretically. We apply results given in Yassaee (1974), (1976) and apply our computer program to show that one can evaluate the probability of type one error for conditional tests in linear model, in general. We show that one can compute the probability integral of Dirichlet distribution by the use of our program which computes the probability integral of inverted Dirichlet distribution.     

Generating random variates from a bicompositional Dirichlet distribution

A composition is a vector of positive components summing to a constant. The sample space of a composition is the simplex, and the sample space of two compositions, a bicomposition, is a Cartesian product of two simplices. We present a way of generating random variates from a bicompositional Dirichlet distribution defined on the Cartesian product of two simplices using the rejection method. We derive a general solution for finding a dominating density function and a rejection constant and also compare this solution to using a uniform dominating density function. Finally, some examples of generated bicompositional random variates, with varying number of components, are presented.     

Modified slashed-Rayleigh distribution

A new family of slash distributions, the modified slashed-Rayleigh distribution, is proposed and studied. This family is an extension of the ordinary Rayleigh distribution, being more flexible in terms of distributional kurtosis. It arises as a quotient of two independent random variables, one being a Rayleigh distribution in the numerator and the other a power of the exponential distribution in denominator. We present properties of the proposed family. In addition, we carry out estimation of the model parameters by moment and maximum likelihood methods. Finally, we conduct a small-scale simulation study to evaluate the performance of the maximum likelihood estimators and apply the results to a real data set, revealing its good performance.     

A General Version of the Short-Tailed Symmetric Distribution

Motivated by Tiku and Vaughan, we consider a type of distribution which is a generalization of the short-tailed symmetric distribution. In particular, we consider the asymptotic behavior of Mills’ ratios for the distribution family. Meanwhile, we obtain the asymptotic distribution of the partial maximum of an independent and identically distributed sequence from the distribution.     

Generalized exponential distributions

The three-parameter gamma and three-parameter Weibull distributions are commonly used for analysing any lifetime data or skewed data. Both distributions have several desirable properties, and nice physical interpretations. Because of the scale and shape parameters, both have quite a bit of flexibility for analysing different types of lifetime data. They have increasing as well as decreasing hazard rate depending on the shape parameter. Unfortunately both distributions also have certain drawbacks. This paper considers a three-parameter distribution which is a particular case of the exponentiated Weibull distribution originally proposed by Mudholkar, Srivastava & Freimer (1995) when the location parameter is not present. The study examines different properties of this model and observes that this family has some interesting features which are quite similar to those of the gamma family and the Weibull family, and certain distinct properties also. It appears this model can be used as an alternative to the gamma model or the Weibull model in many situations. One dataset is provided where the three-parameter generalized exponential distribution fits better than the three-parameter Weibull distribution or the three-parameter gamma distribution.