Given a object with the linear programming problem set up, compute the bounds using the c++ code developed by Alex Balke. Bounds are returned as text but can be converted to R functions using interpret_bounds, or latex code using latex_bounds. This legacy function has been replaced as default with the more efficient optimize_effect_2.
optimize_effect(obj)Object as returned by analyze_graph
An object of class "balkebound" that contains the bounds and logs as character strings
optimize_effect_2 which is the current default
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}")
optimize_effect(obj)
#> lower bound =
#> MAX {
#> - p10_ - p01_
#> }
#>
#> ----------------------------------------
#> upper bound =
#> MIN {
#> - p10_ - p01_ + 1
#> }
#>