Nothing
knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
library(srlTS) library(magrittr) # for pipe
data("LakeHuron") fit_LH <- srlTS(LakeHuron) fit_LH coef(fit_LH)
If you have a univariate time series with suspected trend, such as the EuStockMarkets data set,
data("EuStockMarkets") X <- as.numeric(time(EuStockMarkets)) X_sp <- splines::bs(X-min(X), df = 9) fit_stock <- srlTS(log(EuStockMarkets[,1]), n_lags_max = 400, X = X_sp, w_exo = "unpenalized") tail(coef(fit_stock), 11) # insert plot?
data("nottem") fit_nt <- srlTS(nottem, n_lags_max = 24) coef(fit_nt)
data("UKDriverDeaths") fit_ukdd <- srlTS(UKDriverDeaths, n_lags_max = 24) coef(fit_ukdd)
Adding holiday indicators...
data("sunspot.month") fit_ssm <- srlTS(sunspot.month) fit_ssm
Model summaries
summary(fit_ssm)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.