d <- 5L
lambda <- 8e-2
alpha <- 4e-1
nu <- log2(2 - alpha)
bf <- ScaledBernsteinFunction(
scale = lambda, original = AlphaStableBernsteinFunction(alpha = nu))
ex_qmatrix <- rmo::exQMatrix(bf, d)
ex_intensities <- rmo::exIntensities(bf, d = d)
test_that("`ExtMOParam`-class is correctly initialized", {
parm <- ExtMOParam()
expect_s4_class(parm, "ExtMOParam")
setDimension(parm) <- d
setBernsteinFunction(parm) <- bf
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(parm, ExtMOParam(d, bf))
expect_equal(as(parm, "ExMOParam"), ExMOParam(ex_intensities))
expect_equal(as(parm, "ExMarkovParam"), ExMarkovParam(ex_qmatrix))
})
test_that("`ExtMOParam`-class setters can be used in arbitrary order", {
parm <- ExtMOParam(d, bf)
parm2 <- ExtMOParam()
setBernsteinFunction(parm2) <- bf
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.