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
|
#' @export
model_parameters.draws <- function(model,
centrality = "median",
dispersion = FALSE,
ci = 0.95,
ci_method = "eti",
test = "pd",
rope_range = "default",
rope_ci = 0.95,
exponentiate = FALSE,
keep = NULL,
drop = NULL,
verbose = TRUE,
...) {
out <- .posterior_draws_to_df(model)
# Processing
params <- .extract_parameters_bayesian(
out,
centrality = centrality,
dispersion = dispersion,
ci = ci,
ci_method = ci_method,
test = test,
rope_range = rope_range,
rope_ci = rope_ci,
bf_prior = NULL,
diagnostic = NULL,
priors = FALSE,
keep_parameters = keep,
drop_parameters = drop,
verbose = verbose,
...
)
# exponentiate coefficients and SE/CI, if requested
params <- .exponentiate_parameters(params, exponentiate = exponentiate)
attr(params, "ci") <- ci
attr(params, "object_name") <- insight::safe_deparse_symbol(substitute(model))
class(params) <- c("parameters_model", "see_parameters_model", class(params))
params
}
# Standard Errors ---------------------------------------------
#' @export
standard_error.draws <- function(model, verbose = TRUE, ...) {
params <- .posterior_draws_to_df(model)
.data_frame(
Parameter = colnames(params),
SE = unname(vapply(params, stats::sd, 1, na.rm = TRUE))
)
}
# p-Values ---------------------------------------------
#' @export
p_value.draws <- function(model, ...) {
params <- .posterior_draws_to_df(model)
p <- bayestestR::p_direction(params)
.data_frame(
Parameter = .remove_backticks_from_string(p$Parameter),
p = vapply(p$pd, bayestestR::convert_pd_to_p, 1)
)
}
# helper ------------------------------
.posterior_draws_to_df <- function(x) {
UseMethod(".posterior_draws_to_df")
}
.posterior_draws_to_df.default <- function(x) {
insight::format_error(sprintf("Objects of class `%s` are not yet supported.", class(x)[1]))
}
.posterior_draws_to_df.data.frame <- function(x) {
x
}
.posterior_draws_to_df.draws_df <- function(x) {
insight::check_if_installed("posterior")
datawizard::data_remove(as.data.frame(posterior::as_draws_df(x)), c(".chain", ".iteration", ".draw"))
}
.posterior_draws_to_df.draws_matrix <- .posterior_draws_to_df.draws_df
.posterior_draws_to_df.draws_array <- .posterior_draws_to_df.draws_df
.posterior_draws_to_df.draws_list <- .posterior_draws_to_df.draws_df
.posterior_draws_to_df.draws_rvars <- .posterior_draws_to_df.draws_df
|
