- The
aedseopackage now allows for modelling of binomial data.to_time_series()now acceptssamplesandproportionas arguments. Theproportionargument is the proportion ofcasesin the total testedsamplesat every time point. These variables can be used in theseasonal_onset()function with the family argumentsbinomialorquasibinomialand in theseasonal_burden_levels()function with the family argumentbeta. Attributes;outcome_type,model_outcomeandburden_outcomehave been added to the object classes ofto_time_series(),seasonal_onset()andseasonal_burden_levels()functions to always be able to track if the results are based on count/incidence or binomial data (#100, #101, #103). summary.tsd_onset()now reports proportional observations on the proportion scale and prints an explicitNAreference when no seasonal onset was detected, instead of producing no output (#103).
- The
estimate_disease_threshold()function defaulted to 1 if the selected observation was below 1. Since the package has the possibility to use incidence as an input there is a need for allowing values below 1. This was fixed in #106.
combined_seasonal_output()now includes aseasonal_offsetlogical variables in the output that estimates when the season has ended after the firstseasonal_onset. It can be controlled by the inputs:burden_level_decreaseandsteps_with_decrease(#93).
fit_growth_rate()would incorrectly determine confidence intervals when using ATLAS BLAS/LAPACK (#89).
-
Added new arguments
incidenceandincidence_denominatortoto_time_series()that allows the user to get output as incidence (#84). -
Added new argument
populationtoto_time_series()andfit_growth_rate()that allows the user to add the background population connected to each observation (#83). -
Added new argument
use_offsettoseasonal_onset()that allows the user to add the background population to adjust the growth rate estimations (#83). -
Added new feature to estimate multiple waves in
combined_seasonal_output(#77). -
Added
estimate_disease_threshold()for users to easier estimate the disease specific threshold (#85).
- Observations are now divided into
casesandincidence, which is implemented into all functions in the package. Cases are used as default, but if the user additionally inputspopulationthe output will be incidence (#84).
-
aedseo()is now deprecated. Please useseasonal_onset()instead. A warning is shown when usingaedseo()(#41). -
tsd()is now deprecated. Please useto_time_series()instead. A warning is shown when usingtsd()(#41).
-
Added the
seasonal_burden_levels()function, which calculates burden levels based on data from previous seasons with two different methods; "peak_levels" or "intensity_levels" (#37). -
Added the
fit_percentiles()function, which optimises a user selected distribution and calculates the percentiles based on observations and weights. It is meant to be used within theseasonal_burden_levels()function (#35, #37) - Renamedfit_quantiles()tofit_percentiles()(#60). -
Added
combined_seasonal_output()as the main function to run bothseasonal_onset()andseasonal_burden_levels()to get a combined result for the newest season (#44). -
Added
consecutive_growth_warnings()function to help the user with a method to define the disease-specific threshold (#80). -
Added a new argument
only_current_seasontoseasonal_onset(),seasonal_burden_levels()andcombined_seasonal_output()which gives the possibility to either get output from only the current season or for all available seasons (#45). -
Added
historical_summary()which uses atsd_onsetobject to summarise historical estimations (#75). -
summary()can now summarisetsd_burden_levelobjects (#60). -
plot()andautoplot()can now plottsd_combined_seasonal_outputandtsd_consecutive_growth_warningobjects (#57, #80). -
Added
generate_seasonal_data()to generate synthetic data for testing and documentation purposes (#56). -
Added
seasonal_onset()as a replacement for the deprecatedaedseo()function (#41). -
Added
to_time_series()as a replacement for the deprecatedtsd()function (#41).
-
Enhanced clarity and user guidance in the vignettes:
vignette("generate_seasonal_wave"),vignette("aedseo"),vignette("seasonal_onset")vignette("burden_levels")providing a comprehensive walkthrough of the application of the functions provided by theaedseopackage with detailed explanations and illustrative examples (#56, #57, #58, #59, #60, #61).
-
Improved the
autoplot()function which can now visualise dates as days, weeks and months on the x-axis with thetime_intervalargument (#56). -
Improved the
epi_calendar()function to work for a season spanning new year (#34). -
Using
predict()ontsd_onsetobjects now uses the same time-scale as the given object (#61). That is, thetime_intervalattribute controls if predictions are by "days", "weeks" or "months". -
The
aedseo()function now allows for the choice of adding season as an input argument (#34). -
{checkmate}assertions have been added to enhance user feedback with clearer error messages and to ensure functions operate correctly by validating inputs (#33). -
Improved the
aedseo()function to work withNAvalues. The user now defines how manyNAvalues the function should allow in each window (#32).
-
Added Sofia Myrup Otero as an author of the R package (#55).
-
Added Rasmus Skytte Randløv as a reviewer of the R package (#55).
-
The
disease_thresholdargument now reflects the disease threshold in one time step. If the total number of cases in a window of sizekexceedsdisease_threshold * k, a seasonal onset alarm can be triggered (#32).
- Transferring maintainership of the R package to Lasse Engbo Christiansen.
- Enhanced clarity and user guidance in the introductory vignette, providing a more comprehensive walkthrough of the application of the 'aeddo' algorithm on time series data with detailed explanations and illustrative examples.
-
Updated LICENSE.md to have Statens Serum Institut as a copyright holder.
-
Fixed installation guide for the development version in the README.Rmd and README.md
-
Added Lasse Engbo Christiansen as an author of the R package.
-
Added a new function
epi_calendar()that determines the epidemiological season based on a given date, allowing users to easily categorize dates within or outside specified seasons. -
Introduced additional visualizations in the
autoplot()method, enhancing the capabilities of theplot()method with new displays of observed cases and growth rates.
-
Added the
aedseofunction, which automates the early detection of seasonal epidemic onsets by estimating growth rates for consecutive time intervals and calculating the Sum of Cases (sum_of_cases). -
Introduced
autoplotandplotmethods for visualizingaedseoandaedseo_tsdobjects. These functions allow you to create insightful ggplot2 plots for your data. -
Included the
fit_growth_ratefunction, enabling users to fit growth rate models to time series observations. -
Introduced the
predictmethod foraedseoobjects, which allows you to predict observations for future time steps given the growth rates. -
Added the
summarymethod foraedseoobjects, providing a comprehensive summary of the results. -
Introduced the
tsdfunction, allowing users to create S3aedseo_tsd(time-series data) objects from observed data and corresponding dates.