Estimate the Difference between Profiles — contrast

Given a fitted multimedia model, contrast the mediation and outcome predictions associated wtih two treatment profiles.

Usage

contrast_predictions(model, profile1, profile2, ...)

Arguments

model: An object of class multimedia containing the estimated mediation and outcome models whose mediation and outcome predictions we want to compare.
profile1: An object of class treatment_profile containing the first treatment profile to consider in the difference.
profile2: An object of class treatment_profile containing the second treatment profile to consider in the difference.
...: Additional arguments to pass to predict().

Value

A list with two elements, mediators and outcomes, containing the differences in the predicted M(T') - M(T) and Y(T', M(T')) - Y(T, M(T)) between the two profiles T and T'.

Examples

exper <- demo_joy() |>
    mediation_data("PHQ", "treatment", starts_with("ASV"))
model <- multimedia(exper) |>
    estimate(exper)
t1 <- data.frame(treatment = factor("Treatment"))
t2 <- data.frame(treatment = factor("Control"))
profile1 <- setup_profile(model, t1, t1)
profile2 <- setup_profile(model, t2, t2)
contrast_predictions(model, profile1, profile2)
#> $mediators
#>        ASV1      ASV2        ASV3      ASV4      ASV5
#> 1 0.2370731 0.3572813 -0.03221589 0.2615141 0.2436306
#> 
#> $outcomes
#>          PHQ
#> 1 -0.2694043
#>