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
|
#' @rdname p_value_betwithin
#' @export
dof_betwithin <- function(model) {
if (!insight::is_mixed_model(model)) {
insight::format_error("Model must be a mixed model.")
}
ngrps <- sum(.n_randomeffects(model))
parameters <- insight::find_parameters(model, effects = "fixed")[["conditional"]]
within_effects <- unlist(insight::find_random_slopes(model))
has_intcp <- insight::has_intercept(model)
ddf_within <- ngrps - n_parameters(model)
ddf_between <- insight::n_obs(model, disaggregate = TRUE) - ngrps - n_parameters(model)
if (has_intcp) {
ddf_between <- ddf_between - 1
ddf_within <- ddf_within - 1
}
within_index <- match(within_effects, parameters)
ddf <- stats::setNames(seq_along(parameters), parameters)
if (length(within_index) > 0) {
ddf[match(within_effects, parameters)] <- ddf_within
ddf[-match(within_effects, parameters)] <- ddf_between
} else {
ddf <- ddf_between
}
ddf
}
|
