diff --git a/NAMESPACE b/NAMESPACE
index 0364c3a2..cbc881ac 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -105,6 +105,8 @@ export(util_chisquare_stats_tbl)
 export(util_exponential_aic)
 export(util_exponential_param_estimate)
 export(util_exponential_stats_tbl)
+export(util_f_aic)
+export(util_f_param_estimate)
 export(util_f_stats_tbl)
 export(util_gamma_aic)
 export(util_gamma_param_estimate)
diff --git a/NEWS.md b/NEWS.md
index 40d494fb..c62461cc 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -8,6 +8,7 @@ None
 2. #470 - Add function `util_ztn_binomial_param_estimate()` and `util_rztnbinom_aic()` to estimate the parameters and calculate the AIC for the zero-truncated negative binomial distribution. Also added `util_ztn_binomial_stats_tbl()`
 3. #467 - Add function `util_zero_truncated_poisson_param_estimate()` to estimate
 the parameters of the zero-truncated Poisson distribution.
+4. #472 - Add function `util_f_param_estimate()` and `util_f_aic()` to estimate the parameters and calculate the AIC for the F distribution.
 
 ## Minor Improvements and Fixes
 1. Fix #468 - Update `util_negative_binomial_param_estimate()` to add the use of
diff --git a/R/est-param-f.R b/R/est-param-f.R
new file mode 100644
index 00000000..fed25844
--- /dev/null
+++ b/R/est-param-f.R
@@ -0,0 +1,137 @@
+#' Estimate F Distribution Parameters
+#'
+#' @family Parameter Estimation
+#' @family F Distribution
+#'
+#' @author Steven P. Sanderson II, MPH
+#'
+#' @details This function will attempt to estimate the F distribution parameters
+#' given some vector of values produced by `rf()`. The estimation method
+#' is from the NIST Engineering Statistics Handbook.
+#'
+#' @param .x The vector of data to be passed to the function, where the data
+#' comes from the `rf()` function.
+#' @param .auto_gen_empirical This is a boolean value of TRUE/FALSE with default
+#' set to TRUE. This will automatically create the `tidy_empirical()` output
+#' for the `.x` parameter and use the `tidy_combine_distributions()`. The user
+#' can then plot out the data using `$combined_data_tbl` from the function output.
+#'
+#' @examples
+#' library(dplyr)
+#' library(ggplot2)
+#'
+#' set.seed(123)
+#' x <- rf(100, df1 = 5, df2 = 10, ncp = 1)
+#' output <- util_f_param_estimate(x)
+#'
+#' output$parameter_tbl
+#'
+#' output$combined_data_tbl |>
+#'   tidy_combined_autoplot()
+#'
+#' @return
+#' A tibble/list
+#'
+#' @export
+#'
+
+util_f_param_estimate <- function(.x, .auto_gen_empirical = TRUE) {
+
+  # Tidyeval ----
+  x_term <- as.numeric(.x)
+  n <- length(x_term)
+  minx <- min(as.numeric(x_term))
+  maxx <- max(as.numeric(x_term))
+  m <- mean(as.numeric(x_term))
+  s <- var(x_term)
+
+  # Checks ----
+  if (!inherits(x_term, "numeric")) {
+    rlang::abort(
+      message = "The '.x' parameter must be numeric.",
+      use_cli_format = TRUE
+    )
+  }
+
+  # NIST MME ----
+  m1 <- mean(x_term)
+  m2 <- mean(x_term^2)
+  m3 <- mean(x_term^3)
+  m4 <- mean(x_term^4)
+
+  b1 <- (m3 - m1 * m2) / (m2 - m1^2)
+  b2 <- (m4 - 4 * m1 * m3 + 6 * m1^2 * m2 - 3 * m1^4) / (m2 - m1^2)^2
+
+  df1_mme <- 2 * (2 * b2 - 3 * b1 - 6) / (b1 + 6)
+  df2_mme <- 4 + 2 * b1 * df1_mme / (df1_mme - 2)
+  ncp_mme <- (df1_mme * m1 - df1_mme + 2) / df2_mme
+
+  # Round
+  df1_mme <- ifelse(round(df1_mme, 3) <= 0, abs(round(df1_mme, 3)), round(df1_mme, 3))
+  df2_mme <- ifelse(round(df2_mme, 3) <= 0, abs(round(df2_mme, 3)), round(df2_mme, 3))
+  ncp_mme <- ifelse(round(ncp_mme, 3) <= 0, abs(round(ncp_mme, 3)), round(ncp_mme, 3))
+
+  # Negative Log Likelihood ----
+  # Negative log-likelihood function for the F-distribution
+  nll <- function(params) {
+    df1 <- params[1]
+    df2 <- params[2]
+    ncp <- params[3]
+    if (df1 <= 0 || df2 <= 0 || ncp <= 0) return(Inf) # return Inf if params are not valid
+    -sum(stats::df(x_term, df1, df2, ncp, log = TRUE))
+  }
+
+  # Initial parameter guesses
+  start_params <- c(df1 = 1, df2 = 1, ncp = 0)
+
+  # Use optim to minimize the negative log-likelihood
+  optim_res <- stats::optim(start_params, nll, method = "L-BFGS-B",
+                            lower = c(1e-6, 1e-6, 1e-6))
+
+  # Return the estimated parameters
+  optim_df1 <- round(optim_res$par[[1]], 3)
+  optim_df2 <- round(optim_res$par[[2]], 3)
+  optim_ncp <- round(optim_res$par[[3]], 3)
+
+  # Return Tibble ----
+  if (.auto_gen_empirical) {
+    te <- tidy_empirical(.x = x_term)
+    td_mme <- tidy_f(.n = n,
+                     .df1 = df1_mme, .df2 = df2_mme, .ncp = ncp_mme)
+    td_optim <- tidy_f(.n = n,
+                       .df1 = optim_df1, .df2 = optim_df2, .ncp = optim_ncp)
+    combined_tbl <- tidy_combine_distributions(te, td_mme, td_optim)
+  }
+
+  ret <- dplyr::tibble(
+    dist_type = rep("F Distribution", 2),
+    samp_size = rep(n, 2),
+    min       = minx,
+    max       = maxx,
+    mean      = m,
+    variance  = s,
+    method    = c("MME", "MLE"),
+    df1_est   = c(df1_mme, optim_df1),
+    df2_est   = c(df2_mme, optim_df2),
+    ncp_est   = c(ncp_mme, optim_ncp)
+  )
+
+  # Return ----
+  attr(ret, "tibble_type") <- "parameter_estimation"
+  attr(ret, "family") <- "f_distribution"
+  attr(ret, "x_term") <- .x
+  attr(ret, "n") <- length(x_term)
+
+  if (.auto_gen_empirical) {
+    output <- list(
+      combined_data_tbl = combined_tbl,
+      parameter_tbl     = ret
+    )
+  } else {
+    output <- list(
+      parameter_tbl = ret
+    )
+  }
+
+  return(output)
+}
diff --git a/R/utils-aic-f.R b/R/utils-aic-f.R
new file mode 100644
index 00000000..c70a8a2c
--- /dev/null
+++ b/R/utils-aic-f.R
@@ -0,0 +1,68 @@
+#' Calculate Akaike Information Criterion (AIC) for F Distribution
+#'
+#' This function calculates the Akaike Information Criterion (AIC) for an F
+#' distribution fitted to the provided data.
+#'
+#' @param .x A numeric vector containing the data to be fitted to an F
+#' distribution.
+#'
+#' @return The AIC value calculated based on the fitted F distribution to the
+#' provided data.
+#'
+#' @details
+#' This function fits an F distribution to the input data using maximum
+#' likelihood estimation and then computes the Akaike Information Criterion (AIC)
+#' based on the fitted distribution.
+#'
+#' @examples
+#' # Generate F-distributed data
+#' set.seed(123)
+#' x <- rf(100, df1 = 5, df2 = 10, ncp = 1)
+#'
+#' # Calculate AIC for the generated data
+#' util_f_aic(x)
+#'
+#' @seealso
+#' \code{\link[stats]{rf}} for generating F-distributed data,
+#' \code{\link[stats]{optim}} for optimization.
+#'
+#' @name util_f_aic
+NULL
+
+#' @export
+#' @rdname util_f_aic
+
+util_f_aic <- function(.x) {
+  # Tidyeval
+  x <- as.numeric(.x)
+
+  # Get initial parameter estimates
+  pe <- TidyDensity::util_f_param_estimate(x)$parameter_tbl |>
+    dplyr::filter(method == "MLE")
+
+  # Negative log-likelihood function for the F-distribution
+  nll <- function(params) {
+    df1 <- params[1]
+    df2 <- params[2]
+    ncp <- params[3]
+    if (df1 <= 0 || df2 <= 0 || ncp <= 0) return(Inf) # return Inf if params are not valid
+    -sum(stats::df(x_term, df1, df2, ncp, log = TRUE))
+  }
+
+  # Initial parameter guesses
+  start_params <- c(df1 = 1, df2 = 1, ncp = 0)
+
+  # Use optim to minimize the negative log-likelihood
+  fit_f <- stats::optim(start_params, nll, method = "L-BFGS-B",
+                        lower = c(1e-6, 1e-6, 1e-6))
+
+  # Extract log-likelihood and number of parameters
+  logLik_f <- -fit_f$value
+  k_f <- 3 # Number of parameters for F distribution (df1, df2, ncp)
+
+  # Calculate AIC
+  AIC_f <- 2 * k_f - 2 * logLik_f
+
+  # Return AIC
+  return(AIC_f)
+}
diff --git a/man/tidy_f.Rd b/man/tidy_f.Rd
index 384b9c02..ccd66f74 100644
--- a/man/tidy_f.Rd
+++ b/man/tidy_f.Rd
@@ -90,6 +90,7 @@ Other Continuous Distribution:
 \code{\link{tidy_zero_truncated_geometric}()}
 
 Other F Distribution: 
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_f_stats_tbl}()}
 }
 \author{
diff --git a/man/util_bernoulli_param_estimate.Rd b/man/util_bernoulli_param_estimate.Rd
index 34dfc505..4d5084d3 100644
--- a/man/util_bernoulli_param_estimate.Rd
+++ b/man/util_bernoulli_param_estimate.Rd
@@ -50,6 +50,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_beta_param_estimate.Rd b/man/util_beta_param_estimate.Rd
index 0919670c..2b7bfb39 100644
--- a/man/util_beta_param_estimate.Rd
+++ b/man/util_beta_param_estimate.Rd
@@ -63,6 +63,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_binomial_param_estimate.Rd b/man/util_binomial_param_estimate.Rd
index e890aa9b..b01dd5e5 100644
--- a/man/util_binomial_param_estimate.Rd
+++ b/man/util_binomial_param_estimate.Rd
@@ -53,6 +53,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_burr_param_estimate.Rd b/man/util_burr_param_estimate.Rd
index 421f03af..85f62900 100644
--- a/man/util_burr_param_estimate.Rd
+++ b/man/util_burr_param_estimate.Rd
@@ -50,6 +50,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_cauchy_param_estimate.Rd b/man/util_cauchy_param_estimate.Rd
index 4b5783e7..219511c9 100644
--- a/man/util_cauchy_param_estimate.Rd
+++ b/man/util_cauchy_param_estimate.Rd
@@ -48,6 +48,7 @@ Other Parameter Estimation:
 \code{\link{util_burr_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_chisquare_param_estimate.Rd b/man/util_chisquare_param_estimate.Rd
index f3fd38ca..cb283c90 100644
--- a/man/util_chisquare_param_estimate.Rd
+++ b/man/util_chisquare_param_estimate.Rd
@@ -79,6 +79,7 @@ Other Parameter Estimation:
 \code{\link{util_burr_param_estimate}()},
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_exponential_param_estimate.Rd b/man/util_exponential_param_estimate.Rd
index 1f6809b0..c345161c 100644
--- a/man/util_exponential_param_estimate.Rd
+++ b/man/util_exponential_param_estimate.Rd
@@ -48,6 +48,7 @@ Other Parameter Estimation:
 \code{\link{util_burr_param_estimate}()},
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_f_aic.Rd b/man/util_f_aic.Rd
new file mode 100644
index 00000000..53553af9
--- /dev/null
+++ b/man/util_f_aic.Rd
@@ -0,0 +1,38 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/utils-aic-f.R
+\name{util_f_aic}
+\alias{util_f_aic}
+\title{Calculate Akaike Information Criterion (AIC) for F Distribution}
+\usage{
+util_f_aic(.x)
+}
+\arguments{
+\item{.x}{A numeric vector containing the data to be fitted to an F
+distribution.}
+}
+\value{
+The AIC value calculated based on the fitted F distribution to the
+provided data.
+}
+\description{
+This function calculates the Akaike Information Criterion (AIC) for an F
+distribution fitted to the provided data.
+}
+\details{
+This function fits an F distribution to the input data using maximum
+likelihood estimation and then computes the Akaike Information Criterion (AIC)
+based on the fitted distribution.
+}
+\examples{
+# Generate F-distributed data
+set.seed(123)
+x <- rf(100, df1 = 5, df2 = 10, ncp = 1)
+
+# Calculate AIC for the generated data
+util_f_aic(x)
+
+}
+\seealso{
+\code{\link[stats]{rf}} for generating F-distributed data,
+\code{\link[stats]{optim}} for optimization.
+}
diff --git a/man/util_f_param_estimate.Rd b/man/util_f_param_estimate.Rd
new file mode 100644
index 00000000..c1caaad2
--- /dev/null
+++ b/man/util_f_param_estimate.Rd
@@ -0,0 +1,75 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/est-param-f.R
+\name{util_f_param_estimate}
+\alias{util_f_param_estimate}
+\title{Estimate F Distribution Parameters}
+\usage{
+util_f_param_estimate(.x, .auto_gen_empirical = TRUE)
+}
+\arguments{
+\item{.x}{The vector of data to be passed to the function, where the data
+comes from the \code{rf()} function.}
+
+\item{.auto_gen_empirical}{This is a boolean value of TRUE/FALSE with default
+set to TRUE. This will automatically create the \code{tidy_empirical()} output
+for the \code{.x} parameter and use the \code{tidy_combine_distributions()}. The user
+can then plot out the data using \verb{$combined_data_tbl} from the function output.}
+}
+\value{
+A tibble/list
+}
+\description{
+Estimate F Distribution Parameters
+}
+\details{
+This function will attempt to estimate the F distribution parameters
+given some vector of values produced by \code{rf()}. The estimation method
+is from the NIST Engineering Statistics Handbook.
+}
+\examples{
+library(dplyr)
+library(ggplot2)
+
+set.seed(123)
+x <- rf(100, df1 = 5, df2 = 10, ncp = 1)
+output <- util_f_param_estimate(x)
+
+output$parameter_tbl
+
+output$combined_data_tbl |>
+  tidy_combined_autoplot()
+
+}
+\seealso{
+Other Parameter Estimation: 
+\code{\link{util_bernoulli_param_estimate}()},
+\code{\link{util_beta_param_estimate}()},
+\code{\link{util_binomial_param_estimate}()},
+\code{\link{util_burr_param_estimate}()},
+\code{\link{util_cauchy_param_estimate}()},
+\code{\link{util_chisquare_param_estimate}()},
+\code{\link{util_exponential_param_estimate}()},
+\code{\link{util_gamma_param_estimate}()},
+\code{\link{util_geometric_param_estimate}()},
+\code{\link{util_hypergeometric_param_estimate}()},
+\code{\link{util_logistic_param_estimate}()},
+\code{\link{util_lognormal_param_estimate}()},
+\code{\link{util_negative_binomial_param_estimate}()},
+\code{\link{util_normal_param_estimate}()},
+\code{\link{util_pareto_param_estimate}()},
+\code{\link{util_poisson_param_estimate}()},
+\code{\link{util_triangular_param_estimate}()},
+\code{\link{util_uniform_param_estimate}()},
+\code{\link{util_weibull_param_estimate}()},
+\code{\link{util_zero_truncated_poisson_param_estimate}()},
+\code{\link{util_ztn_binomial_param_estimate}()}
+
+Other F Distribution: 
+\code{\link{tidy_f}()},
+\code{\link{util_f_stats_tbl}()}
+}
+\author{
+Steven P. Sanderson II, MPH
+}
+\concept{F Distribution}
+\concept{Parameter Estimation}
diff --git a/man/util_f_stats_tbl.Rd b/man/util_f_stats_tbl.Rd
index 2e4d4dbd..0320ce67 100644
--- a/man/util_f_stats_tbl.Rd
+++ b/man/util_f_stats_tbl.Rd
@@ -30,7 +30,8 @@ tidy_f() |>
 }
 \seealso{
 Other F Distribution: 
-\code{\link{tidy_f}()}
+\code{\link{tidy_f}()},
+\code{\link{util_f_param_estimate}()}
 
 Other Distribution Statistics: 
 \code{\link{util_bernoulli_stats_tbl}()},
diff --git a/man/util_gamma_param_estimate.Rd b/man/util_gamma_param_estimate.Rd
index d312ad43..b9155eb6 100644
--- a/man/util_gamma_param_estimate.Rd
+++ b/man/util_gamma_param_estimate.Rd
@@ -49,6 +49,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
 \code{\link{util_logistic_param_estimate}()},
diff --git a/man/util_geometric_param_estimate.Rd b/man/util_geometric_param_estimate.Rd
index b27d5840..fda7e4e5 100644
--- a/man/util_geometric_param_estimate.Rd
+++ b/man/util_geometric_param_estimate.Rd
@@ -51,6 +51,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
 \code{\link{util_logistic_param_estimate}()},
diff --git a/man/util_hypergeometric_param_estimate.Rd b/man/util_hypergeometric_param_estimate.Rd
index 305dc352..175a208b 100644
--- a/man/util_hypergeometric_param_estimate.Rd
+++ b/man/util_hypergeometric_param_estimate.Rd
@@ -74,6 +74,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_logistic_param_estimate}()},
diff --git a/man/util_logistic_param_estimate.Rd b/man/util_logistic_param_estimate.Rd
index aa84b2fe..f6151695 100644
--- a/man/util_logistic_param_estimate.Rd
+++ b/man/util_logistic_param_estimate.Rd
@@ -59,6 +59,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_lognormal_param_estimate.Rd b/man/util_lognormal_param_estimate.Rd
index c4e6393f..013a1bb1 100644
--- a/man/util_lognormal_param_estimate.Rd
+++ b/man/util_lognormal_param_estimate.Rd
@@ -58,6 +58,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_negative_binomial_param_estimate.Rd b/man/util_negative_binomial_param_estimate.Rd
index 59d012e6..fc482958 100644
--- a/man/util_negative_binomial_param_estimate.Rd
+++ b/man/util_negative_binomial_param_estimate.Rd
@@ -65,6 +65,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_normal_param_estimate.Rd b/man/util_normal_param_estimate.Rd
index f0d2d382..33cd8e87 100644
--- a/man/util_normal_param_estimate.Rd
+++ b/man/util_normal_param_estimate.Rd
@@ -58,6 +58,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_pareto_param_estimate.Rd b/man/util_pareto_param_estimate.Rd
index 07b3dac9..efd74e72 100644
--- a/man/util_pareto_param_estimate.Rd
+++ b/man/util_pareto_param_estimate.Rd
@@ -58,6 +58,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_poisson_param_estimate.Rd b/man/util_poisson_param_estimate.Rd
index ecf8ea72..25a09c9f 100644
--- a/man/util_poisson_param_estimate.Rd
+++ b/man/util_poisson_param_estimate.Rd
@@ -52,6 +52,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_triangular_param_estimate.Rd b/man/util_triangular_param_estimate.Rd
index 576ae5f0..5c250e9b 100644
--- a/man/util_triangular_param_estimate.Rd
+++ b/man/util_triangular_param_estimate.Rd
@@ -57,6 +57,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_uniform_param_estimate.Rd b/man/util_uniform_param_estimate.Rd
index c20954c5..3abf36c2 100644
--- a/man/util_uniform_param_estimate.Rd
+++ b/man/util_uniform_param_estimate.Rd
@@ -49,6 +49,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_weibull_param_estimate.Rd b/man/util_weibull_param_estimate.Rd
index c4c03337..47d6a490 100644
--- a/man/util_weibull_param_estimate.Rd
+++ b/man/util_weibull_param_estimate.Rd
@@ -49,6 +49,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_zero_truncated_poisson_param_estimate.Rd b/man/util_zero_truncated_poisson_param_estimate.Rd
index 7286ed40..9e9c25c3 100644
--- a/man/util_zero_truncated_poisson_param_estimate.Rd
+++ b/man/util_zero_truncated_poisson_param_estimate.Rd
@@ -76,6 +76,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},
diff --git a/man/util_ztn_binomial_param_estimate.Rd b/man/util_ztn_binomial_param_estimate.Rd
index a1df1c99..e5db1139 100644
--- a/man/util_ztn_binomial_param_estimate.Rd
+++ b/man/util_ztn_binomial_param_estimate.Rd
@@ -59,6 +59,7 @@ Other Parameter Estimation:
 \code{\link{util_cauchy_param_estimate}()},
 \code{\link{util_chisquare_param_estimate}()},
 \code{\link{util_exponential_param_estimate}()},
+\code{\link{util_f_param_estimate}()},
 \code{\link{util_gamma_param_estimate}()},
 \code{\link{util_geometric_param_estimate}()},
 \code{\link{util_hypergeometric_param_estimate}()},