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 106 107 108 109 110 111 112
|
# simulate_model -----------------------------------------
#' @export
simulate_model.svyglm.nb <- simulate_model.default
#' @export
simulate_model.svyglm.zip <- simulate_model.default
# standard erors -----------------------------------------
#' @export
standard_error.svyglm.nb <- function(model, ...) {
if (!isNamespaceLoaded("survey")) {
requireNamespace("survey", quietly = TRUE)
}
se <- sqrt(diag(stats::vcov(model, stderr = "robust")))
.data_frame(
Parameter = .remove_backticks_from_string(names(se)),
SE = as.vector(se)
)
}
#' @export
standard_error.svyglm.zip <- standard_error.svyglm.nb
#' @export
standard_error.svyglm <- function(model, ...) {
cs <- stats::coef(summary(model))
se <- cs[, 2]
.data_frame(
Parameter = .remove_backticks_from_string(names(se)),
SE = as.vector(se)
)
}
# confidence intervals -----------------------------------
#' @export
ci.svyglm.nb <- ci.tobit
#' @export
ci.svyglm.glimML <- ci.tobit
#' @export
ci.svyglm.zip <- ci.tobit
# p values -----------------------------------------------
#' @export
p_value.svyglm <- function(model, verbose = TRUE, ...) {
cs <- stats::coef(summary(model))
if (ncol(cs) < 4) {
if (isTRUE(verbose)) {
warning("Could not retrieve p-values.", call. = FALSE)
}
return(NULL)
}
p <- cs[, 4]
.data_frame(
Parameter = .remove_backticks_from_string(rownames(cs)),
p = as.vector(p)
)
}
#' @export
p_value.svyolr <- function(model, verbose = TRUE, ...) {
cs <- stats::coef(summary(model))
if (ncol(cs) < 3) {
if (isTRUE(verbose)) {
warning("Could not retrieve p-values.", call. = FALSE)
}
return(NULL)
}
p <- 2 * stats::pt(abs(cs[, 3]), df = degrees_of_freedom(model, method = "any"), lower.tail = FALSE)
.data_frame(
Parameter = .remove_backticks_from_string(rownames(cs)),
p = as.vector(p)
)
}
#' @export
p_value.svyglm.nb <- function(model, ...) {
if (!isNamespaceLoaded("survey")) {
requireNamespace("survey", quietly = TRUE)
}
est <- stats::coef(model)
se <- sqrt(diag(stats::vcov(model, stderr = "robust")))
p <- 2 * stats::pt(abs(est / se), df = degrees_of_freedom(model, method = "any"), lower.tail = FALSE)
.data_frame(
Parameter = .remove_backticks_from_string(names(p)),
p = as.vector(p)
)
}
#' @export
p_value.svyglm.zip <- p_value.svyglm.nb
|