Nothing
# Check that the calculation of the expected information for the GEV
# distribution agrees with mev::gev.infomat()
sigma <- 1
mevMat <- cbind(
c(1/sigma^2, -0.422784335098467/sigma^2, 0.41184033042644/sigma),
c(-0.422784335098467/sigma^2, 1.82368066085288/sigma^2, 0.332484907160274/sigma),
c(0.41184033042644/sigma, 0.332484907160274/sigma, 2.42360605517703)
)
# Based on the constants
evilsMat1 <- cbind(c(gev11e(1, 0), gev12e0Constant, gev13e0Constant),
c(gev12e0Constant, gev22e0Constant, gev23e0Constant),
c(gev13e0Constant, gev23e0Constant, gev33e0Constant))
# Based on the functions that create the constants
evilsMat1f <- cbind(c(gev11e(1, 0), gev12e0Fn(), gev13e0Fn()),
c(gev12e0Fn(), gev22e0Fn(), gev23e0Fn()),
c(gev13e0Fn(), gev23e0Fn(), gev33e0Fn()))
# Based on calling the functions with xi = 0
evilsMat2 <- cbind(c(gev11e(1, 0), gev12e(1, 0), gev13e(1, 0)),
c(gev12e(1, 0), gev22e(1, 0), gev23e(1, 0)),
c(gev13e(1, 0), gev23e(1, 0), gev33e(0)))
test_that("Expected Information constants agree with mev::gev.infomat()", {
testthat::expect_equal(mevMat, evilsMat1)
})
test_that("Expected Information constant funs agree with mev::gev.infomat()", {
testthat::expect_equal(mevMat, evilsMat1f)
})
test_that("Expected Information constants agree with mev::gev.infomat()", {
testthat::expect_equal(mevMat, evilsMat2)
})
test_that("gevExpInfo(shape = 0) agrees with mev::gev.infomat()", {
testthat::expect_equal(mevMat, gevExpInfo(scale = 1, shape = 0),
ignore_attr = TRUE)
})
# Check that gevExpInfo() throws an error when shape <= -0.5
test_that("gevExpInfo() errors for shape = -1/2", {
testthat::expect_error(gevExpInfo(scale = 1, shape = -1/2))
})
test_that("gevExpInfo() errors for shape = -1", {
testthat::expect_error(gevExpInfo(scale = 1, shape = -1))
})
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.