data-raw/scan-Data.R

In jazznbass/scan: Single-Case Data Analyses for Single and Multiple Baseline Designs

### EXAMPLE DATA for scan ###---------------------------------

library(scan)

filenames <- c()

# example A24 -------------------------------------------------------------

example_A24 <- scdf(
  injuries = c(
    A = 239, 263, 264, 277, 283, 296, 228,
    B = 136, 133, 106, 123, 96, 97, 103, 138, 129, 136, 142, 113, 
        177, 151, 106, 71
  ),
  year = 1996:2018, dvar = "injuries", mvar = "year", name = "A24"
)

scdf_attr(example_A24, "info") <- "Number of injuries on a German autobahn before and after implementation of a speedlimit (130km/h)."
scdf_attr(example_A24, "author") <- "Ministerium fuer Infrastruktur und Landesplanung. Land Brandenburg."

filenames <- c(filenames, "example_A24")

# exampleAB_ABC ------

case1 <- scdf(
  values = c(58, 58, 53, 55, 45, 52, 46, 51, 50, 48, 63, 67, 56, 55, 64, 62, 61, 54, 64, 59, 63, 44, 60, 64, 65), 
  phase_design = c(A = 10, B = 15)
)

case2 <- scdf(
  values = c(54, 50, 50, 47, 57, 44, 55, 54, 66, 55, 67, 59, 55, 66, 53, 46, 45, 55, 48, 55, 48, 50, 48, 45, 53, 46), 
  phase_design = c(A = 8, B = 7, C = 11)
)

case3 <- scdf(
  values = c(56, 57, 52, 54, 43, 48, 48, 56, 51, 66, 63, 59, 57, 56, 74, 75, 74, 68, 73, 67, 62, 52, 59, 63, 54, 64, 65, 61, 55, 63, 55), 
  phase_design = c(A = 9, B = 5, C = 7, D = 10)
)

exampleAB_mpd <- c(case1, case2, case3)

scdf_attr(exampleAB_mpd, "info") <- "A multiple phase design study."
scdf_attr(exampleAB_mpd, "author") <- "Juergen Wilbert"

filenames <- c(filenames, "exampleAB_mpd")

# exampleAB_simple from the scan Manual -----

exampleAB_simple <- c(
  scdf(score = c(5, 7, 8, 5, 7, 12, 16, 18, 15, 14, 19), 
       B_start = 6, name = "Charlotte", dvar = "score"),
  scdf(score = c(3, 4, 3, 5, 7, 8, 7, 9, 8, 10, 12), 
       B_start = 5, name = "Theresa", dvar = "score"),
  scdf(score = c(9, 8, 8, 7, 5, 7, 13, 14, 15, 12, 16), 
       B_start = 7, name = "Antonia", dvar = "score")
)

scdf_attr(exampleAB_simple, "info") <- "A simple multiple baseline AB Design."
scdf_attr(exampleAB_simple, "author") <- "Juergen Wilbert"

filenames <- c(filenames, "exampleAB_simple")

# exampleA1B1A2B2_zvt ------------

p1 <- scdf(
  zvt = c(47, 58, 76, 63, 71, 59, 64, 69, 72, 77, 76, 73),
  d2 = c(131, 134, 141, 141, 140, 140, 138, 140, 141, 140, 138, 140),
  phase_design = c(A1 = 3, B1 = 3, A2 = 3, B2 = 3),
  day = 1:12, dvar = "zvt", pvar = "part", mvar = "day", name = "Tick"
)
p2 <- scdf(
  zvt = c(51, 58, 70, 65, 67, 63, 64, 70, 70, 68, 69, 70),
  d2 = c(100, 126, 130, 130, 137, 133, 136, 137, 135, 128, 137, 138),
  phase_design = c(A1 = 3, B1 = 3, A2 = 3, B2 = 3),
  day = 1:12, dvar = "zvt", pvar = "part", mvar = "day", name = "Trick"
)
p3 <- scdf(
  zvt = c(54, 57, 51, 61, 57, 53, 58, 57, 60, 55, 58, 58),
  d2 = c(89, 116, 114, 131, 132, 130, 128, 131, 130, 129, 118, 131),
  phase_design = c(A1 = 3, B1 = 3, A2 = 3, B2 = 3),
  day = 1:12, dvar = "zvt", pvar = "part", mvar = "day", name = "Track"
)

exampleA1B1A2B2_zvt <- c(p1, p2, p3)

filenames <- c(filenames, "exampleA1B1A2B2_zvt")
rm(list = c("p1", "p2", "p3"), envir = .GlobalEnv)


# Grosche, M. (2011). Effekte einer direkt-instruktiven Förderung  --------

Grosche2011 <- c(
  scdf(c(3.39, 2.09, 1.46, 3.23, 3.06, 2.83, 2.44, 1.86, 2.17, 2.59, 
         3.92, 4.14, 3.31, 2.69, 4.27, 4.66, 3.20, 4.98, 4.42),
    B_start = 38,
    mt = c(0, 1, 7, 8, 13, 15, 38, 42, 45, 49, 50, 56, 63, 64, 70, 
           77, 78, 91, 92)
  ),
  scdf(c(13.00, 17.40, 13.00, 7.74, 8.82, 7.50, 5.82, 6.26, 5.60, 
         12.79, 9.17, 10.59, 9.23, 10.71, 7.61, 9.78, 8.33, 12.00, 
         12.19, 10.70, 11.36, 7.89, 10.82, 11.56),
    B_start = 38,
    mt = c(0, 1, 6, 10, 13, 15, 36, 38, 41, 43, 48, 50, 55, 59, 62, 
           64, 69, 71, 77, 78, 85, 91, 94, 101)
  ),
  scdf(c(7.69, 7.27, 7.02, 5.37, 4.57, 3.02, 4.68, 4.97, 5.73, 7.63, 
         6.98, 7.50, 6.91, 4.86, 6.71, 5.67, 5.00, 6.14, 8.14, 4.53),
    B_start = 59,
    mt = c(0, 1, 7, 8, 14, 15, 38, 42, 43, 49, 50, 56, 59, 63, 64, 70, 
           71, 80, 91, 97)
  )
)

names(Grosche2011) <- c("Eva", "Georg", "Olaf")

scdf_attr(Grosche2011, "info") <- "Direct instruction intervention on reading accuracy."
scdf_attr(Grosche2011, "author") <- "Grosche, M. (2011). Effekte einer direkt-instruktiven F\u00f6rderung der Lesegenauigkeit. Empirische Sonderp\u00E4dagogik, 3(2), 147-161."

filenames <- c(filenames, "Grosche2011")


# Data from students learning vocabulary by heart and checking the --------
case1 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(0, 1, 0, 0, 1, 2, 5, 6, 7, 10, 10, 15, 16, 14, 17, 18, 
             17, 19, 18, 18),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 
         18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Lisa (Turkish)"
)

case2 <- scdf(
  phase_design = c(A = 5, B = 20),
  values = c(2, 1, 3, 0, 1, 4, 8, 8, 8, 12, 13, 13, 15, 14, 15, 
             16, 18, 18, 20, 20, 20, 19, 20, 19, 20),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 
         17, 18, 19, 20, 21, 22, 23, 24, 25),
  dvar = "values",
  mvar = "mt",
  name = "Patrick (Spanish)"
)

case3 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(0, 1, 0, 1, 0, 2, 4, 3, 0, 1, 3, 6, 8, 8, 10, 4, 6, 9, 11, 10),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Anna (Twi)"
)

case4 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(0, 0, 0, 0, 0, 1, 1, 4, 2, 4, 3, 5, 5, 6, 6, 6, 7, 10, 14, 12),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Lena (Turkish)"
)

case5 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(0, 0, 0, 1, 0, 2, 0, 2, 2, 3, 2, 3, 1, 5, 4, 3, 6, 8, 6, 7),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Nicola (Turkish)"
)

case6 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(0, 0, 0, 0, 0, 0, 2, 3, 2, 3, 4, 0, 3, 5, 4, 6, 5, 7, 7, 8),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Melanie (Swedish)"
)

case7 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(1, 1, 2, 1, 2, 1, 4, 2, 3, 3, 2, 3, 1, 3, 3, 4, 2, 4, 3, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Kira (Italian)"
)

case8 <- scdf(
  phase_design = c(A = 5, B = 21),
  values = c(3, 4, 7, 6, 8, 10, 11, 8, 9, 13, 17, 11, 14, 9, 13, 14, 
             16, 19, 18, 20, 20, 20, 19, 20, 18, 20),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 
         18, 19, 20, 21, 22, 23, 24, 25, 26),
  dvar = "values",
  mvar = "mt",
  name = "Lynn (Spanisch)"
)

case9 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(2, 0, 1, 1, 2, 5, 4, 4, 5, 8, 6, 7, 6, 6, 8, 9, 7, 8, 8, 7),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 
         18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Lars (Turkish)"
)

case10 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(0, 1, 1, 0, 0, 2, 1, 3, 2, 3, 4, 5, 4, 4, 5, 6, 5, 3, 6, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 
         18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Katharina (Turkish)"
)

case11 <- scdf(
  phase_design = c(A = 5, B = 15),
  values = c(1, 0, 2, 1, 0, 2, 3, 5, 4, 6, 5, 7, 4, 5, 5, 8, 7, 9, 10, 9),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 
         18, 19, 20),
  dvar = "values",
  mvar = "mt",
  name = "Anja (Italian)"
)

byHeart2011 <- c(case1, case2, case3, case4, case5, case6, case7, case8, case9, case10, case11)
scdf_attr(byHeart2011, "info") <- "Data from university students learning vocabulary by heart and checking their progress with 20 flashcards each session."
scdf_attr(byHeart2011, "author") <- "Juergen Wilbert, 2011"
filenames <- c(filenames, "byHeart2011")


# Grünke, M., Wilbert, J., & Stegemann-Calder, K. (2013). Analyzin --------

case1 <- scdf(
  phase_design = c(A = 4, B = 14),
  score = c(6, 5, 4, 5, 8, 10, 10, 8, 8, 10, 9, 9, 10, 9, 9, 10, 9, 9),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18),
  dvar = "score",
  mvar = "mt",
  name = ""
)

case2 <- scdf(
  phase_design = c(A = 7, B = 11),
  score = c(3, 4, 2, 2, 3, 4, 3, 9, 9, 8, 9, 8, 10, 9, 7, 9, 10, 9),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18),
  dvar = "score",
  mvar = "mt",
  name = ""
)

case3 <- scdf(
  phase_design = c(A = 6, B = 12),
  score = c(4, 5, 4, 3, 4, 3, 8, 9, 9, 9, 8, 10, 9, 8, 10, 9, 8, 9),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18),
  dvar = "score",
  mvar = "mt",
  name = ""
)

case4 <- scdf(
  phase_design = c(A = 8, B = 10),
  score = c(4, 5, 5, 3, 4, 4, 3, 5, 10, 8, 9, 10, 9, 10, 9, 7, 9, 8),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18),
  dvar = "score",
  mvar = "mt",
  name = ""
)

