test_that("Test basic association test.", {
# Type 1 error.
n <- 1e3
withr::local_seed(101)
x <- cbind(1, rnorm(n))
g <- replicate(1e3, rbinom(n = n, size = 2, prob = 0.25))
storage.mode(g) <- "numeric"
y <- as.numeric(x %*% c(1, 1)) + stats::rnorm(n)
p_score <- DINT(y = y, G = g, X = x, test = "Score", simple = TRUE)
p_wald <- DINT(y = y, G = g, X = x, test = "Wald", simple = TRUE)
# Approximately correct type 1 error.
t1e_score <- mean(p_score <= 0.05)
t1e_wald <- mean(p_wald <= 0.05)
expect_equal(t1e_score, 0.05, tolerance = 0.05)
expect_equal(t1e_wald, 0.05, tolerance = 0.05)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.