d <- 5L
lambda <- 8e-2
alpha <- 4e-1
rho <- 3 * alpha / (4 - alpha)
tau <- alpha / (2 - alpha)
nu <- alpha
bf <- ScaledBernsteinFunction(
scale = lambda,
original = SumOfBernsteinFunctions(
first = LinearBernsteinFunction(scale = 1 - nu),
second = ConstantBernsteinFunction(constant = nu)
))
ex_qmatrix <- rmo::exQMatrix(bf, d)
ex_intensities <- rmo::exIntensities(bf, d = d)
test_that("`ArmageddonExtMO2FParam`-class is correctly initialized", {
parm <- ArmageddonExtMO2FParam()
expect_s4_class(parm, "ArmageddonExtMO2FParam")
setDimension(parm) <- d
setBernsteinFunction(parm) <- bf
setLambda(parm) <- lambda
setNu(parm) <- nu
expect_error(validObject(parm), NA)
expect_equal(getDimension(parm), d)
expect_equal(getExQMatrix(parm), ex_qmatrix)
expect_equal(getExIntensities(parm), ex_intensities)
expect_equal(getBernsteinFunction(parm), bf)
expect_equal(getLambda(parm), lambda)
expect_equal(getNu(parm), nu)
expect_equal(getRho(parm), rho)
expect_equal(getTau(parm), tau)
expect_equal(getAlpha(parm), alpha)
expect_equal(parm, ArmageddonExtMO2FParam(d, lambda, nu))
expect_equal(parm, ArmageddonExtMO2FParam(d, lambda, rho = rho))
expect_equal(parm, ArmageddonExtMO2FParam(d, lambda, tau = tau))
expect_equal(parm, ArmageddonExtMO2FParam(d, lambda, alpha = alpha))
expect_equal(as(parm, "ExtMOParam"), ExtMOParam(d, bf))
expect_equal(as(parm, "ExMOParam"), ExMOParam(ex_intensities))
expect_equal(as(parm, "ExMarkovParam"), ExMarkovParam(ex_qmatrix))
})
test_that("`ArmageddonExtMO2FParam`-class setters can be used in arbitrary order", { # nolint
parm <- ArmageddonExtMO2FParam(d, lambda, nu)
parm2 <- ArmageddonExtMO2FParam()
setDimension(parm2) <- d
setLambda(parm2) <- lambda
setRho(parm2) <- rho
expect_equal(parm, parm2)
parm2 <- ArmageddonExtMO2FParam()
setDimension(parm2) <- d
setLambda(parm2) <- lambda
setTau(parm2) <- tau
expect_equal(parm, parm2)
parm2 <- ArmageddonExtMO2FParam()
setDimension(parm2) <- d
setNu(parm2) <- nu
setLambda(parm2) <- lambda
expect_equal(parm, parm2)
parm2 <- ArmageddonExtMO2FParam()
setLambda(parm2) <- lambda
setDimension(parm2) <- d
setNu(parm2) <- nu
expect_equal(parm, parm2)
parm2 <- ArmageddonExtMO2FParam()
setNu(parm2) <- nu
setDimension(parm2) <- d
setLambda(parm2) <- lambda
expect_equal(parm, parm2)
parm2 <- ArmageddonExtMO2FParam()
setLambda(parm2) <- lambda
setNu(parm2) <- nu
setDimension(parm2) <- d
expect_equal(parm, parm2)
parm2 <- ArmageddonExtMO2FParam()
setNu(parm2) <- nu
setLambda(parm2) <- lambda
setDimension(parm2) <- d
expect_equal(parm, parm2)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.