case5 <- scdf(
  phase_design = c(A = 5, B = 13),
  score = c(4, 5, 5, 4, 5, 8, 9, 9, 8, 10, 10, 10, 8, 10, 9, 9, 10, 8),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18),
  dvar = "score",
  mvar = "mt",
  name = ""
)

case6 <- scdf(
  phase_design = c(A = 4, B = 14),
  score = c(4, 3, 3, 2, 7, 8, 10, 9, 10, 10, 10, 10, 9, 9, 9, 8, 9, 8),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18),
  dvar = "score",
  mvar = "mt",
  name = ""
)

GruenkeWilbert2014 <- c(case1, case2, case3, case4, case5, case6)
names(GruenkeWilbert2014) <- c("Anton", "Bob", "Paul", "Robert", "Sam", "Tim")
scdf_attr(GruenkeWilbert2014, "info") <- "Data from an intervention study on text comprehension."
scdf_attr(GruenkeWilbert2014, "author") <- "Gruenke, M., Wilbert, J., & Stegemann-Calder, K. (2013). Analyzing the effects of story mapping on the reading comprehension of children with low intellectual abilities. Learning Disabilities: A Contemporary Journal, 11(2), 51-64."

filenames <- c(filenames, "GruenkeWilbert2014")


# Grosche, Lüke, Wilbert (in preparation) ---------------------------------

Grosche2014 <-
  list(
    "Case1" =
      list(
        "Material 1" = scdf(
          c(5, 11, 12, 13, 20, 21, 33, 34, 24, 41, 45, 24, 29, 
            34, 30, 41, 30, 31, 39, 37, 39), 
          mt = c(2, 4, 8, 10, 15, 16, 17, 22, 23, 24, 29, 31, 
                 37, 43, 45, 50, 52, 57, 58, 64, 66), 
          B_start = 15
        ),
        "Material 2" = scdf(
          c(9, 12, 16, 14, 12, 12, 14, 13, 15, 16, 16, 16, 19, 
            25, 27, 30, 30), 
          mt = c(15, 16, 17, 22, 23, 24, 29, 31, 37, 43, 45, 
                 50, 52, 57, 58, 64, 66), 
          B_start = 45
        )
      ),
    "Case2" =
      list(
        "Material 1" = scdf(
          c(17, 16, 19, 18, 17, 16, 22, 32, 27, 24, 27, 28, 26, 32, 
            37, 37, 40, 39, 48, 51), 
          mt = c(2, 4, 8, 10, 15, 16, 17, 22, 23, 30, 32, 37, 43, 45, 
                 50, 52, 57, 58, 64, 66), 
          B_start = 15
        ),
        "Material 2" = scdf(
          c(18, 16, 19, 18, 14, 18, 17, 17, 21, 29, 31, 24, 36, 43, 45, 59), 
          mt = c(15, 16, 17, 22, 23, 30, 32, 37, 43, 45, 50, 52, 57, 
                 58, 64, 66), 
          B_start = 45
        ),
        "Material 3" = scdf(
          c(6, 4, 5, 3, 11, 7, 7, 4, 8, 9, 14, 17, 15, 17, 33, 32), 
          mt = c(15, 16, 17, 22, 23, 30, 32, 37, 43, 45, 
                 50, 52, 57, 58, 64, 66), 
          B_start = 58
        )
      ),
    "Case3" =
      list(
        "Material 1" = scdf(
          c(24, 23, 20, 26, 27, 32, 35, 38, 35, 30, 40, 35, 37, 33, 27, 
            42, 37, 43, 52, 47, 36), 
          mt = c(2, 6, 8, 10, 14, 15, 16, 21, 22, 23, 28, 35, 38, 43, 45, 
                 50, 52, 57, 59, 64, 66), 
          B_start = 14),
        "Material 2" = scdf(
          c(31, 22, 31, 34, 26, 37, 39, 32, 24, 28, 35, 41, 36, 41, 42, 42, 43), 
          mt = c(14, 15, 16, 21, 22, 23, 28, 35, 38, 43, 45, 50, 52, 
                 57, 59, 64, 66), 
          B_start = 35),
        "Material 3" = scdf(
          c(23, 20, 16, 21, 23, 16, 29, 29, 18, 28, 21, 32, 34, 30, 35, 42, 36), 
          mt = c(14, 15, 16, 21, 22, 23, 28, 35, 38, 43, 45, 50, 
                 52, 57, 59, 64, 66), 
          B_start = 52)
      ),
    "Case4" =
      list(
        "Material 1" = scdf(
          c(16, 11, 17, 17, 23, 33, 41, 25, 25, 27, 33, 32, 27, 30, 
            38, 17, 38, 29, 34, 39, 44), 
          mt = c(1, 3, 7, 9, 14, 16, 21, 22, 23, 29, 31, 35, 38, 43, 
                 45, 50, 59, 63, 64, 65, 66), 
          B_start = 14),
        "Material 2" = scdf(
          c(9, 20, 21, 18, 18, 16, 22, 27, 28, 21, 30, 30, 33, 30, 31, 45, 39), 
          mt = c(14, 16, 21, 22, 23, 29, 31, 35, 38, 43, 45, 50, 
                 59, 63, 64, 65, 66), 
          B_start = 35),
        "Material 3" = scdf(
          c(2, 8, 4, 8, 3, 4, 8, 5, 4, 6, 8, 11, 21, 21, 24, 25, 30), 
          mt = c(14, 16, 21, 22, 23, 29, 31, 35, 38, 43, 45, 
                 50, 59, 63, 64, 65, 66), 
          B_start = 59)
      )
  )


scdf_attr(Grosche2014, "info") <- "Data from a multiple material multi person intervention study on reading."
scdf_attr(Grosche2014, "author") <- "Michael Grosche, Timo Lueke, and Juergen Wilbert"

filenames <- c(filenames, "Grosche2014")

# Behavioral data (compliance in percent) from a study by Christia --------

case1 <- scdf(
  phase_design = c(A = 10, B = 27),
  compliance = c(25, 20.8, 39.6, 75, 45, 39.6, 54.2, 50, 28.1, 40, 52.1, 
                 31.3, 15.6, 29.2, 43.8, 25, 12.5, 40, 29.16666667, 25, 
                 31.25, 22.91666667, 45.83333333, 37.5, 50, 25, 25, 20, 
                 14.58333333, 43.75, 25, 43.75, 29.16666667, 12.5, 18.75, 
                 27.5, 25),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 
         19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 
         35, 36, 37),
  dvar = "compliance",
  mvar = "mt",
  name = "Adam"
)

case2 <- scdf(
  phase_design = c(A = 10, B = 19),
  compliance = c(25, 20.8, 39.6, 75, 45, 14.6, 45.8, 33.3, 31.3, 32.5, 4.2, 
                 45.8, 31.3, 9.4, 6.3, 4.2, 10, 31.25, 25, 12.5, 22.91666667, 
                 37.5, 6.25, 4.166666667, 8.333333333, 8.333333333, 37.5, 
                 17.5, 10.41666667),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 
         20, 21, 22, 23, 24, 25, 26, 27, 28, 29),
  dvar = "compliance",
  mvar = "mt",
  name = "Berta"
)

case3 <- scdf(
  phase_design = c(A = 10, B = 66),
  compliance = c(58.3, 14.6, 29.2, 25, 54.2, 91.7, 22.9, 33.3, 25, 50, 33.3, 
                 14.6, 50, 45.8, 45.8, 33.3, 85.4, 91.66666667, 60.41666667, 
                 66.66666667, 83.33333333, 83.33333333, 100, 81.25, 90, 66.66666667, 75, 100, 52.08333333, 58.33333333, 56.25, 66.66666667, 70.83333333, 
                 66.66666667, 83.33333333, 72.91666667, 89.58333333, 81.25, 
                 87.5, 85.41666667, 77.08333333, 85.41666667, 85.41666667, 
                 91.66666667, 85.41666667, 91.66666667, 89.58333333, 
                 91.66666667, 85.41666667, 81.25, 85.41666667, 79.16666667, 
                 89.58333333, 87.5, 85.41666667, 91.66666667, 95.83333333, 
                 87.5, 91.66666667, 47.91666667, 81.25, 68.75, 95.83333333, 
                 85.41666667, 79.16666667, 95.83333333, 100, 81.25, 93.75, 
                 95.83333333, 83.33333333, 83.33333333, 87.5, 100, 90, 100),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 
         20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 
         37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 
         54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 
         71, 72, 73, 74, 75, 76),
  dvar = "compliance",
  mvar = "mt",
  name = "Christian"
)

case4 <- scdf(
  phase_design = c(A = 10, B = 66),
  compliance = c(65.6, 37.5, 58.3, 72.9, 33.3, 59.4, 77.1, 54.2, 68.8, 43.8, 62.5, 64.6, 60.4, 81.3, 79.2, 77.1, 85.4, 79.16666667, 66.66666667, 75, 70.83333333, 58.33333333, 77.08333333, 77.08333333, 68.75, 75, 79.16666667, 79.16666667, 68.75, 91.66666667, 64.58333333, 87.5, 58.33333333, 64.58333333, 87.5, 85.41666667, 75, 83.33333333, 87.5, 83.33333333, 87.5, 89.58333333, 93.75, 91.66666667, 91.66666667, 70.83333333, 87.5, 77.08333333, 85.41666667, 70.83333333, 93.75, 93.75, 87.5, 85.41666667, 93.75, 90, 100, 75, 83.33333333, 83.33333333, 93.75, 91.66666667, 87.5, 79.16666667, 87.5, 87.5, 75, 85.41666667, 87.5, 87.5, 75, 91.66666667, 83.33333333, 85.41666667, 79.16666667, 87.5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76),
  dvar = "compliance",
  mvar = "mt",
  name = "David"
)

Huber2014 <- c(case1, case2, case3, case4)
scdf_attr(Huber2014, "info") <- "Behavioral data (compliance in percent)."
scdf_attr(Huber2014, "author") <- "Christian Huber"

filenames <- c(filenames, "Huber2014")

# Leidig 2018 Good Behavior Game ------------------------------------------

case1 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(4, 1, 2, NA, 2, 3, 1, 1, 3, 4, 3, 4, 4, 4, 4, 4, NA, NA, NA, NA, 4, 3, 5, 4, 5, 0, 2, NA, 3, 4, 3, 4, NA, 4, 2, NA, NA, 3, 4, 4, 3, 3, 4, 4, 4, 3, 3, 4, 5, NA, 3, NA, 3, 1, 1, 4, 1, 4, NA, NA, 5, 3, 3, 4, 3, 3, 4, 3, 5, 5, NA, NA, NA, 3, 4, 4, 4, 3, NA, NA, 3, 3, 5, 5, 5, 3, NA, NA, 5, 2, 2, NA, NA, 3, NA, 2, 4, 4, 5, 4, 5, 3, 4, NA, NA, 4, 4, 3),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 1, 1, NA, 1, 1, 1, 2, 0, 0, 0, 1, 0, 0, 0, 0, NA, NA, NA, NA, 1, 0, 0, 0, 0, 1, 0, NA, 1, 0, 0, 0, NA, 0, 1, NA, NA, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, NA, 0, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, NA, NA, NA, 0, 0, 0, 0, 1, NA, NA, 0, 1, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, NA, 2, NA, 1, 0, 1, 1, 0, 0, 0, 0, NA, NA, 0, 0, 1),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "1a1"
)

