# TEST plot.mulTree
# Testing get.ylim
test_that("get.ylim works", {
# Error
expect_error(
get.ylim("bla")
)
# Output is 0.99 + 1.01 (max and min of 1 +/- 1%)
expect_equal(
get.ylim(1), c(0.99, 1.01)
)
# Output is 0.99 + 10.01
expect_equal(
get.ylim(matrix(seq(from=0, to=9), 5)), c(0, 9.09)
)
})
# Testing get.width
test_that("get.width works", {
# Error
expect_error(
get.width("bla")
)
# Output is 0,0,2,2
expect_equal(
get.width(1,1,1), c(0,0,2,2)
)
# Output is 9.7, 9.7, 10.3, 10.3
expect_equal(
get.width(c(0.10, 0.20, 0.30), 10, 3), c(9.7, 9.7, 10.3, 10.3)
)
})
# Testing get.height
test_that("get.height works", {
# Error
expect_error(
get.height("bla")
)
expect_error(
get.height(1,1,1)
)
expect_error(
get.height(matrix(seq(1:15), nrow=3), 4, 1)
)
# Output works
expect_equal(
get.height(matrix(seq(1:15), nrow=3), 1, 1), c(4,13,13,4)
)
expect_equal(
get.height(matrix(seq(1:15), nrow=3), 3, 2), c(9,12,12,9)
)
})
# Testing plot examples
test_that("plot.mulTree examples work", {
## read in the data
data(lifespan.mcmc)
## summarising the results
summarized_data <- summary(lifespan.mcmc)
## plotting the results
expect_null(plot(summarized_data))
## Same plot using more options
expect_null(
plot(summarized_data, horizontal = TRUE, ylab = "", ylim = c(-2,2),
main = "Posterior distributions", cex.terms = 0.5, cex.coeff = 0.8,
terms = c("Intercept", "BodyMass", "Volancy", "Phylogeny", "Residuals"),
col = c("red"), cex.main = 0.8)
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.