Nothing
context("Preprocess_SCD function")
# MB design
MB_dat <- data.frame(
Student = rep(LETTERS[1:4], each = 10),
Phase = rep(rep(c("Base","Trt"), each = 5), 4),
Time = rep(1:10, 4),
Behavior = rnorm(40)
)
dat_dup_session <- MB_dat
dat_dup_session$Time <- rep(1:8, 5)
MB_dat_clean <- preprocess_SCD(case = Student, phase = Phase,
session = Time, outcome = Behavior,
design = "MBP", data = MB_dat)
test_that("The returned dataset is consistent with the input dataset for MBP designs.", {
expect_error(preprocess_SCD(case = Student, phase = Phase,
session = Time, outcome = Behavior,
design = "MB", data = dat_dup_session))
expect_equal(names(MB_dat), names(MB_dat_clean)[1:4])
data("Laski")
dat1_Laski <- preprocess_SCD(case = case, phase = treatment,
session = time, outcome = outcome,
design = "MBP", data = Laski)
dat2_Laski <- preprocess_SCD(case = case, phase = treatment,
session = time, outcome = outcome,
design = "MBP", center = 4, data = Laski)
dat3_Laski <- preprocess_SCD(case = Laski$case, phase = Laski$treatment,
session = Laski$time, outcome = Laski$outcome,
design = "MB")
expect_error(preprocess_SCD(data = Laski))
expect_equivalent(dat1_Laski, dat3_Laski)
expect_equal(Laski$case, dat1_Laski$case)
expect_equal(Laski$treatment, dat1_Laski$treatment)
expect_equal(Laski$time, dat1_Laski$time)
expect_equal(Laski$outcome, dat1_Laski$outcome)
expect_equal(Laski$time, dat2_Laski$time + 4)
})
test_that("The returned dataset is consistent with the input dataset for TR designs.", {
data("Anglesea")
dat1_Ang <- preprocess_SCD(case = case, phase = condition,
session = session, outcome = outcome,
design = "TR", data = Anglesea)
dat2_Ang <- preprocess_SCD(case = Anglesea$case, phase = Anglesea$condition,
session = Anglesea$session, outcome = Anglesea$outcome,
design = "TR")
expect_equivalent(dat1_Ang, dat2_Ang)
expect_equal(Anglesea$case, dat1_Ang$case)
expect_equal(Anglesea$condition, dat1_Ang$condition)
expect_equal(Anglesea$session, dat1_Ang$session)
expect_equal(Anglesea$outcome, dat1_Ang$outcome)
Ang_case <- Anglesea$case
Ang_condition <- Anglesea$condition
Ang_session <- Anglesea$session
Ang_outcome <- Anglesea$outcome
dat3_Ang <- preprocess_SCD(case = Ang_case, phase = Ang_condition,
session = Ang_session, outcome = Ang_outcome,
design = "TR")
expect_equivalent(dat1_Ang, dat3_Ang)
})
test_that("The returned dataset is consistent with the input dataset for RMBB designs.", {
data("Thiemann2001")
dat1_Thiemann <- preprocess_SCD(design = "RMBB",
case = case, series = series,
phase = treatment, session = time,
outcome = outcome, data = Thiemann2001)
expect_equal(Thiemann2001$case, dat1_Thiemann$case)
expect_equal(Thiemann2001$series, dat1_Thiemann$series)
expect_equal(Thiemann2001$treatment, dat1_Thiemann$treatment)
expect_equal(Thiemann2001$time, dat1_Thiemann$time)
expect_equal(Thiemann2001$outcome, dat1_Thiemann$outcome)
dat2_Thiemann <- preprocess_SCD(design = "RMBB",
case = Thiemann2001$case, series = Thiemann2001$series,
phase = Thiemann2001$treatment, session = Thiemann2001$time,
outcome = Thiemann2001$outcome)
expect_equivalent(dat1_Thiemann, dat2_Thiemann)
Thi_case <- Thiemann2001$case
Thi_series <- Thiemann2001$series
Thi_phase <- Thiemann2001$treatment
Thi_session <- Thiemann2001$time
Thi_outcome <- Thiemann2001$outcome
dat3_Thiemann <- preprocess_SCD(design = "RMBB",
case = Thi_case, series = Thi_series,
phase = Thi_phase, session = Thi_session,
outcome = Thi_outcome)
expect_equivalent(dat1_Thiemann, dat3_Thiemann)
scramble <- sample(1:nrow(Thiemann2001))
dat4_Thiemann <- preprocess_SCD(design = "RMBB",
case = case, series = series,
phase = treatment, session = time,
outcome = outcome, data = Thiemann2001[scramble,])
expect_equivalent(dat1_Thiemann[scramble,], dat4_Thiemann)
})
test_that("The returned dataset is consistent with the input dataset for CMB designs.", {
data("Bryant2018")
dat1_Bry <- preprocess_SCD(design = "CMB",
cluster = group, case = case,
phase = treatment, session = session,
outcome = outcome, data = Bryant2018)
Bryant <- na.omit(Bryant2018)
expect_equal(Bryant$group, dat1_Bry$group)
expect_equal(Bryant$case, dat1_Bry$case)
expect_equal(Bryant$treatment, dat1_Bry$treatment)
expect_equal(Bryant$session, dat1_Bry$session)
expect_equal(Bryant$outcome, dat1_Bry$outcome)
dat2_Bry <- preprocess_SCD(design = "CMB",
cluster = Bryant2018$group, case = Bryant2018$case,
phase = Bryant2018$treatment, session = Bryant2018$session,
outcome = Bryant2018$outcome)
expect_equivalent(dat1_Bry, dat2_Bry)
Bry_group <- Bryant2018$group
Bry_case <- Bryant2018$case
Bry_treatment <- Bryant2018$treatment
Bry_session <- Bryant2018$session
Bry_outcome <- Bryant2018$outcome
dat3_Bry <- preprocess_SCD(design = "CMB",
cluster = Bry_group, case = Bry_case,
phase = Bry_treatment, session = Bry_session,
outcome = Bry_outcome)
expect_equivalent(dat1_Bry, dat3_Bry)
scramble <- sample(1:nrow(Bryant))
dat4_Bry <- preprocess_SCD(design = "CMB",
cluster = group, case = case,
phase = treatment, session = session,
outcome = outcome, data = Bryant[scramble,])
expect_equivalent(dat1_Bry[, 1:6], dat4_Bry[order(scramble), 1:6])
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.