Nothing
#context("qad and zeta1")
test_that("zeta1 returns exact values", {
#Example 1
x <- c(1,2,3,4,5)
y <- c(1,2,3,4,5)
expect_equal(qad::zeta1(x, y, resolution = 5), 0.9)
expect_equal(qad::zeta1(x, y, resolution = 2), 0.6)
expect_equal(qad::zeta1(x, y, resolution = 1), 0)
expect_equal(qad::zeta1(x, y, resolution = 10), 0.9)
expect_equal(qad::zeta1(x,y), qad::zeta1(y,x))
#Example 2
x <- c(1,2,3,4,5,6,7,8,9,10)
y <- c(1,1,1,1,3,3,1,1,1,1)
expect_equal(qad::zeta1(x, y, resolution = 10), 0.48)
expect_equal(qad::zeta1(y, x, resolution = 10), 0.24)
#Example (checkmin copula)
NN <- c(2,10,22,44,122,200)
for(N in NN){
x <- 1:N
y <- 1:N
expect_equal(qad::zeta1(x,y, resolution = N), 1-1/(2*N))
expect_equal(qad::zeta1(y,x, resolution = N), 1-1/(2*N))
}
})
test_that("class of qad output and correct data output", {
n <- 20
x <- runif(n)
y <- runif(n)
expect_output(qad::qad(x,y))
fit <- qad::qad(x,y, print = FALSE)
expect_identical(class(fit), "qad")
expect_equal(fit$n, n)
expect_equal(fit$resolution, floor(sqrt(n)))
expect_equal(colSums(fit$mass_matrix), rowSums(fit$mass_matrix))
expect_equal(NROW(fit$data), n)
expect_equal(NCOL(fit$data), 2)
})
n <- 100
x <- runif(n)
y <- runif(n)
z <- runif(n)
df <- data.frame(x,y,z)
pw_fit <- pairwise.qad(df)
test_that("pairwise qad", {
expect_equal(class(pw_fit), "list")
expect_equal(names(pw_fit), c("q",
"max.dependence",
"asymmetry",
"q_p.values",
"q_p.values.adjusted",
"max.dependence_p.values",
"max.dependence_p.values.adjusted",
"asymmetry_p.values",
"asymmetry_p.values.adjusted",
"resolution",
"n_removed_00",
"p.adjust.method")
)
expect_equal(class(pw_fit$q), "data.frame")
expect_equal(class(pw_fit$max.dependence), "data.frame")
expect_equal(class(pw_fit$asymmetry), "data.frame")
expect_equal(class(pw_fit$q_p.values), "data.frame")
expect_equal(class(pw_fit$max.dependence_p.values), "data.frame")
expect_equal(class(pw_fit$asymmetry_p.values), "data.frame")
expect_equal(class(pw_fit$resolution), "data.frame")
})
test_that("heatmap.qad", {
expect_equal(class(qad::heatmap.qad(pw_fit, select = "dependence")), c("gg", "ggplot"))
expect_equal(class(qad::heatmap.qad(pw_fit, select = "dependence",
fontsize = 5, significance = T, scale = "rel", color = "viridis", title = "test")), c("gg", "ggplot"))
})
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.