Fixes #138

NAMESPACE

@@ -59,6 +59,7 @@ export(util_beta_param_estimate)
 export(util_binomial_param_estimate)
 export(util_cauchy_param_estimate)
 export(util_cauchy_stats_tbl)
+export(util_chisquare_stats_tbl)
 export(util_exponential_param_estimate)
 export(util_f_stats_tbl)
 export(util_gamma_param_estimate)

NEWS.md

@@ -23,6 +23,7 @@ qq, quantile and probability plots to a single graph.
 16. Fix #141 - Add function `util_cauchy_stats_tbl()`
 17. Fix #140 - Add function `util_t_stats_tbl()`
 18. Fix #139 - Add function `util_f_stats_tbl()`
+19. Fix #138 - Add function `util_chisquare_stats_tbl()`
 
 ## Minor Fixes and Improvements
 1. Fix #110 - Bug fix, correct the `p` calculation in `tidy_poisson()` that will

R/stats-chisquare-tbl.R

@@ -0,0 +1,85 @@
+#' Distribution Statistics
+#'
+#' @family Chisquare
+#' @family Distribution Statistics
+#'
+#' @author Steven P. Sanderson II, MPH
+#'
+#' @details This function will take in a tibble and returns the statistics
+#' of the given type of `tidy_` distribution. It is required that data be
+#' passed from a `tidy_` distribution function.
+#'
+#' @description Returns distribution statistics in a tibble.
+#'
+#' @param .data The data being passed from a `tidy_` distribution function.
+#'
+#' @examples
+#' tidy_chisquare() %>%
+#'   util_chisquare_stats_tbl()
+#'
+#' @return
+#' A tibble
+#'
+#' @export
+#'
+
+util_chisquare_stats_tbl <- function(.data){
+
+    # Immediate check for tidy_ distribution function
+    if (!"tibble_type" %in% names(attributes(.data))){
+        rlang::abort(
+            message = "You must pass data from the 'tidy_dist' function.",
+            use_cli_format = TRUE
+        )
+    }
+
+    if (attributes(.data)$tibble_type != "tidy_chisquare"){
+        rlang::abort(
+            message = "You must use 'tidy_chisquare()'",
+            use_cli_format = TRUE
+        )
+    }
+
+    # Data
+    data_tbl <- tibble::as_tibble(.data)
+
+    atb <- attributes(data_tbl)
+
+    stat_mean   <- atb$.df
+    stat_median <- (atb$.df - 2/3)
+    stat_mode   <- ifelse(
+        atb$.df > 2,
+        atb$.df - 2,
+        "undefined"
+    )
+    stat_range = "0 to Inf"
+    stat_sd     <- sqrt(2 * atb$.df)
+    stat_coef_var <- sqrt(2/atb$.df)
+    stat_skewness <- ((2^1.5) / sqrt(atb$.df))
+    stat_kurtosis <- 3 + (12/atb$.df)
+
+    # Data Tibble
+    ret <- tibble::tibble(
+        tidy_function = atb$tibble_type,
+        function_call = atb$dist_with_params,
+        distribution = atb$tibble_type %>%
+            stringr::str_remove("tidy_") %>%
+            stringr::str_to_title(),
+        distribution_type = atb$distribution_family_type,
+        points = atb$.n,
+        simulations = atb$.num_sims,
+        mean = stat_mean,
+        median = stat_median,
+        mode = stat_mode,
+        std_dv = stat_sd,
+        coeff_var = stat_coef_var,
+        skewness = stat_skewness,
+        kurtosis = stat_kurtosis,
+        computed_std_skew = tidy_skewness_vec(data_tbl$y),
+        computed_std_kurt = tidy_kurtosis_vec(data_tbl$y)
+    )
+
+    # Return
+    return(ret)
+
+}