partition <- list(1L:2L, 3L:6L, 7L:8L)
composition <- purrr::map_int(partition, length)
lambda <- 8e-2
alpha <- c(0.3, 0.8)
d <- sum(composition)
rho <- 3 * alpha / (4 - alpha)
tau <- alpha / (2 - alpha)
fraction <- (2 * alpha[[1]] + 1 - alpha[[2]]) / 2
models <- purrr::map(composition, ~{
ArmageddonExtMO2FParam(.x, lambda, alpha = (alpha[[2]] - alpha[[1]]) / (1 - fraction))
})
models <- c(
list(ArmageddonExtMO2FParam(d, lambda, alpha = alpha[[1]] / fraction)),
models)
nu <- purrr::map_dbl(models[1:2], getNu)
test_that("`ArmageddonH2ExtMO3FParam`-class is correctly initialized", {
parm <- ArmageddonH2ExtMO3FParam()
expect_s4_class(parm, "ArmageddonH2ExtMO3FParam")
setComposition(parm) <- composition
setFraction(parm) <- fraction
setLambda(parm) <- lambda
setNu(parm) <- nu
expect_error(validObject(parm), NA)
expect_equal(getDimension(parm), d)
expect_equal(getComposition(parm), composition)
expect_equal(getPartition(parm), partition)
expect_equal(getModels(parm), models)
expect_equal(getFraction(parm), fraction)
expect_equal(getLambda(parm), lambda)
expect_equal(getNu(parm), nu)
expect_equal(getRho(parm), rho)
expect_equal(getTau(parm), tau)
expect_equal(parm, ArmageddonH2ExtMO3FParam(composition, lambda, nu, fraction))
expect_equal(parm, ArmageddonH2ExtMO3FParam(composition, lambda, fraction = fraction, rho = rho))
expect_equal(parm, ArmageddonH2ExtMO3FParam(composition, lambda, fraction = fraction, tau = tau))
expect_equal(as(parm, "H2ExMarkovParam"), H2ExMarkovParam(fraction, models))
expect_equal(as(parm, "H2ExMOParam"), H2ExMOParam(fraction, models))
expect_equal(as(parm, "H2ExtMOParam"), H2ExtMOParam(fraction, models))
})
test_that("`ArmageddonH2ExtMO3FParam`-class setters can be used in arbitrary order", { # nolint
parm <- ArmageddonH2ExtMO3FParam(composition, lambda, nu, fraction)
parm2 <- ArmageddonH2ExtMO3FParam()
setComposition(parm2) <- composition
setLambda(parm2) <- lambda
setRho(parm2) <- rho
expect_equal(parm, parm2)
parm2 <- ArmageddonH2ExtMO3FParam()
setComposition(parm2) <- composition
setLambda(parm2) <- lambda
setTau(parm2) <- tau
expect_equal(parm, parm2)
parm2 <- ArmageddonH2ExtMO3FParam()
setComposition(parm2) <- composition
setRho(parm2) <- rho
setLambda(parm2) <- lambda
expect_equal(parm, parm2)
parm2 <- ArmageddonH2ExtMO3FParam()
setLambda(parm2) <- lambda
setComposition(parm2) <- composition
setRho(parm2) <- rho
expect_equal(parm, parm2)
parm2 <- ArmageddonH2ExtMO3FParam()
setRho(parm2) <- rho
setComposition(parm2) <- composition
setLambda(parm2) <- lambda
expect_equal(parm, parm2)
parm2 <- ArmageddonH2ExtMO3FParam()
setLambda(parm2) <- lambda
setRho(parm2) <- rho
setComposition(parm2) <- composition
expect_equal(parm, parm2)
parm2 <- ArmageddonH2ExtMO3FParam()
setRho(parm2) <- rho
setLambda(parm2) <- lambda
setComposition(parm2) <- composition
expect_equal(parm, parm2)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.