tests/testthat/test-Variance.R

In RcppTN: Rcpp-Based Truncated Normal Distribution RNG and Family

context("Variance Calculation")

test_that("Variances  are identical for symmetric cases", {
lows <- c(-1, 5, -Inf, 4, 4, -Inf, 7)
highs <- c(1, Inf, 10, 7, 4.1, Inf, 100)

for(case in 1:length(lows)) {
    out1 <- vtn(.mean = 0,
                .low = lows[case],
                .high = highs[case]
                )
    out2 <- vtn(.mean = 0,
                .high = -lows[case],
                .low = -highs[case]
                )
    expect_true(abs(out1 - out2) < .01)
}
}
          )

test_that("Variances are positive.", {
    means <- seq(0, -37, by = -1)
    for (case in 1:length(means)) {
        out1 <- vtn(.mean = means[case],
                    .low = 0,
                    .high = Inf
                    )
        expect_true(out1 > 0)
    }
    means <- seq(0, 37, by = 1)
    for (case in 1:length(means)) {
        out1 <- vtn(.mean = means[case],
                    .low = -Inf,
                    .high = 0
                    )
        expect_true(out1 > 0)
    }
}
          )

test_that("Extreme variances are defined.", {
    means <- seq(0, -37, by = -1)
    for (case in 1:length(means)) {
        out1 <- vtn(.mean = means[case],
                    .low = 0,
                    .high = Inf
                    )
        expect_true(is.finite(out1))
    }
    means <- seq(0, 37, by = 1)
    for (case in 1:length(means)) {
        out1 <- vtn(.mean = means[case],
                    .low = -Inf,
                    .high = 0
                    )
        expect_true(is.finite(out1))
    }
}
          )