Causal Mediation Analysis for Multiple Outcome/Treatment/Mediator Combinations
'mediations' can be used to process a set of outcome/treatment/mediator
combinations through the mediate function to produce a series
of causal mediation analysis results.
mediations(datasets, treatment, mediators, outcome, covariates = NULL,
families = c("gaussian", "gaussian"), tau.m = 0.5, tau.y = 0.5,
LowerY = NULL, UpperY = NULL, interaction = FALSE,
conf.level = 0.95, sims = 500, boot = FALSE, weights = NULL, ...)datasets |
a named list of data frames. Each data frame has a separate treatment variable. The names of each data frame must begin with the exact name of the treatment variable that is contained in that dataset (see example below). |
treatment |
a vector of character strings indicating the names of the treatment variables, with length equal to the length of 'datasets'. Each treatment variable must be included in the data frame listed in the same position of list 'datasets' and its name must match the first part of the corresponding data frame. |
mediators |
a vector of character strings indicating the names of the mediators contained within each data frame. All of the mediators will be used with each treatment variable and hence must be included in each data frame of 'datasets'. |
outcome |
a vector of character strings indicating the names of the outcome variables contained within each data frame. All of the outcomes will be used with each treatment variable and must be in each data frame. |
covariates |
a character string representing the set of pre-treatment covariate names (as they appear in the data frame) to be included in each model. The value must take the form of standard model formula, with each additive component separated by "+", etc. (see example below). All covariates must be in each data frame. Default is 'NULL'. |
families |
a vector of length two specifying the types of the mediator and outcome models. Currently only supports "gaussian" (for linear regression), "binomial" (for binary probit), "oprobit" (for ordered probit) and "quantile" (for quantile regression, see 'tau'). For the outcome the tobit model ("tobit") is also available in addition to the mediator model options. |
tau.m |
a numeric value specifying the quantile to be used for a
quantile regression for the mediator model. Only relevant if the first
element of 'families' is "quantile". See |
tau.y |
a numeric value specifying the quantile to be used for a
quantile regression for the outcome model. Only relevant if the second
element of 'families' is "quantile". See |
LowerY |
a numeric value indicating the lower bound for the tobit
outcome model. See |
UpperY |
a numeric value indicating the upper bound for the tobit
outcome model. See |
interaction |
a logical value indicating whether the treatment and
mediator variables should be interacted. This will apply to applications of
|
conf.level |
confidence level used in each application of the
|
sims |
an integer indicating the desired number of simulations for inference. This will apply to all applications of 'mediate' to all the treatment/mediator/outcome combinations. |
boot |
a logical value, indicating whether or not nonparametric
bootstrap should be used in each |
weights |
a single valued vector of a character string indicating a weight variable to be used in all model fitting. |
... |
other arguments passed to |
This function processes multiple treatment/mediators/outcome
variable combinations to produce a collected set of output ready for
analysis or graphing. In principle, this is a function designed to
facilitate running causal mediation analyses on multiple models that share
the same basic specification (i.e. the types of parametric models and the
set of pre-treatment covariates) except the treatment, mediator and outcome
variables can differ across specifications. The function works by looping
over a set of data frames that are pre-loaded into the workspace. Each one
of these data frames has a specific treatment variable that is used for
analysis with that data frame. Then the code runs causal mediation analysis
via mediate on every combination of the treatment, mediator,
and outcomes specified in these arguments. This allows the users to explore
whether different mediators transmit the effect of the treatment variable
on a variety of outcome variables. A single set of pre-treatment control
variables can be specified in 'covariates', which will be used throughout.
The 'mediations' function can be used with either multiple mediators and a single outcome, a single mediator and multiple outcomes, or multiple mediators and outcomes. For example, with three different treatments, user will create three different data frames, each containing a treatment variable. In addition, if there are also four different mediators, each of these will be contained in each data frame, along with the outcome variable. The function will estimate all of the combinations of treatment variables and mediators instead of separate lines of code being written for each one.
The default value of 'covariates' is 'NULL' and no covariate will be included in either mediator or outcome models without a custom value. It should be noted that users typically should have pre-treatment covariates to make the sequential ignorability assumption more plausible.
There are several limitations to the code. First, it works only with a subset of the model types that will be accommodated if 'mediate' is used individually (see the 'families' argument above for details). Second, one cannot specify separate sets of covariates for different treatment/mediator/outcome combinations. Users should use 'mediate' separately for individual models if more flexibility is required in their specific applications.
An object of class "mediations" (or "mediations.order" if the outcome
model is ordered probit), a list of "mediate" ("mediate.order") objects
produced by applications of mediate for the specified
treatment/mediator/outcome combinations. The elements are named based on
the names of the outcome, treatment, and mediator variables, each separated
by a "." (see example below).
Dustin Tingley, Harvard University, dtingley@gov.harvard.edu; Teppei Yamamoto, Massachusetts Institute of Technology, teppei@mit.edu.
mediate, summary.mediations,
plot.mediations, rq, tobit.
## Not run:
# Hypothetical example
datasets <- list(T1 = T1, T2 = T2)
# List of data frames corresponding to the two different treatment variables
#"T1vsCont" and "T2vsCont".
# Each data set has its respective treatment variable.
mediators <- c("M1", "M2")
# Vector of mediator names, all included in each data frame.
outcome <- c("Ycont1","Ycont2")
# Vector of outcome variable names, again all included in each data frame.
treatment <- c("T1vsCont", "T2vsCont")
# Vector of treatment variables names; must begin with identical strings with dataset
# names in 'datasets'.
covariates <- c("X1 + X2")
# Set of covariates (in each data set), entered using the standard model formula format.
x <- mediations(datasets, treatment, mediators, outcome, covariates,
families=c("gaussian","gaussian"), interaction=FALSE,
conf.level=.90, sims=50)
# Runs 'mediate' iteratively for each variable combinations, with 'lm' on both mediator
# and outcome model.
summary(x) # tabular summary of results for all model combinations
plot(x) # graphical summary of results for all model combinations at once
plot(x$Ycont1.T1vsCont.M1)
# Individual 'mediate' outputs are stored as list elements and can be
# accessed using the usual "$" operator.
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.