simulate_bounds: Simulate bounds

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

Simulate bounds

Description

Run a simple simulation based on the bounds. For each simulation, sample the set of counterfactual probabilities from a uniform distribution, translate into a multinomial distribution, and then compute the objective and the bounds in terms of the observable variables.

Usage

simulate_bounds(obj, bounds, nsim = 1000)

Arguments

obj	Object as returned by analyze_graph
bounds	Object as returned by optimize_effect_2
nsim	Number of simulation replicates

Value

A data frame with columns: objective, bound.lower, bound.upper

Examples

b <- initialize_graph(graph_from_literal(X -+ Y, Ur -+ X, Ur -+ Y))
obj <- analyze_graph(b, constraints = NULL, effectt = "p{Y(X = 1) = 1} - p{Y(X = 0) = 1}")
bounds <- optimize_effect_2(obj)
simulate_bounds(obj, bounds, nsim = 5)

causaloptim documentation built on Oct. 17, 2024, 9:08 a.m.