---
title: Dynamic Bayesian Predictive Synthesis
url: https://www.emergentmind.com/topics/dynamic-bayesian-predictive-synthesis
type: topic
---

# Dynamic Bayesian Predictive Synthesis

Dynamic Bayesian Predictive Synthesis (BPS) is a formally Bayesian framework for combining multiple predictive distributions from models, forecasters, or other “agents” into a single predictive distribution whose calibration parameters evolve over time. In its core representation, the reported agent densities are treated as data, latent agent states are drawn from those densities, and a user-specified synthesis function maps the latent states to the outcome distribution. The static synthesis identity was given a rigorous proof by Masuda and Irie, while the dynamic time-series framework was developed by McAlinn and West through state-space models that allow time-varying biases, miscalibration, inter-dependencies, and forecast-combination weights [2409.09660] [1601.07463].

## 1. Foundational representation

The fundamental BPS equation starts with a decision maker (DM) facing $K$ agent forecast densities $h_k(x_k)$ and an information set $\mathcal H=\{h_1,\dots,h_K\}$. The DM’s posterior predictive density for outcome $y$ is

$$
p(y\mid \mathcal H)=\int \alpha(y\mid x_1,\dots,x_K)\prod_{k=1}^K h_k(x_k)\,dx_1\cdots dx_K,
$$

where $\alpha(y\mid x_{1:K})$ is a proper density in $y$ for each $x_{1:K}$ [2409.09660]. In the agent-opinion analysis perspective, the densities $h_j$ are themselves data, the latent variables $x_{1:J}$ are “agent states,” and the synthesis function $\alpha(y\mid x_{1:J})$ encodes how the DM believes these latent forecasts relate to the realized outcome [1601.07463].

Masuda and Irie show that the static formula follows from a partially specified prior in which only marginal means of agent-reported probabilities are fixed, together with a consistency condition requiring that marginalization reproduces the DM’s unconditioned prior and that extreme forecasts match Bayes-rule values. They give two proofs: a proof by marginalization using induction on $K$, and a proof via an artificial independent prior that shares the same marginal means [2409.09660]. The first proof moves from the linear-in-deviations form for discretized forecast probabilities to a discrete mixture over joint-extreme conditional probabilities, and then takes a fine-grid limit. The second proof uses iterated expectations under a convenient independent prior to recover the same discrete mixture, and again takes the grid limit to the integral form [2409.09660].

A central structural point is the equivalence between two representations: a linear-deviation form and a discrete-mixture form. In the continuous limit, the mixture form shows that $\alpha(y\mid x)$ is the DM’s implied conditional density of $y$ when each agent is degenerate at $x_k$ [2409.09660]. The graphical interpretation is a two-stage generative model in which $x_k\sim h_k(\cdot)$ independently and then $y\sim \alpha(\cdot\mid x_{1:K})$. Masuda and Irie emphasize that the product $\prod h_k$ in this representation does **not** require agent-forecast independence in the DM’s actual prior; independence enters only as a proof device in one argument and as a conditional feature of the generative representation [2409.09660].

## 2. Dynamic state-space constructions

Dynamic BPS embeds the static kernel in a time-indexed state-space model. At time $t$, agents deliver predictive densities $h_{k,t}(x_{k,t}\mid \mathcal F_{t-1})$, the synthesis function may depend on a latent state $\theta_t$, and $\theta_t$ evolves by a Markovian law such as a random walk or a dynamic linear model [2409.09660]. In general form, one writes

$$
p(y_{t+h}\mid \mathcal D_t)=\int \alpha(y_{t+h}\mid x_{1:J,t+h\mid t},\Psi_{t+h})
\prod_{j=1}^J \pi_{j,t}(x_{j,t+h\mid t})\,dx_{1:J,t+h\mid t},
$$

with latent states and parameters $\Psi_{t+h}$ governing the synthesis function [2311.12671].

The canonical specification is linear-Gaussian. McAlinn and West formulate the synthesis density as

