Learn R Programming

⚠️There's a newer version (1.0.0) of this package.Take me there.

causaloptim: An Interface to Specify Causal Graphs and Compute Bounds on Causal Effects

When causal quantities are not identifiable from the observed data, it still may be possible to bound these quantities using the observed data. We outline a class of problems for which the derivation of tight bounds is always a linear programming problem and can therefore, at least theoretically, be solved using a symbolic linear optimizer. We provide a user friendly graphical interface for setting up such problems via DAGs, which only allow for problems within this class to be depicted. The user can then define linear constraints to further refine their assumptions to meet their specific problem, and then specify a causal query using a text interface. The program converts this user defined DAG, query, and constraints, and returns tight bounds. The bounds can be converted to R functions to evaluate them for specific datasets, and to latex code for publication.

Development status

This package is in stable development. The interface is unlikely to have major changes at this time. New features may be added over time.

Installation

install.packages("causaloptim")
# or
remotes::install_github("sachsmc/causaloptim")

Or use the web application: https://sachsmc.shinyapps.io/causaloptimweb/

Usage

Launch the shiny app to get started, results are saved in the results object:

results <- specify_graph()

References

M.C. Sachs, G. Jonzon, E.E. Gabriel, and A. Sjölander, "A General Method for Deriving Tight Symbolic Bounds on Causal Effects", 2022. Journal of Computational and Graphical Statistics, https://www.tandfonline.com/doi/full/10.1080/10618600.2022.2071905 .

A. Balke and J. Pearl, "Counterfactual Probabilities: Computational Methods,Bounds, and Applications" UCLA Cognitive Systems Laboratory, Technical Report (R-213-B). In R. Lopez de Mantaras and D. Poole (Eds.), Proceedings of the Conference on Uncertainty in Artificial Intelligence (UAI-94), Morgan Kaufmann, San Mateo, CA, 46-54, July 29-31, 1994. https://ftp.cs.ucla.edu/pub/stat_ser/R213-B.pdf .

Copy Link

Version

Install

install.packages('causaloptim')

Monthly Downloads

360

Version

0.9.7

License

MIT + file LICENSE

Issues

Pull Requests

Stars

Forks

Maintainer

Michael Sachs

Last Published

March 24th, 2023

Functions in causaloptim (0.9.7)

latex_bounds

Latex bounds equations
linear_expression

Compute the scalar product of a vector of numbers and a vector of strings
numberOfValues

Get the number of values of a given variable in the graph
plot_graphres

Plot the analyzed graph object
optimize_effect_2

Run the optimizer
optimize_effect

Run the Balke optimizer
print.linearcausalproblem

Print the causal problem
pastestar

Paste with asterisk sep
expand_cond

Expand potential outcome conditions
symb.subtract

Symbolic subtraction
find_cycles

Find cycles in a graph
specify_graph

Shiny interface to specify network structure and compute bounds
plot.linearcausalproblem

Plot the graph from the causal problem
list_to_path

Recursive function to translate an effect list to a path sequence
reduce.sets

Algebraically reduce sets
shortentxt

Shorten strings to 80 characters wide
get_default_effect

Define default effect for a given graph
print_nvals

Print the number of values of each variable/vertex of the analyzed graph object
linear_term

Compute the product of a single numeric scalar and a single string
interpret_bounds

Convert bounds string to a function
update_effect

Update the effect in a linearcausalproblem object
simulate_bounds

Simulate bounds
opt_effect

Compute a bound on the average causal effect
parse_constraints

Parse text that defines a the constraints
parse_effect

Parse text that defines a causal effect
btm_var

Recursive function to get the last name in a list
analyze_graph

Analyze the causal graph to determine constraints and objective
create_q_matrix

Translate response functions into matrix of counterfactuals
create_effect_vector

Translate target effect to vector of response variables
evaluate_objective

Compute the scalar product of a vector of numbers and a vector of both numbers and strings
constant_term

Compute the scalar product of two numeric vectors of the same length
causaloptim-package

An Interface to Specify Causal Graphs and Compute Bounds on Causal Effects
create_R_matrix

Create constraint matrix
create_response_function

Translate regular DAG to response functions
const.to.sets

Translate lists of constraints to lists of vectors