case2 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(4, 2, 1, NA, 3, 3, 1, 1, 2, 2, 4, 3, 4, 4, 4, 4, NA, NA, NA, NA, 3, 5, 5, 4, 4, 4, 2, 4, 4, 3, 2, 3, 3, 4, 2, 4, 4, 4, 2, 4, 3, 1, 5, 4, 3, 3, 3, 4, 2, NA, 3, 4, 3, 1, 3, 3, 4, 4, NA, NA, 3, 4, 4, 5, NA, 4, 5, 5, 4, 4, NA, NA, NA, 4, 3, 4, 4, 3, NA, NA, 3, 4, 5, 5, 5, 4, NA, NA, 3, 2, 5, 5, 5, 4, 0, 2, 5, 5, 5, 3, 4, 4, 5, NA, NA, 4, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(0, 1, 2, NA, 1, 2, 1, 1, 3, 4, 0, 0, 0, 1, 0, 0, NA, NA, NA, NA, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 1, 2, 0, 0, 1, 2, 0, 0, 2, NA, 0, 1, 2, 0, 0, 0, 0, 0, NA, NA, 0, 1, 1, 0, NA, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 1, NA, NA, 1, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 2, NA, 2, 0, 1, 1, 0, 0, 0, 0, NA, NA, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "1a2"
)

case3 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(3, 4, 4, NA, 4, 4, 4, 2, 3, 4, 5, 4, 4, 5, 4, 5, NA, NA, NA, NA, 4, 4, 5, NA, 5, 4, 4, 5, 4, 4, 5, 5, 5, 5, 5, 4, 4, 4, 4, 3, 4, 3, 4, 5, 5, 4, 5, 4, 5, NA, 5, NA, 4, 5, 5, 5, 5, 5, NA, NA, 3, 3, 5, 5, 3, 4, NA, 3, 5, 5, NA, NA, NA, 3, 4, 4, 4, 3, NA, NA, 5, 4, 5, 5, 4, 3, NA, NA, 3, 4, 2, NA, 5, 2, 3, 4, 4, 3, 4, 4, 4, 3, 4, NA, NA, 4, 4, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(2, 2, 1, NA, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, 0, NA, 0, 0, 0, 0, 0, 0, NA, NA, 1, 1, 0, 0, 0, 1, NA, 1, 0, 0, NA, NA, NA, 0, 0, 1, 0, 1, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, 0, 3, 3, 1, 0, 1, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "1a3"
)

case4 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(4, 3, 4, NA, 3, 4, 3, 3, 4, 4, 5, 4, 4, 5, 4, 5, NA, NA, NA, NA, 5, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 3, 4, 4, 5, 4, 5, 5, 5, NA, NA, NA, NA, 5, NA, 5, 5, 4, 5, 5, 5, 5, 5, NA, NA, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, NA, NA, NA, 5, 5, 5, 5, 5, NA, NA, 5, 5, 5, 5, 5, 5, NA, NA, 4, 5, 4, 5, 5, 5, 0, 4, 5, 5, 5, 5, 5, 4, 5, NA, NA, 5, 5, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 2, 3, NA, 4, 1, 2, 3, 1, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 2, NA, 1, 0, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "1a4"
)

case5 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(3, 3, 4, NA, 1, 4, 4, 3, 3, 4, 4, 2, 4, 5, 4, 4, NA, NA, NA, NA, 3, 4, 5, 4, 5, 4, 5, 5, 5, 4, 5, 5, 5, 5, 5, 4, 4, 5, 4, 3, 4, 5, 4, 5, 4, 5, 4, 4, 5, NA, 3, NA, 4, 2, 4, 4, 4, 5, NA, NA, 4, 5, 5, 5, 4, 5, NA, 5, 5, 4, NA, NA, NA, 3, 4, 4, 4, 3, NA, NA, 5, 4, 5, 5, 5, 5, NA, NA, 3, 4, 4, NA, 5, 4, 3, 4, 5, 5, 5, 5, 5, 4, 5, NA, NA, 5, 5, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a", "1a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(2, 2, 2, NA, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 0, 0, 0, NA, 0, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, NA, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 1, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, 0, 1, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "1a5"
)

case6 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(1, 1, 1, NA, 3, 3, 2, 2, 4, 4, 4, 2, 3, 4, 3, 3, NA, NA, NA, NA, 3, 4, 2, 3, 3, 3, 3, 4, 3, 4, 3, 2, 3, 3, NA, NA, NA, 4, 4, 5, 4, 4, NA, 3, NA, 5, 4, NA, NA, NA, 4, 4, NA, NA, 4, 4, 2, NA, NA, 3, 3, NA, NA, NA, 4, NA, 4, NA, NA, NA, 2, NA, NA, NA, 4, 4, 5, 5, NA, NA, 5, 2, 3, 4, 4, 3, NA, NA, 4, NA, 2, NA, 4, 4, 3, 3, NA, NA, 3, 3, 3, 3, 5, NA, 3, NA, 3, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(4, 2, 2, NA, 2, 1, 3, 1, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, 1, NA, 0, 0, NA, NA, NA, 0, 0, NA, NA, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, 0, NA, 0, NA, NA, NA, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, NA, 1, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, NA, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2a1"
)

case7 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(4, 1, NA, NA, 3, 3, 1, 2, 4, 4, 3, 1, 2, NA, 2, 1, NA, NA, NA, NA, 2, 3, 3, 1, 3, 3, 3, 3, 1, 0, 3, 1, 2, 1, NA, NA, NA, 3, 3, 1, 2, 3, NA, NA, 1, 1, NA, 4, NA, 3, 1, 2, NA, NA, 3, 3, 1, NA, NA, 3, 4, NA, NA, NA, 2, NA, 2, NA, NA, NA, 2, NA, NA, NA, NA, NA, 3, 4, NA, NA, 5, 2, 2, 5, 4, 1, NA, NA, 5, NA, 2, NA, 1, NA, 2, NA, NA, NA, 3, 2, 3, 3, NA, NA, 3, NA, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 3, NA, NA, 2, 2, 2, 1, 0, 0, 0, 0, 0, NA, 0, 0, NA, NA, NA, NA, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 2, 0, NA, NA, NA, 0, 0, 1, 0, 0, NA, NA, 0, 0, NA, 0, NA, 0, 0, 0, NA, NA, 0, 0, 1, NA, NA, 0, 0, NA, NA, NA, 0, NA, 1, NA, NA, NA, 0, NA, NA, NA, NA, NA, 0, 0, NA, NA, 0, 0, 1, 0, 0, 0, NA, NA, 0, NA, 0, NA, 1, NA, 0, NA, NA, NA, 1, 0, 0, 0, NA, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2a2"
)

case8 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(4, 0, 1, NA, 2, 3, 2, 1, 4, 3, 3, 1, 2, NA, 3, 2, NA, NA, NA, NA, 3, 3, 2, 2, 2, 2, 3, 2, 2, 3, 3, 2, 3, 3, NA, NA, NA, 4, 4, 4, 4, 4, NA, 3, 4, 4, 3, NA, NA, 2, 4, 4, NA, NA, 3, 3, 2, NA, NA, 3, 4, NA, NA, NA, 3, 3, 2, NA, NA, NA, NA, NA, NA, NA, 3, 4, 4, 4, NA, NA, 5, 3, 4, 5, 3, 3, NA, NA, 5, NA, 2, NA, 4, 4, 4, 3, NA, NA, 4, 4, 4, 3, NA, NA, NA, NA, NA, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 3, 4, NA, 2, 2, 2, 1, 0, 0, 0, 0, 0, NA, 0, 0, NA, NA, NA, NA, 0, 1, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, NA, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, 0, 0, 0, NA, NA, NA, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, NA, 0, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, NA, NA, NA, NA, NA, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2a3"
)

case9 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(3, 1, 1, NA, 2, 1, 1, 1, 4, 3, 4, 3, 3, 4, 2, 3, NA, NA, NA, NA, 3, 4, 3, 3, 3, 3, 3, 3, 3, 3, 4, 3, 4, 3, NA, NA, NA, 4, 4, 5, 4, 5, NA, 4, 4, 5, 2, NA, NA, 4, 4, 4, NA, NA, 4, 4, 2, NA, NA, 3, 4, NA, NA, NA, 3, 4, 3, NA, NA, NA, 2, NA, NA, NA, 4, 3, 4, 4, NA, NA, 4, 3, 4, 5, 4, 3, NA, NA, 5, NA, 4, NA, 4, 4, 4, 3, NA, NA, 4, 3, 5, 4, NA, NA, 3, NA, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(4, 4, 4, NA, 3, 4, 4, 3, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 2, 0, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, NA, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, 0, 0, 0, NA, NA, NA, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 0, NA, 0, NA, 0, 0, 1, 0, NA, NA, 0, 0, 0, 0, NA, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2a4"
)

case10 <- scdf(
  phase_design = c(A = 8, B = 100),
  academic_engagement = c(4, 1, 4, NA, 1, 2, 1, NA, 4, 4, NA, 2, 2, 4, 3, 3, NA, NA, NA, NA, 2, 3, 3, 3, 2, 3, 3, 4, 2, 3, 4, 2, 4, 3, NA, NA, NA, NA, NA, NA, NA, 5, NA, 3, 4, 5, 3, NA, NA, 4, 4, 4, NA, NA, NA, 4, 2, NA, NA, 3, 4, NA, NA, NA, NA, 3, 3, NA, NA, NA, 2, NA, NA, NA, 3, 2, 4, 4, NA, NA, 4, 2, 4, 5, 4, 2, NA, NA, 4, NA, 3, NA, 4, 4, 3, 3, NA, NA, 2, 4, 4, 3, NA, NA, 4, NA, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a", "2a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 1, 0, NA, 1, 2, 1, NA, 0, 0, NA, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, NA, NA, 0, 0, NA, NA, 0, 0, NA, NA, NA, NA, 0, 0, NA, NA, NA, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, 1, NA, 0, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, 0, NA, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2a5"
)

case11 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(NA, 1, 2, NA, 2, 2, 2, NA, 2, 2, 2, 2, 2, 2, 2, 2, NA, NA, NA, NA, NA, 3, 3, 4, 4, 4, 3, 2, NA, 2, 4, 4, 2, 2, NA, NA, NA, 3, 2, 3, NA, 3, 3, 3, 3, 3, 3, 3, 3, 3, NA, 3, 3, 3, 3, 3, NA, NA, NA, 2, 2, 2, 3, 2, 2, 3, NA, NA, 2, 2, NA, NA, NA, NA, 2, 2, 3, NA, NA, NA, 3, 2, 2, 2, 3, 2, NA, NA, NA, 2, 2, 3, 2, NA, NA, 2, NA, NA, NA, NA, NA, NA, 2, NA, 2, 2, 2, 2),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 4, 4, NA, 3, 3, 3, NA, 3, 2, 3, 1, 3, 2, 3, 4, NA, NA, NA, NA, NA, 0, 0, 0, 0, 0, 0, 0, NA, 1, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, 0, 1, NA, NA, NA, 0, 3, 3, 0, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, NA, 1, 3, 0, NA, NA, NA, 0, 0, 0, 1, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, 0, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2b1"
)

case12 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(NA, 3, 4, NA, 4, 3, 3, NA, 4, 5, 4, 4, 4, 4, 4, 4, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 4, 4, 4, 4, NA, NA, NA, 2, 3, 4, NA, 4, 4, 4, NA, NA, NA, 5, 5, 5, 4, 4, 5, 4, 4, 5, NA, NA, NA, 5, 3, 3, 3, 4, 4, 5, NA, NA, 5, 5, NA, NA, NA, NA, 4, 5, 5, NA, NA, NA, 5, 5, 5, 5, 5, 5, NA, NA, NA, 5, 4, 5, 5, 5, NA, 5, NA, NA, NA, NA, NA, NA, NA, NA, 5, 4, 5, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 4, 3, NA, 3, 3, 3, NA, 3, 1, 2, 1, 3, 2, 3, 4, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 1, 0, NA, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 2, 2, 1, 0, 1, 0, NA, NA, 0, 0, NA, NA, NA, NA, 0, 1, 0, NA, NA, NA, 0, 0, 0, 1, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2b2"
)

case13 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(NA, 2, NA, NA, 3, 2, 2, NA, 2, 3, 3, NA, 3, 2, 2, 2, NA, NA, NA, NA, NA, 2, 3, 4, 2, 2, NA, 3, NA, 3, 3, 2, 3, 2, NA, NA, NA, 2, NA, NA, NA, NA, NA, 2, 3, 3, 3, 2, 3, 3, 3, 3, 3, 3, 3, 3, NA, NA, NA, 4, 3, 3, 3, 3, 2, 3, NA, NA, NA, 3, NA, NA, NA, NA, 2, 2, 2, NA, NA, NA, 3, 3, 2, 3, 3, 3, NA, NA, NA, 3, 3, 2, 3, 3, NA, 3, NA, NA, NA, NA, NA, NA, 3, NA, 2, 2, 3, 3),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 3, NA, NA, 1, 2, 2, NA, 2, 2, 1, NA, 2, 1, 2, 2, NA, NA, NA, NA, NA, 2, 0, 0, 1, 1, NA, 0, NA, 0, 0, 1, 0, 0, NA, NA, NA, 0, NA, NA, NA, NA, NA, 0, 0, 0, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, 1, 0, 1, 2, 0, NA, NA, 0, 0, NA, NA, NA, NA, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2b3"
)

case14 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(NA, 2, 1, NA, 3, 2, 3, NA, 3, 3, 3, 2, 3, 2, 2, 2, NA, NA, NA, NA, NA, 4, 3, 4, 4, 4, 3, 3, NA, 2, 3, 4, 2, 2, NA, NA, NA, 2, 3, 3, NA, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, NA, NA, NA, 4, 3, 3, 3, 3, 2, 3, NA, NA, 3, 3, NA, NA, NA, NA, 3, 4, 3, NA, NA, NA, 3, 3, 3, 3, 3, 3, NA, NA, NA, 3, 3, 2, 3, 3, NA, 3, NA, NA, NA, NA, NA, NA, 3, NA, 3, 3, 3, 3),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 2, 2, NA, 1, 1, 1, NA, 2, 1, 1, 3, 2, 2, 3, 2, NA, NA, NA, NA, NA, 0, 0, 0, 0, 0, 2, 0, NA, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, 1, 0, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, NA, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2b4"
)

case15 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(NA, 2, 3, NA, 3, 3, 3, NA, 3, 3, 3, 3, 3, 3, 3, 3, NA, NA, NA, NA, NA, 4, 3, 4, 3, 3, 4, 3, NA, NA, NA, 4, 3, NA, NA, NA, NA, 4, 4, 3, NA, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, NA, NA, NA, 4, 3, 3, 3, 3, 3, 4, NA, NA, 3, 3, NA, NA, NA, NA, 3, 4, 3, NA, NA, NA, 4, 4, 4, 4, 4, 4, NA, NA, NA, 4, 4, 4, 4, 4, NA, 4, NA, NA, NA, NA, NA, NA, 4, NA, 4, 4, 4, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b", "2b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 3, 2, NA, 2, 2, 1, NA, 2, 1, 2, 1, 1, 2, 3, 2, NA, NA, NA, NA, NA, 0, 0, 0, 1, 1, 0, 1, NA, NA, NA, 0, 0, NA, NA, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 1, 2, 0, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, NA, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, 0, 0, 1, 0, NA, 0, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2b5"
)

case16 <- scdf(
  phase_design = c(A = 38, B = 70),
  academic_engagement = c(4, 4, 4, NA, 4, 3, 3, 4, 4, 4, NA, 4, 4, 4, 3, 4, NA, NA, NA, NA, 3, 4, 4, 4, 4, 3, 4, 3, 4, 4, NA, NA, 2, 4, NA, NA, NA, 4, 4, 3, 5, 4, 5, 4, 4, 4, 4, 5, NA, 4, 4, 4, NA, 4, 5, 5, 5, 4, NA, 4, NA, 4, 3, 4, 4, 4, 4, NA, NA, NA, 4, NA, NA, 3, 4, 4, 3, 5, NA, NA, 4, 4, 4, 4, 4, 4, NA, NA, NA, NA, NA, 4, NA, 5, 5, 4, NA, NA, NA, 5, 4, 4, 4, NA, NA, NA, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(2, 1, 1, NA, 2, 2, 1, 1, 1, 1, NA, 1, 1, 2, 2, 0, NA, NA, NA, NA, 2, 1, 1, 2, 1, 2, 1, 2, 1, 4, NA, NA, 3, 1, NA, NA, NA, 1, 1, 2, 0, 0, 1, 0, 0, 0, 1, 0, NA, 0, 0, 0, NA, 0, 1, 0, 0, 0, NA, 1, NA, 1, 2, 1, 2, 2, 1, NA, NA, NA, 1, NA, NA, 2, 1, 1, 1, 1, NA, NA, 2, 1, 2, 1, 1, 1, NA, NA, NA, NA, NA, 1, NA, 0, 0, 1, NA, NA, NA, 0, 1, 1, 1, NA, NA, NA, 1, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2c1"
)

case17 <- scdf(
  phase_design = c(A = 38, B = 70),
  academic_engagement = c(4, 5, 4, NA, 4, 4, 4, 4, 3, 3, NA, 3, 4, 4, 2, 4, NA, NA, NA, NA, 3, 5, 4, 4, 3, 3, 3, 4, 4, 3, NA, NA, 4, 3, NA, NA, NA, 4, 4, 4, 5, 4, 1, 5, 5, 5, 3, 4, 3, 4, 4, 5, 5, 4, 4, 5, 5, 4, NA, 4, NA, 5, 4, 4, 4, 4, 4, NA, NA, NA, 4, NA, NA, 4, 4, 4, 4, 5, NA, NA, 3, 4, NA, 4, 4, 4, NA, NA, NA, NA, NA, 5, NA, 5, 5, 4, NA, 5, NA, 5, 5, 4, 4, NA, 5, NA, 5, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 0, 1, NA, 1, 2, 0, 1, 1, 2, NA, 1, 1, 1, 2, 0, NA, NA, NA, NA, 1, 0, 0, 1, 2, 1, 2, 1, 1, 0, NA, NA, 1, 2, NA, NA, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, NA, 1, NA, 0, 1, 0, 1, 1, 1, NA, NA, NA, 0, NA, NA, 1, 0, 1, NA, 1, NA, NA, 2, 0, NA, 1, NA, 1, NA, NA, NA, NA, NA, 0, NA, 0, 0, 1, NA, 0, NA, 0, 0, 1, 1, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2c2"
)

case18 <- scdf(
  phase_design = c(A = 38, B = 70),
  academic_engagement = c(NA, 5, 3, NA, 4, 4, 4, 4, 4, 4, NA, 4, 3, 5, 4, NA, NA, NA, NA, NA, 4, 5, 4, 4, 4, 4, 4, 5, 4, 4, NA, NA, 3, 4, NA, NA, NA, 4, 5, 5, 5, 5, 5, 5, 5, 5, 4, 4, 4, 4, 5, 5, 5, 4, 4, NA, NA, NA, NA, 5, NA, 5, 4, 4, 5, 4, 4, NA, NA, NA, 5, NA, NA, NA, NA, NA, NA, NA, NA, NA, 4, 5, 4, 5, 4, 5, NA, NA, NA, NA, NA, 4, NA, 5, 5, 5, NA, NA, NA, 5, 5, 5, 5, NA, 5, NA, 5, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 1, 1, NA, 1, 1, 0, 1, 1, 1, NA, 0, 2, 0, 1, NA, NA, NA, NA, NA, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0, NA, NA, 2, 0, NA, NA, NA, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, NA, NA, NA, 0, NA, 0, 1, 0, 0, 1, 0, NA, NA, NA, 0, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1, 0, 2, 0, NA, 0, NA, NA, NA, NA, NA, 1, NA, 0, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2c3"
)

case19 <- scdf(
  phase_design = c(A = 38, B = 70),
  academic_engagement = c(2, 2, 3, NA, 3, 3, NA, 2, 2, 2, NA, 3, 3, 3, 3, 4, NA, NA, NA, NA, 4, 4, 3, 4, 4, 4, NA, 4, 3, 3, NA, NA, 3, 4, NA, NA, NA, 2, 3, 3, 5, 4, 4, 5, 3, 5, 3, 2, 4, 3, 3, NA, NA, 4, 4, 4, 3, 4, NA, NA, NA, 5, 4, 4, 5, 5, 4, NA, NA, NA, 3, NA, NA, 3, NA, 4, NA, NA, NA, NA, 4, 4, 4, 4, 4, 5, NA, NA, NA, NA, NA, 5, NA, 5, 5, 5, NA, NA, NA, 5, 4, 5, 4, NA, 4, NA, 5, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 1, 0, NA, 2, 0, NA, 1, 1, 1, NA, 0, 1, 1, 1, 0, NA, NA, NA, NA, 0, 0, 3, 1, 0, 1, NA, 1, 0, 0, NA, NA, 0, 0, NA, NA, NA, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, NA, NA, 0, 1, 0, 0, 0, NA, NA, NA, 0, 1, 0, 0, 1, 0, NA, NA, NA, 0, NA, NA, 1, NA, 0, NA, NA, NA, NA, 1, 0, 1, 1, NA, 0, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, NA, NA, NA, 0, 1, 0, 1, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2c4"
)

case20 <- scdf(
  phase_design = c(A = 38, B = 70),
  academic_engagement = c(3, 3, 4, NA, 4, 3, 2, 3, 4, 4, NA, 3, 4, 3, 4, 3, NA, NA, NA, NA, 3, 4, 3, 4, 4, 4, 4, 5, 3, 4, NA, NA, 4, 4, NA, NA, NA, 2, 4, 5, 5, 4, 5, NA, 5, 5, 4, 3, 5, 5, 5, 5, 4, 4, 5, 5, 4, 5, NA, 5, NA, 5, 5, 4, 5, 5, 5, NA, NA, NA, 3, NA, NA, 4, 4, 5, 5, 4, NA, NA, 4, 5, 5, 5, 5, 5, NA, NA, NA, NA, NA, 5, NA, 5, 5, 5, NA, 5, NA, 5, 5, 5, 5, NA, 5, NA, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c", "2c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(0, 0, 0, NA, 0, 0, 0, 0, 0, 0, NA, 0, 0, 1, 0, 0, NA, NA, NA, NA, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, NA, NA, 0, 0, NA, NA, NA, 0, 0, 0, 0, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, NA, 0, NA, 0, 0, 0, 0, 0, 0, NA, NA, NA, 0, NA, NA, 0, 0, 0, NA, 0, NA, NA, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, NA, 0, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "2c5"
)

case21 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(4, 4, 5, NA, 3, 4, 3, 4, 4, 3, 3, NA, NA, NA, NA, 4, NA, NA, NA, NA, 3, 5, 4, 4, 3, 3, 4, NA, 4, NA, 3, 4, 4, 5, 3, 3, 2, NA, 4, 4, 4, NA, 3, NA, 2, 2, 4, 3, NA, 3, 4, 3, 4, 3, 4, 4, 3, 3, NA, 3, 3, 3, NA, NA, 3, 4, 4, 4, NA, 3, NA, NA, NA, 3, NA, 3, 4, 4, NA, NA, 3, 5, NA, 3, 5, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 4, NA, NA, NA, NA, NA, 4, NA, 4, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(2, 0, 0, NA, 1, 1, 1, 1, 0, 1, 1, NA, NA, NA, NA, 2, NA, NA, NA, NA, 0, 0, 0, 0, 0, 0, 1, NA, 0, NA, 0, 0, 0, 0, 0, 0, 2, NA, 0, 0, 1, NA, 0, NA, 1, 1, 0, 1, NA, 1, 1, 1, 1, 2, 1, 1, 2, 2, NA, 2, 2, 2, NA, NA, 2, 1, 0, 0, NA, 2, NA, NA, NA, 1, NA, 1, NA, 1, NA, NA, 0, NA, NA, 1, 0, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1, NA, NA, NA, NA, NA, 1, NA, 0, 1),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3a1"
)

case22 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(5, 5, 3, NA, 3, 4, 3, 4, 4, NA, 3, 3, 3, 3, 4, 4, NA, NA, NA, NA, 4, 5, 4, 5, 3, 5, 4, NA, 4, 5, 5, 5, 4, 5, 5, 5, 4, NA, 5, 5, 5, NA, 5, NA, 5, 4, 4, 4, 5, 3, 4, 4, 5, 5, 3, 4, 4, 4, NA, 4, 4, 5, NA, NA, 4, 4, 4, 4, NA, 4, NA, NA, NA, 4, NA, 4, 5, 5, NA, NA, 4, 5, NA, 4, 5, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 5, NA, NA, NA, NA, NA, 4, NA, 4, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 0, 1, NA, 0, NA, 2, 1, 1, NA, 2, 3, 1, 2, 2, 2, NA, NA, NA, NA, 0, 0, 0, 0, 1, 0, 1, NA, 0, 0, 0, 2, 0, 0, 0, 0, 1, NA, 0, 0, 1, NA, 0, NA, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 2, 1, 1, 1, NA, 1, 1, 1, NA, NA, 1, 1, 0, 0, NA, 1, NA, NA, NA, 1, NA, 1, NA, 0, NA, NA, 0, NA, NA, 1, 0, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1, NA, NA, NA, NA, NA, 1, NA, 1, 1),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3a2"
)

case23 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(3, 5, 3, NA, NA, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, NA, NA, NA, NA, 3, 5, 3, 4, 3, 3, 4, NA, 4, 3, 4, 4, 4, 5, 2, 3, 2, NA, 4, 5, 4, NA, 3, NA, 2, 2, 4, 2, 3, 3, 4, 3, 5, 5, 4, 4, 4, 4, NA, 4, 4, 4, NA, NA, 5, 4, 4, 4, NA, 4, NA, NA, NA, 3, NA, 3, 4, 5, NA, NA, 3, 4, NA, 3, 5, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 4, NA, 5, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 0, 1, NA, NA, 2, 2, 1, 2, 1, 2, 2, 1, 1, 1, 1, NA, NA, NA, NA, 0, 0, 1, 0, 0, 1, 1, NA, 0, 2, 0, 0, 0, 0, 1, 0, 2, NA, 0, 0, 1, NA, 0, NA, 1, 1, 0, 1, 2, 1, 1, 1, 0, 0, 1, 1, 1, 0, NA, 0, 1, 1, NA, NA, 0, 1, 0, 0, NA, 0, NA, NA, NA, 1, NA, 1, NA, 0, NA, NA, 0, NA, NA, 1, 0, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0, NA, NA, NA, NA, NA, 1, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3a3"
)

case24 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(4, 5, 3, NA, 5, 5, 4, 4, 5, NA, 3, 3, 3, 4, 4, 4, NA, NA, NA, NA, 4, 5, 4, 5, 4, 4, 4, NA, 4, NA, NA, NA, NA, 5, NA, NA, 4, NA, 5, 5, 5, NA, 5, NA, 5, 3, 4, 4, NA, 3, 4, 4, 5, 5, 5, 5, 5, 4, NA, 3, 4, 5, NA, NA, 5, 4, 4, 4, NA, 4, NA, NA, NA, 4, NA, 4, 5, 5, NA, NA, 4, 5, NA, 4, 5, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 5, NA, NA, NA, NA, NA, 4, NA, 4, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 1, 1, NA, 0, 1, 2, 1, 2, NA, 2, 2, 2, 1, 1, 1, NA, NA, NA, NA, 0, 0, 2, 0, 1, 0, 1, NA, 0, NA, NA, NA, NA, 0, NA, NA, 0, NA, 0, 0, 1, NA, 1, NA, 0, 1, 1, 1, NA, 1, 1, 1, 0, 0, 1, 0, 0, 0, NA, 2, 1, 1, NA, NA, 0, 1, 0, 0, NA, 2, NA, NA, NA, 1, NA, 1, NA, 0, NA, NA, 0, NA, NA, 1, 0, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1, NA, NA, NA, NA, NA, 1, NA, 1, 1),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3a4"
)

case25 <- scdf(
  phase_design = c(A = 18, B = 90),
  academic_engagement = c(3, 3, 3, NA, 3, 4, 3, 4, 4, 4, 3, 3, 3, 4, 1, 4, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 5, 3, NA, 4, NA, NA, 3, 4, 3, NA, 3, 4, 4, 5, 5, 5, 5, 4, 4, NA, 4, 4, 5, NA, NA, 4, 4, 4, 4, NA, 4, NA, NA, NA, 4, NA, 4, 5, 4, NA, NA, 4, 5, NA, 4, 5, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 5, NA, NA, NA, NA, NA, 4, NA, 4, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a", "3a"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(2, 2, 2, NA, 2, 2, 1, 1, 2, 2, 3, 3, 2, 1, 3, 2, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0, 2, NA, 1, NA, NA, 1, 1, 2, NA, 1, 1, 1, 0, 0, 1, 0, 1, 1, NA, 1, 1, 1, NA, NA, 1, 1, 0, 0, NA, 0, NA, NA, NA, 1, NA, 1, NA, 0, NA, NA, 0, NA, NA, 1, 0, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1, NA, NA, NA, NA, NA, 1, NA, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3a5"
)

case26 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(4, 4, 5, NA, NA, NA, NA, NA, 3, 3, 3, 2, 4, 3, NA, 4, NA, NA, NA, NA, 3, 4, 3, 3, 2, 3, NA, NA, 2, 3, 4, 3, NA, 3, 3, 4, 4, 3, 3, 4, NA, 2, NA, 3, 4, 3, NA, NA, 3, 4, NA, 3, 3, NA, 3, 4, 2, 4, NA, 2, NA, NA, NA, NA, 4, 4, 4, 3, NA, 3, 4, NA, NA, 4, NA, 3, 3, 3, NA, NA, 3, NA, 3, 3, 4, NA, NA, NA, 3, 4, NA, 3, 3, 4, 4, NA, 4, NA, 3, 4, 4, NA, NA, NA, NA, 3, 2, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(1, 0, 1, NA, NA, NA, NA, NA, 2, 3, 2, 4, 1, 2, NA, 1, NA, NA, NA, NA, 2, 2, 3, 2, 4, 2, NA, NA, 2, 2, 1, 2, NA, 1, 0, 0, 1, 2, 2, 1, NA, 3, NA, 2, 1, 1, NA, NA, 2, 1, NA, 2, 2, NA, 2, 0, 3, 1, NA, 3, NA, NA, NA, NA, 1, 2, 1, 3, NA, 2, 1, NA, NA, 1, NA, 2, 3, 2, NA, NA, 2, NA, 2, 2, 1, NA, NA, NA, 2, 1, NA, 2, 2, 2, 1, NA, 1, NA, 2, 2, 1, NA, NA, NA, NA, 3, 3, 2),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3b1"
)

case27 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(4, 2, 4, NA, NA, NA, NA, NA, 4, 3, 4, 3, 4, 4, NA, 4, NA, NA, NA, NA, 3, 4, 4, 3, 3, 3, NA, NA, 4, 4, 4, 4, NA, 3, 2, 3, 3, 3, 2, 3, NA, 3, NA, NA, NA, 4, NA, 1, 4, 3, NA, 4, 3, NA, 3, 4, 3, 3, NA, 4, NA, NA, NA, NA, 4, 4, 3, 4, 4, 3, 4, NA, NA, 4, 4, 4, 3, 4, NA, NA, 4, 4, 3, 3, 4, NA, NA, NA, 3, 4, NA, 2, 4, 4, 3, NA, 4, NA, 3, 2, 3, NA, NA, NA, 4, 3, 4, 3),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(4, 3, 2, NA, NA, NA, NA, NA, 1, 3, 2, 3, 2, 1, NA, 2, NA, NA, NA, NA, 2, 1, 2, 3, 2, 0, NA, NA, 1, 1, 2, 2, NA, 1, 2, 2, 2, 2, 3, 3, NA, 2, NA, NA, NA, 1, NA, 4, 0, 2, NA, 2, 3, NA, 3, 1, 2, 2, NA, 1, NA, NA, NA, NA, 0, 1, 2, 1, 2, 2, 1, NA, NA, 1, 1, 1, 4, 2, NA, NA, 1, 2, 2, 3, 1, NA, NA, NA, 2, 1, NA, 4, 1, 2, 2, NA, 0, NA, 2, 3, 2, NA, NA, NA, 2, 3, 2, 2),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3b2"
)

case28 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(2, 2, 3, NA, NA, NA, NA, NA, 2, 2, NA, 1, 3, 2, NA, 3, NA, NA, NA, NA, 2, 2, 2, 2, 2, 2, NA, NA, 1, NA, NA, 3, NA, 2, 2, 3, 4, 3, 3, 4, NA, 3, NA, 2, 3, 2, NA, NA, 1, 2, NA, 2, 3, NA, 2, 3, NA, 3, NA, 4, NA, NA, NA, NA, 3, 2, 4, 1, 3, 2, 4, NA, NA, 3, 3, 3, 3, 3, NA, NA, 5, NA, 3, 2, 2, NA, NA, NA, 2, 3, NA, 3, 2, 2, 3, NA, 4, NA, 2, 2, 4, NA, NA, NA, 3, 3, 3, 3),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(3, 4, 3, NA, NA, NA, NA, NA, 4, 3, NA, 3, 1, 3, NA, 3, NA, NA, NA, NA, 4, 3, 4, 3, 3, 3, NA, NA, 4, NA, NA, 3, NA, 2, 2, 2, 1, 1, 2, 2, NA, 1, NA, 2, 2, 2, NA, NA, 0, 3, NA, 3, 3, NA, 3, 1, NA, 2, NA, 1, NA, NA, NA, NA, 1, 3, 1, 3, 2, 3, 0, NA, NA, 1, 1, 2, 2, 2, NA, NA, 0, NA, 1, 3, 3, NA, NA, NA, 3, 2, NA, 2, 3, 3, 2, NA, 1, NA, 3, 3, 1, NA, NA, NA, 2, 3, 2, 2),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3b3"
)

case29 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(3, 3, 2, NA, NA, NA, NA, NA, NA, 2, 2, 2, 3, NA, NA, 3, NA, NA, NA, NA, 2, 3, 3, 2, 2, 1, NA, NA, 3, NA, NA, 2, NA, 2, 3, NA, 3, 4, 3, 3, NA, 3, NA, 2, 1, 2, NA, NA, 1, 2, NA, NA, 2, NA, 2, 4, NA, 3, NA, 2, NA, NA, NA, NA, 3, 3, NA, 2, 4, 3, 4, NA, NA, 3, 4, 3, 4, 3, NA, NA, 3, 2, 4, 2, 3, NA, NA, NA, 3, 3, NA, 3, 3, NA, 3, NA, 4, NA, 3, 4, 5, NA, NA, NA, 5, 4, 3, 4),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(3, 3, 3, NA, NA, NA, NA, NA, NA, 3, 3, 4, 2, NA, NA, 1, NA, NA, NA, NA, 4, 2, 2, 3, 3, 3, NA, NA, 2, NA, NA, 3, NA, 2, 2, NA, 2, 1, 2, 2, NA, 1, NA, 3, 3, 2, NA, NA, 3, 3, NA, NA, 3, NA, 2, 1, NA, 1, NA, 2, NA, NA, NA, NA, 1, 1, NA, 2, 1, 2, 0, NA, NA, 1, 0, 2, 1, 2, NA, NA, 3, 3, 1, 2, 1, NA, NA, NA, 1, 2, NA, 2, 2, NA, 1, NA, 0, NA, 1, 0, 0, NA, NA, NA, 0, 1, 2, 1),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3b4"
)

case30 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, 3, 3, 3, 4, 3, NA, 4, NA, NA, NA, NA, 3, 4, 3, 3, 2, 3, NA, NA, 3, 3, 3, 3, NA, 3, 3, 3, 4, 4, NA, 4, NA, 4, NA, 4, 4, 3, NA, 5, 4, 4, NA, 5, 4, NA, 4, 4, 4, 5, NA, NA, NA, NA, NA, NA, 4, 4, NA, NA, 3, 4, NA, NA, NA, NA, 4, 3, 5, 4, NA, NA, 4, 3, 4, 3, 4, NA, NA, NA, 3, 3, NA, 4, 4, 5, 3, NA, 4, NA, 4, 4, 5, NA, NA, NA, 5, 5, 4, 5),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b", "3b"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, 3, 2, 3, 1, 1, NA, 2, NA, NA, NA, NA, 1, 1, 3, 3, 3, 2, NA, NA, 1, 2, 1, 2, NA, 0, 1, 1, 1, 0, NA, 1, NA, 1, NA, 1, 0, 0, NA, 0, 0, 0, NA, 0, 0, NA, 1, 0, 1, 0, NA, NA, NA, NA, NA, NA, 0, 1, NA, NA, 2, 0, NA, NA, NA, NA, 0, 1, 1, 0, NA, NA, 0, 2, 0, 1, 1, NA, NA, NA, 1, 2, NA, 1, 1, 1, 2, NA, 0, NA, 0, 0, 0, NA, NA, NA, 0, 1, 0, 0),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3b5"
)

case31 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(NA, 1, 3, NA, 3, NA, 3, 1, 1, 4, NA, 4, 5, 3, NA, NA, NA, NA, NA, NA, NA, 1, 4, 1, 3, NA, 3, 3, 3, 4, NA, 5, 5, 5, 5, NA, 5, 4, 5, NA, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, 5, 3, NA, 5, NA, NA, NA, NA, NA, NA, NA, 5, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, NA, NA, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, 5, 5, NA, 2, 5, NA, NA, NA, NA, NA),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 2, 0, NA, 0, NA, 1, 3, 1, 0, NA, 0, 0, 0, NA, NA, NA, NA, NA, NA, NA, 4, 0, 4, 1, NA, 1, 2, 2, 0, NA, 0, 0, 0, 0, NA, 0, 0, 0, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0, NA, 2, 0, NA, NA, NA, NA, NA),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3c1"
)

case32 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(NA, 2, 4, NA, 4, NA, 2, 3, 1, 5, NA, 5, 3, 4, NA, NA, NA, NA, NA, NA, NA, 4, 4, 2, 3, NA, 2, 3, 4, 5, NA, 5, 4, 5, 5, NA, 5, 5, 5, NA, NA, 4, NA, 5, 5, NA, NA, 5, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, NA, NA, NA, NA, NA, NA, 5, NA, 4, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, NA, NA, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, 5, 5, NA, 5, 5, NA, NA, NA, NA, NA),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 1, 1, NA, 1, NA, 2, 1, 2, 0, NA, 0, 1, 1, NA, NA, NA, NA, NA, NA, NA, 0, 1, 2, 1, NA, 1, 2, 1, 0, NA, 0, 1, 0, 0, NA, 0, 0, 0, NA, NA, 0, NA, 0, 0, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0, NA, 0, 0, NA, NA, NA, NA, NA),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3c2"
)

