# Illustrative normal-approximation summaries; no participant data or fitted study.
# These synthetic coefficient/SE inputs demonstrate a partly updated paragraph.
adjusted_fit <- list(estimate = 0.24, se = 0.136)
unadjusted_fit <- list(estimate = 0.27, se = 0.123)
summarize_fit <- function(fit) {
  c(estimate = fit$estimate,
    ci_lower = fit$estimate - qnorm(0.975) * fit$se,
    ci_upper = fit$estimate + qnorm(0.975) * fit$se,
    p_value = 2 * pnorm(-abs(fit$estimate / fit$se)))
}
adjusted_result <- summarize_fit(adjusted_fit)
unadjusted_result <- summarize_fit(unadjusted_fit)
print(adjusted_result, digits = 10)
print(unadjusted_result, digits = 10)