Results are currently different for `calculate_power(g, marginal_power = c(0, 0, 0), sim_seed = 51123, force_closure = TRUE)` and `calculate_power(g, marginal_power = c(0, 0, 0), sim_seed = 51123)`. Should validation exclude 0/1 for marginal power?