A mixed-effect model for positive responses augmented by zeros

Stat Med. 2015 May 10;34(10):1761-78. doi: 10.1002/sim.6450. Epub 2015 Feb 11.

Abstract

In this research article, we propose a class of models for positive and zero responses by means of a zero-augmented mixed regression model. Under this class, we are particularly interested in studying positive responses whose distribution accommodates skewness. At the same time, responses can be zero, and therefore, we justify the use of a zero-augmented mixture model. We model the mean of the positive response in a logarithmic scale and the mixture probability in a logit scale, both as a function of fixed and random effects. Moreover, the random effects link the two random components through their joint distribution and incorporate within-subject correlation because of the repeated measurements and between-subject heterogeneity. A Markov chain Monte Carlo algorithm is tailored to obtain Bayesian posterior distributions of the unknown quantities of interest, and Bayesian case-deletion influence diagnostics based on the q-divergence measure is performed. We apply the proposed method to a dataset from a 24 hour dietary recall study conducted in the city of São Paulo and present a simulation study to evaluate the performance of the proposed methods.

Keywords: Bayesian inference; gamma distribution; log-normal distribution; mixed models; random effects; usual intake; zero-augmented distributions.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms
  • Bayes Theorem
  • Brazil
  • Computer Simulation
  • Diet / statistics & numerical data*
  • Humans
  • Likelihood Functions
  • Linear Models
  • Markov Chains
  • Mental Recall
  • Models, Statistical*
  • Monte Carlo Method
  • Poisson Distribution