The mixture of left–right truncated normal distributions

This paper introduces the mixture of left–right truncated normal distributions, from the spreads between bid and ask prices, as a statistical model for handle non-normality of asset price returns. It has been proved that there is only one maximum for the likelihood function of the new model.     

Testing homogeneity in semiparametric mixture case–control models

Parametric and semiparametric mixture models have been widely used in applications from many areas, and it is often of interest to test the homogeneity in these models. However, hypothesis testing is non standard due to the fact that several regularity conditions do not hold under the null hypothesis. We consider a semiparametric mixture case–control model, in the sense that the density ratio of two distributions is assumed to be of an exponential form, while the baseline density is unspecified. This model was first considered by Qin and Liang (2011 Qin, J., Liang, K.Y. (2011). Hypothesis testing in a mixture case–control model. Biometrics 67(1):182–198.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar], biometrics), and they proposed a modified score statistic for testing homogeneity. In this article, we consider alternative testing procedures based on supremum statistics, which could improve power against certain types of alternatives. We demonstrate the connection and comparison among the proposed and existing approaches. In addition, we provide a unified theoretical justification of the supremum test and other existing test statistics from an empirical likelihood perspective. The finite-sample performance of the supremum test statistics was evaluated in simulation studies.     

Value-at-risk forecasting based on Gaussian mixture ARMA–GARCH model

In this paper, we develop a new forecasting algorithm for value-at-risk (VaR) based on ARMA–GARCH (autoregressive moving average–generalized autoregressive conditional heteroskedastic) models whose innovations follow a Gaussian mixture distribution. For the parameter estimation, we employ the conditional least squares and quasi-maximum-likelihood estimator (QMLE) for ARMA and GARCH parameters, respectively. In particular, Gaussian mixture parameters are estimated based on the residuals obtained from the QMLE of GARCH parameters. Our algorithm provides a handy methodology, spending much less time in calculation than the existing resampling and bias-correction method developed in Hartz et al. [Accurate value-at-risk forecasting based on the normal-GARCH model, Comput. Stat. Data Anal. 50 (2006), pp. 3032–3052]. Through a simulation study and a real-data analysis, it is shown that our method provides an accurate VaR prediction.     

Improved auxiliary mixture sampling for hierarchical models of non-Gaussian data

The article considers Bayesian analysis of hierarchical models for count, binomial and multinomial data using efficient MCMC sampling procedures. To this end, an improved method of auxiliary mixture sampling is proposed. In contrast to previously proposed samplers the method uses a bounded number of latent variables per observation, independent of the intensity of the underlying Poisson process in the case of count data, or of the number of experiments in the case of binomial and multinomial data. The bounded number of latent variables results in a more general error distribution, which is a negative log-Gamma distribution with arbitrary integer shape parameter. The required approximations of these distributions by Gaussian mixtures have been computed. Overall, the improvement leads to a substantial increase in efficiency of auxiliary mixture sampling for highly structured models. The method is illustrated for finite mixtures of generalized linear models and an epidemiological case study.     

Bayesian covariance matrix estimation using a mixture of decomposable graphical models

We present a Bayesian approach to estimating a covariance matrix by using a prior that is a mixture over all decomposable graphs, with the probability of each graph size specified by the user and graphs of equal size assigned equal probability. Most previous approaches assume that all graphs are equally probable. We show empirically that the prior that assigns equal probability over graph sizes outperforms the prior that assigns equal probability over all graphs in more efficiently estimating the covariance matrix. The prior requires knowing the number of decomposable graphs for each graph size and we give a simulation method for estimating these counts. We also present a Markov chain Monte Carlo method for estimating the posterior distribution of the covariance matrix that is much more efficient than current methods. Both the prior and the simulation method to evaluate the prior apply generally to any decomposable graphical model.     

Higher order C(α) tests with applications to mixture models

A class of second-order C(α) tests is proposed for testing composite hypotheses when no optimal test among the (first-order) C(α) tests, as defined by Neyman (Probability and Statistics, The Harald Cramer Volume, Almqvist and Wiksell, Uppsala, Sweden, 1959, p. 213), exists. The form of an optimal test of the new subclass is also presented. The procedure is seen to be easily extendable to C(α) tests of third or higher orders. Two examples of univariate normal mixtures are considered where the procedure is successfully applied. For the corresponding multivariate models, union–intersection tests of Roy (Ann. Math. Statist. 24 (1953) 220) are derived by combining the above optimal tests of the univariate problems.     

