# Create Model
model <- (
outcome('Y')
+ within_predictor('X')
+ effect_size(icc = 0.1)
)
# Induce MCAR data on outcome
set.seed(19723)
model |> power_analysis(50, 5, 50, mechanism = MCAR(0.25)) -> powersim_mcar
# Induce MAR data on outcome
set.seed(19723)
model |> power_analysis(
50, 5, 50,
mechanism = MAR(0.25, 'X', 0.6)
) -> powersim_mar
Run the code above in your browser using DataLab