Decomposable additive forecast splits a series into a piecewise-linear trend with automatically placed rate changepoints plus a Fourier-series seasonal pattern, fits them jointly, and projects h steps ahead with a prediction band.
The series is modelled as y(t) = trend(t) + seasonality(t) + noise. The trend is piecewise-linear: a base growth rate plus a set of candidate 'changepoints' (evenly spaced over the first 80% of history) where the rate is allowed to shift. A sparsity (Laplace) prior on those rate shifts keeps most of them at exactly zero, so only the changepoints the data really need stay active — this is the automatic changepoint selection. Seasonality is a sum of sine/cosine harmonics at the chosen period, entered as ordinary regressors. Everything is fit at once by penalised least squares (a MAP estimate).
The forecast interval reflects trend uncertainty: future rate changes are assumed to arrive at the same frequency and typical size as the changepoints seen in history, so the band widens with the horizon. With no changepoints and no seasonality the model reduces to an ordinary least-squares linear trend; with a flat trend it reduces to harmonic (Fourier) regression.
One numeric series (row order = time) + forecast horizon h + a seasonal period (or automatic detection) + the number of candidate trend changepoints + interval level.
A trend + seasonality component decomposition chart, a history-plus-forecast fan chart with the prediction band, a table of the active trend changepoints and their rate changes, and a table of point forecasts with lower/upper bounds.
The trend component shows the underlying level with season and noise removed; a changepoint with a large rate change marks where growth sped up or slowed. The seasonal component shows the repeating within-cycle pattern and integrates to about zero over a full period.
The point forecast drifts along the extrapolated trend with the seasonal pattern continued; the widening band shows how quickly trend uncertainty accumulates. Fit is by penalised least squares (no closed-form oracle), validated by synthetic recovery, exact component coherence, and its OLS reductions.