need to sort changepoint times in forecasts

nicholasjclark · Nov 22, 2023 · 4b7d708 · 4b7d708
1 parent 8c33790
commit 4b7d708
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/R/trends.R b/R/trends.R
@@ -1408,7 +1408,7 @@ forecast_trend = function(trend_model, use_lv, trend_pars,
 
         # Combine with changepoints from the history
         deltas <- c(trend_pars$delta_trend[[x]], deltas_new)
-        changepoint_ts <- c(trend_pars$changepoints, t_change_new)
+        changepoint_ts <- sort(c(trend_pars$changepoints, t_change_new))
 
         # Generate a trend draw
         draw <- piecewise_linear(t = 1:max(time),
@@ -1441,7 +1441,7 @@ forecast_trend = function(trend_model, use_lv, trend_pars,
 
         # Combine with changepoints from the history
         deltas <- c(trend_pars$delta_trend[[x]], deltas_new)
-        changepoint_ts <- c(trend_pars$changepoints, t_change_new)
+        changepoint_ts <- sort(c(trend_pars$changepoints, t_change_new))
 
         # Get historical capacities
         oldcaps <- trend_pars$cap[[x]]

diff --git a/src/mvgam.dll b/src/mvgam.dll