Description Usage Arguments Details Value Author(s) Examples
View source: R/BANOVA.multi.mediation.R
BANOVA.multi.mediation
is a function for analysis of multiple possibly correlated mediators.
These mediators are assumed to have no causal influence on each other.
Both singlelevel and multilevel models can be analyzed.
1  BANOVA.multi.mediation(sol_1, sol_2, xvar, mediators, individual = FALSE)

sol_1 
an object of class "BANOVA" returned by BANOVA.run function with a fitted model for an outcome variable regressed on a causal variable, a mediator, and, possibly, moderators and control variables. The outcome variable can follow Normal, T, Poisson, Bernoulli, Binomial, Truncated Normal and ordered Multinomial distributions. 
sol_2 
an object of class "BANOVA" returned by BANOVA.run function, which contains an outcome of the analysis for multiple Multivariate Normal mediators regressed on a casual variable and other possible moderators and control variables. 
xvar 
a character string that specifies the name of the causal variable used in both models. 
mediators 
a vector with character strings, which specifies the names of the mediator variables used in the models. 
individual 
logical indicator of whether to output effects for individual units in the
analysis (TRUE or FALSE). This analysis requires a multilevel 
The function extends BANOVA.mediation
to the case with multiple possibly
correlated mediators. For details about mediation analysis performed in BANOVA see
the help page for the BANOVA.mediation.
BANOVA.multi.mediation
estimates and tests specific indirect effects of the causal
variable conveyed through each mediator. Furthermore, the total indirect effect of the causal
variables are computed as a sum of the specific indirect effects.
The function prints multiple tables with mediated effects. Tables with direct effects of the causal variable and mediators on the outcome variable, as well as direct effects of the causal variable on the mediators include a posterior mean and 95% credible intervals of the effects. Next, the function displays on the console tables with specific indirect effects and effect sizes of the mediators, followed by the TIE of the causal variable. These tables include the mean, 95% credible intervals, and twosided Bayesian pvalues.
Returns an object of class "BANOVA.multi.mediation"
. The returned object is a list
containing:

table or tables with the direct effect. 

is returned if 

a list with tables of the effects of the mediator on the outcome 

a list with tables of the effect of the causal variable on the mediator 
indir_effects 
tables of the indirect effect 
individual_indirect 
is returned if 
effect_sizes 
a list with effect sizes on individual mediators 
total_indir_effects 
table or tables with the total indirect effect of the causal variable 
xvar 
the name of the causal variable 
mediators 
the names of the mediating variables 
individual 
the value of the argument individual (TRUE or FALSE) 
Anna Kopyakova
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23  # Use the colorad data set
data(colorad)
# Add a second mediator to the data set
colorad$blur_squared < (colorad$blur)^2
# Prepare mediators to be analyzed in the Multivariate Normal model
mediators < cbind(colorad$blur, colorad$blur_squared)
colnames(mediators) < c("blur", "blur_squared")
colorad$mediators < mediators
# Build and analyze the model for the outcome variable
model < BANOVA.model('Binomial')
banova_binom_model < BANOVA.build(model)
res_1 < BANOVA.run(y ~ typic, ~ color + blur + blur_squared, fit = banova_binom_model,
data = colorad, id = 'id', num_trials = as.integer(16),
iter = 2000, thin = 1, chains = 2)
# Build and analyze the model for the mediators
model < BANOVA.model('multiNormal')
banova_multi_norm_model < BANOVA.build(model)
res_2 < BANOVA.run(mediators ~ typic, ~ color, fit = banova_multi_norm_model,
data = colorad, id = 'id', iter = 2000, thin = 1, chains = 2)
# Calculate (moderated) effects of "typic" mediated by "blur" and "blur_squared"
results < BANOVA.multi.mediation(res_1, res_2, xvar='typic', mediators=c("blur", "blur_squared"))

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.