summary.medsens {mediation} | R Documentation |
Function to perform sensitivity analysis on mediation effect for violations of sequential ignorability assumption. The procedure estimates the mediation effect allowing for a correlation between the error terms of the outcome model and the mediator model. The extent of this correlation is expressed in terms of the parameter rho. Sensitivity analysis is possible with continuous mediator and continuous outcome and binary outcome and continuous mediator. Future versions of the function will also permit sensitivity analysis for a continuous outcome and binary mediator. Output from the function can be passed through summary or plot functions which display estimated mediation effects for given values of rho.
## S3 method for class 'medsens': summary(object, ...) ## S3 method for class 'summary.medsens': print(x, ...)
x |
Output from medsens function. |
object |
Output from medsens function. |
... |
Additional arguments to be passed. |
Luke Keele, Ohio State University, keele.4@osu.edu , Dustin Tingley, Princeton University, dtingley@princeton.edu, Teppei Yamamoto, Princeton University, tyamamot@princeton.edu, Kosuke Imai, Princeton University, kimai@princeton.edu
Imai, Kosuke, Luke Keele and Dustin Tingley (2009) A General Approach to Causal Mediation Analysis. Imai, Kosuke, Luke Keele and Teppei Yamamoto (2009) Identification, Inference, and Sensitivity Analysis for Causal Mediation Effects. Imai, Kosuke, Luke Keele, Dustin Tingley, and Teppei Yamamoto. (2009) "Causal Mediation Analysis Using R" in Advances in Social Science Research Using R, ed. H. D. Vinod New York: Springer. Forthcoming.
See also medsens
## Not run: ######################################### #Continuous mediator and continuous outcome ######################################### #Fit parametric models model.m <- lm(job_seek ~ treat + depress1, data=jobs) model.y <- lm(depress2 ~ treat + job_seek + depress1, data=jobs) #Pass model objects through mediate function med.cont <- mediate(model.m, model.y, treat="treat", mediator="job_seek", sims=1000) #Pass mediate output through medsens function sens.cont <- medsens(med.cont, sims=1000, rho.by=.1) #Use summary function to display values of rho where 95 summary(sens.cont) ######################################### #binary outcome and continuous mediator ######################################### model.m <- lm(job_seek ~ treat + depress1 + econ_hard + sex + age + occp + marital + nonwhite + educ + income, data=jobs) model.y <- glm(work1 ~ treat + job_seek + depress1 + econ_hard + sex + age + occp + marital + nonwhite + educ + income, family=binomial(link="probit"), data=jobs) med.bin <- mediate(model.m, model.y, treat="treat", mediator="job_seek", sims=1000) sens.bin <- medsens(med.bin, sims=1000, rho.by=.1) summary(sens.bin) ## End(Not run)