diff --git a/CRAN-SUBMISSION b/CRAN-SUBMISSION
deleted file mode 100644
index f976aad..0000000
--- a/CRAN-SUBMISSION
+++ /dev/null
@@ -1,3 +0,0 @@
-Version: 0.1.2
-Date: 2023-11-27 07:58:31 UTC
-SHA: b1c8a4f556e73869245052017d7cb49f7ba21742
diff --git a/R/seasonal_burden_levels.R b/R/seasonal_burden_levels.R
index 3ccfef4..1c373f0 100644
--- a/R/seasonal_burden_levels.R
+++ b/R/seasonal_burden_levels.R
@@ -78,6 +78,7 @@ seasonal_burden_levels <- function(
   disease_threshold = 20,
   n_peak = 6,
   only_current_season = TRUE,
+  family = NULL,
   ...
 ) {
   # Check input arguments
@@ -127,7 +128,7 @@ seasonal_burden_levels <- function(
     # Run quantiles_fit function
     quantiles_fit <- weighted_seasonal_tsd |>
       dplyr::select("observation", "weight") |>
-      fit_quantiles(weighted_observations = _, conf_levels = conf_levels, ...)
+      fit_quantiles(weighted_observations = _, conf_levels = conf_levels, family = family, ...)
 
     # If method intensity_levels was chosen; use the high level from the `fit_quantiles` function as the high
     # level and the disease_threshold as the very low level. The low and medium levels are defined as the relative