Description Usage Arguments Details Value References Examples

View source: R/medweightcont.R

Causal mediation analysis (evaluation of natural direct and indirect effects) of a continuous treatment based on weighting by the inverse of generalized propensity scores as suggested in Hsu, Huber, Lee, and Lettry (2020).

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
medweightcont(
y,
d,
m,
x,
d0,
d1,
ATET = FALSE,
trim = 0.1,
lognorm = FALSE,
bw = NULL,
boot = 1999,
cluster = NULL
)
``` |

`y` |
Dependent variable, must not contain missings. |

`d` |
Continuous treatment, must not contain missings. |

`m` |
Mediator(s), may be a scalar or a vector, must not contain missings. |

`x` |
Pre-treatment confounders of the treatment, mediator, and/or outcome, must not contain missings. |

`d0` |
Value of |

`d1` |
Value of |

`ATET` |
If FALSE, the average treatment effect (ATE) and the corresponding direct and indirect effects are estimated. If TRUE, the average treatment effect on the treated (ATET) and the corresponding direct and indirect effects are estimated. Default is |

`trim` |
Trimming rule for discarding observations with too large weights in the estimation of any mean potential outcome. That is, observations with a weight> |

`lognorm` |
If FALSE, a linear model with normally distributed errors is assumed for generalized propensity score estimation. If TRUE, a lognormal model is assumed. Default is |

`bw` |
Bandwith for the second order Epanechnikov kernel functions of the treatment. If set to NULL, bandwidth computation is based on the rule of thumb for Epanechnikov kernels, determining the bandwidth as the standard deviation of the treatment times 2.34/( |

`boot` |
Number of bootstrap replications for estimating standard errors. Default is 1999. |

`cluster` |
A cluster ID for block or cluster bootstrapping when units are clustered rather than iid. Must be numerical. Default is NULL (standard bootstrap without clustering). |

Estimation of causal mechanisms (natural direct and indirect effects) of a continuous treatment under a selection on observables assumption assuming that all confounders of the treatment and the mediator, the treatment and the outcome, or the mediator and the outcome are observed. Units are weighted by the inverse of their conditional treatment densities (known as generalized propensity scores) given the mediator and/or observed confounders, which are estimated by linear or loglinear regression. Standard errors are obtained by bootstrapping the effects.

A medweightcont object contains two components, `results`

and `ntrimmed`

:

`results`

: a 3X5 matrix containing the effect estimates in the first row ("effects"), standard errors in the second row ("se"), and p-values in the third row ("p-value").
The first column provides the total effect, namely the average treatment effect (ATE) if `ATET=FALSE`

or the average treatment effect on the treated (ATET), i.e. those with D=`d1`

, if `ATET=TRUE`

.
The second and third columns provide the direct effects under treatment and control, respectively ("dir.treat", "dir.control"). The fourth and fifth columns provide the indirect effects under treatment and control, respectively ("indir.treat", "indir.control").

`ntrimmed`

: number of discarded (trimmed) observations due to extreme propensity score values.

Hsu, Y.-C., Huber, M., Lee, Y.-Y., Lettry, L. (2020): "Direct and indirect effects of continuous treatments based on generalized propensity score weighting", Journal of Applied Econometrics, forthcoming.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
# A little example with simulated data (10000 observations)
## Not run:
n=10000
x=runif(n=n,min=-1,max=1)
d=0.25*x+runif(n=n,min=-2,max=2)
d=d-min(d)
m=0.5*d+0.25*x+runif(n=n,min=-2,max=2)
y=0.5*d+m+0.25*x+runif(n=n,min=-2,max=2)
# The true direct and indirect effects are all equal to 0.5
output=medweightcont(y,d,m,x,d0=2,d1=3,ATET=FALSE,trim=0.1,
lognorm=FALSE,bw=NULL,boot=19)
round(output$results,3)
output$ntrimmed
## End(Not run)
``` |

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.