# Create a `tsd` object
tsd_data <- generate_seasonal_data()
# Create a `tsd_onset` object
tsd_onset <- seasonal_onset(
tsd = tsd_data,
k = 3,
disease_threshold = 100,
season_start = 21,
season_end = 20,
level = 0.95,
family = "poisson",
only_current_season = TRUE
)
# Print the summary
summary(tsd_onset)
Run the code above in your browser using DataLab