Nothing
#' @templateVar class nlrq
#' @template title_desc_tidy
#'
#' @param x A `nlrq` object returned from [quantreg::nlrq()].
#' @template param_confint
#' @template param_unused_dots
#'
#' @evalRd return_tidy(regression = TRUE)
#'
#' @aliases nlrq_tidiers
#' @export
#' @seealso [tidy()], [quantreg::nlrq()]
#' @family quantreg tidiers
#' @examplesIf rlang::is_installed("quantreg")
#'
#' # load modeling library
#' library(quantreg)
#'
#' # build artificial data with multiplicative error
#' set.seed(1)
#' dat <- NULL
#' dat$x <- rep(1:25, 20)
#' dat$y <- SSlogis(dat$x, 10, 12, 2) * rnorm(500, 1, 0.1)
#'
#' # fit the median using nlrq
#' mod <- nlrq(y ~ SSlogis(x, Asym, mid, scal),
#' data = dat, tau = 0.5, trace = TRUE
#' )
#'
#' # summarize model fit with tidiers
#' tidy(mod)
#' glance(mod)
#' augment(mod)
#'
tidy.nlrq <- function(x, conf.int = FALSE, conf.level = 0.95, ...) {
check_ellipses("exponentiate", "tidy", "nlrq", ...)
ret <- as_tidy_tibble(
coef(summary(x)),
new_names = c("estimate", "std.error", "statistic", "p.value")
)
if (conf.int) {
x_summary <- summary(x)
a <- (1 - conf.level) / 2
cv <- qt(c(a, 1 - a), df = x_summary[["rdf"]])
ret[["conf.low"]] <- ret[["estimate"]] + (cv[1] * ret[["std.error"]])
ret[["conf.high"]] <- ret[["estimate"]] + (cv[2] * ret[["std.error"]])
}
as_tibble(ret)
}
#' @templateVar class nlrq
#' @template title_desc_glance
#'
#' @inherit tidy.nlrq params examples
#'
#' @evalRd return_glance(
#' "tau",
#' "logLik",
#' "AIC",
#' "BIC",
#' "df.residual"
#' )
#'
#' @export
#' @seealso [glance()], [quantreg::nlrq()]
#' @family quantreg tidiers
#' @inherit tidy.nlrq examples
glance.nlrq <- function(x, ...) {
n <- length(x[["m"]]$fitted())
s <- summary(x)
as_glance_tibble(
tau = x[["m"]]$tau(),
logLik = logLik(x),
AIC = AIC(x),
BIC = AIC(x, k = log(n)),
df.residual = s[["rdf"]],
na_types = "rrrri"
)
}
#' @templateVar class nlrq
#' @template title_desc_tidy
#'
#' @param x A `nlrq` object returned from [quantreg::nlrq()].
#' @inherit augment.nls params examples return
#'
#' @export
#' @seealso [augment()], [quantreg::nlrq()]
#' @family quantreg tidiers
#' @include stats-nls-tidiers.R
#'
augment.nlrq <- augment.nls
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.