Added the fit_growth_rate function

DESCRIPTION

@@ -2,11 +2,10 @@ Package: aedseo
 Title: What the Package Does (One Line, Title Case)
 Version: 0.0.0.9000
 Authors@R: 
-    person("First", "Last", , "first.last@example.com", role = c("aut", "cre"),
-           comment = c(ORCID = "YOUR-ORCID-ID"))
+    person("Kasper", "Schou Telkamp", ,"ksst@ssi.dk", role = c("aut", "cre"),
+           comment = c(ORCID = "0009-0001-5126-0190"))
 Description: What the package does (one paragraph).
 License: MIT + file LICENSE
-Author@R: person("Kasper", "Schou Telkamp", email = "KSST@ssi.dk")
 Encoding: UTF-8
 Roxygen: list(markdown = TRUE)
 RoxygenNote: 7.2.3
@@ -19,5 +18,6 @@ Imports:
     base,
     magrittr,
     MASS,
+    rlang,
     stats,
     tibble

R/aedseo-package.R

@@ -1,4 +1,5 @@
 #' @keywords internal
+#' @aliases aedseo-package
 "_PACKAGE"
 
 ## usethis namespace: start

R/fit_growth_rate.R

@@ -0,0 +1,42 @@
+
+
+
+fit_growth_rate <- function(
+    observations,
+    level = 0.95,
+    family = c(
+      stats::poisson(link = "log"),
+      stats::quasipoisson(link = "log")))
+  {
+
+  # Calculate the length of the series
+  n <- base::length(observations)
+
+  # Construct a data.frame wit growth rate data
+  growth_data <- stats::aggregate(
+    observations,
+    by = list(growth_rate = rep(1:n, each=1)),
+    FUN = sum)
+
+  # Fit the model using the specified family
+  growth_fit <- stats::glm(
+    formula = x ~ growth_rate,
+    data = growth_data,
+    family = family)
+
+  # Profile the model
+  growth_profile <- stats::profile(fitted = growth_fit)
+
+  # Collect the estimates
+  ans <- c(
+    stats::coef(object = growth_fit)["growth_rate"],
+    stats::confint(object = growth_profile,
+                   parm = "growth_rate",
+                   level = level))
+
+  return(list(
+    estimate = ans,
+    level = level,
+    family = family))
+
+}