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 Sun Oct 30 18:00:27 2022 =-=-=-=-=
=-=-=-=-= Iteration 10 Sun Oct 30 18:00:27 2022 =-=-=-=-=
=-=-=-=-= Iteration 20 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 30 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 40 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 50 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 60 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 70 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 80 Sun Oct 30 18:00:28 2022 =-=-=-=-=
=-=-=-=-= Iteration 90 Sun Oct 30 18:00:29 2022 =-=-=-=-=
> 
> ## Plot all model state components.
> plot(model, "comp")
> 
> ## Plot the specific holiday state component.
> plot(ss[[2]], model)
> 
> ## Try again with some shrinkage.  With only 3 holidays there won't be much
> ## shrinkage.  
> ss2 <- AddLocalLevel(list(), y)
> 
> ## Plot the specific holiday state component.
> ss2 <- AddHierarchicalRegressionHoliday(ss2, y, holiday.list = holiday.list)
> model2 <- bsts(y, state.specification = ss2, niter = niter)
=-=-=-=-= Iteration 0 Sun Oct 30 18:00:30 2022 =-=-=-=-=
=-=-=-=-= Iteration 10 Sun Oct 30 18:00:30 2022 =-=-=-=-=
=-=-=-=-= Iteration 20 Sun Oct 30 18:00:30 2022 =-=-=-=-=
=-=-=-=-= Iteration 30 Sun Oct 30 18:00:30 2022 =-=-=-=-=
=-=-=-=-= Iteration 40 Sun Oct 30 18:00:30 2022 =-=-=-=-=
=-=-=-=-= Iteration 50 Sun Oct 30 18:00:30 2022 =-=-=-=-=
=-=-=-=-= Iteration 60 Sun Oct 30 18:00:31 2022 =-=-=-=-=
=-=-=-=-= Iteration 70 Sun Oct 30 18:00:31 2022 =-=-=-=-=
=-=-=-=-= Iteration 80 Sun Oct 30 18:00:31 2022 =-=-=-=-=
=-=-=-=-= Iteration 90 Sun Oct 30 18:00:31 2022 =-=-=-=-=
> 
> plot(model2, "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