---
title: 'MACROCAST: Vintage-Consistent TSFM'
url: https://www.emergentmind.com/topics/macrocast
type: topic
---

# MACROCAST: Vintage-Consistent TSFM

Searching arXiv for MACROCAST and closely related TSFM papers to ground the article in the cited literature.
First, locating the MACROCAST paper itself.
MACROCAST is a lightweight, vintage-consistent time series foundation model designed specifically for real-time macroeconomic forecasting. It is defined by an end-to-end training and evaluation protocol that eliminates two leakage mechanisms that are especially consequential in macroeconomic applications: temporal contamination, in which a model has seen the realized values of the same series it later forecasts, and revision bias, in which a model is trained on fully revised data rather than the preliminary releases available to a forecaster at the time [2606.28670]. The model is pretrained only on purely synthetic time series and fine-tuned only on synthetic panels generated from econometric models estimated on vintage-specific ALFRED data, so no observed future or revised value enters the model at any stage [2606.28670].

## 1. Real-time forecasting setting and leakage problem

Real-time macro forecasting is formulated around the information set actually available at the forecast date. In this setting, statistical agencies release preliminary estimates that are later revised, sometimes substantially. MACROCAST is motivated by the observation that a model trained on fully revised data has implicitly learned from information unavailable to a real-time forecaster [2606.28670].

