A new form of non-linear autoregressive time series is proposed to model so
lar radiation data, by specifying joint marginal distributions at low lags
to be multivariate Gaussian mixtures. The model is also a type of multiproc
ess dynamic linear model, but with the advantage that the likelihood has a
closed form.