This package utilises the Šesták–Berggren equation alongside the Arrhenius equation to make a simple and consistent way for a user to carry out the calculations and predictions required by accelerated stability studies. Currently the package works with decreasing variables, you may choose to transform your increasing variable into a decreasing one but note that your choice of transformation can impact the outcome.
The available functions within the package are as follows:
step1_down() Fit the one-step Šesták–Berggren kinetic
model.step1_plot_desc() Plot the stability data.step1_plot_pred() Plot the stability data and visualise
the predictions.step1_plot_CI() Plot the stability data and visualise
the predictions with confidence intervals.step1_plot_PI() Plot the stability data and visualise
the predictions with prediction intervals.step1_plot_T() Plot the stability data and visualise
the predictions with focus on one temperature.excursion() Predict a temperature excursion for a
product.step1_down_rmse() Calculate Root Mean Square Error
(RMSE) for the one-step Šesták–Berggren kinetic model.step1_plot_diagnostic() Generate residual diagnostic
plots from a step1_down fit.step1_sample_mvt() Take a selected number of samples
from the multivariate t distribution.Install AccelStab the following way -
install.packages("AccelStab")
library(AccelStab)Log an issue here or contact a moderator.