Nothing
context("check precog.test accuracy")
set.seed(112)
library(smerc)
data("neast")
data("neastw")
# setup up arguments
coords <- neast[, c("x", "y"), drop = TRUE]
pop <- neast$population
cases <- neast$cases
# precog
out <- precog.test(coords = coords, cases = cases,
pop = pop, w = neastw, nsim = 999)
# correct clusters
cl1 <- c(88, 96, 83, 78, 89, 91, 157, 127, 77, 141, 84,
210, 182, 205, 198, 213, 199, 194, 172, 216)
cl2 <- c(161, 163, 196, 202, 179, 221, 201, 183, 171,
102, 112)
cl3 <- c(24, 23, 230, 13, 227, 16)
test_that("check accuracy for precog.test with neast", {
expect_equal(out$clusters[[1]]$locids, cl1)
expect_equal(round(out$clusters[[1]]$r, 2), 22412.31)
expect_equal(out$clusters[[1]]$pop, 5504246)
expect_equal(out$clusters[[1]]$cases, 12646)
expect_equal(round(out$clusters[[1]]$exp, 2), 10984.75)
expect_equal(round(out$clusters[[1]]$smr, 2), 1.15)
expect_equal(round(out$clusters[[1]]$rr, 2), 1.19)
expect_equal(round(out$clusters[[1]]$loglik, 4), 148.8356)
expect_equal(round(out$clusters[[1]]$test_statistic, 4), 148.8356)
expect_equal(out$clusters[[1]]$pvalue, 0.001)
expect_equal(out$clusters[[2]]$locids, cl2)
expect_equal(round(out$clusters[[2]]$r, 2), 29469.36)
expect_equal(out$clusters[[2]]$pop, 1583988)
expect_equal(out$clusters[[2]]$cases, 3778)
expect_equal(round(out$clusters[[2]]$exp, 2), 3161.14)
expect_equal(round(out$clusters[[2]]$smr, 2), 1.20)
expect_equal(round(out$clusters[[2]]$rr, 2), 1.21)
expect_equal(round(out$clusters[[2]]$loglik, 5), 60.03691)
expect_equal(round(out$clusters[[2]]$test_statistic, 5), 60.03691)
expect_equal(out$clusters[[2]]$pvalue, 0.001)
expect_equal(out$clusters[[3]]$locids, cl3)
expect_equal(round(out$clusters[[3]]$r, 2), 7464.83)
expect_equal(out$clusters[[3]]$pop, 987552)
expect_equal(out$clusters[[3]]$cases, 2255)
expect_equal(round(out$clusters[[3]]$exp, 2), 1970.84)
expect_equal(round(out$clusters[[3]]$smr, 2), 1.14)
expect_equal(round(out$clusters[[3]]$rr, 2), 1.15)
expect_equal(round(out$clusters[[3]]$loglik, 5), 20.27554)
expect_equal(round(out$clusters[[3]]$test_statistic, 5), 20.27554)
expect_equal(out$clusters[[3]]$pvalue, 0.002)
})
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.