parameter_estimates: A column matrix of β, γ, and θ parameter values obtained from a COMMA analysis function. Parameter estimates should be supplied in the following order: 1) β
(intercept, slope), 2) γ (intercept and slope from the M = 1 mechanism, intercept and slope from the M = 2 mechanism), and 3) θ
(intercept, slope, coefficient for x, slope coefficient for m, slope coefficient for c, and, optionally, slope coefficient for xm if using).
sigma_estimate: A numeric value specifying the estimated standard deviation. This value is only required if outcome_distribution
is "Normal". Default is 1. For non-Normal outcome distributions, the value should be NULL.
outcome_distribution: A character string specifying the distribution of the outcome variable. Options are "Bernoulli", "Normal", or "Poisson".
interaction_indicator: A logical value indicating if an interaction between x and m should be used to generate the outcome variable, y.
x_matrix: A numeric matrix of predictors in the true mediator and outcome mechanisms. x_matrix should not contain an intercept and no values should be NA.
z_matrix: A numeric matrix of covariates in the observation mechanism. z_matrix should not contain an intercept and no values should be NA.
c_matrix: A numeric matrix of covariates in the true mediator and outcome mechanisms. c_matrix should not contain an intercept and no values should be NA.
Returns
COMMA_boot_sample returns a list with the bootstrap sample data: - obs_mediator: A vector of observed mediator values.
true_mediator: A vector of true mediator values.
outcome: A vector of outcome values.
x_matrix: A matrix of predictor values in the true mediator mechanism. Identical to that supplied by the user.
z_matrix: A matrix of predictor values in the observed mediator mechanism. Identical to that supplied by the user.
c_matrix: A matrix of covariates. Identical to that supplied by the user.