case33 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(NA, 1, 2, NA, 3, NA, 1, 4, 2, 3, NA, 3, 2, 2, NA, NA, NA, NA, NA, NA, NA, 1, 1, 1, 3, NA, 0, 2, 4, 5, NA, 4, 5, 5, 5, NA, 5, 5, 5, NA, NA, NA, 5, 5, 5, NA, 5, 5, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, NA, NA, NA, NA, NA, NA, 5, NA, 5, 5, 5, 4, NA, NA, NA, 4, 5, NA, 5, 5, NA, NA, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, 5, 5, NA, 3, 5, NA, NA, NA, NA, NA),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 3, 0, NA, 1, NA, 2, 1, 1, 2, NA, 1, 0, 0, NA, NA, NA, NA, NA, NA, NA, 4, 3, 3, 1, NA, 2, 3, NA, 1, NA, 0, 0, 0, 0, NA, 0, 0, 1, NA, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 1, NA, NA, NA, 1, 0, NA, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0, NA, 1, 0, NA, NA, NA, NA, NA),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3c3"
)

case34 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(NA, 4, 4, NA, 4, NA, 4, 4, 5, 5, NA, 5, 5, 4, NA, NA, NA, NA, NA, NA, NA, 4, 5, 4, 4, NA, 3, 4, 5, 5, NA, 4, 5, 5, 5, NA, 5, NA, 3, NA, NA, 5, 5, 5, 5, NA, NA, 5, 5, 5, NA, 5, 4, 5, 5, NA, 5, NA, NA, NA, NA, NA, NA, NA, 5, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, NA, NA, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, 5, 5, NA, 5, 5, NA, NA, NA, NA, NA),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 1, 1, NA, 1, NA, 1, 1, 0, 0, NA, 0, 0, 1, NA, NA, NA, NA, NA, NA, NA, 0, 0, 0, 1, NA, 2, 1, 0, 1, NA, 1, 0, 0, 0, NA, 0, NA, 0, NA, NA, 0, 0, 0, 0, NA, NA, 0, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0, NA, 0, 0, NA, NA, NA, NA, NA),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3c4"
)