$$
\alpha_t(y_t\mid \theta_t,x_t)=N(y_t\mid F_t'\theta_t,V_t),
$$

with

$$
F_t=(1,x_{t1},x_{t2},\dots,x_{tJ})', \qquad
\theta_t=(\theta_{t0},\theta_{t1},\dots,\theta_{tJ})',
$$

and state evolution

$$
\theta_t=G_t\theta_{t-1}+\omega_t,\qquad \omega_t\sim N(0,W_t),
$$

often with $G_t=I$ so that the coefficients follow a random walk [1601.07463]. In this form, the intercept absorbs residual common bias, the coefficients on the latent agent states encode dynamic calibration, and $W_t$ controls the rate at which synthesis weights adapt [1601.07463]. A closely related DLM formulation writes $y_t\mid z_t,\theta_t,V_t\sim N(F_t'\theta_t,V_t)$, where $z_t=(z_{t1},\dots,z_{tJ})'$ are auxiliary draws from the component densities $f_{j,t}$ and the synthesized one-step-ahead density is obtained by integrating over both $z_t$ and $\theta_t$ [2510.07180].

Several named special cases appear in the recent literature. In the tree-based synthesis paper, the common parametric choice is a time-varying-parameter linear regression
$$
\alpha(y_{t+h}\mid x,\Psi_{t+h})=
N\!\left(y_{t+h}\mid c_{t+h}+\sum_{j=1}^J(\gamma_j+\beta_{j,t+h})x_{j,t+h\mid t},\sigma_{t+h}^2\right),
$$
where $c_{t+h}$ is a time-varying intercept, $\gamma_j$ are static combination weights, and $\beta_{j,t+h}$ are dynamic weight adjustments. The variants BPS-CONST and BPS-RW respectively set $\beta_{j,t+h}=0$ for all $t$ or evolve $\beta$ as a Gaussian random walk [2311.12671].

Dynamic BPS also has a multivariate state-space generalization. In the multivariate macroeconomic formulation, each agent supplies a $p$-variate density forecast $h_{k,t}(x_{t,k})$ for a latent vector, the synthesis density is $N_p(y_t\mid B_t x_t,V_t)$, and a structured dynamic latent-factor model yields a DFSUR representation $y_t=F_t\theta_t+\nu_t$, $\nu_t\sim N_q(0,\Sigma_t)$ with random-walk evolution for $\theta_t$ [1711.01667]. This retains the basic synthesis logic while allowing dynamic inter-series dependence through $\Sigma_t$ and through the joint posterior of the synthesis coefficients [1711.01667].

## 3. Sequential learning and posterior computation

Posterior computation in dynamic BPS is driven by the latent agent states and the time-varying synthesis parameters. In the linear-Gaussian case, the standard posterior target is
$$
p(\theta_{1:T},V_{1:T},x_{1:T}\mid y_{1:T},H_{1:T}),
$$
and McAlinn and West use a two-block Gibbs scheme: conditional on $x_{1:T}$, the model is a conditionally Gaussian DLM and $\{\theta_t,V_t\}_{t=1}^T$ are sampled by forward filtering–backward sampling (FFBS); conditional on $\theta_t,V_t,y_t$, each $x_t$ is sampled from
$$
p(x_t\mid \theta_t,V_t,y_t,H_t)\propto
N(y_t\mid F_t'\theta_t,V_t)\prod_{j=1}^J h_{tj}(x_{tj}),
$$
with closed-form Gaussian updates when the agent densities are normal and scale-mixture augmentation when they are Student-$t$ [1601.07463].

A conjugate filtering form is available for the DLM synthesis. With a normal–inverse–gamma prior on $(\theta_0,V_0)$, forward filtering updates the prior moments $(a_t,R_t)$, constructs the one-step predictive distribution, and then applies Kalman-style posterior recursions for $(m_t,C_t,n_t,d_t)$ [2510.07180]. This yields a fully Bayesian sequential updating scheme in which the synthesis coefficients are literally time-varying intercept and weights on the component forecasts [2510.07180].

Outcome-dependent dynamic BPS requires richer latent structure. In the mixture-based formulation with dynamic biases, covariance parameters, and Dirichlet base weights, posterior inference targets $p(\theta_t,x_t\mid y_{1:t})$ using a Gibbs sampler augmented by a discrete indicator selecting which mixture component “fired.” To maintain conjugate parametric priors for the next time point, the Monte Carlo samples are re-projected to a Normal-Inverse-Wishart distribution for $(\beta_t,\Sigma_t)$ and a Dirichlet distribution for $q_t$ by minimizing Kullback–Leibler divergence, a step described as “variational Bayes” in the paper [1803.01984].

Sequential computation has also been specialized to avoid repeated full MCMC. For linear and Gaussian synthesis, a Rao-Blackwellized particle filter marginalizes out $(\theta_t,\nu_t)$ and tracks only a particle approximation to $p(x_{1:t}\mid y_{1:t},H_{1:t})$, using the closed-form Student-$t$ one-step predictive density as the particle weight. When the effective sample size falls below a threshold, the algorithm discards the particle set and invokes an MCMC approximation to the posterior, then restarts the filter with equal-weight particles [2308.15910]. In the reported U.S. quarterly inflation example, this RBCPF plus timely MCMC intervention maintained synthesized-posterior accuracy at a small fraction of the cost of MCMC-only computation, reported as approximately $0.6$ seconds versus $360$ seconds per time point [2308.15910].

The same paper proposes adaptive discount-factor estimation through power-discounted likelihoods. A grid of discount-factor pairs is evaluated by discounted log predictive likelihoods, and softmax or argmax weights over the grid define a second-layer mixture predictive density. In the 2020–2022 inflation burst, fixed-discount DBPS adapted slowly, whereas discount-factor estimation through the Loss Discounting Framework produced wider intervals peaking at the burst and improved log-predictive-density ratios [2308.15910].

## 4. Synthesis-function families and major variants

Dynamic BPS is defined by the synthesis function $\alpha_t$, and much of the methodology consists of specifying that function so as to capture calibration structure that static linear pools cannot represent.

| Variant | Synthesis function | Distinctive feature |
|---|---|---|
| Outcome-dependent pool | Mixture with $w_{k,t}(y_t)$ and optional baseline $h_{0t}(y_t)$ | Weights depend on forecast outcome |
| Linear-Gaussian dynamic BPS | $N(y_t\mid F_t'\theta_t,V_t)$ | Dynamic intercepts and weights in a DLM |
| Tree-based BPS-RT | Regression-tree prior means for $\gamma_j,\beta_{j,t+h}$ | Nonparametric dependence on weight modifiers |
| DRQS / FDRQS | Asymmetric Laplace synthesis | Direct quantile targeting |
| Network BPS | Affine-logit pool $\sigma(\beta_t' z_{t,e})$ | Dynamic mechanism weights on edges |

In outcome-dependent BPS, the synthesis density can be written in discrete-mixture form as
$$
\alpha_t(y_t\mid x_t)=\omega_{0t}(x_t)h_{0t}(y_t)+\sum_{k=1}^K \omega_{kt}(x_t)\delta_{x_{kt}-\beta_{kt}}(y_t),
$$
which implies
$$
p(y_t\mid \mathcal F_{t-1})=w_{0t}(y_t)h_{0t}(y_t)+\sum_{k=1}^K w_{kt}(y_t)p_k(y_t\mid \mathcal F_{t-1}).
$$
The baseline density $h_{0t}$ is included to address model-set incompleteness, and the functions $\alpha_{kt}(x_t)$ defining $\omega_{kt}(x_t)=q_{kt}\alpha_{kt}(x_t)$ can take Gaussian “well,” softmax, “consensus,” or “herding” forms [1803.01984]. Because $w_{kt}(y_t)$ depends on the realized value of $y_t$, this formulation permits regional calibration in the tails or in other outcome regimes [1803.01984].

Tree-based BPS-RT replaces the static or random-walk prior on synthesis weights by nonparametric functions of user-chosen “weight modifiers” such as ex-post forecast scores, moments of the agent densities, cross-sectional dispersion, financial conditions, uncertainty, or a time trend. The prior means $\mu^\gamma_j(\cdot)$ and $\mu^\beta_j(\cdot)$ are learned by regression trees or sums of trees, while local scaling parameters follow horseshoe priors [2311.12671]. Although the tree structures are fixed through time, the modifiers $z^\beta_{j,t+h\mid t}$ vary with $t$, so the implied combination weights vary with time; when an indicator crosses a tree threshold, the prior mean for a dynamic coefficient can jump from one leaf value to another [2311.12671].

Dynamic Bayesian regression quantile synthesis (DRQS) changes the synthesis target from conditional means to conditional quantiles by using the asymmetric Laplace distribution in the observation equation,
$$
y_t=F_t'\theta_t+\epsilon_t,\qquad \epsilon_t\sim AL(\tau,\sigma_t),
$$
so that $F_t'\theta_t$ is the $\tau$th conditional quantile of $y_t$ [2603.11474]. Its multivariate extension, factor DRQS (FDRQS), imposes a time-varying latent factor structure on the synthesis weights to leverage cross-sectional dependencies across multiple time series [2603.11474].

A further extension adapts BPS to dynamic networks. In that setting, each structural mechanism—such as a stochastic block model, generalized random dot-product graph, Chung–Lu model, or Adamic–Adar score—is treated as an agent returning edge-probability forecasts. The synthesis layer stacks the agent logits into features $z_{t,e}$ and models edge probabilities through
$$
P_{t,e}=\sigma(\beta_t' z_{t,e}),\qquad \beta_t=\beta_{t-1}+\omega_t,
$$
with $\omega_t\sim N(0,W_t)$ [2606.26136]. This is a dynamic affine-logit pool rather than a Gaussian regression, but it preserves the core BPS architecture of latent agent forecasts plus dynamic synthesis weights [2606.26136].

## 5. Empirical domains and observed behavior

The original dynamic BPS time-series study examined quarterly U.S. inflation, short-term interest rates, and unemployment from 1961Q1 to 2014Q4 using four agent DLMs. In the 1990Q1–2014Q4 evaluation period, one-step BPS achieved $\mathrm{MSFE}_1=0.0512$, compared with $0.0598$ for the next best agent, and the horizon-customized four-step model BPS(4) achieved $\mathrm{MSFE}_4=0.3686$ compared with $0.4001$ for direct BPS, with LPDR gains of order $+70$ to $+85$ nats over single agents at horizon $4$ [1601.07463]. The same study reports that BPS and especially BPS(4) dominate Bayesian Model Averaging, linear and log pools, and the DeCo approach, with particularly marked differences during crisis periods [1601.07463].

In the outcome-dependent pooling paper, dynamic BPS is illustrated on Euro–USD daily log-prices over July–December 2016, forecasting five trading days ahead with time-varying autoregressions and DLMs as agents. Using “consensus” weights, posterior trajectories reveal evolving model biases, time-varying cross-model correlations that hit lows around the U.S.-election volatility and recover thereafter, and adaptive base weights that translate into effective outcome-dependent weights in the MCMC draws. In that setting, BPS outperforms each individual model, a static equal-weight pool, classical BMA scoring one day ahead, and an equal-weight pool augmented with baseline density [1803.01984].

Chernis, McAlinn, and collaborators evaluate tree-based dynamic BPS in two macroeconomic applications. For euro-area GDP growth using the ECB Survey of Professional Forecasters at the two-quarter horizon, a single-tree BPS-RT with stochastic volatility and AVG-SCORES weight modifiers achieves approximately $10$–$15\%$ lower CRPS than BPS-RW, with statistical significance by Diebold–Mariano, similar RMSE to benchmarks, well-calibrated PIT histograms, and sparse tree structures with one to two splits capturing crisis versus tranquil periods. For U.S. inflation from FRED-QD indicators, single-tree BPS-RT cuts CRPS by approximately $5$–$12\%$ relative to BPS-RW at horizon $1$, with larger improvements at horizon $4$, and adapts quickly to the 2021–2022 inflation spike [2311.12671].

Quantile-targeted dynamic BPS has been tested on U.S. inflation-at-risk and global GDP growth-at-risk. In the inflation application, DRQS achieves lower relative cumulative scores especially in the right tail and at the four-quarter horizon, broadening intervals appropriately and avoiding extreme mis-forecasting seen in some agents. In the global GDP application, FDRQS attains the lowest total relative cumulative score at both horizons, outperforms all agents and univariate DRQS, and during the COVID-19 shock in 2020 its RTCS “barely rises” while DRQS and the agent models exhibit large score deteriorations [2603.11474].

Recent work also carries dynamic BPS into domains beyond univariate macroeconomic density combination. In network forecasting, the synthesis layer returns calibrated edge forecasts and inference on mechanism weights, and the theory proves that a single snapshot identifies and estimates the weights at rate $(n^2\rho_n)^{-1/2}$, with a sharp distinguishability threshold $\Delta_n\asymp (n^2\rho_n)^{-1/2}$ and a tracking bound that pays only at regime switches [2606.26136]. In portfolio analysis, predictive-synthesis ideas have been extended in two related directions: Bayesian Predictive Decision Synthesis tilts model-specific predictive densities by realized utility scores in sequential portfolio rebalancing [2405.01598], and a later study uses BPS with dynamic linear models to combine multiple asset-return prediction models and then construct mean-variance and quantile-based portfolios from the synthesized posterior [2510.07180].

## 6. Conceptual relations, misconceptions, and active directions

Dynamic BPS is broader than classical Bayesian Model Averaging. The outcome-dependent pooling paper states explicitly that classical BMA is a special case in which the outcome-dependent synthesis weights reduce to outcome-independent BMA weights and no baseline density is used [1803.01984]. By contrast, BPS can incorporate a baseline or “safe-haven” density to address model-set incompleteness, can calibrate models differently across outcome regions, and can represent dependence structures among forecasts rather than treating models as competing, mutually exclusive explanations [1803.01984].

A common misconception is that the factorized term $\prod_k h_k(x_k)$ implies that the DM assumes the agents are independent. Masuda and Irie clarify that no assumption of agent-forecast independence is required in the DM’s actual prior; independence appears only in one proof as a convenient artificial prior sharing the same marginal means, and in the graphical interpretation conditional on the reported forecasts [2409.09660]. Another misconception is that dynamic BPS is synonymous with random-walk Gaussian regression. The literature now includes outcome-dependent pools, regression-tree synthesis, multivariate latent-factor forms, asymmetric-Laplace quantile synthesis, affine-logit network synthesis, and sequential particle-filter implementations [1803.01984] [2311.12671] [2603.11474] [2606.26136] [2308.15910].

The current methodological frontier is heterogeneous. In BPS-RT, the reported best empirical performance came from $S=1$ rather than large BART-style ensembles, which the authors interpret as suggesting that combination weights are largely stable with occasional regime shifts; they propose learning $S$ automatically or using spike-and-slab priors to control complexity [2311.12671]. The same paper notes that weight modifiers must be chosen by the user, that multivariate BPS-RT is a natural extension, and that sequential Monte Carlo may be preferable in real-time or very high-frequency settings [2311.12671]. DRQS and FDRQS push BPS toward direct quantile synthesis and cross-series factor structures [2603.11474], while the dynamic-network paper adds explicit identification theory and per-switch localization results [2606.26136]. Taken together, these developments indicate that dynamic BPS functions less as a single model class than as a coherent synthesis principle: latent agent forecasts are combined through a Bayesian calibration layer whose structure, dynamics, and inferential machinery can be adapted to the forecasting target at hand.

Source: https://www.emergentmind.com/topics/dynamic-bayesian-predictive-synthesis