Nothing
#' Pseudo Likelihood Ratio Test under the
#' zero-inflated Poisson model with relative reporting rate parametrization
#' @inheritParams pvlrt
#' @param ... additional arguments passed to pvlrt
#'
#' @seealso \link{pvlrt}
#'
#' @note
#' \code{lrt_zi_poisson()} is a wrapper for \code{pvlrt()} with
#' \code{parametrization = "rrr"}.
#'
#' @returns
#'
#' Returns a `pvlrt` object. See \link{pvlrt} for more details.
#'
#'
#' @examples
#'
#' data("statin46")
#'
#' # 500 bootstrap iterations (nsim) in the example below
#' # are for quick demonstration only --
#' # we recommended setting nsim to 10000 (default) or bigger
#' test1 <- lrt_zi_poisson(statin46, nsim = 500)
#' test1
#'
#' @md
#' @export
lrt_zi_poisson <- function(contin_table,
nsim = 1e4,
...) {
dots <- list(...)
if (!is.null(dots$parametrization)) {
msg <- glue::glue(
"parametrization' is automatically set to 'rrr' in \\
`lrt_zi_poisson()`. For 'rr' parametrization use \\
pvlrt() or lrt_poisson()"
)
warning(msg)
}
dots$parametrization <- "rrr"
inargs <- list(
contin_table = contin_table,
nsim = nsim
) %>%
c(dots)
out <- tryCatch(
do.call(pvlrt, inargs),
error = function(e) e
)
if (is(out, "error")) {
stop(out$message)
}
out
}
#' Likelihood Ratio Test under the (vanilla, non-zero-inflated) Poisson model
#' @inheritParams pvlrt
#' @inheritParams lrt_zi_poisson
#'
#' @seealso \link{pvlrt}
#' @note
#' \code{lrt_poisson()} and \code{lrt_vanilla_poisson()}
#' are both wrappers for \code{pvlrt()} with
#' \code{omega_vec = rep(0, ncol(contin_table))}
#'
#'
#' @examples
#'
#' data("statin46")
#'
#' # 500 bootstrap iterations (nsim) in the example below
#' # are for quick demonstration only --
#' # we recommended setting nsim to 10000 (default) or bigger
#'
#' # no grouping -- each drug forms its own class
#' test1 <- lrt_poisson(lovastatin, nsim = 500)
#'
#'
#' @returns
#'
#' Returns a `pvlrt` object. See \link{pvlrt} for more details.
#'
#' @md
#' @export
lrt_poisson <- function(contin_table,
nsim = 1e4,
parametrization = "rrr",
...) {
dots <- list(...)
for (nm in c("zi_prob", "omega_vec", "omega_est_vec")) {
if (!is.null(dots[[nm]])) {
msg <- glue::glue(
"{nm}' is automatically set to 0 in `lrt_poisson()`. \\
For non-zero zi probability use `pvlrt()` or `lrt_zi_poisson()` \\
and supply as '{nm}'"
)
warning(msg)
dots$omega_vec <- dots$omega_est_vec <- dots$zi_prob <- NULL
}
}
if (parametrization %in% c("rrr", "lambda")) {
dots$omega_vec <- rep(0, ncol(contin_table))
}
inargs <- list(
contin_table = contin_table,
nsim = nsim,
parametrization = parametrization
) %>%
c(dots)
out <- tryCatch(
do.call(pvlrt, inargs),
error = function(e) e
)
if (is(out, "error")) {
stop(out$message)
}
out
}
#' @rdname lrt_poisson
#' @export
lrt_vanilla_poisson <- lrt_poisson
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.