Application of marginalized zero-inflated models when mediators have excess zeroes

Mediation analysis has become increasingly popular over the last decade as researchers are interested in assessing mechanistic pathways for intervention. Although available methods have increased, there are still limited options for mediation analysis with zero-inflated count variables where the distribution of response has a "cluster" of data at the zero value (i.e. distribution of number of cigarettes smoked per day, where nonsmokers cluster at zero cigarettes). The currently available methods do not obtain unbiased population average effects of mediation effects. In this paper, we propose an extension of the counterfactual approach to mediation with direct and indirect effects to scenarios where the mediator is a count variable with excess zeroes by utilizing the Marginalized Zero-Inflated Poisson Model (MZIP) for the mediator model. We derive direct and indirect effects for continuous, binary, and count outcomes, as well as adapt to allow mediator-exposure interactions. Our proposed work allows straightforward calculation of direct and indirect effects for the overall population mean values of the mediator, for scenarios in which researchers are interested in generalizing direct and indirect effects to the population. We apply this novel methodology to an application observing how alcohol consumption may explain sex differences in cholesterol and assess model performance via a simulation study comparing the proposed MZIP mediator framework to existing methods for marginal mediator effects.