case35 <- scdf(
  phase_design = c(A = 28, B = 80),
  academic_engagement = c(NA, 2, 5, NA, 4, NA, 4, 5, 4, 4, NA, 5, 5, 4, NA, NA, NA, NA, NA, NA, NA, 4, 5, 4, 4, NA, 3, 4, 5, 5, NA, 3, 5, 5, 5, NA, 5, 5, 2, NA, NA, 5, 5, 5, 5, NA, NA, 5, NA, 5, NA, 5, 5, 5, 5, NA, 5, NA, NA, NA, NA, NA, NA, NA, 5, NA, 5, 5, 5, 3, NA, NA, NA, 5, 5, NA, 5, 5, NA, NA, NA, 5, 5, 5, 5, NA, NA, NA, 5, 5, NA, 5, 5, 5, 5, NA, 5, NA, 5, 5, NA, 4, 5, NA, NA, NA, NA, NA),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108),
  classID = c("3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c", "3c"),
  weekday = c(3, 4, 5, 1, 2, 3, 4, 4, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5),
  disruptive_behavior = c(NA, 2, 0, NA, 0, NA, 1, 0, 1, 1, NA, 0, 0, 1, NA, NA, NA, NA, NA, NA, NA, 0, 0, 0, 0, NA, 1, 0, 0, 0, NA, 1, 0, 0, 0, NA, 0, 0, 2, NA, NA, 0, 0, 0, 0, NA, NA, 0, NA, 0, NA, 0, 0, 0, 0, NA, 0, NA, NA, NA, NA, NA, NA, NA, 0, NA, 0, 0, 0, 2, NA, NA, NA, 0, 0, NA, 0, 0, NA, NA, NA, 0, 0, 0, 0, NA, NA, NA, 0, 0, NA, 0, 0, 0, 0, NA, 0, NA, 0, 0, NA, 1, 0, NA, NA, NA, NA, NA),
  dvar = "academic_engagement",
  mvar = "mt",
  name = "3c5"
)

