Convert bounds string to a function
interpret_bounds(bounds, parameters)The bounds element as returned by optimize_effect
Character vector defining parameters, as returned by analyze_graph
A function that takes arguments for the parameters, i.e., the observed probabilities and returns a vector of length 2: the lower bound and the upper bound.
b <- graph_from_literal(X -+ Y, Ur -+ X, Ur -+ Y)
V(b)$leftside <- c(0,0,0)
V(b)$latent <- c(0,0,1)
V(b)$nvals <- c(2,2,2)
E(b)$rlconnect <- E(b)$edge.monotone <- c(0, 0, 0)
obj <- analyze_graph(b, constraints = NULL, effectt = "p{Y(X = 1) = 1} - p{Y(X = 0) = 1}")
bounds <- optimize_effect(obj)
bounds_func <- interpret_bounds(bounds$bounds, obj$parameters)
bounds_func(.1, .1, .4, .3)
#> lower upper
#> 1 -0.5 0.5
# vectorized
do.call(bounds_func, lapply(1:4, function(i) runif(5)))
#> lower upper
#> 1 -0.8947922 0.10520781
#> 2 -0.8787121 0.12128795
#> 3 -1.0501947 -0.05019472
#> 4 -1.2155072 -0.21550717
#> 5 -1.0243593 -0.02435934