Helper function to evaluate the direct effect

evaluate_direct_effect(Qbarbar, use_conditional = TRUE)