test_that("output is correct class", {
outp <- TD(-2, rnorm(15))
expect_is(outp, "htest")
})
test_that("all TD functionality works", {
controls <- rnorm(20)
expect_error(
TD(-2, controls, alternative = "g"),
NA
)
expect_error(
TD(-2, controls, alternative = "l"),
NA
)
expect_error(
TD(-2, controls, alternative = "t"),
NA
)
controls[1] <- NA
expect_error(
TD(-2, controls, na.rm = TRUE),
NA
)
expect_error(
TD(-2, controls, conf_int = FALSE, na.rm = TRUE),
NA
)
})
test_that("errors are occuring as they should", {
controls <- rnorm(15)
expect_error(TD(NA, controls), "Case is NA")
controls[1] <- NA
expect_error(TD(-2, controls), "Controls contains NA, set na.rm = TRUE to proceed")
expect_error(TD(-2, 1), "Not enough obs. Set sd and sample size for input of controls to be treated as mean")
expect_error(TD(-2, controls = 0, sd = 1), "Input sample size")
expect_error(TD(rnorm(2), controls, na.rm = TRUE), "Case should only have 1 observation")
expect_error(TD(2, controls, na.rm = TRUE, conf_level = 2), "Confidence level must be between 0 and < 1")
expect_error(TD(2, controls, na.rm = TRUE, conf_level = -2), "Confidence level must be between 0 and < 1")
})
test_that("summary input works as expected", {
t1 <- TD(-2, controls = 0, sd = 1, sample_size = 15)[["statistic"]][["t"]]
expect_equal(round(t1, 3), -1.936)
t2 <- TD(2, controls = 0, sd = 1, sample_size = 15)[["statistic"]][["t"]]
expect_equal(round(t2, 3), 1.936)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.