Leidig2018 <- c(case1, case2, case3, case4, case5, case6, case7, case8, case9, case10, case11, case12, case13, case14, case15, case16, case17, case18, case19, case20, case21, case22, case23, case24, case25, case26, case27, case28, case29, case30, case31, case32, case33, case34, case35)
#scdf_attr(Leidig2018, "info") <- NULL
scdf_attr(Leidig2018, "author") <- "Leidig, T., Casale, G., Wilbert, J., Hennemann, T., Volpe, R. J., Briesch, A., & Grosche, M. (2022). Individual, generalized, and moderated effects of the good behavior game on at-risk primary school students: A multilevel multiple baseline study using behavioral progress monitoring. Frontiers in Education, 7. https://www.frontiersin.org/articles/10.3389/feduc.2022.917138"
filenames <- c(filenames, "Leidig2018")

Leidig2018_l2 <- data.frame(
  class = c("1a", "1a", "1a", "1a", "1a", "2a", "2a", "2a", "2a", "2a", "2b", "2b", "2b", "2b", "2b", "2c", "2c", "2c", "2c", "2c", "3a", "3a", "3a", "3a", "3a", "3b", "3b", "3b", "3b", "3b", "3c", "3c", "3c", "3c", "3c"),
  case = c("1a1", "1a2", "1a3", "1a4", "1a5", "2a1", "2a2", "2a3", "2a4", "2a5", "2b1", "2b2", "2b3", "2b4", "2b5", "2c1", "2c2", "2c3", "2c4", "2c5", "3a1", "3a2", "3a3", "3a4", "3a5", "3b1", "3b2", "3b3", "3b4", "3b5", "3c1", "3c2", "3c3", "3c4", "3c5"),
  gender = c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0),
  migration = c(0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0),
  first_language_german = c(1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1),
  SDQ_TOTAL = c(10, 11, 6, 8, 8, 9, 7, 16, 21, 22, 13, 14, 16, 14, 9, 4, 9, 13, 11, 16, 4, 6, 13, 7, 15, 15, 13, 10, 11, 19, 8, 9, 12, 12, 10),
  SDQ_EXTERNALIZING = c(9, 11, 6, 5, 7, 8, 4, 13, 15, 18, 10, 13, 13, 10, 8, 4, 6, 9, 4, 6, 3, 6, 8, 6, 12, 13, 11, 7, 10, 10, 6, 9, 9, 11, 10),
  SDQ_INTERNALIZING = c(1, 0, 0, 3, 1, 1, 3, 3, 6, 4, 3, 1, 3, 4, 1, 0, 3, 4, 7, 10, 1, 0, 5, 1, 3, 2, 2, 3, 1, 9, 2, 0, 3, 1, 0),
  ITRF_TOTAL = c(11, 21, 18, 14, 13, 9, 16, 23, 32, 33, 31, 15, 33, 16, 6, 4, 6, 16, 15, 23, 5, 12, 10, 3, 15, 35, 27, 23, 23, 21, 14, 15, 15, 17, 19),
  ITRF_ACADEMIC = c(7, 10, 1, 7, 4, 2, 16, 14, 16, 13, 24, 8, 20, 11, 2, 1, 0, 2, 15, 17, 0, 12, 0, 0, 8, 22, 17, 16, 6, 17, 1, 4, 9, 13, 10),
  ITRF_BEHAVIOR = c(4, 11, 17, 7, 9, 7, 0, 9, 16, 20, 7, 7, 13, 5, 4, 3, 6, 14, 0, 6, 5, 0, 10, 3, 7, 13, 10, 7, 17, 4, 13, 11, 6, 4, 9)
)

