An imputation approach for fitting two-part mixed effects models for longitudinal semi-continuous data

Two-part mixed effects models are often used for analyzing longitudinal data with many zeros. Typically, these models are formulated with binary and continuous components separately with random effects that are correlated between the two components. Researchers have developed maximum-likelihood and Bayesian approaches for fitting these models that often require using particular software packages or very specialized software. We propose an imputation approach that will allow practitioners to separately use standard linear and generalized linear mixed models to estimate the fixed effects for two-part mixed effects models with complex random effects structures. An approximation to the conditional distribution of positive measurements given an individual's pattern of non-zero measurements is proposed that can be easily estimated and then imputed from. We show that for a wide range of parameter values, the imputation approach results in nearly unbiased estimation and can be implemented with standard software. We illustrate the proposed imputation approach for the analysis of longitudinal clinical trial data with many zeros.