Nothing
context("Test logic of trunc_constant")
test_that("LRRi, LRRd, SMD, Tau and their standard errors are positive for count data.", {
set.seed(20210323)
A_dat <- rpois(10, 3)
B_dat <- rep(0, 5)
C_dat <- rep(10, 5)
D_dat <- rpois(5, 3)
E_dat <- rep(10, 10)
F_dat <- rep(2, 5)
G_dat <- rep(10, 5)
ABd <- calc_ES(A_data = A_dat, B_data = B_dat,
improvement = "decrease", scale = "count",
ES = c("LRRd","LRRi","SMD","Tau")) %>% dplyr::select(-baseline_SD)
ABd <- rbind(ABd, Tau(A_dat, B_dat, improvement = "decrease", SE = "Hanley"))
ABi <- calc_ES(A_data = A_dat, B_data = B_dat,
improvement = "increase", scale = "count",
ES = c("LRRd","LRRi","SMD","Tau")) %>% dplyr::select(-baseline_SD)
ABi <- rbind(ABi, Tau(A_dat, B_dat, improvement = "increase", SE = "Hanley"))
expect_equal(ABd$Est, -ABi$Est)
expect_equal(ABd$SE, ABi$SE)
expect_equal(ABd$CI_lower, -ABi$CI_upper)
expect_equal(ABd$CI_upper, -ABi$CI_lower)
expect_true(all(ABd$SE > 0))
CDd <- calc_ES(A_data = C_dat, B_data = D_dat,
improvement = "decrease", scale = "count", std_dev = "pool",
ES = c("LRRd","LRRi","SMD","Tau")) %>% dplyr::select(-pooled_SD)
CDd <- rbind(CDd, Tau(C_dat, D_dat, improvement = "decrease", SE = "Hanley"))
CDi <- calc_ES(A_data = C_dat, B_data = D_dat,
improvement = "increase", scale = "count", std_dev = "pool",
ES = c("LRRd","LRRi","SMD","Tau")) %>% dplyr::select(-pooled_SD)
CDi <- rbind(CDi, Tau(C_dat, D_dat, improvement = "increase", SE = "Hanley"))
expect_equal(CDd$Est, -CDi$Est)
expect_equal(CDd$SE, CDi$SE)
expect_equal(CDd$CI_lower, -CDi$CI_upper)
expect_equal(CDd$CI_upper, -CDi$CI_lower)
expect_true(all(CDd$SE > 0))
expect_warning(SMD(A_data = C_dat, B_data = D_dat, std_dev = "baseline"))
EFd <- calc_ES(A_data = E_dat, B_data = F_dat,
improvement = "decrease", scale = "count",
ES = c("LRRd","LRRi","Tau"))
EFd <- rbind(EFd, Tau(E_dat, F_dat, improvement = "decrease", SE = "Hanley"))
EFi <- calc_ES(A_data = E_dat, B_data = F_dat,
improvement = "increase", scale = "count",
ES = c("LRRd","LRRi","Tau"))
EFi <- rbind(EFi, Tau(E_dat, F_dat, improvement = "increase", SE = "Hanley"))
expect_equal(EFd$Est, -EFi$Est)
expect_equal(EFd$SE, EFi$SE)
expect_equal(EFd$CI_lower, -EFi$CI_upper)
expect_equal(EFd$CI_upper, -EFi$CI_lower)
expect_true(all(EFd$SE > 0))
expect_warning(SMD(A_data = E_dat, B_data = F_dat, std_dev = "baseline"))
expect_warning(SMD(A_data = E_dat, B_data = F_dat, std_dev = "pool"))
EGd <- calc_ES(A_data = E_dat, B_data = G_dat,
improvement = "decrease", scale = "count",
ES = c("LRRd","LRRi","Tau"))
EGd <- rbind(EGd, Tau(E_dat, G_dat, improvement = "decrease", SE = "Hanley"))
EGd <- rbind(EGd, Tau(E_dat, G_dat, improvement = "decrease", SE = "null"))
EGi <- calc_ES(A_data = E_dat, B_data = G_dat,
improvement = "increase", scale = "count",
ES = c("LRRd","LRRi","Tau"))
EGi <- rbind(EGi, Tau(E_dat, G_dat, improvement = "increase", SE = "Hanley"))
EGi <- rbind(EGi, Tau(E_dat, G_dat, improvement = "increase", SE = "null"))
expect_equal(EGd$Est, -EGi$Est)
expect_equal(EGd$SE, EGi$SE)
expect_equal(EGd$CI_lower, -EGi$CI_upper)
expect_equal(EGd$CI_upper, -EGi$CI_lower)
expect_true(all(EGd$SE > 0))
expect_warning(SMD(A_data = E_dat, B_data = G_dat, std_dev = "baseline"))
expect_warning(SMD(A_data = E_dat, B_data = G_dat, std_dev = "pool"))
})
test_that("trunc_const() works for batched datasets.", {
data("Schmidt2007")
Schmidt2007$Metric <- ifelse(Schmidt2007$Metric == "count","Natural Count", "Percentage")
expect_error(
output_pkg <-
batch_calc_ES(dat = Schmidt2007,
grouping = c(Behavior_type, Case_pseudonym, Phase_num),
condition = Condition,
outcome = Outcome,
session_number = Session_number,
baseline_phase = "A",
intervention_phase = "B",
ES = c("LRRd", "LRRi"),
improvement = direction,
pct_change = FALSE,
scale = Metric,
intervals = NA,
observation_length = NA,
D_const = NA,
std_dev = "baseline",
confidence = 0.95,
Kendall = FALSE,
pretest_trend = FALSE,
format = "long"
)
)
})
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.