## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")

## ----psu-table, echo = FALSE--------------------------------------------------
psu_example <- data.frame(
  date     = rep(as.Date(c("2024-06-03", "2024-06-04", "2024-06-05")), each = 2),
  day_type = rep(c("weekday", "weekend"), times = 3),
  psu_id   = paste0("PSU-", seq_len(6))
)
knitr::kable(
  psu_example,
  col.names = c("Date", "Day type", "PSU"),
  caption   = "Six PSUs formed by crossing 3 days × 2 strata."
)

## ----single-count, message = FALSE--------------------------------------------
library(tidycreel)

calendar <- data.frame(
  date = as.Date(c("2024-06-03", "2024-06-04", "2024-06-05")),
  day_type = rep("weekday", 3)
)
design <- creel_design(calendar, date = date, strata = day_type)

# n_anglers holds angler-hours (count × open hours already combined)
counts <- data.frame(
  date      = as.Date(c("2024-06-03", "2024-06-04", "2024-06-05")),
  day_type  = rep("weekday", 3),
  n_anglers = c(150L, 230L, 190L)
)
design <- add_counts(design, counts)
result <- estimate_effort(design)
result$estimates

## ----multi-count, message = FALSE---------------------------------------------
calendar_m <- data.frame(
  date     = as.Date(c("2024-06-03", "2024-06-04", "2024-06-05", "2024-06-06")),
  day_type = rep("weekday", 4)
)
design_m <- creel_design(calendar_m, date = date, strata = day_type)

multi_counts <- data.frame(
  date = as.Date(rep(
    c("2024-06-03", "2024-06-04", "2024-06-05", "2024-06-06"),
    each = 2
  )),
  day_type = rep("weekday", 8),
  count_time = rep(c("am", "pm"), 4),
  n_anglers = c(10L, 20L, 20L, 30L, 30L, 40L, 40L, 50L)
)
design_m <- add_counts(design_m, multi_counts, count_time_col = count_time)
result_m <- estimate_effort(design_m)
result_m$estimates

## ----progressive, message = FALSE---------------------------------------------
calendar_p <- data.frame(
  date     = as.Date(c("2024-06-03", "2024-06-04", "2024-06-05")),
  day_type = rep("weekday", 3)
)
design_p <- creel_design(calendar_p, date = date, strata = day_type)

prog_counts <- data.frame(
  date        = as.Date(c("2024-06-03", "2024-06-04", "2024-06-05")),
  day_type    = rep("weekday", 3),
  n_anglers   = c(15L, 20L, 25L),
  shift_hours = rep(8, 3)
)
design_p <- add_counts(
  design_p, prog_counts,
  count_type = "progressive",
  circuit_time = 2,
  period_length_col = shift_hours
)

# Per-day Ê_d values stored after C × T expansion
design_p$counts

result_p <- estimate_effort(design_p)
result_p$estimates