An Approximate Bayesian Marginal Likelihood Approach for Estimating Finite Mixtures

Estimation of finite mixture models when the mixing distribution support is unknown is an important problem. This article gives a new approach based on a marginal likelihood for the unknown support. Motivated by a Bayesian Dirichlet prior model, a computationally efficient stochastic approximation version of the marginal likelihood is proposed and large-sample theory is presented. By restricting the support to a finite grid, a simulated annealing method is employed to maximize the marginal likelihood and estimate the support. Real and simulated data examples show that this novel stochastic approximation and simulated annealing procedure compares favorably with existing methods.     

Mixtures of Mixtures Based on Multiresolution Analysis Theory

In this article, the approaches for exploiting mixtures of mixtures are expanded by using the Multiresolution family of probability density functions (MR pdf). The flexibility and the properties of local analysis of the MR pdf facilitate the location of subpopulations into a given population. In order to do this, two algorithms are provided.

The MR model is more flexible in adapting to the different subpopulations than the traditional mixtures. In addition, the problems of identification of mixtures distributions and the label-switching do not appear in the MR pdf context.     

Estimators for the Finite Mixture of Rayleigh Model Based on Progressively Censored Data

In this article, based on progressively Type-II censored samples from a heterogeneous population that can be represented by a finite mixture of two-component Rayleigh lifetime model, the problem of estimating the parameters and some lifetime parameters (reliability and hazard functions) are considered. Both Bayesian and maximum likelihood estimators are of interest. A class of natural conjugate prior densities is considered in the Bayesian setting. The Bayes estimators are obtained using both the symmetric (squared error) loss function, and the asymmetric (LINEX and General Entropy) loss functions. It has been seen that the estimators obtained can be easily evaluated for this type of censoring by using suitable numerical methods. Finally, the performance of the estimates have been compared on the basis of their simulated maximum square error via a Monte Carlo simulation study.     

Optimal designs for experiments with mixtures: a survey

This is a survey article on known results about analytic solutions and numerical solutions of optimal designs for various regression models for experiments with mixtures. The regression models include polynomial models, models containing homogeneous functions, models containing inverse terms and ratios, log contrast models, models with quantitative variables, and mod els containing the amount of mixture, Optimality criteria considered include D-, A-, E-,φ_p- and I_λ-Optimalities. Uniform design and uniform optimal design for mixture components, and efficiencies of the {q,2} simplex-controid design are briefly discussed.     

High-order extensions of the Double Chain Markov Model

The Double Chain Markov Model is a fully Markovian model for the representation of time-series in random environments. In this article, we show that it can handle transitions of high-order between both a set of observations and a set of hidden states. In order to reduce the number of parameters, each transition matrix can be replaced by a Mixture Transition Distribution model. We provide a complete derivation of the algorithms needed to compute the model. Three applications, the analysis of a sequence of DNA, the song of the wood pewee, and the behavior of young monkeys show that this model is of great interest for the representation of data that can be decomposed into a finite set of patterns.     

On the type-B integral equation and the distribution of wilks statist for testing independence of several groups of variables

In this paper, WILKS'type-B integral equation is solved in the general form of a series of beta functions and a series of weighted gamma functions as proposed by WALD and BROOKNER 1941. The coefficients in both representations can be obtained by explicit recurrence relartions, therefore the results solve many distributional problems and have the fewest computational difficulties of any representation that has surfaced to date. The radius of convergence of the second series representation is given, whereas the convergence property of the first series representation is given, whereas the convergence property of the first series representation was studied by WALD and Brookner. The exact null distributions of WILKS' statistic A for testing the independence of several groups of variables and of V = -log A are given. The coefficients in all the series representration can be computed recursilvely and hence can be obtained easily with the help of modern computatinal facilities     

The Parshvanath yantram yields a constant-sum partially balanced incomplete block design

Yantrams have been used to generate mixture designs in the interior of a simplex. In this note, we show a connection between Parshvanath yantram and a particular partially balanced incomplete block design. This block design is rather special and somewhat unexpected due to the feature that sum of the treatment symbols in any block is constant.     

The zero-inflated Poisson model and the decayed, missing and filled teeth index in dental epidemiology

For frequency counts, the situation of extra zeros often arises in biomedical applications. This is demonstrated with count data from a dental epidemiological study in Belo Horizonte (the Belo Horizonte caries prevention study) which evaluated various programmes for reducing caries. Extra zeros, however, violate the variance–mean relationship of the Poisson error structure. This extra-Poisson variation can easily be explained by a special mixture model, the zero-inflated Poisson (ZIP) model. On the basis of the ZIP model, a graphical device is presented which not only summarizes the mixing distribution but also provides visual information about the overall mean. This device can be exploited to evaluate and compare various groups. Ways are discussed to include covariates and to develop an extension of the conventional Poisson regression. Finally, a method to evaluate intervention effects on the basis of the ZIP regression model is described and applied to the data of the Belo Horizonte caries prevention study.     

On the three-dimensional negative binomial distribution

We consider a particular generalization of the negative binomial distribution to the multivariate case obtained through a specification of the probability generating function as the negative power of a certain polynomial. The probability function itself has previously been derived for the two-dimensional case only, and inference in the multivariate negative binomial distribution has been restricted to the use of composite likelihood based on one- or two-dimensional marginals. In this article, we derive the three-dimensional probability function as a sum with all terms positive and study the range of possible parameter values. We illustrate the use of the three-dimensional distribution for modeling three correlated SAR images.     

A unified approach to stochastic comparisons of multivariate mixture models

In this paper, we consider a unified approach to stochastic comparisons of random vectors corresponding to two general multivariate mixture models. These stochastic comparisons are made with respect to multivariate hazard rate, reversed hazard rate and likelihood ratio orders. As an application, results are presented for stochastic comparisons of generalized multivariate frailty models.