## To test gsHR() we use zn2hr() which is already tested.
# For comparing floating-point numbers, an exact match cannot be expected.
# For such test cases, the tolerance is set to 1e-6 (= 0.000001),
# a sufficiently low value
#-------------------------------------------------------------------------------
test_that(desc = "test: output validation",
code = {
xgs <- gsSurv(lambdaC = .2, hr = .5, eta = .1, T = 2, minfup = 1.5)
obs <- gsHR(z = xgs$lower$bound, i = 1:3, x = xgs)
expected <- zn2hr(z = xgs$lower$bound, n = as.vector(xgs$eDC + xgs$eDE))
expect_lte(abs(obs[1] - expected[1]), 1e-6)
expect_lte(abs(obs[2] - expected[2]), 1e-6)
expect_lte(abs(obs[3] - expected[3]), 1e-6)
})
test_that(desc = "output validation",
code = {
xgs <- gsSurv(lambdaC = .2, hr = .5, eta = .1, T = 2, minfup = 1.5)
obs <- gsHR(z = xgs$upper$bound, i = 1:3, x = xgs)
expected <- zn2hr(z = xgs$upper$bound, n = as.vector(xgs$eDC + xgs$eDE))
expect_lte(abs(obs[1] - expected[1]), 1e-6)
expect_lte(abs(obs[2] - expected[2]), 1e-6)
expect_lte(abs(obs[3] - expected[3]), 1e-6)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.