This Article Full Text (PDF) References Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Alert me to new issues of the journal Add to Saved Citations Download to citation manager Add to My Marked Citations Citing Articles Citing Articles via Google Scholar Citing Articles via Scopus Google Scholar Articles by Liu, J. Articles by Dey, D. K Search for Related Content Social Bookmarking                         What's this?

Skew random effects in multilevel binomial models

an alternative to nonparametric approach

Biometrics Division, The Cancer Institute of New Jersey, USA and Department of Biostatistics, School of Public Health, University of Medicine and Dentistry of New Jersey, USA

Dipak K Dey

Department of Statistics, University of Connecticut, Storrs, USA

Compared to modelling observable data, it is more difficult to choose a suitable distribution to describe latent variables since no prior knowledge or observable information can be used and only normal or nonparametric distributions are mainly applied to random effects for generalized linear mixed models (GLMMs) in the literature. To enhance the modelling toolkit, this article investigates a class of parametric skew elliptical random effects in multilevel binomial regression models using a Bayesian approach; the class includes skew normal, skew Students’ t-distributions and others. Skewness mechanism is considered through multiplying skewness parameter by standardized folded elliptical random variables, and the posterior sampling is realized by working on a binary skewness indicator (BSI) instead of continuous for parameter identifiability. Simulation study shows that the original continuous skewness parameter and the posterior mean of BSI may have dichotomous signs to describe the directional (right/left) skewness; thus we address the importance of assuming specific random effects distribution and interpreting the skewness carefully. The methodology is exemplified through reanalyzing a teratogenic activity study of two niacin analogs published in the biological literature, and sampling-based model comparison shows that the parametric skew normal random effects model works largely better than nonparametric Dirichlet process mixture models for this data set.

Key Words: Binomial • binary skewness indicator (BSI) • generalized linear mixed model (GLMM) • multilevel • random effects • sign switching • skewness

Statistical Modelling, Vol. 8, No. 3, 221-241 (2008)
DOI: 10.1177/1471082X0800800301


CiteULike    Complore    Connotea    Del.icio.us    Digg    Reddit    Technorati    Twitter    What's this?