Nothing
context("Utilities")
test_that("eigenCpp", {
mat <- kern_linear(1:3)
mat <- mat %*% mat
res1 <- eigenCpp(mat)
res2 <- eigen(mat)
expect_equal(res1$values, sort(res2$values), tolerance = 1e-6)
})
test_that("fastSquare", {
mat <- kern_linear(1:3)
res1 <- mat %*% mat
res2 <- fastSquare(mat)
expect_equal(res1, res2, tolerance = 1e-6)
})
test_that("fastSquareRoot", {
mat <- kern_linear(1:3)
mat2 <- mat %*% mat
res1 <- fastSquareRoot(mat2)
res2 <- fastSquareRoot2(mat2)
tmp <- eigen(mat2)
res3 <- tmp$vec %*% (t(tmp$vec) * sqrt(tmp$val))
expect_equal(res1, res2, tolerance = 1e-6)
expect_equal(res1, res3, tolerance = 1e-6)
expect_equal(res2, res3, tolerance = 1e-6)
})
test_that("fastSquareRoot", {
mat <- kern_linear(1:3)
mat2 <- mat %*% mat
res1 <- fastSquareRoot(mat2)
res2 <- fastSquareRoot2(mat2)
tmp <- eigen(mat2)
res3 <- tmp$vec %*% (t(tmp$vec) * sqrt(tmp$val))
expect_equal(res1, res2, tolerance = 1e-6)
expect_equal(res1, res3, tolerance = 1e-6)
expect_equal(res2, res3, tolerance = 1e-6)
})
test_that("fastVDiag", {
# for some reason it overwrites to environment...
mat <- kern_linear(1:3)
mat <- mat %*% mat
tmp <- eigenCpp(mat)
res <- fastVDiag(tmp$vec, tmp$val)
expect_equal(mat, res, tolerance = 1e-6)
})
test_that("check_theta", {
mod <- kernL(y ~ ., gen_smooth(10), kernel = "poly", est.offset = TRUE)
res <- capture.output(check_theta(mod))
expect_equal(res[1], "theta consists of 3:")
expect_equal(res[2], "log(lambda), log(offset), log(psi)")
})
test_that("gg_col_hue", {
res1 <- gg_col_hue(10)
expect_warning(res2 <- ipriorColPal(10))
expect_warning(res3 <- ggColPal(10))
expect_equal(res1, res2)
})
test_that("is.ipriorKernel", {
mod <- kernL(y ~ ., gen_smooth(10), kernel = "poly", est.offset = TRUE)
expect_true(is.ipriorKernel(mod))
})
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.