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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
skip_on_cran()
skip_if_not_installed("ordinal")
d <- data.frame(
Stim = c(
"New", "New", "New", "New", "New", "New",
"Old", "Old", "Old", "Old", "Old", "Old"
),
Response = c(
"Confidence1", "Confidence2", "Confidence3", "Confidence4",
"Confidence5", "Confidence6", "Confidence1", "Confidence2",
"Confidence3", "Confidence4", "Confidence5", "Confidence6"
),
w = c(320, 295, 243, 206, 174, 159, 136, 188, 208, 256, 302, 333),
stringsAsFactors = FALSE
)
m1 <- ordinal::clm(ordered(Response) ~ Stim,
scale = ~Stim,
link = "probit",
data = d, weights = w
)
m2 <- ordinal::clm2(ordered(Response) ~ Stim,
scale = ~Stim,
link = "probit",
data = d, weights = w
)
test_that("model_parameters.clm", {
mp <- model_parameters(m1)
expect_equal(
mp$Parameter,
c(
"Confidence1|Confidence2", "Confidence2|Confidence3", "Confidence3|Confidence4",
"Confidence4|Confidence5", "Confidence5|Confidence6", "StimOld",
"StimOld"
),
tolerance = 1e-4
)
expect_equal(
mp$Component,
c("intercept", "intercept", "intercept", "intercept", "intercept", "location", "scale"),
tolerance = 1e-4
)
expect_equal(
mp$Coefficient,
c(-0.72845, -0.15862, 0.26583, 0.69614, 1.23477, 0.55237, -0.04069),
tolerance = 1e-4
)
mp <- model_parameters(m1, exponentiate = TRUE)
expect_equal(
mp$Coefficient,
c(0.48266, 0.85332, 1.30451, 2.006, 3.4376, 0.55237, -0.04069),
tolerance = 1e-4
)
expect_snapshot(print(mp))
})
test_that("model_parameters.clm2", {
mp <- model_parameters(m2)
expect_equal(
mp$Parameter,
c(
"Confidence1|Confidence2", "Confidence2|Confidence3", "Confidence3|Confidence4",
"Confidence4|Confidence5", "Confidence5|Confidence6", "StimOld",
"StimOld"
),
tolerance = 1e-4
)
expect_equal(
mp$Component,
c("intercept", "intercept", "intercept", "intercept", "intercept", "location", "scale"),
tolerance = 1e-4
)
expect_equal(
mp$Coefficient,
c(-0.72845, -0.15862, 0.26583, 0.69614, 1.23477, 0.55237, -0.04069),
tolerance = 1e-4
)
mp <- model_parameters(m2, exponentiate = TRUE)
expect_equal(
mp$Coefficient,
c(0.48266, 0.85332, 1.30451, 2.006, 3.4376, 0.55237, -0.04069),
tolerance = 1e-4
)
expect_snapshot(print(mp))
})
test_that("model_parameters.clmm, exponentiate works w/o component column", {
data(wine, package = "ordinal")
mox <- ordinal::clmm(rating ~ temp + contact + (1 | judge), data = wine)
out1 <- model_parameters(mox, exponentiate = FALSE)
out2 <- model_parameters(mox, exponentiate = TRUE)
expect_equal(out1$Coefficient, c(-1.62367, 1.51337, 4.22853, 6.08877, 3.063, 1.83488, 1.13113), tolerance = 1e-4)
expect_equal(out2$Coefficient, c(0.19717, 4.54199, 68.61606, 440.87991, 21.39156, 6.26441, 1.13113), tolerance = 1e-4)
expect_identical(attributes(out1)$coefficient_name, "Log-Odds")
expect_identical(attributes(out2)$coefficient_name, "Odds Ratio")
})
|