filenames <- c(filenames, "Leidig2018_l2")

# Example from Beretvas, S., & Chung, H. (2008). An evaluation of  --------

Beretvas2008 <- scdf(c(0.7, 1.60, 1.40, 1.60, 1.90, 1.20, 1.30, 1.60, 10.00, 10.80, 11.90, 11.00, 13.00, 12.70, 14.00), B_start = 9)
scdf_attr(Beretvas2008, "author") <- "Beretvas, S., & Chung, H. (2008). An evaluation of modified R2-change effect size indices for single-subject experimental designs. Evidence-Based Communication Assessment and Intervention, 2, 120-128."

filenames <- c(filenames, "Beretvas2008")


# Example from Huitema et al., 2000 ---------------------------------------

Huitema2000 <- scdf(score = c(58, 70, 79, 85, 93, 94, 105, 110, 122, 130, 122, 110, 103, 88, 81, 68, 53, 50, 36, 31), B_start = 11, dvar = "score")
scdf_attr(Huitema2000, "info") <- "Example from Huitema, B. E., & Mckean, J. W. (2000). Design specification issues in time-series intervention models. Educational and Psychological Measurement, 60(1), 38-58."

filenames <- c(filenames, "Huitema2000")

# Example from Waddell et al., 2011 ---------------------------------------

Waddell2011 <- scdf(score = c(5.2, 5, 6, 6.1, 4.8, 5.3, 7, 6.2, 6.1, 7.2, 6.2, 7, 8.4, 11, 10, 11), B_start = 6, dvar = "score")
scdf_attr(Waddell2011, "info") <- "Example from Waddell, D. E., Nassar, S. L., & Gustafson, S. A. (2011). Single-Case Design in Psychophysiological Research: Part II: Statistical Analytic Approaches. Journal of Neurotherapy, 15(2), 160 - 169."

filenames <- c(filenames, "Waddell2011")

# Parker, R. I., Vannest, K. J., Davis, J. L., & Sauber, S. B. (2011) --------

Parker2011 <- scdf(c(2, 3, 5, 3, 4, 5, 5, 7, 6), phase_design = c(A = 4, B = 5))
scdf_attr(Parker2011, "author") <- "Parker, R. I., Vannest, K. J., Davis, J. L., & Sauber, S. B. (2011). Combining Nonoverlap and Trend for Single-Case Research: Tau-U. Behavior Therapy, 42(2), 284-299."

filenames <- c(filenames, "Parker2011")

# Parker 2009 ------
Parker2009 <- c(
  scdf(c(A1 = 0,0,0,0,0,0,0, 
         B1 = 0, 50, 30, 0, 0,40, 50, 70, 40, 60, 60, 70, 70, 60, 60, 60, 
         A2 = 30, 50, 0, 0, 
         B2 = 30, 80, 70, 80, 80, 100, 70, 100, 90,70, 70, 80)),
  author = "Parker, R. I., Vannest, K. J., & Brown, L. (2009). The improvement rate difference for single-case research. Exceptional Children, 75(2), 135-150."
)

filenames <- c(filenames, "Parker2009")

Parker2009b <- c(
  scdf(c(A=4, 3, 4, 3, 4, 7, 5, 2, 3, 2, B=5, 9, 7, 9, 7, 5, 9, 11, 11, 10, 9), name = "Parker2009b"),
  author = "Parker, R. I., & Vannest, K. (2009). An improved effect size for single-case research: Nonoverlap of all pairs. Behavior Therapy, 40(4), 357-367."
)

filenames <- c(filenames, "Parker2009b")

# Parker 2007 -------

Parker2007 <- c(
  scdf(c(A = 20, 18, 20, B = 20, 21, 23), name = "Adam"),
  scdf(c(A = 18, 17, 16, 20, B = 19, 22, 19, 20, 20, 24), name = "Bob"),
  scdf(c(A = 19, 18, 24, 22, 21, 19, B = 30, 21, 23, 28, 32, 34), name = "Carol"),
  author = "Parker, R. I., Hagan-Burke, S., & Vannest, K. (2007). Percentage of All Non-Overlapping Data (PAND) An Alternative to PND. The Journal of Special Education, 40(4), 194-204.
"
)

filenames <- c(filenames, "Parker2007")


# Borckardt, J. J., & Nash, M. R. (2014). Simulation modelling ana --------

Borckardt2014 <- scdf(c(7, 9, 8, 7, 7, 8, 8, 3, 4, 2, 4, 2, 1, 4, 5, 3, 2, 1, 2, 4, 4), B_start = 8)
scdf_attr(Borckardt2014, "author") <- "Borckardt, J. J., & Nash, M. R. (2014). Simulation modelling analysis for small sets of single-subject data collected over time. Neuropsychological Rehabilitation, 24(3-4), 492-506."

filenames <- c(filenames, "Borckardt2014")

# Lenz, A. S. (2013). Calculating Effect Size in Single-Case Resea --------

Lenz2013 <- scdf(c(22, 25, 23, 26, 20, 22, 25, 26, 25, 22, 22, 18, 21, 18, 20, 21, 18, 16, 16, 14, 14, 12), phase_design = c("A" = 7, "B" = 15))
scdf_attr(Lenz2013, "author") <- "Lenz, A. S. (2013). Calculating Effect Size in Single-Case Research: A Comparison of Nonoverlap Methods. Measurement and Evaluation in Counseling and Development, 46(1), 64-73."

filenames <- c(filenames, "Lenz2013")

# example SSD package -----------------------------------------------------

SSDforR2017 <- scdf(
  esteem = c(3, 4, 2, 5, 3, 4, 2, 2, 3, 2, 1, 2, 2, 2, 1, 2, 1, 0, 0, 0),
  phase_design = c(A = 6, B = 6, B1 = 8), dvar = "esteem", pvar = "pesteem"
)

scdf_attr(SSDforR2017, "info") <- "Example from the SSDforR package."
scdf_attr(SSDforR2017, "author") <- "Charles Auerbach, PhD & Wendy Zeitlin, PhD; Yeshiva University, Wurzweiler school of social work."

filenames <- c(filenames, "SSDforR2017")

## examples from Tarlow 2017

Tarlow2017 <- c(
  scdf(c(A = 9, 6, 11, 5, B = 4, 7, 8, 7, 3, 7, 1), name = "example_a"),
  scdf(c(A = 1, 2, 3, 4, 5, B = 2, 1, 0, 0, 0), name = "example_b"),
  scdf(c(A = 33, 25, 17, 25, 14, 13,14, B = 14, 15, 15, 4, 6, 9, 5 ,4 ,2 ,2 ,8, 11 ,7), name = "example_c"),
  author = "Tarlow, K. R. (2017). An Improved Rank Correlation Effect Size Statistic for Single-Case Designs: Baseline Corrected Tau. Behavior Modification, 41(4), 427-467. https://doi.org/10.1177/0145445516676750"
)

filenames <- c(filenames,"Tarlow2017")

# Parker 2011b

Parker2011b <- c(
scdf(c(A = 20, 20, 26, 25, 22, 23, B = 28, 25, 24, 27, 30, 30, 29)),
author = "Parker, R. I., Vannest, K. J., & Davis, J. L. (2011). Effect Size in Single-Case Research: A Review of Nine Nonoverlap Techniques. Behavior Modification, 35(4), 303-322. https://doi.org/10.1177/0145445511399147"
)

filenames <- c(filenames, "Parker2011b")

# exampleABAB-Design --------------------------------------------------

case1 <- scdf(
  phase_design = c(A = 10, B = 10, A = 10, B = 10),
  points = c(7, 14, 14, 14, 18, 15, 5, 8, 15, 13, 34, 31, 26, 38, 26, 37, 26, 
             25, 24, 25, 10, 10, 7, 6, 8, 17, 13, 19, 17, 6, 29, 25, 26, 30, 
             24, 35, 24, 25, 40, 36),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20,
         21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 
         38, 39, 40),
  dvar = "points",
  mvar = "mt",
  name = "Howard"
)

