CRAN Package Check Results for Package bsts

Last updated on 2026-05-21 07:49:17 CEST.

Additional issues

rchk

Check Details

Version: 0.9.11
Check: examples
Result: ERROR Running examples in ‘bsts-Ex.R’ failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: regression.holiday > ### Title: Regression Based Holiday Models > ### Aliases: AddRegressionHoliday RegressionHolidayStateModel > ### AddHierarchicalRegressionHoliday > ### HierarchicalRegressionHolidayStateModel > ### Keywords: models > > ### ** Examples > > trend <- cumsum(rnorm(730, 0, .1)) > dates <- seq.Date(from = as.Date("2014-01-01"), length = length(trend), by = "day") > y <- zoo(trend + rnorm(length(trend), 0, .2), dates) > > AddHolidayEffect <- function(y, dates, effect) { + ## Adds a holiday effect to simulated data. + ## Args: + ## y: A zoo time series, with Dates for indices. + ## dates: The dates of the holidays. + ## effect: A vector of holiday effects of odd length. The central effect is + ## the main holiday, with a symmetric influence window on either side. + ## Returns: + ## y, with the holiday effects added. + time <- dates - (length(effect) - 1) / 2 + for (i in 1:length(effect)) { + y[time] <- y[time] + effect[i] + time <- time + 1 + } + return(y) + } > > ## Define some holidays. > memorial.day <- NamedHoliday("MemorialDay") > memorial.day.effect <- c(.3, 3, .5) > memorial.day.dates <- as.Date(c("2014-05-26", "2015-05-25")) > y <- AddHolidayEffect(y, memorial.day.dates, memorial.day.effect) > > presidents.day <- NamedHoliday("PresidentsDay") > presidents.day.effect <- c(.5, 2, .25) > presidents.day.dates <- as.Date(c("2014-02-17", "2015-02-16")) > y <- AddHolidayEffect(y, presidents.day.dates, presidents.day.effect) > > labor.day <- NamedHoliday("LaborDay") > labor.day.effect <- c(1, 2, 1) > labor.day.dates <- as.Date(c("2014-09-01", "2015-09-07")) > y <- AddHolidayEffect(y, labor.day.dates, labor.day.effect) > > ## The holidays can be in any order. > holiday.list <- list(memorial.day, labor.day, presidents.day) > > ## In a real example you'd want more than 100 MCMC iterations. > niter <- 100 > > ## Fit the model > ss <- AddLocalLevel(list(), y) > ss <- AddRegressionHoliday(ss, y, holiday.list = holiday.list) > model <- bsts(y, state.specification = ss, niter = niter) =-=-=-=-= Iteration 0 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 10 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 20 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 30 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 40 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 50 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 60 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 70 Wed May 20 17:07:52 2026 =-=-=-=-= =-=-=-=-= Iteration 80 Wed May 20 17:07:53 2026 =-=-=-=-= =-=-=-=-= Iteration 90 Wed May 20 17:07:53 2026 =-=-=-=-= > > ## Plot all model state components. > plot(model, "comp") Error in apply(curves, 2, quantile, probs = qtl, na.rm = TRUE) : dim(X) must have a positive length Calls: plot ... .PlotStateContribution -> PlotDynamicDistribution -> t -> apply Execution halted Examples with CPU (user + system) or elapsed time > 5s user system elapsed add.monthly.annual.cycle 3.126 0.02 5.134 Flavor: r-devel-linux-x86_64-debian-gcc