Nothing
context("w2n")
test_that("The output is in the right format", {
nbStates <- 2
nbCovs <- 2
par <- c(0.001, 0.999, 0.5, 0.7, 1500.3, 0.001)
parSize <- c(2, 1)
bounds <- matrix(c(0, 1, 0, 1, 0, 1,
0, 1, 0, Inf, 0, Inf),
byrow = TRUE, ncol = 2)
beta <- matrix(rnorm(6), ncol = 2, nrow = 3)
delta <- c(0.6, 0.4)
wpar <- n2w(par, bounds, beta, delta, nbStates, FALSE)
p <- w2n(wpar, bounds, parSize, nbStates, nbCovs, FALSE, FALSE)
expect_equal(length(p$stepPar), parSize[1]*nbStates)
expect_equal(length(p$anglePar), parSize[2]*nbStates)
expect_equal(dim(p$beta), dim(beta))
expect_equal(length(p$delta), length(delta))
})
test_that("w2n and n2w are inverse", {
nbStates <- 2
nbCovs <- 2
par <- c(0.001, 0.999, 0.7, 0.7, 1500.3, 0.001)
parSize <- c(2, 1)
bounds <- matrix(c(0, 1, 0, 1, 0, 1,
0, 1, -Inf, Inf, -Inf, Inf),
byrow = TRUE, ncol = 2)
beta <- matrix(rnorm(6), ncol = 2, nrow = 3)
delta <- c(0.6, 0.4)
wpar <- n2w(par, bounds, beta, delta, nbStates, FALSE)
p <- w2n(wpar, bounds, parSize, nbStates, nbCovs, FALSE, FALSE)
expect_equal(c(t(p$stepPar)), par[1:(parSize[1]*nbStates)], tolerance = 1e-10)
expect_equal(c(t(p$anglePar)), par[(parSize[1]*nbStates+1):(sum(parSize)*nbStates)], tolerance = 1e-10)
expect_equal(p$beta, beta, tolerance = 1e-10)
expect_equal(p$delta, delta, tolerance = 1e-10)
})
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.