Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/mediate.designs.R

Estimating Average Causal Mediation Effects under the Single Experiment Design

1 2 | ```
mediate.sed(outcome, mediator, treat, data, SI = FALSE, sims = 1000,
conf.level = 0.95, boot = FALSE)
``` |

`outcome` |
name of the outcome variable in 'data'. The variable must be binary (factor or numeric 0/1) if 'SI' is FALSE. |

`mediator` |
name of the mediator in 'data'. The variable must be binary (factor or numeric 0/1) if 'SI' is FALSE and discrete if TRUE. |

`treat` |
name of the treatment variable in 'data'. Must be binary (factor or numeric 0/1). |

`data` |
a data frame containing all the above variables. |

`SI` |
whether the sequential ignorability assumption is made. |

`sims` |
number of bootstrap simulations. Only relevant when 'SI' is TRUE. |

`conf.level` |
level of the returned two-sided confidence intervals. Only relevant when 'SI' is TRUE. |

`boot` |
a logical value. if 'FALSE' a large sample Delta method approximation is used for confidence intervals; if 'TRUE' nonparametric bootstrap will be used. Default is 'FALSE'. Only relevant if 'SI' is TRUE. |

'mediate.sed' estimates average causal mediation effects for the single experiment design. The two options are to use either the sequential ignorability (SI) assumption in which nonparametric estimates of the average causal mediation effect are produced, or, to relax the SI assumption and to calculate the nonparametric bounds on the average causal mediation effect.

This function calculates average causal mediation effects (ACME) for the single experiment design, where the treatment is randomized and the mediator/outcome variables are measured. The user specifies whether they want non-parametric point estimates based on the sequential ignorability (SI) assumption, or nonparametric bounds without the SI assumption.

`mediate.sed`

returns an object of class
"`mediate.design`

", a list that contains the components listed below.

The `summary`

function can be used to obtain a table of the results.

`d0, d1` |
point estimates or lower/upper bounds for causal mediation effects under the control and treatment conditions, respectively. |

`d0.ci, d1.ci` |
confidence intervals for average causal mediation effects for the nonparametric estimates. The confidence level is set at the value specified in 'conf.level'. The value exists only when 'SI' is TRUE. |

`boot` |
logical, the 'boot' argument used. |

`conf.level` |
the confidence level used. |

`sims` |
number of bootstrap simulations used for confidence interval calculation. |

`nobs` |
number of observations used. |

`design` |
indicates the design. Equals either "SED.NP.SI" or "SED.NP.NOSI". |

Dustin Tingley, Harvard University, dtingley@gov.harvard.edu; Teppei Yamamoto, Massachusetts Institute of Technology, teppei@mit.edu.

Tingley, D., Yamamoto, T., Hirose, K., Imai, K. and Keele, L. (2014). "mediation: R package for Causal Mediation Analysis", Journal of Statistical Software, Vol. 59, No. 5, pp. 1-38.

Imai, K., Tingley, D. and Yamamoto, T. (2012) Experimental Designs for Identifying Causal Mechanisms. Journal of the Royal Statistical Society, Series A (Statistics in Society)"

Imai, K., Keele, L., Tingley, D. and Yamamoto, T. (2011). Unpacking the Black Box of Causality: Learning about Causal Mechanisms from Experimental and Observational Studies, American Political Science Review, Vol. 105, No. 4 (November), pp. 765-789.

Imai, K., Keele, L. and Yamamoto, T. (2010) Identification, Inference, and Sensitivity Analysis for Causal Mediation Effects, Statistical Science, Vol. 25, No. 1 (February), pp. 51-71.

Imai, K., Keele, L., Tingley, D. and Yamamoto, T. (2009) "Causal Mediation Analysis Using R" in Advances in Social Science Research Using R, ed. H. D. Vinod New York: Springer.

`mediate`

, `summary.mediate.design`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ```
# Example 1: Bounds without SI assumption
data(boundsdata)
data.SED <- subset(boundsdata, manip == 0)
bound1 <- mediate.sed("out", "med", "ttt", data.SED, SI=FALSE)
summary(bound1)
# Example 2: Nonparametric estimate of ACME under SI assumption
# Example with JOBS II Field Experiment
data(jobs)
foo.1 <- mediate.sed("depress2", "job_disc", "treat", jobs, SI=TRUE)
summary(foo.1)
foo.2 <- mediate.sed("depress2", "job_disc", "treat", jobs, SI=TRUE, boot=TRUE)
summary(foo.2)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.