MAZE: Mediation Analysis for ZEro-inflated mediators

In MAZE: Mediation Analysis for Zero-Inflated Mediators

Mediation Analysis for ZEro-inflated mediators

Description

A novel mediation modeling approach to address zero-inflated mediators containing true zeros and false zeros.

Usage

MAZE(
  data,
  distM = c("zilonm", "zinbm", "zipm"),
  K = 1,
  selection = "AIC",
  X,
  M,
  Y,
  Z = NULL,
  XMint = c(TRUE, FALSE),
  x1,
  x2,
  zval = NULL,
  mval = 0,
  B = 20,
  seed = 1,
  ncore = 1
)

Arguments

data	a data frame containing variables: an independent variable X, a mediator M, an outcome Y, and confounder variables Z (if any). See example dataset: data(zinb10)
distM	a vector with choices of the distribution of mediator to try with. One or more of 'zilonm', 'zinbm', and 'zipm' for zero-inflated log-normal, negative binomial, and Poisson mediators respectively. Default is c('zilonm', 'zinbm', 'zipm') where all three distributions are fitted and the final mediation model is selected by model selection criterion selection
K	a vector with choices of the number of component K in the zero-inflated mixture mediators to try with. Default is K=1 for zero-inflated (non-mixture) mediators
selection	model selection criterion when more than one model (combination of different values in distM and K) is fitted. Either 'AIC' or 'BIC'. Default is 'AIC'
X	name of the independent variable. Can be continuous or discrete
M	name of the mediator variable. Non-negative values
Y	name of the outcome variable. Continuous values
Z	name(s) of confounder variables (if any)
XMint	a logical vector of length 2 indicating whether to include the two exposure-mediator interaction terms between (i) X and 1_{(M>0)} and (ii) X and M. Default is c(TRUE, FALSE), which only includes the first
x1	the first value of independent variable of interest
x2	the second value of independent variable of interest
zval	a vector of value(s) of confounders to be conditional on when estimating effects
mval	the fixed value of mediator to be conditional on when estimating CDE
B	the upper bound value B to be used in the probability mechanism of observing false zeros
seed	an optional seed number to control randomness
ncore	number of cores available for parallel computing

Details

For an independent variable X, a zero-inflated mediator M and a continuous outcome variable Y, the following regression equation is used to model the association between Y and (X,M):

Y_{xm1_{(m>0)}}=β_0+β_1m+β_2 1_{(m>0)}+β_3x+β_4x1_{(m>0)}+β_5xm+ε

Users can choose to include either one, both, or none of the two exposure-mediator interaction terms between (i) X and 1_{(M>0)} and (ii) X and M using the argument XMint.

For mediators, zero-inflated log-normal, zero-inflated negative binomial, and zero-inflated Poisson distributions are considered and can be specified through the argument distM.

The indirect and direct effects (NIE1, NIE2, NIE, NDE, and CDE) are estimated for X changing from x1 to x2. When confounders are present, the conditional effects are estimated given the fixed value zval.

Value

a list containing:

• results_effects: a data frame for the results of estimated effects (NIE1, NIE2, NIE, NDE, and CDE)

• results_parameters: a data frame for the results of model parameters

• selected_model_name: a string for the distribution of M and number of components K selected in the final mediation model

• BIC: a numeric value for the BIC of the final mediation model

• AIC: a numeric value for the AIC of the final mediation model

• models: a list with all fitted models

• analysis2_out: a list with output from analysis2() function (used for internal check)

Author(s)

Meilin Jiang [e-mail] and Zhigang Li [e-mail]

Examples

data(zinb10)

maze_out <- MAZE(data = zinb10,
                 distM = c('zilonm', 'zinbm', 'zipm'),  K = 1,
                 selection = 'AIC',
                 X = 'X', M = 'Mobs', Y = 'Y', Z = NULL,
                 XMint = c(TRUE, FALSE),
                 x1 = 0, x2 = 1, zval = NULL, mval = 0,
                 B = 20, seed = 1)
## results of selected mediation model
maze_out$results_effects # indirect and direct effects
maze_out$selected_model_name # selected distribution of the mediator and number of components K
maze_out$results_parameters # model parameters
maze_out$BIC; maze_out$AIC # BIC and AIC of the selected mediation model

MAZE documentation built on Feb. 16, 2023, 5:21 p.m.