From 192bd149aae88f56874eddc4e7c68b221c181c51 Mon Sep 17 00:00:00 2001 From: Bill Denney Date: Thu, 20 May 2021 09:39:42 -0400 Subject: [PATCH] Remove boot_fit_nb() --- R/add_ci_negbin.R | 15 ++++----------- man/add_ci.negbin.Rd | 4 +++- 2 files changed, 7 insertions(+), 12 deletions(-) diff --git a/R/add_ci_negbin.R b/R/add_ci_negbin.R index 1ff188e..674e7d0 100644 --- a/R/add_ci_negbin.R +++ b/R/add_ci_negbin.R @@ -69,7 +69,9 @@ #' fit <- MASS::glm.nb(y ~ x1, data = df) #' df <- df[sample(100),] #' add_ci(df, fit, names = c("lcb", "ucb")) -#' +#' \dontrun{ +#' add_ci(df, fit, names = c("lcb", "ucb"), type = "boot") +#' } #' @export add_ci.negbin <- function(df, fit, alpha = 0.05, names = NULL, yhatName = "pred", @@ -123,15 +125,6 @@ parametric_ci_negbin <- function(df, fit, alpha, names, yhatName, response){ data.frame(df) } - -boot_fit_nb<- function(data, df, fit, lvl, indices){ - data_temp <- data[indices,] - form <- fit$call[2] - ##temp_fit <- glm.nb(form, data = data_temp) - temp_fit <- update(fit, data = data_temp) - predict(temp_fit, newdata = df, type = lvl) -} - boot_ci_negbin <- function(df, fit, alpha, names, yhatName, response, nSims){ if (response){ @@ -144,7 +137,7 @@ boot_ci_negbin <- function(df, fit, alpha, names, yhatName, response, nSims){ out <- predict(fit, df, type = lvl) boot_obj <- boot(data = fit$model, - statistic = boot_fit_nb, + statistic = boot_fit, R = nSims, fit = fit, df = df, diff --git a/man/add_ci.negbin.Rd b/man/add_ci.negbin.Rd index 5114928..fce3c74 100644 --- a/man/add_ci.negbin.Rd +++ b/man/add_ci.negbin.Rd @@ -74,7 +74,9 @@ df <- data.frame(x1 = x1, y = y) fit <- MASS::glm.nb(y ~ x1, data = df) df <- df[sample(100),] add_ci(df, fit, names = c("lcb", "ucb")) - +\dontrun{ +add_ci(df, fit, names = c("lcb", "ucb"), type = "boot") +} } \seealso{ \code{\link{add_pi.negbin}} for prediction intervals for