A mixture of beta–Dirichlet processes prior for Bayesian analysis of event history data

In this paper, we propose a mixture of beta–Dirichlet processes as a nonparametric prior for the cumulative intensity functions of a Markov process. This family of priors is a natural extension of a mixture of Dirichlet processes or a mixture of beta processes which are devised to compromise advantages of parametric and nonparametric approaches. They give most of their prior mass to the small neighborhood of a specific parametric model. We show that a mixture of beta–Dirichlet processes prior is conjugate with Markov processes. Formulas for computing the posterior distribution are derived. Finally, results of analyzing credit history data are given.     

Order statistics and their concomitants from multivariate normal mean–variance mixture distributions with application to Swiss Markets Data

In this article, by considering a multivariate normal mean–variance mixture distribution, we derive the exact joint distribution of linear combinations of order statistics and their concomitants. From this general result, we then deduce the exact marginal and conditional distributions of order statistics and their concomitants arising from this distribution. We finally illustrate the usefulness of these results by using a Swiss markets dataset.     

Robust mixture modeling using multivariate skew <Emphasis Type="Italic">t</Emphasis> distributions

This paper presents a robust mixture modeling framework using the multivariate skew t distributions, an extension of the multivariate Student’s t family with additional shape parameters to regulate skewness. The proposed model results in a very complicated likelihood. Two variants of Monte Carlo EM algorithms are developed to carry out maximum likelihood estimation of mixture parameters. In addition, we offer a general information-based method for obtaining the asymptotic covariance matrix of maximum likelihood estimates. Some practical issues including the selection of starting values as well as the stopping criterion are also discussed. The proposed methodology is applied to a subset of the Australian Institute of Sport data for illustration.     

Inference for serological surveys investigating past exposures to infections resulting in long-lasting immunity – an approach using finite mixture models with concomitant information

This paper is concerned with developing a latent class mixture modelling technique which efficiently exploits data from serological surveys aiming to investigate past exposures to infections resulting in long-term or life-lasting immunity. Mixture components featured by antibody assays’ distribution are associated with the serological groups in the population, whilst the probability mixture that an individual belongs to the positive serological group is regarded as an age-dependent prevalence. The latter embeds a mechanistic model which explains the infection process, accounting for heterogeneities, contact patterns in the population and incorporating elements of study design. A Bayesian framework for statistical inference using Markov chain Monte Carlo estimation methods naturally accommodates missing responses in the data and allows straightforward assessement of uncertainties in nonlinear models. The applicability of the method is illustrated by investigating past exposure to varicella zoster virus infection in pre-school children, using data from a large scale UK cohort study which included a cross-sectional serological survey based on oral fluid samples.     

Flexible mixture modeling via the multivariate t distribution with?the?Box-Cox transformation: an?alternative to?the?skew-t distribution

Cluster analysis is the automated search for groups of homogeneous observations in a data set. A popular modeling approach for clustering is based on finite normal mixture models, which assume that each cluster is modeled as a multivariate normal distribution. However, the normality assumption that each component is symmetric is often unrealistic. Furthermore, normal mixture models are not robust against outliers; they often require extra components for modeling outliers and/or give a poor representation of the data. To address these issues, we propose a new class of distributions, multivariate t distributions with the Box-Cox transformation, for mixture modeling. This class of distributions generalizes the normal distribution with the more heavy-tailed t distribution, and introduces skewness via the Box-Cox transformation. As a result, this provides a unified framework to simultaneously handle outlier identification and data transformation, two interrelated issues. We describe an Expectation-Maximization algorithm for parameter estimation along with transformation selection. We demonstrate the proposed methodology with three real data sets and simulation studies. Compared with a wealth of approaches including the skew-t mixture model, the proposed t mixture model with the Box-Cox transformation performs favorably in terms of accuracy in the assignment of observations, robustness against model misspecification, and selection of the number of components.