case2 <- scdf(
  phase_design = c(A = 15, B = 8, A = 7, B = 10),
  points = c(13, 15, 10, 14, 10, 13, 15, 12, 9, 16, 6, 10, 16, 13, 7, 30, 30, 
             35, 23, 37, 37, 21, 26, 5, 9, 16, 10, 13, 6, 13, 7, 18, 5, 17, 6, 
             13, 11, 6, 10, 15),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 
         21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 
         38, 39, 40),
  dvar = "points",
  mvar = "mt",
  name = "Sheldon"
)

case3 <- scdf(
  phase_design = c(A = 20, B = 7, A = 6, B = 7),
  points = c(19, 12, 7, 13, 8, 18, 11, 10, 7, 18, 7, 19, 7, 7, 7, 13, 10, 5, 
             10, 16, 21, 31, 26, 24, 23, 27, 23, 7, 12, 6, 16, 7, 19, 22, 24, 
             38, 39, 26, 22, 36),
  mt = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20,
         21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 
         38, 39, 40),
  dvar = "points",
  mvar = "mt",
  name = "Leonard"
)

exampleABAB <- c(case1, case2, case3)
scdf_attr(exampleABAB, "info") <- "Randomly created data with uniform distribution."
filenames <- c(filenames, "exampleABAB")



### exampleA1B1A2B2-Design---------
set.seed(4321)
dat <- round(
  c(runif(10, 5, 20), runif(10, 15, 40), runif(10, 8, 25), runif(10, 20, 35))
)
A <- scdf(dat, phase_design = c(A1 = 10, B1 = 10, A2 = 10, B2 = 10))
dat <- round(
  c(runif(15, 5, 20), runif(8, 15, 40), runif(7, 8, 25), runif(10, 20, 35))
)
B <- scdf(dat, phase_design = c(A1 = 15, B1 = 8, A2 = 7, B2 = 10))
dat <- round(
  c(runif(20, 5, 20), runif(7, 15, 40), runif(6, 8, 25), runif(7, 20, 35))
)
C <- scdf(dat, phase_design = c(A1 = 20, B1 = 7, A2 = 6, B2 = 7))
exampleA1B1A2B2 <- c(A, B, C)
names(exampleA1B1A2B2) <- c("Pawel", "Moritz", "Jannis")

filenames <- c(filenames, "exampleA1B1A2B2")

### example of ABC-Design--------------
set.seed(421)
A <- scdf(
  round(c(rnorm(10, 50, 10), rnorm(20, 70, 10))), 
  phase_design = c(A = 10, B = 10, C = 10), name = "Heintje"
)
B <- scdf(
  round(c(rnorm(15, 50, 10), rnorm(15, 70, 10))), 
  phase_design = c(A = 15, B = 8, C = 7), name = "Heino"
)
C <- scdf(
  round(c(rnorm(20, 50, 10), rnorm(10, 70, 10))), 
  phase_design = c(A = 20, B = 7, C = 3), name = "Hannelore"
)
exampleABC <- c(A, B, C)
names(exampleABC) <- c("Marie", "Rosalind", "Lise")
filenames <- c(filenames, "exampleABC")

# example of AB-Design -----------------

case1 <- scdf(
  values = c(54, 53, 56, 58, 52, 61, 62, 71, 66, 64, 78, 70, 74, 82, 77, 86, 
             68, 80, 86, 87), 
  phase_design = c(A = 5, B = 15),
  name = "Johanna"
)

case2 <- scdf(
  values = c(41, 59, 56, 51, 52, 57, 56, 67, 75, 66, 69, 68, 73, 77, 79, 86, 
             82, 75, 83, 89), 
  phase_design = c(A = 5, B = 15),
  name = "Karolina"
)

case3 <- scdf(
  values = c(55, 58, 53, 50, 52, 55, 68, 68, 81, 67, 78, 73, 72, 78, 81, 78, 
             71, 85, 80, 76), 
  phase_design = c(A = 5, B = 15),
  name = "Anja"
)

exampleAB <- c(case1, case2, case3)

scdf_attr(exampleAB, "info") <- "Randomly created data with normal distributed dependent variable."

filenames <- c(filenames, "exampleAB")

# example of AB-Design decreasing -------------------
case1 <- scdf(
  values = c(108, 105, 109, 112, 100, 94, 105, 97, 73, 73, 82, 93, 81, 93, 59, 
             74, 83, 61, 84, 78), 
  phase_design = c(A = 7, B = 13),
  name = "Peter"
)

case2 <- scdf(
  values = c(92, 118, 117, 105, 109, 102, 113, 102, 83, 74, 86, 69, 78, 89, 93,
             83, 75, 77, 69, 79), 
  phase_design = c(A = 8, B = 12),
  name = "Tony"
)

case3 <- scdf(
  values = c(126, 94, 90, 95, 106, 96, 67, 84, 79, 82, 76, 80, 91, 74, 65, 73,
             89, 74, 62, 63), 
  phase_design = c(A = 6, B = 14),
  name = "Bruce"
)

exampleAB_decreasing <- c(case1, case2, case3)

scdf_attr(exampleAB_decreasing, "info") <- "Random data-set from a poisson distribution. Level effect is negative."

filenames <- c(filenames, "exampleAB_decreasing")

# example for AB multilevel analyzes -----
set.seed(123456)
n <- 25

design <- design(
  n = n, rtt = 0.8,
  phase_design = list(A = round(runif(n, 5, 10)), B = round(runif(n, 15, 25))),
  level = list(A = 0, B = 1),
  slope = list(A = 0, B = 0.1),
  trend = 0.05,
  start_value = round(rnorm(n, 50, 10)),
  s = 10
)

male <- random_scdf(design = design, random_names = "male")

design <- design(
  n = n, rtt = 0.8,
  phase_design = list(A = round(runif(n, 5, 10)), B = round(runif(n, 15, 25))),
  level = list(A = 0, B = 1.7),
  slope = list(A = 0, B = 0.1),
  trend = list(0.05),
  start_value = round(rnorm(n, 50, 10)),
  s = 10
)

female <- random_scdf(design = design, random_names = "female")

exampleAB_50 <- c(male, female)

exampleAB_50.l2 <- data.frame(
  case = names(exampleAB_50),
  sex = factor(c(rep("m", n), rep("f", n))),
  age = round(runif(n * 2, 9, 13), 0)
)

filenames <- c(filenames, "exampleAB_50", "exampleAB_50.l2")

# example for ABC multilevel analyzes -----
set.seed(123456)
n <- 150

design <- design(
  n = n, rtt = 0.9,
  phase_design = list(A = round(runif(n, 5, 10)), 
                      B = round(runif(n, 15, 25)), 
                      C = round(runif(n, 5, 10))),
  level = list(A = 0, B = rnorm(n, 0.5, 0.5), C = 0),
  slope = list(A = 0, B = rnorm(n, 0.1, 0.1), C = 0),
  trend = 0.05,
  start_value = round(rnorm(n, 50, 10)),
  s = 10
)

exampleABC_150 <- random_scdf(design = design, random_names = TRUE, round = 0)

scdf_attr(exampleABC_150, "info") <- "Random data-set for testing out hplm. Level and slope effects vary."

filenames <- c(filenames, "exampleABC_150")


### exampleAB_add for additional variables ----------------
case1 <- scdf(
  phase_design = c(Base = 15, Medication = 25),
  day = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 
          20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 
          37, 38, 39, 40),
  wellbeing = c(46, 49, 49, 49, 50, 47, 45, 59, 58, 59, 59, 43, 46, 52, 55, 53, 
                63, 65, 55, 62, 62, 64, 47, 54, 60, 64, 57, 55, 57, 57, 56, 63, 
                55, 53, 58, 57, 65, 61, 63, 67),
  cigarrets = c(2, 5, 4, 1, 2, 4, 4, 0, 2, 3, 2, 1, 4, 5, 5, 5, 4, 1, 3, 1, 4, 
                1, 1, 1, 4, 3, 5, 4, 3, 5, 0, 0, 4, 1, 0, 1, 2, 5, 2, 4),
  depression = c(7, 6, 1, 4, 7, 2, 6, 0, 3, 6, 8, 7, 9, 6, 3, 5, 10, 10, 2, 5, 
                 1, 3, 5, 3, 1, 1, 9, 2, 10, 3, 0, 0, 9, 2, 5, 7, 3, 3, 3, 7),
  dvar = "wellbeing",
  mvar = "day",
  name = "Rolf"
)

exampleAB_add <- c(case1)
scdf_attr(exampleAB_add, "info") <- "Random data-set for testing out plm with additional variables."

filenames <- c(filenames, "exampleAB_add")

### exampleABC_outlier  ----------
exampleABC_outlier <- exampleABC
exampleABC_outlier[[1]]$values[c(3, 10, 12)] <- 120
exampleABC_outlier[[2]]$values[c(7, 12, 18)] <- 120
exampleABC_outlier[[3]]$values[c(9, 17, 20)] <- 120
names(exampleABC_outlier) <- c("Bernadette", "Penny", "Amy")
scdf_attr(exampleABC_outlier, "info") <- "Random data-set based on exampleABC but with outliers."

filenames <- c(filenames, "exampleABC_outlier")

### example binomial  ----------

case1 <- scdf(
  values = c(1, 3, 3, 3, 5, 3, 0, 2, 4, 3, 12, 13, 15, 11, 15, 12, 16, 16, 17, 
             16, 16, 16, 17, 19, 17, 15, 16, 14, 15, 19),
  trials = c(20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 
             20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20), 
  phase_design = c(A = 10, B = 20),
  name = "Christiano"
)

case2 <- scdf(
  values = c(3, 2, 2, 3, 2, 4, 2, 2, 7, 4, 13, 13, 15, 13, 15, 14, 14, 15, 16, 
             14, 18, 16, 15, 16, 18, 16, 17, 16, 18, 15),
  trials = c(20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 
             20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20), 
  phase_design = c(A = 10, B = 20),
  name = "Lionel"
)

case3 <- scdf(
  values = c(5, 1, 2, 0, 2, 4, 2, 3, 1, 3, 16, 11, 18, 12, 17, 14, 15, 17, 16, 
             14, 12, 16, 18, 16, 18, 14, 17, 18, 18, 15),
  trials = c(20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 
             20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20), 
  phase_design = c(A = 10, B = 20),
  name = "Neymar"
)

exampleAB_score <- c(case1, case2, case3)

scdf_attr(exampleAB_score, "info") <- "Random data-set for binomial data."
filenames <- c(filenames, "exampleAB_score")


# write data -----

.createDataFileScan <- function(filenames) {
  library(usethis)
  for(i in 1:length(filenames)) {
    
    do.call("use_data", list(as.name(filenames[i]), overwrite = TRUE))
    #usethis::use_data(str2lang(filenames[i]), overwrite = TRUE)    
  }

  #save(
  #  list = filenames, 
  #  file = file.path("..", "data", "scan-Data.RData"), 
  #  version = 2
  #)
  rm(list = filenames, envir = .GlobalEnv)
}

.createDataFileScan(filenames)

rm(list = "design", envir = .GlobalEnv)
rm(list = paste0("case", 1:35), envir = .GlobalEnv)
rm(list = c("A", "B", "C", "dat", "n", "filenames"), envir = .GlobalEnv)
rm(list = c("male", "female"), envir = .GlobalEnv)