Estimate Causal Effects of Population Treatment Policies Assuming Clustered Interference

Source: R/Estimator.R

This function implements the estimators from Barkley et al. (2017) for estimating causal effects ("FX") of treatment policies from an observational study when clustered interference is assumed. Clustered interference is also often known as "partial" interference. For the manuscript introducing the methods in clusteredinterference, see: URL https://arxiv.org/abs/1711.04834

policyFX(data, formula, alphas, k_samps = NULL, ..., verbose = FALSE,
  root_options = NULL, nAGQ = 2, return_matrices = FALSE)

Arguments

data

A data.frame (not a tibble). Columns of factor types are not recommended and will sometimes throw (defensive) errors.
formula

The formula defines the different components of the method. The components are specified by outcome | treatment ~ f(covariates) + (1|cluster_id) | cluster_id. The middle component is passed to glmer, so treatment ~ f(covariates) + (1|cluster_id) specifies the model form for the propensity score (i.e., treatment) model. See Details.
alphas

A numeric vector for the probabilities corresponding to the policies of interest. Each entry must be between 0 and 1.
k_samps

The maximum number of vectors to evaluate to estimate the counterfactual probabilities (i.e., \(\omega(A,N,\alpha)\)). Setting to 0 avoids approximation at the cost of increased computation time. Recommended to set <= 5.
...

The dots argument. The user may supply their own target_grid through the dots argument. The target_grid can be made through exported function makeTargetGrid
verbose

A Boolean on whether to print output to stderr. Defaults to FALSE.
root_options

These are passed to multiroot function.
nAGQ

This is the number of Adaptive Gaussian Quadrature points used in the glmer model fitting computation. Defaults to 2. It is recommended to use more than 1.
return_matrices

A Boolean on whether to return the "bread" and "meat" matrices in the sandwich variance. Defaults to FALSE.

Value

A list object including:

  1. estimates: A tidy data.frame with columns estimand, estimate, var, se, LCI and UCI for 95% CI's, and more information.

  2. parameters: An untidy list of the point estimates of all (target and nuisance) parameters.

  3. variance_matrices: When return_matrices is TRUE this is a list object for the "bread" and "meat" matrices in the sandwich variance calculations for each estimand. Otherwise, it is a list object with length 0.

  4. propensity_scores: The estimated propensity scores for each cluster.

  5. model: The treatment model object.

  6. formula: The full formula argument provided, after coercion to a Formula object

Details

These estimators are based on inverse probability-weighting by the propensity score for treatment (IPW) to estimate causal effects of counterfactual policies of interest (i.e., "policy effects") when clustered interference is assumed. The policies of interest correspond to counterfactual scenarios in which treatment may be correlated within clusters.

This method estimates causal contrasts of these policies by estimating the counterfactual treatment probabilities; taking the correlation structures into account requires heavy computational resources, so the user should be patient.

The modeling formula for the propensity score (i.e., treatment) model is specified via the formula formal argument. An example of a model logit-linear fixed effects would be Y | A ~ X1 + X2 + (1 | cluster_ID) | cluster_ID. A similar model that also includes an interaction term is Y | A ~ X1 + X2 + X1:X2 + (1 | cluster_ID) | cluster_ID.

References

Barkley, B. G., Hudgens, M. G., Clemens, J. D., Ali, M., and Emch, M. E. (2017). Causal Inference from Observational Studies with Clustered Interference. arXiv preprint arXiv:1711.04834. (URL: https://arxiv.org/abs/1711.04834.)

Bradley C Saul and Michael G Hudgens (2017). A Recipe for inferference: Start with Causal Inference. Add Interference. Mix Well with R. Journal of Statistical Software 82(2), pp. 1-21. doi: <10.18637/jss.v082.i02> (URL: http://doi.org/10.18637/jss.v082.i02). https://cran.r-project.org/package=inferference. https://github.com/bsaul/inferference.

Bradley Saul (2017). geex: An API for M-Estimation. https://cran.r-project.org/package=geex. https://github.com/bsaul/geex.

See also

Please see the main package vignette at vignette("estimate-policyFX"). It describes the necessary arguments, as well as some extra functionality.

Examples

# NOT RUN {
toy_data <- clusteredinterference::toy_data
causal_fx <- policyFX(
  data = toy_data,
  formula = Outcome | Treatment ~ Age + Distance + (1 | Cluster_ID) | Cluster_ID,
  alphas = c(.3, .5),
  k_samps = 1,
  verbose = FALSE
)
# }