1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
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
)
})
|
