swh:1:snp:2c68a6c5a8af2f06ac2c0225927f25b54fd1f9d0
Tip revision: 428249f43a9c6fd0c425b28deb5fee51a9525d69 authored by Dominique Makowski on 18 September 2022, 01:46:03 UTC
version 0.13.0
version 0.13.0
Tip revision: 428249f
test-contr.R
test_that("contr.equalprior | gen", {
skip_on_cran()
set.seed(1234)
k <- 3
g <- 4.1
contr1 <- contr.equalprior(k, contrasts = TRUE)
contr2 <- contr.equalprior(k, contrasts = FALSE)
samps1 <- replicate(ncol(contr1), {
rnorm(4e3, 0, g)
})
samps2 <- replicate(ncol(contr2), {
rnorm(4e3, 0, g)
})
means1 <- t(contr1 %*% t(samps1))
means2 <- t(contr2 %*% t(samps2))
expect_equal(mean(apply(means1, 2, sd)), mean(apply(means2, 2, sd)), tolerance = 0.1)
})
test_that("contr.equalprior | pairs", {
skip_on_cran()
set.seed(1234)
k <- 3
g <- 4.1
contr1 <- contr.equalprior_pairs(k, contrasts = TRUE)
contr2 <- contr.equalprior_pairs(k, contrasts = FALSE)
samps1 <- replicate(ncol(contr1), {
rnorm(4e3, 0, g)
})
samps2 <- replicate(ncol(contr2), {
rnorm(4e3, 0, g)
})
means1 <- t(contr1 %*% t(samps1))
means2 <- t(contr2 %*% t(samps2))
w <- matrix(c(
-1, 1, 0,
1, 0, -1,
0, -1, 1
), 3, 3)
pairs1 <- t(w %*% t(means1))
pairs2 <- t(w %*% t(means2))
expect_equal(mean(apply(pairs1, 2, sd)), g, tolerance = 0.1)
expect_equal(mean(apply(pairs1, 2, sd)), mean(apply(pairs2, 2, sd)), tolerance = 0.1)
})
test_that("contr.equalprior | dev", {
skip_on_cran()
set.seed(1234)
k <- 3
g <- 4.1
contr1 <- contr.equalprior_deviations(k, contrasts = TRUE)
contr2 <- contr.equalprior_deviations(k, contrasts = FALSE)
samps1 <- replicate(ncol(contr1), {
rnorm(4e3, 0, g)
})
samps2 <- replicate(ncol(contr2), {
rnorm(4e3, 0, g)
})
means1 <- t(contr1 %*% t(samps1))
means2 <- t(contr2 %*% t(samps2))
expect_equal(mean(apply(means1, 2, sd)), g, tolerance = 0.1)
expect_equal(mean(apply(means1, 2, sd)), mean(apply(means2, 2, sd)), tolerance = 0.1)
})