The formal setup is vintage-specific. Let \(v\) denote the ALFRED vintage, interpreted as the month of release, and let the panel at vintage \(v\) be \(Y^{(v)}_{1:T_v} \in \mathbb{R}^{N\times T_v}\). Because revisions occur across releases, an entry \(Y^{(v)}_{i,t}\) can differ across vintages \(v \neq v'\). At forecast origin \(t\), training uses the expanding window through \(t-1\) from the vintage-\(v\) panel only, and the forecast for series \(i\) at horizon \(h\) is

\[
\hat{y}_{i,t+h\mid t,v} \;=\; f\bigl(\,\bm{y}^{(v)}_{1:t-1}\,\bigr),
\]

with no observation from vintages \(v'>v\) used in either training or evaluation [2606.28670].

Within this framework, MACROCAST distinguishes two leakage types. **Temporal contamination** refers to training on observed values of the same target series that will later be forecast out-of-sample. **Revision bias** refers to training on final, revised values when only preliminary releases were available in real time [2606.28670]. The model’s central methodological claim is that it eliminates both forms of leakage end-to-end.

Forecast errors are evaluated against “six-month delayed actuals,” defined as the FRED-MD vintage released approximately six months after the target date. This choice is intended to avoid conflating forecast errors with subsequent revisions while also not relying on especially noisy first estimates [2606.28670].

## 2. Architecture and forecasting mechanism

MACROCAST builds on TempoPFN’s linear recurrent neural network design rather than a transformer architecture [2606.28670]. In the description given for the model, this choice reduces quadratic attention costs, preserves a state-space interpretation, and enables fast per-vintage recalibration [2606.28670]. The recurrence is

\[
\bm{y}_t = g(\bm{H}_t,\bm{x}_t), \qquad
\bm{H}_t = \bm{A}(\bm{x}_t)\bm{H}_{t-1} + \bm{B}(\bm{x}_t)\bm{x}_t,
\]

with input-dependent state transition matrices \(A(\cdot)\) and \(B(\cdot)\), linear in the state and nonlinear in the inputs, together with a nonlinear output map \(g(\cdot)\) [2606.28670].

The instantiated recurrence is the GatedDeltaProduct update. The state update is stacked in three layers, and each layer is augmented with a short moving-average filter, normalization, and a nonlinear projection [2606.28670]. The hidden state dimension is \(128\), and the total parameter count is approximately \(1\) million [2606.28670].

For multi-step prediction, MACROCAST uses **state weaving**, which alternates observed-history and forecast stages. In this procedure, the model feeds its own 1-step predictions back as inputs for subsequent steps, yielding multi-horizon forecasts in the econometric analogue of iterating a VAR forward [2606.28670]. This mechanism is paired with linear-recurrence \(O(T)\) training complexity and single-pass multi-horizon inference [2606.28670].

The model is designed for unbalanced macro panels with ragged edges. Each series is standardized, and missing observations are replaced by a learnable missing-value embedding rather than being imputed [2606.28670]. The training and evaluation context spans up to about ten years of monthly data, while fine-tuning operates on slices of length \(L=240\) months [2606.28670].

MACROCAST also produces probabilistic forecasts. It predicts nine conditional quantiles per series and horizon, with \(\tau \in \{0.1,0.2,\dots,0.9\}\), and uses the median \((\tau=0.5)\) as the point forecast [2606.28670]. The training objective is the joint pinball loss over series, horizons, and quantiles:

\[
\mathcal{L}
= \sum_{i=1}^{N}\sum_{\tau \in \{0.1,\dots,0.9\}}\sum_{h=1}^{H}
\Bigl[\tau\,\max\bigl(y_{i,T+h} - \hat{y}^{(\tau)}_{i,T+h},0\bigr) + (1-\tau)\,\max\bigl(\hat{y}^{(\tau)}_{i,T+h} - y_{i,T+h},0\bigr)\Bigr].
\]

The quantiles are not post-processed to enforce monotonicity, and the reported comparisons use the median only, so quantile crossing does not affect RMSFE [2606.28670].

## 3. Leakage-free training pipeline

MACROCAST’s training pipeline has two stages, both structured to preserve vintage consistency and exclude exposure to future or revised observations [2606.28670].

The first stage is pretraining on purely synthetic time series. This stage uses approximately ten million distinct synthetic series, \(150{,}000\) optimization steps, and AdamW with linear warm-up and cosine decay; the reported cost is \(28\) hours on a single NVIDIA A6000, approximately one GPU-day [2606.28670]. The pretrained weights are then frozen as initialization for all fine-tuning runs [2606.28670].

The pretraining generators combine the Chronos synthetic library with TempoPFN augmentations. The covered process classes include trends, seasonality, periodic oscillations, smooth kernels, regime changes, outliers, and mean-reverting dynamics [2606.28670]. The generator family includes trend-times-seasonality with multiplicative Weibull noise, sine-wave families, Gaussian-process kernels including SE, periodic, rational quadratic, and white noise, sawtooth and step-function processes, anomaly and spike processes, CauKer structural causal processes over a random DAG, and Ornstein–Uhlenbeck dynamics with regime-switching parameters [2606.28670]. Post-generation transforms include reversals, sign flips, regime insertion, amplitude modulation, calendar seasonality, resampling artifacts, TSMix convex mixing, and missing-value masks [2606.28670]. At each optimization step, a series and horizon from \(1\) to \(60\) steps are drawn, and the pinball loss is evaluated on held-out future points [2606.28670].

The second stage is fine-tuning on synthetic panels calibrated to vintage-specific ALFRED data. For each vintage \(v\), the cleaned stationary panel \(y^{(v)}_{1:T_v}\) is the sole input to econometric model estimation [2606.28670]. MACROCAST then calibrates four generators and simulates \(400\)-month synthetic panels with a \(60\)-step burn-in [2606.28670]. The four generators are a Dynamic Factor Model, a Bayesian VAR, a univariate autoregressive generator, and a block bootstrap [2606.28670].

The Dynamic Factor Model is specified as

\[
\bm{y}_t^{(v)} = \bm{\Lambda}^{(v)} \bm{f}_t^{(v)} + \bm{\varepsilon}_t^{(v)},\quad
\bm{f}_t^{(v)} = \sum_{\ell=1}^{p} \bm{A}_\ell^{(v)} \bm{f}_{t-\ell}^{(v)} + \bm{u}_t^{(v)},
\]

with idiosyncratic noise \(\mathrm{diag}(\sigma^{(v)}_\varepsilon)\), factor innovations \(\Sigma_u^{(v)}\), EM estimation in state-space form via DynamicFactorMQ, factor count capped at \(k=3\) with fallback to \(2\) or \(1\), loadings estimated by OLS, and factor dynamics estimated with VAR(2) on smoothed factors [2606.28670].

The Bayesian VAR is

\[
\bm{y}_t^{(v)} = \bm{c}^{(v)} + \sum_{\ell=1}^{p} \bm{B}_\ell^{(v)} \bm{y}_{t-\ell}^{(v)} + \bm{\eta}_t^{(v)},\quad
\bm{\eta}_t^{(v)}\sim\mathcal{N}(\bm{0},\bm{\Sigma}^{(v)}),
\]

with \(p=2\) and Minnesota-type priors [2606.28670]. Own-lag coefficients shrink toward \(1\) for unit-root series or \(0\) for stationary series, while cross-lags shrink toward \(0\) with strength proportional to \(\lambda_1\lambda_2\hat{\sigma}_i/(\hat{\sigma}_j \ell^d)\). The baseline hyperparameters are \(\lambda_1=0.10\), \(\lambda_2=0.50\), and \(d=1\); for the full \(126\)-variable BVAR, tighter shrinkage \(\lambda_1=0.05\), \(\lambda_2=0.30\) is used [2606.28670]. To accommodate early vintages with limited sample length, both clustered BVARs and a full BVAR are estimated [2606.28670].

The univariate autoregressive generator estimates

\[
y_{i,t}^{(v)} = \phi^{(v)}_{i,1} y_{i,t-1}^{(v)} + \phi^{(v)}_{i,2} y_{i,t-2}^{(v)} + e_{i,t}^{(v)},\quad
e_{i,t}^{(v)}\sim \mathcal{N}(0,\sigma^{2\,(v)}_i),
\]

by OLS for each series and simulates independently as

\[
y^{(s)}_{i,t} = \sum_{\ell=1}^{p_i} \hat\phi^{(v)}_{i,\ell}\, y^{(s)}_{i,t-\ell} + \hat\sigma^{(v)}_i\, z_{i,t},\quad
z_{i,t}\sim\mathcal{N}(0,1).
\]

If data are too short or coefficients are explosive \((>|0.99|)\), the generator falls back to small-variance white noise \(\mathcal{N}(0,0.01^2)\) [2606.28670]. The more general ARIMA form is also noted,

\[
(1 - \textstyle\sum_{i=1}^p \phi_i L^i)\,(1-L)^d\,y_t
= (1 + \sum_{j=1}^q \theta_j L^j)\,\epsilon_t,
\]

with MACROCAST’s AR generator corresponding to \(d=0, q=0\) [2606.28670].

The block bootstrap samples contiguous blocks from \(y^{(v)}_{1:T_v}\), with random starts and random lengths, concatenating them to fill \(400\) months and clipping extreme outliers greater than \(25\times\mathrm{MAD}\) [2606.28670]. Two regimes are used: \(24\)–\(60\) months and \(36\)–\(96\) months [2606.28670].

Per vintage, the fine-tuning corpus contains \(2{,}100\) stored panels: \(500\) standard bootstrap, \(400\) long bootstrap, \(400\) AR, \(300\) DFM, \(300\) clustered BVAR, and \(200\) full BVAR [2606.28670]. Training examples are then constructed with stochastic augmentation by drawing a random \(240\)-month window, selecting a random subset of \(n_c \in \{32,\dots,126\}\) series, permuting columns, adding Gaussian noise with \(\kappa = 0.05\times \mathrm{MAD}\) per series, and forcing inclusion of outlier-prone series such as help-wanted with probability \(0.3\) [2606.28670].

Fine-tuning uses \(6{,}000\) AdamW steps per vintage, a short linear warm-up of approximately \(250\) steps, cosine decay, gradient clipping with \(L^2\) norm \(1.0\), and mixed precision in bfloat16 [2606.28670]. The gradient-based fine-tuning stage completes in approximately nine minutes on a modern GPU, and the full monthly refresh including econometric estimation and simulation runs in well under \(45\) minutes [2606.28670].

## 4. Evaluation protocol and benchmark design

The evaluation uses U.S. monthly FRED-MD data retrieved through ALFRED/FRED so that both monthly vintages and revision histories are available [2606.28670]. FRED-MD contains \(137\) series, of which \(14\) are excluded because they are discontinued, irregular, or have publication lags greater than two months, leaving \(123\) series [2606.28670].

The evaluation window runs from August 1999 through December 2024. At each vintage month \(t\), models are trained on data through \(t-1\) using an expanding window [2606.28670]. For the ragged edge, series lagging two months are carried forward so that they have a \(t-1\) value, while series with lags of at least three months are excluded [2606.28670]. Forecasts are made for horizons \(h=0\) through \(12\), where \(h=0\) is a nowcast for month \(t\), and targets from January through June 2020 are excluded to avoid the extreme pandemic shock [2606.28670].

Forecast comparison is based on RMSFE ratios relative to an AR(1) benchmark:

\[
\mathrm{RMSFE}_{i,j,h} =
\sqrt{\frac{\sum_{\tau} \bigl(y_{j,\tau+h}-\hat{y}_{i,j,\tau+h}\bigr)^2}
{\sum_{\tau} \bigl(y_{j,\tau+h}-\hat{y}_{\text{AR(1)},j,\tau+h}\bigr)^2}},
\]

with reporting at the median across series and also within category subsets [2606.28670]. Statistical significance is assessed with the Diebold–Mariano test using the Harvey–Leybourne–Newbold finite-sample correction, a one-sided alternative that the model beats AR(1), and Bartlett kernels with \(h-1\) lags for the long-run variance [2606.28670]. Tables use majority-vote stars when more than half of the series in a category satisfy \(p<0.10\), \(0.05\), or \(0.01\) [2606.28670].

The benchmark set comprises AR(1), a BVAR with conjugate Normal–inverse-Wishart priors and hyperparameters chosen by marginal likelihood, a DFM implemented through factor-augmented direct projections with factors obtained by principal components and EM and the number of factors selected by Bai–Ng information criteria up to eight, Chronos-2 and Moirai2-Small evaluated zero-shot with public pretrained weights, and MACROCAST itself, which is pretrained once on synthetic data and fine-tuned monthly per vintage on synthetic panels estimated from that vintage’s ALFRED data [2606.28670].

A notable design feature is that MACROCAST is never tuned to any single target series; it is evaluated in a series-level zero-shot mode [2606.28670]. This suggests that the model is intended to serve as a panel-level forecaster rather than a bank of separately specialized univariate systems.

## 5. Empirical results

Across the full sample, MACROCAST is reported to improve on AR(1), match or beat Chronos-2, and outperform the econometric BVAR and DFM benchmarks, while preserving the stated leakage-free guarantees [2606.28670].

For the \(123\)-series full-sample evaluation against six-month delayed actuals, MACROCAST’s median RMSFE ratios relative to AR(1) are \(0.959^\*\) at \(h=0\), \(0.971^\*\) at \(h=1\), \(0.987\) at \(h=3\), \(0.992\) at \(h=6\), \(0.995\) at \(h=9\), and \(0.995\) at \(h=12\) [2606.28670]. Chronos-2 records \(0.952\), \(0.977\), \(0.993\), \(0.992\), \(0.998\), and \(0.996\) at the same horizons, while BVAR and DFM are centered around one and edge above one from medium horizons onward; Moirai2-Small tracks AR(1) [2606.28670].

Breadth is emphasized as a separate performance criterion. MACROCAST beats AR(1) for the majority of series at all horizons \(h=1,\dots,12\), and at \(h=1\) approximately \(86\%\) of series fall below AR(1), with a dense mass just below one and a highly compressed upper tail, compared with heavier right tails for BVAR and DFM at longer horizons [2606.28670]. The broader summary is that MACROCAST improves on AR(1) for roughly \(80\%\) of series-horizon pairs [2606.28670].

On curated subsets, the model shows mixed but generally strong behavior. On the “Medium” set of \(18\) headline indicators, MACROCAST achieves \(0.974^\*\) at \(h=1\), while Chronos-2 is best at \(h=3\) and \(h=6\) with \(0.965\) and \(0.967\) [2606.28670]. On the “Large” set of \(30\) indicators, MACROCAST is best at four of six horizons, including \(h=1: 0.970^\*\), \(h=3: 0.975\), and \(h=12: 0.996\) [2606.28670].

The category breakdown indicates heterogeneity across macro domains. At \(h=1\), MACROCAST is best in Output \((0.984^\*)\), Labor \((0.958^\*)\), Money \((0.970^\*)\), Interest/FX \((0.983^\*)\), and is closest to one in Stocks; DFM is best for Consumption \((0.979)\), and Chronos-2 leads Housing \((0.884)\) and Prices \((0.953)\) [2606.28670]. At \(h=6\), MACROCAST leads Output, Labor, Money, Prices, and Stocks, while Chronos-2 leads Housing and Interest/FX [2606.28670]. At \(h=12\), MACROCAST leads Output, Labor, Consumption, Interest/FX, and Stocks, while Chronos-2 remains strongest in Housing and marginally in Money [2606.28670].

Several individual series are highlighted. For nonfarm payrolls (PAYEMS), MACROCAST is best across models at every horizon from \(h=1\) onward, with \(0.747^{***}\) at \(h=1\), \(0.851^{**}\) at \(h=3\), \(0.849^{**}\) at \(h=6\), \(0.859^{**}\) at \(h=9\), and \(0.900^\*\) at \(h=12\) [2606.28670]. For industrial production (INDPRO), MACROCAST stays at or below AR(1) throughout, with examples including \(0.974\) at \(h=3\) and \(0.992\) at \(h=6\), whereas BVAR and DFM deteriorate beyond \(h=3\) [2606.28670]. For the unemployment rate (UNRATE), it is described as the most stable across horizons, with \(0.980^{**}\) at \(h=3\), \(0.992\) at \(h=6\), and approximately one at long horizons [2606.28670]. For the federal funds rate (FEDFUNDS), MACROCAST is strong at short and medium horizons, including \(0.863^{***}\) at \(h=1\), \(0.893^{***}\) at \(h=3\), and \(0.952^{**}\) at \(h=6\), though Chronos-2 is marginally ahead at \(h=6\) to \(9\) [2606.28670]. For capacity utilization (CUMFNS), MACROCAST is best from \(h=3\) to \(h=12\), including \(0.975\) at \(h=3\) and \(0.996\) at \(h=12\) [2606.28670].

The qualitative robustness narrative is that MACROCAST’s gains are concentrated around turning points and stress episodes, especially the 2001 recession and the 2008–2009 crisis, whereas AR(1) is difficult to beat during calm expansions [2606.28670]. The long-horizon distributions remain anchored below one for MACROCAST, while BVAR shifts rightward, with an example median rising from \(0.99\) at \(h=1\) to \(1.03\) at \(h=12\) [2606.28670].

## 6. Interpretation, practical use, and limitations

MACROCAST’s practical significance lies in combining leakage-free real-time training with a computational budget compatible with monthly recalibration. The pretrained model requires approximately one GPU-day once, and each vintage-specific fine-tune takes approximately nine minutes, with an end-to-end refresh time below \(45\) minutes including econometric estimation and synthetic simulation [2606.28670]. This makes genuine month-by-month deployment operationally feasible.

In applied use, the inputs are monthly ALFRED/FRED-MD vintages transformed to stationarity using FRED-MD transformation codes [2606.28670]. The ragged edge is handled by carrying forward one month for series with two-month publication lags and excluding longer-lag series [2606.28670]. Because missing values are handled through a learned embedding, the model supports unbalanced panels without explicit imputation [2606.28670]. Optimization uses AdamW with warm-up and cosine decay, gradient clipping, bfloat16 mixed precision, slices of length \(L=240\), and \(n_c \in [32,126]\) series per batch, with small noise augmentation at \(\kappa=0.05\times\mathrm{MAD}\) [2606.28670].

The model’s central methodological contribution is not merely architectural. It is the coupling of a lightweight linear-RNN TSFM to a vintage-consistent synthetic-data regime in which pretraining uses only simulated series and fine-tuning uses only synthetic panels generated from econometric models estimated on the contemporaneous vintage [2606.28670]. This suggests a template for domains in which revision histories and temporal leakage are structurally important, although that broader generalization is an inference rather than a directly demonstrated result.

Several limitations are explicit. First, the macro calibration depends on synthetic panels drawn from vintage-estimated econometric models; this removes leakage and provides diverse training trajectories, but it may miss fine-grained features not captured by the parametric suite [2606.28670]. Second, the evaluation focuses on U.S. monthly FRED-MD data, and extension to mixed-frequency nowcasting with weekly or daily indicators is presented as a natural direction [2606.28670]. Third, despite regime-switching generators and bootstrap preservation of real-time co-movements, large structural breaks remain difficult; January through June 2020 are excluded for this reason [2606.28670]. Fourth, quantile outputs are unconstrained and may cross, so future work could impose monotone quantile layers or density calibration for improved interval forecasts [2606.28670]. Additional proposed extensions include structural identification, scenario analysis, expansion to other economies or domains, and hybrid fine-tuning that blends synthetic vintages with carefully curated leakage-safe observed signals such as survey nowcasts [2606.28670].

## 7. Position within macroeconomic forecasting research

MACROCAST is situated at the intersection of foundation-model methodology and classical real-time macroeconometrics. It retains familiar macroeconometric components—Dynamic Factor Models, Bayesian VARs, autoregressive persistence, and bootstrap resampling—but uses them to generate synthetic training corpora rather than as the deployed forecasters themselves [2606.28670]. In that sense, the model can be viewed as a TSFM whose domain adaptation is mediated by vintage-specific econometric simulation.

Its comparative claims are correspondingly specific. Relative to Chronos-2, described as the strongest currently available TSFM in the reported comparisons, MACROCAST matches or surpasses it in aggregate and category-level evaluations while avoiding both temporal contamination and revision bias [2606.28670]. Relative to BVAR and DFM benchmarks, it outperforms them under the same leakage-free real-time protocol [2606.28670]. Relative to naive autoregressive forecasting, it improves on AR(1) for roughly \(80\%\) of series-horizon pairs and for approximately \(86\%\) of one-month-ahead forecasts [2606.28670].

The broader significance of these results lies in the formulation of “vintage consistency” as a first-class design constraint for foundation models in macroeconomics [2606.28670]. Existing TSFMs may score well under conventional retrospective evaluation but still violate the information set available to a forecaster. MACROCAST’s contribution is to define a training, fine-tuning, and evaluation regime in which no future or revised observation is ever ingested. A plausible implication is that its empirical results should be interpreted not only as forecast-accuracy outcomes, but also as evidence that leakage-free TSFMs need not forfeit competitiveness with stronger but less strictly real-time alternatives.

Source: https://www.emergentmind.com/topics/macrocast