Description Usage Arguments Details Value References Examples

Dynamic treatment effect estimation for assessing the average effects of sequences of treatments (consisting of two sequential treatments). Combines estimation based on (doubly robust) efficient score functions with double machine learning to control for confounders in a data-driven way.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ```
dyntreatDML(
y2,
d1,
d2,
x0,
x1,
s = NULL,
d1treat = 1,
d2treat = 1,
d1control = 0,
d2control = 0,
trim = 0.01,
MLmethod = "lasso",
fewsplits = FALSE,
normalized = TRUE
)
``` |

`y2` |
Dependent variable in the second period (=outcome period), must not contain missings. |

`d1` |
Treatment in the first period, must be discrete, must not contain missings. |

`d2` |
Treatment in the second period, must be discrete, must not contain missings. |

`x0` |
Covariates in the baseline period (prior to the treatment in the first period), must not contain missings. |

`x1` |
Covariates in the first period (prior to the treatment in the second period), must not contain missings. |

`s` |
Indicator function for defining a subpopulation for whom the treatment effect is estimated as a function of the subpopulation's distribution of |

`d1treat` |
Value of the first treatment in the treatment sequence. Default is 1. |

`d2treat` |
Value of the second treatment in the treatment sequence. Default is 1. |

`d1control` |
Value of the first treatment in the control sequence. Default is 0. |

`d2control` |
Value of the second treatment in the control sequence. Default is 0. |

`trim` |
Trimming rule for discarding observations with products of treatment propensity scores in the first and second period that are smaller than |

`MLmethod` |
Machine learning method for estimating the nuisance parameters based on the |

`fewsplits` |
If set to |

`normalized` |
If set to |

Estimation of the causal effects of sequences of two treatments under sequential conditional independence, assuming that all confounders of the treatment in either period and the outcome of interest are observed. Estimation is based on the (doubly robust) efficient score functions for potential outcomes, see e.g. Bodory, Huber, and Laffers (2020), in combination with double machine learning with cross-fitting, see Chernozhukov et al (2018). To this end, one part of the data is used for estimating the model parameters of the treatment and outcome equations based machine learning. The other part of the data is used for predicting the efficient score functions. The roles of the data parts are swapped (using 3-fold cross-fitting) and the average dynamic treatment effect is estimated based on averaging the predicted efficient score functions in the total sample. Standard errors are based on asymptotic approximations using the estimated variance of the (estimated) efficient score functions.

A `dyntreatDML`

object contains ten components, `effect`

, `se`

, `pval`

, `ntrimmed`

, `meantreat`

, `meancontrol`

, `psd1treat`

, `psd2treat`

, `psd1control`

, and `psd2control`

:

`effect`

: estimate of the average effect of the treatment sequence.

`se`

: standard error of the effect estimate.

`pval`

: p-value of the effect estimate.

`ntrimmed`

: number of discarded (trimmed) observations due to low products of propensity scores.

`meantreat`

: Estimate of the mean potential outcome under the treatment sequence.

`meancontrol`

: Estimate of the mean potential outcome under the control sequence.

`psd1treat`

: P-score estimates for first treatment in treatment sequence.

`psd2treat`

: P-score estimates for second treatment in treatment sequence.

`psd1control`

: P-score estimates for first treatment in control sequence.

`psd2control`

: P-score estimates for second treatment in control sequence.

Bodory, H., Huber, M., Laffers, L. (2020): "Evaluating (weighted) dynamic treatment effects by double machine learning", working paper, arXiv preprint arXiv:2012.00370.

Chernozhukov, V., Chetverikov, D., Demirer, M., Duflo, E., Hansen, C., Newey, W., Robins, J. (2018): "Double/debiased machine learning for treatment and structural parameters", The Econometrics Journal, 21, C1-C68.

van der Laan, M., Polley, E., Hubbard, A. (2007): "Super Learner", Statistical Applications in Genetics and Molecular Biology, 6.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 | ```
# A little example with simulated data (2000 observations)
## Not run:
n=2000
# sample size
p0=10
# number of covariates at baseline
s0=5
# number of covariates that are confounders at baseline
p1=10
# number of additional covariates in period 1
s1=5
# number of additional covariates that are confounders in period 1
x0=matrix(rnorm(n*p0),ncol=p0)
# covariate matrix at baseline
beta0=c(rep(0.25,s0), rep(0,p0-s0))
# coefficients determining degree of confounding for baseline covariates
d1=(x0%*%beta0+rnorm(n)>0)*1
# equation of first treatment in period 1
x1=matrix(rnorm(n*p1),ncol=p1)
# covariate matrix for covariates of period 1
beta1=c(rep(0.25,s1), rep(0,p1-s1))
# coefficients determining degree of confounding for additonal covariates of period 1
d2=(x0%*%beta0+x1%*%beta1+0.5*d1+rnorm(n)>0)*1
# equation of second treatment in period 2
y2=x0%*%beta0+x1%*%beta1+1*d1+0.5*d2+rnorm(n)
# outcome equation in period 2
output=dyntreatDML(y2=y2,d1=d1,d2=d2,x0=x0,x1=x1,
d1treat=1,d2treat=1,d1control=0,d2control=0)
cat("dynamic ATE: ",round(c(output$effect),3),", standard error: ",
round(c(output$se),3), ", p-value: ",round(c(output$pval),3))
output$ntrimmed
# The true effect of the treatment sequence is 1.5
## 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.