## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")

## ----eval=FALSE---------------------------------------------------------------
# priors <- create_mediation_prior_specification(
#   coefficient_sd = 0.75,
#   group_sd_scale = 1
# )
# 
# spec <- specify_multilevel_gaze_mediation(
#   prepared,
#   mediator_family = "gaussian",
#   outcome_family = "bernoulli",
#   priors = priors,
#   random_slopes = "mediator_x",
#   missingness_policy = "error"
# )

## ----eval=FALSE---------------------------------------------------------------
# fit <- fit_multilevel_gaze_mediation(
#   prepared,
#   mediator_family = "gaussian",
#   outcome_family = "bernoulli",
#   priors = priors,
#   random_slopes = "mediator_x",
#   missingness_policy = "error"
# )
# check_mediation_convergence(fit)
# estimate_within_indirect_effect(fit)
# posterior_predictive_check_mediation(fit)

## ----eval=FALSE---------------------------------------------------------------
# narrow <- create_mediation_prior_specification(coefficient_sd = 0.5)
# wide <- create_mediation_prior_specification(coefficient_sd = 1.5)
# 
# fit_narrow <- fit_multilevel_gaze_mediation(
#   prepared, priors = narrow, random_slopes = "mediator_x"
# )
# fit_wide <- fit_multilevel_gaze_mediation(
#   prepared, priors = wide, random_slopes = "mediator_x"
# )