The inputs are M individual forecasts y^1,…,y^M for a common target variable y at a common horizon h. Each forecast comes from a distinct model (VAR, DSGE, ARIMA, factor model, judgmental survey, etc.) or from the same model class with different specifications. The forecasts may be point forecasts (scalars) or density forecasts (full predictive distributions). Historical forecast errors are needed for weight estimation: a holdout sample of T0 periods where both the forecasts and the realized values are observed.
The combined forecast is y^c=∑m=1Mwmy^m, where wm≥0 and ∑mwm=1. The weights wm are chosen to minimize the forecast-error variance of the combination. Under the assumption that forecast errors have constant covariance Σe (M×M matrix), the optimal weights are w∗=Σe−1ι/(ι′Σe−1ι), where ι is a vector of ones. This is the minimum-variance portfolio from finance. In practice, Σe must be estimated from the holdout sample, introducing estimation error that can make the estimated optimal weights worse than equal weights.
Weight estimation methods range from simple to complex. Equal weights (wm=1/M) ignore the covariance structure entirely but avoid estimation error. Inverse-MSE weights set wm∝1/MSEm, using only individual model performance without cross-model covariance. Regression-based weights run the regression yt=β0+∑mβmy^m,t+εt and use β^m as weights (Granger-Ramanathan, 1984), but this requires enough holdout data to estimate M+1 parameters reliably. Bayesian model averaging sets wm proportional to the posterior model probability p(Mm∣data), computed from the marginal likelihood of each model. Time-varying weights use exponential decay: wm,t∝exp(−α∑s=1t(ys−y^m,s)2) with a forgetting factor α.