Generates simulated single-case design data using a piecewise multilevel model.
Usage
simulate_scd_data(
n_cases = 5,
n_phase_a = 5,
n_phase_b = 5,
effect_size = 0,
trend_effect = 0,
autocorr = 0,
baseline_mean = 50,
baseline_sd = 10,
baseline_trend = 0
)Arguments
- n_cases
Number of cases to simulate.
- n_phase_a
Number of time points in phase A.
- n_phase_b
Number of time points in phase B.
- effect_size
Standardized effect size for the phase change.
- trend_effect
Standardized trend effect after phase change.
- autocorr
Lag-1 autocorrelation coefficient.
- baseline_mean
Mean of the baseline phase.
- baseline_sd
Standard deviation of the baseline phase.
- baseline_trend
Linear trend in the baseline phase.
Value
A data frame with columns case, time, phase,
phase_num, time_center, time_since_change, and
outcome.
Examples
simulate_scd_data(n_cases = 3, n_phase_a = 5, n_phase_b = 5,
effect_size = 1.0)
#> case time phase phase_num time_center time_since_change outcome
#> 1 C1 1 A 0 0 0 67.33578
#> 2 C1 2 A 0 1 0 64.48657
#> 3 C1 3 A 0 2 0 65.18193
#> 4 C1 4 A 0 3 0 46.15993
#> 5 C1 5 A 0 4 0 68.27125
#> 6 C1 6 B 1 5 1 54.48508
#> 7 C1 7 B 1 6 2 51.34246
#> 8 C1 8 B 1 7 3 56.56169
#> 9 C1 9 B 1 8 4 70.62876
#> 10 C1 10 B 1 9 5 68.13058
#> 11 C2 1 A 0 0 0 68.03483
#> 12 C2 2 A 0 1 0 48.94931
#> 13 C2 3 A 0 2 0 59.82453
#> 14 C2 4 A 0 3 0 32.86697
#> 15 C2 5 A 0 4 0 41.67980
#> 16 C2 6 B 1 5 1 71.00492
#> 17 C2 7 B 1 6 2 58.26180
#> 18 C2 8 B 1 7 3 61.78812
#> 19 C2 9 B 1 8 4 53.01571
#> 20 C2 10 B 1 9 5 50.39551
#> 21 C3 1 A 0 0 0 40.24577
#> 22 C3 2 A 0 1 0 46.61423
#> 23 C3 3 A 0 2 0 61.52347
#> 24 C3 4 A 0 3 0 54.05101
#> 25 C3 5 A 0 4 0 45.29078
#> 26 C3 6 B 1 5 1 58.66749
#> 27 C3 7 B 1 6 2 72.26682
#> 28 C3 8 B 1 7 3 63.32944
#> 29 C3 9 B 1 8 4 56.52912
#> 30 C3 10 B 1 9 5 59.01449