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
105
|
skip_if_not_installed("tinytable")
test_that("include_reference, on-the-fly factors", {
data(mtcars)
d <- as.data.frame(mtcars)
d$gear <- as.factor(d$gear)
d$am <- as.factor(d$am)
m1 <- lm(mpg ~ as.factor(gear) + factor(am) + hp, data = mtcars)
m2 <- lm(mpg ~ gear + am + hp, data = d)
out1 <- model_parameters(m1, include_reference = TRUE)
out2 <- model_parameters(m2, include_reference = TRUE)
expect_snapshot(print(out1))
expect_snapshot(print(out2))
expect_equal(attributes(out1)$pretty_names, attributes(out2)$pretty_names, ignore_attr = TRUE)
expect_equal(out1$Coefficient, out2$Coefficient, tolerance = 1e-4)
out <- compare_parameters(m1, m2, include_reference = TRUE)
expect_snapshot(print_md(out, engine = "tt"))
})
skip_if(getRversion() < "4.3.3")
skip_if_not_installed("datawizard")
test_that("include_reference, on-the-fly factors", {
data(mtcars)
d <- as.data.frame(mtcars)
d$gear <- as.factor(d$gear)
d$am <- as.factor(d$am)
m1 <- lm(mpg ~ as.factor(gear) + factor(am) + hp, data = mtcars)
m2 <- lm(mpg ~ gear + am + hp, data = d)
out1 <- model_parameters(m1, include_reference = TRUE)
out3 <- mtcars |>
datawizard::data_modify(gear = factor(gear), am = as.factor(am)) |>
lm(formula = mpg ~ gear + am + hp) |>
model_parameters(include_reference = TRUE)
expect_equal(attributes(out1)$pretty_names, attributes(out3)$pretty_names, ignore_attr = TRUE)
})
test_that("include_reference, with pretty formatted cut", {
data(mtcars)
mtcars$mpg_cut <- cut(mtcars$mpg, breaks = c(0, 20, 30, 100))
m <- lm(wt ~ mpg_cut, data = mtcars)
out <- parameters(m, include_reference = TRUE)
expect_identical(
attributes(out)$pretty_names,
c(
`(Intercept)` = "(Intercept)", `mpg_cut(0,20]` = "mpg cut [>0-20]",
`mpg_cut(20,30]` = "mpg cut [>20-30]", `mpg_cut(30,100]` = "mpg cut [>30-100]"
)
)
})
test_that("include_reference, different contrasts", {
data("mtcars")
mtcars$cyl <- factor(mtcars$cyl)
mtcars$gear <- factor(mtcars$gear)
m <- lm(mpg ~ cyl + gear, data = mtcars, contrasts = list(cyl = datawizard::contr.deviation))
out <- model_parameters(m, include_reference = TRUE)
expect_snapshot(print(out))
m <- lm(mpg ~ cyl + gear, data = mtcars)
out <- model_parameters(m, include_reference = TRUE)
expect_snapshot(print(out))
m <- lm(
mpg ~ cyl + gear,
data = mtcars,
contrasts = list(
cyl = datawizard::contr.deviation,
gear = contr.sum
)
)
out <- model_parameters(m, include_reference = TRUE)
expect_snapshot(print(out))
m <- lm(
mpg ~ cyl + gear,
data = mtcars,
contrasts = list(
cyl = contr.SAS,
gear = contr.sum
)
)
out <- model_parameters(m, include_reference = TRUE)
expect_snapshot(print(out))
m <- lm(
mpg ~ cyl + gear,
data = mtcars,
contrasts = list(
cyl = contr.SAS,
gear = contr.treatment
)
)
out <- model_parameters(m, include_reference = TRUE)
expect_snapshot(print(out))
})
|
