set.seed(345)
# Create an example `tsd` object
time_series <- generate_seasonal_data()
autoplot(time_series)
# Create an `tsd_onset` object
time_series_with_onset <- seasonal_onset(
tsd = time_series,
k = 3,
level = 0.95,
family = "quasipoisson"
)
autoplot(time_series_with_onset)
# Define `disease_threshold`
disease_threshold <- 150
# Create a `tsd_onset_and_burden` object
tsd_onset_burden <- combined_seasonal_output(
tsd = time_series,
disease_threshold = disease_threshold
)
autoplot(tsd_onset_burden)
# Create an `tsd_onset` object
tsd_onset <- seasonal_onset(
tsd = time_series,
k = 5,
family = "quasipoisson",
season_start = 21,
only_current_season = FALSE
)
tsd_growth_warning <- consecutive_growth_warnings(tsd_onset)
autoplot(tsd_growth_warning)
Run the code above in your browser using DataLab