context("WMWssp")
# data
x <- c(315,375,356,374,412,418,445,403,431,410,391,475,379)
y <- x - 20
# calculate sample size, true result
result_WMWsspN <- 111.3719511
result_WMWsspp <- 0.3491124
result_power <- 0.9 # 0.7
suppressWarnings(RNGversion("3.5.0"))
set.seed(1)
result_simpower <- WMWssp::WMWssp(x, y, alpha = 0.05, power = 0.8, t = 1/2, simulation = TRUE, nsim = 10)$result[3, ]
test_that("function WMWssp", {
expect_equivalent(WMWssp::WMWssp(x, y, alpha = 0.05, power = 0.8, t = 1/2)$result[4,], result_WMWsspN, tolerance=1e-4)
expect_equivalent(WMWssp::WMWssp(x, y, alpha = 0.05, power = 0.8, t = 1/2)$result[3,], result_WMWsspp, tolerance=1e-4)
expect_output(print(WMWssp::WMWssp(x, y, alpha = 0.05, power = 0.8, t = 1/2)))
expect_output(summary(WMWssp::WMWssp(x, y, alpha = 0.05, power = 0.8, t = 1/2)))
expect_equivalent(result_power, result_simpower)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.