Nothing
test_that("previous values can be reproduced", {
set.seed(20191204)
cpos <- suppressWarnings(
find_critical_pos(rho = c(.1, .7), sample_size_max = 1e3,
n_studies = 40e3)
)
# https://github.com/nicebread/corEvol
nicebread <- matrix(c(252, 66, 360, 96, 474, 129), nrow = 2)
cpos_old <- read.csv("cpos4.csv")
cpos_old$precision_rel <- as.numeric(c(NA, NA))
diff_rel <- (nicebread - cpos[, 2:4]) / nicebread
expect_equal(cpos, cpos_old, check.attributes = FALSE)
# average relative deviation within 4%
expect_true(abs(mean(unlist(diff_rel))) < 0.04,
info = print(round(abs(mean(unlist(diff_rel))), 2)))
# individual relative deviation within 10%
expect_true(all(abs(diff_rel) < .10), info = print(diff_rel))
cpos_mc <- suppressWarnings(
find_critical_pos(rho = c(.1, .7), sample_size_max = 1e3,
n_studies = 40e3,
n_cores = 2)
)
diff_rel_mc <- (nicebread - cpos_mc[, 2:4]) / nicebread
# average relative deviation within 4%
expect_true(abs(mean(unlist(diff_rel_mc))) < 0.04,
info = print(round(abs(mean(unlist(diff_rel_mc))), 2)))
# individual relative deviation within 10%
expect_true(all(abs(diff_rel_mc) < .10), info = print(diff_rel_mc))
diff_rel_mc_replace <- (nicebread - cpos_mc[, 2:4]) / nicebread
cpos_mc_replace <- suppressWarnings(
find_critical_pos(rho = c(.1, .7), sample_size_max = 1e3,
n_studies = 40e3,
n_cores = 2,
replace = FALSE)
)
# average relative deviation within 4%
expect_true(abs(mean(unlist(diff_rel_mc_replace))) < 0.04,
info = round(abs(mean(unlist(diff_rel_mc_replace))), 2))
# individual relative deviation within 10%
expect_true(all(abs(diff_rel_mc_replace) < .10),
info = print(diff_rel_mc_replace))
})
test_that("unloading package works", {
expect_null(detach(package:fastpos, unload = TRUE))
})
test_that("create_pop_inexact (not used atm) works", {
expect_equal(round(cor(create_pop_inexact(0.5, 1e6))[1, 2], 2), 0.50)
})
test_that("create_pop works", {
expect_equal(round(cor(create_pop(0.5, 1e6))[1, 2], 2), 0.50)
})
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.