skip_if_not_installed("MCMCglmm")
data(PlodiaPO, package = "MCMCglmm")
set.seed(123)
m1 <- MCMCglmm::MCMCglmm(
PO ~ plate,
random = ~FSfamily,
data = PlodiaPO,
verbose = FALSE,
nitt = 1300,
burnin = 300,
thin = 1
)
test_that("ci", {
expect_equal(
ci(m1)$CI_low,
c(0.97495, 0.03407),
tolerance = 0.01
)
})
test_that("se", {
expect_equal(
standard_error(m1)$SE,
c(0.02309, 0.00509),
tolerance = 0.01
)
})
test_that("p_value", {
expect_equal(
p_value(m1)$p,
c(0, 0),
tolerance = 0.01
)
})
test_that("model_parameters", {
expect_equal(
model_parameters(m1, centrality = "mean", verbose = FALSE)$Mean,
c(1.0132, 0.04232),
tolerance = 0.01
)
})
test_that("model_parameters", {
expect_equal(
model_parameters(m1, centrality = "median", verbose = FALSE)$Median,
c(1.01382, 0.04207),
tolerance = 0.01
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.