Structural Panel Bayesian VAR with Multivariate Time-Varying Volatility to Jointly Deal with Structural Changes, Policy Regime Shifts, and Endogeneity Issues

: This paper improves a standard Structural Panel Bayesian Vector Autoregression model in order to jointly deal with issues of endogeneity, because of omitted factors and unobserved heterogeneity, and volatility, because of policy regime shifts and structural changes. Bayesian methods are used to select the best model solution for examining if international spillovers come from multivariate volatility, time variation, or contemporaneous relationship. An empirical application among Central-Eastern and Western Europe economies is conducted to describe the performance of the methodology, with particular emphasis on the Great Recession and post-crisis periods. A simulated example is also addressed to highlight the performance of the estimating procedure. Findings from evidence-based forecasting are also addressed to evaluate the impact of an ongoing pandemic crisis on the global economy.

Introduction
Vector Autoregressions (VARs) are widely used when studying macroeconomicfinancial linkages to detect interdependencies and co-movements among multiple economic time-series. In the simplest form, error terms in the VAR models are assumed to have constant variances. While convenient, assuming time-invariant coefficients and variances, it turns out to be highly restrictive in capturing the evolution and thus the dynamics of multiple economic time-series. When time-varying series are introduced in a VAR to highlight the evolving relationship between multiple economic-financial variables, its state space structure need to be modeled and then used in the empirical analysis to estimate unobserved time variations and volatility. Since allowing time-varying coefficients and/or volatility introduces too many parameters than data points, the literature have proposed random processes to time variations or volatility to deal with the curse of dimensionality (see, for instance, Koop and Korobilis (2016)). The randomness in these time-varying parameters fits sufficiently well with Bayesian methods because the probability that a random sequence of events occurs can be interpreted through a reasonable expectationrepresenting a state of knowledge-or a quantification of a personal belief. More precisely, in Bayesian statistics, one assumes to know something about the possible values of the parameters of interest prior to seeing the data. Thus, the inputs required for estimating that quantity are prior knowledge (or belief) plus data; contrary to frequentist statistics, where these inputs are represented by the data only.
In this context, multicountry Bayesian VAR (BVAR) models have given a new impulse to the literature to evaluate macroeconomic-financial linkages, to test specification hypotheses (see, for instance, Korobilis (2016)), and to conduct policy exercises (see, e.g., Ciccarelli et al. (2018); Canova and Ciccarelli (2009) ;Canova et al. (2007Canova et al. ( , 2012; and Koop (1996)). Nevertheless, although estimation of time-varying structures is feasible with a large homogeneous cross-section, heterogeneous dynamics due to an unexpected shock combined with not directly observed or measured factors make it difficult to exploit cross-sectional information estimating time-series variations in multicountry setups (see, for instance, Koop and Korobilis (2016)). More precisely, these empirical models tend to be non-structural and constrained because of time-invariant or exogenous factors in the system. Thus, when formulating policies or forecasting, it is not possible to identify-for example-the reasons underlying different cross-country reactions given an unexpected shock, the causality between real and financial variables, how additional transmission channels allow shocks to spill over, and how economic and institutional implications matter in driving shock transmission. Furthermore, many economic-financial issues show consistent and dramatic structural breaks associated with events such as economic recessions, currency crises, and epidemiology (such as the current COVID-19 pandemic). Such changes in regime are-even now-the core to some of the most important questions in macroeconomics. Recent empirical studies on the conduct of monetary policy along with fiscal consolidation periods have pointed out significant policy paradigm shifts and their implications over the past decades (see, e.g., Cette et al. (2016); Foerster and Matthes (2020); Stock and Watson (2012); Coibion and Gorodnichenko (2012); and Daly et al. ( 2016)). These shifts seem to have concerned strategy and policy development, monetary policy shocks, and policy interactions. Parallel to these shifts, both the volatility linkages between economic-financial data and the uncertainty in financial markets have changed substantially over time as well. For instance, the Great Recession prompted a renewed debate concerning objectives and conduct of policy issues due to structural changes that have taken place. In particular, slower growth led to questions about how the systematic conduct of policy tools should respond to these shifts. An another issue for changing fundamental policy instruments due to structural changes is also related to the economy that tends to shift between periods of high and low growth. These repeated regime changes may not occur very frequently, but needs to be dealt with when performing forecasting and policy-making. More precisely, dynamic models should account for the possibility of future shifts in the growth rate when setting policy. Indeed, if inter-linkages between real and financial sectors are changing how they react as the economy undergoes structural shifts, cross-country interdependencies and feedback effects matter and then need to be accounted for.
Recently, Pacifico (2019b) developed a structural version of the BVAR-labeled as Structural Panel Bayesian VAR (SPBVAR)-in order to deal with model misspecification and unobserved heterogeneity problems when jointly modeling and quantifying multicountry data using the information contained in a large set of endogenous and economicfinancial variables 1 . The SPBVAR model focuses on a state-space factorization structure, where the factors driving the coefficients of the multicountry VAR are restricted to evolve over time as random walks so as to: (i) reduce the number of parameters; (ii) allow for the evaluation of permanent shifts; (iii) investigate any type of coefficient factors via their interactions; and (iv) replace volatility changes by coefficient changes. This latter turns out to be highly restrictive to evaluate multiple time-varying change-points (or structural breaks) when studying macroeconomic and financial time-series. For example, international business cycle dynamics, policy interactions, and interdependencies and co-movements among different sectors and countries have changed substantially during the recent global crisis and successive consolidation periods. In addition, the increasing volatility and uncertainty in financial markets have confirmed the close volatility linkage between economic-financial data and thus the need to investigate shifts in either coefficients or volatility when describing these changes in a time-varying multicountry framework (see, e.g., Primiceri (2005); Canova and Gambetti (2009) ;Clark (2009) ;Cogley et al. (2010); and Sims and Zha (2006)).
My approach and empirical application aim to contribute to this debate. More precisely, this paper develops a modified version of the SPBVAR model allowing for multi-variate time-varying volatility to deal with three issues jointly: structural changes, policy regime shifts, and endogeneity issues. Specifically, this model takes the name of multicountry SPBVAR with Multivariate Time-varying Volatility (SPBVAR-MTV) and includes two additional components: a set of lagged endogenous variables to investigate fiscal and monetary policy implications and interactions, and-particularly-time-varying log-volatilities to deal with multiple structural breaks and policy regime shifts over time. Concerning the latter, I imply Autoregressive Conditional Heteroskedasticity in Mean (ARCH-M) model effects in the representation of the outcomes, and it would be an useful way of modeling time-varying conditional second moments providing an alternative to the stochastic volatility specification (see, e.g., Cogley and Sargent (2005); Carriero et al. (2019); and Clark (2011)). In contrast to SPBVAR, in the state-space factorization structure, the variance in error terms is allowed to be time-variant and volatility changes are not more replaced by coefficient changes. More precisely, let the framework be multidimensional (panel data analysis), I model the covariance matrix in time-varying log-volatilities interpreting it as having been broken into a collection of smaller matrices (blocks). Every collection is then estimated using a hierarchical (conjugate) prior specification strategy to allow cross-unit interdependent relationships. The computational costs involved in using that specification are moderate since the high dimensionality is avoided via Bayesian inference and Monte Carlo Markov Chain (MCMC) implementations. For instance, Kalman-Filter technique is used to get appropriate posterior distributions for time-varying coefficients and Metropolis-Hastings algorithm is used to draw posteriors for log-volatilities evolving over time. A Structural Normal Linear Regression (SNLR) model is obtained via Bayesian methods to work with smaller systems in which all the regressors are endogenous, observable, directly measured, and time-varying linear combinations of the right-hand variables of the SPBVAR-MTV model. I also account for three more indices 2 in order to quantify international spillover effects and thus evaluate their size (or intensity in terms of volatility) and dynamics (or spreading) among countries and sectors over time. (i) The Bilateral Net Spillover Effect (BNSE) is used to account for cross-unit interdependencies, feedback effects from the impulse variables, and temporary or persistent long-run effects of a potential shock (or excess spillover effects). (ii) The Systemic Contribution (SC) index is used to evaluate sequential features associated with systemic events. (iii) The Total Contagion Index (TCI) is addressed to investigate contagion measures in real economy and financial markets when dealing with both issues of endogeneity (because of omitted variables and unobserved heterogeneity) and volatility (because of policy regime shifts and structural changes).
In this paper, the SPBVAR-MTV model incorporates the econometric literature on standard Time-Varying Parameter Vector Autoregressions (TVP-VARs) with stochastic volatility, become a benchmark model for analysing and forecasting the evolving interrelationships between multiple macroeconomic variables (see, e.g., Koop et al. (2009);Koop (1996); Liu and Morley (2014);D'Agostino et al. (2013); and Clark and Ravazzolo (2015a)). Despite the empirical success of these flexible time-varying models, they show a relevant limit about their potential and feasible over-parametrization. More precisely, on the methodological side, the literature makes out two popular Bayesian methods for TVP-VARs with stochastic volatility: Marginal Likelihood (ML), evaluating how likely the observed data are occurred within the system, and Deviance Information Criterion (DIC), trading off between model fit and model complexity. As regards ML estimates, they are usually obtained by using the harmonic mean 3 of a conditional likelihood 4 that tends to have a substantial bias selecting the wrong model (see, for instance, Chan and Grant (2015) and Frühwirth-Schnatter and Wagner (2008)). Concerning DIC procedure, the MCMC integration based on the conditional likelihood tends to associate higher probability to the most complex models (overfitting 5 ).
The methodology proposed in this paper overtakes these limits by using analytical integrations for integrating out the time-varying volatilities. More precisely, integrated likelihood evaluation is achieved by integrating out the time-varying parameters analytically (e.g., Kalman-Filter technique), whereas the log-volatilities are integrated out numerically via importance sampling. The latter consists of two steps: (i) the Metropolis-Hastings algorithm is used to draw posteriors for time-varying log-volatilities from the proposal density distribution, and then (ii) the Newton-Raphson (N-R) method 6 algorithm is involved to find the maximum of the (log) conditional density. In this way, the computational costs are further reduced focusing on band and sparse matrix algorithms instead of the conventional Kalman filter.
An empirical application is developed by accounting for the Central, Eastern, and Western European (CEWE) countries in order to include a large pool of advanced and emerging economies, with particular emphasis to the Great Recession and successive postcrisis periods. The United States (US) are included in the analysis to assess international spillover effects and possible contagion measures among financial markets. In this study, I focus on the latest two alternative monetary policy regimes that have been in place since the 1990 (see, for instance, Kallianiotis (2019)): (1) the Inflation Stabilization Era (ISE) from 1994 to 2008 and (2) the Zero Interest Rate Era (ZIRE) from 2008 to 2015. I also consider two more additional periods: (1) 2006q1-2009q4 to investigate possible commonality between financial markets and real economy during the Great Recession and (2) 2010q1-2018q4 to evaluate fiscal implications and policy perspectives during post-crisis consolidation. Moreover, since most of countries joined in with Euro Area (EA), one is also able to investigate how policy regime shifts and endogeneity issues matter when studying macroeconomic-financial linkages. The analysis focuses on five main questions. First, I investigate how different economic-institutional characteristics affect the transmission of fiscal and monetary policy shocks among countries and sectors. Second, I investigate how policy interactions affect the benefits of consolidation among countries. Third, I evaluate how endogeneity and volatility issues affect inter-sector and inter-country linkages in panel setups. Fourth, I evaluate the role of policy regime shifts and their interactions when structural changes and contagion effects matter.
A simulated experiment-compared to related works-is also addressed to highlight the performance of the estimating procedure developed in this study using some Monte Carlo simulations. The findings prove that the hierarchical structural framework with timevarying log-volatilities perform better conditional forecasts when studying macroeconomicfinancial linkages with structural breaks and volatility changes.
The remainder of this paper is organized as follows. Section 2 introduces the econometric model, its features, and the estimation procedure. Section 3 describes the dynamic analysis focusing on prior assumptions strategy, posterior distributions, and MCMC implementations. Section 4 presents the data and the empirical analysis. Section 5 addresses a counterfactual assessment on macroeconomic-financial linkages in multicountry dynamic setups by investigating in depth how structural changes and policy regime shifts affect the spreading and the evolution of international spillover effects, with particular attention on triggering events and policy recommendations for decision makers. Section 6 provides a simulated example through Monte Carlo simulations to demonstrate the performance of the estimation method compared to some related existing approaches. The final section contains some concluding remarks.

Econometric Model
Compared to Pacifico (2019b), I extend and improve the standard version of the multicountry SPBVAR model in order to jointly account for time-varying parameters and multivariate volatility evolving over time.
Thus, the SPBVAR-MTV model developed in this study includes two additional components: (i) a set of lagged endogenous variables in order to assess different policy regimes and their interactions and (ii) time-varying log-volatilities to capture further evolving inter-relationships between multiple economic-financial data. The model has the form: where the subscripts i, j = 1, 2, . . . , N are country indices, t = 1, 2, . . . , T denotes time, L stands for the lag operator, a im,0 is an N M · 1 vector of intercepts for each i, the subscripts (m, q,q, ξ) denote the endogenous variables in each equation within the system observed for i, with ( , ,¨ ,φ) referring to the ones observed for j and independent of i, and ε im,t ∼ i.i.d.N(0, Σ t ) is an N M · 1 vector of disturbance terms. Given the SPBVAR-MTV in (1), all variables within the system are so defined. (i) Y im,t is an N M · 1 vector of outcomes to be predicted for each i for a given m. (ii) A im,j is an N M · N M matrix of coefficients for each pair of countries (i, j) for a given m, with m = 1, . . . , M, and Y im,t−λ is an N M · 1 vector of observed lagged variables for each i for a given m to address economic-financial issues. More precisely, stacking for m, I decompose it i,t−λ denoting lagged outcomes to capture the persistence and Y c i,t−λ including lagged control variables such as general economic-financial conditions. (iii) B iq,j is an NQ · NQ matrix of coefficients for each pair of countries (i, j) for a given q, with q = 1, . . . , Q, and W iq,t−λ is an NQ · 1 vector including a set of (directly) observed lagged variables for each i for a given q to evaluate economic-institutional interdependencies (such as additional transmission channels). (iv)B iq,j¨ is an NQ · NQ matrix of coefficients for each pair of countries (i, j) for a givenq, withq = 1, . . . ,Q, and W iq,t−λ is an NQ · 1 vector including a set of additional (directly) observed lagged variables for each i for a givenq to evaluate policy implications and interactions (such as policy tools). (v) C iξ,jφ is an NΞ · NΞ matrix of coefficients for each pair of countries (i, j) for a given ξ, with ξ = 1, . . . , Ξ, and Z iξ,t−λ is an NΞ · 1 vector including a set of observed lagged proxy 7 variables for each i for a given ξ to address further economic-financial linkages (e.g., territorial competitiveness and infrastructural system) and economic-institutional implications (e.g., competitiveness developments and macroeconomic imbalances). Here, all variables in the system are endogenous and time-varying.
The two main differences with respect to a standard SPBVAR lie in the additional componentB iq,j¨ Ẅ iq,t−λ and -particularly -in the variance-covariance matrix of the vector of innovations (ε im,t ). More precisely, Σ t = diag exp(h 1t ), exp(h 2t ), . . . , exp(h Nt ) , where h t = (h 1t , h 2t , . . . , h Nt ) denotes the time-varying log-volatilities stacked for i and following a random walk process: where Σ ht = diag(σ 2 h,1t , σ 2 h,2t , . . . , σ 2 h,Nt ) is a block diagonal covariance matrix and h 0 denotes the initial conditions to be estimated. More precisely, the matrix Σ ht is broken into a collection of smaller matrices (blocks). Such a specification implies Autoregressive Conditional Heteroskedasticity in Mean (ARCH-M) model effects in the representation of Y t , and it would be an useful way of modeling time-varying conditional second moments providing an alternative to the stochastic volatility specification. The main novelty in this study is that, in the state-transition Equation (2), the variance in v t and then volatility changes are allowed to be time-variant. The random-walk assumption in (2) is very common in the time-varying VAR literature, having the advantage of focusing on permanent shifts and reducing the number of parameters in the estimation procedure. In addition, the block diagonality of Σ ht guarantees the identifiability of the log-volatilities (h t ). The computational costs involved in using that specification are then moderate using hierarchical (conjugate) priors through Bayesian inference and obtain posterior distributions via MCMC algorithms.
In Equation (1), the dynamic relationships are allowed to be unit-specific and all the (potential) structural changes are allowed to vary over time. In addition, whenever the matrices A im,j (L), B iq,j (L),B iq,j¨ (L), and C iξ,jφ (L) differ 8 for some L, cross-unit lagged interdependencies and structural time variations matter, and then dynamic feedback and interactions among countries and variables are possible. Thus, the framework of the model (1) makes it able to connect the empirical results to the existing literature and contemporaneous developments when quantifying international business cycles, evaluating policy interactions, and performing conditional forecasting. Nevertheless, even if these features add flexibility to the specification, it is very costly making the number of coefficients of (1) very large. Indeed, they are increased by N(M + Q +Q + Ξ)l factors.
Let k = N[M + Q +Q + Ξ]l be the number of all matrix coefficients in each equation of the SPBVAR-MTV model for each pair of countries (i, j), a 1 · k vector X t = (I, Y im,t−1 , Y im,t−2 , . . . , Y im,t−l , W iq,t−1 , W iq,t−2 , . . . , W iq,t−l ,Ẅ iq,t−1 ,Ẅ iq,t−2 , . . . ,Ẅ iq,t−l , Z iξ,t−1 , Z iξ,t−2 , . . . , Z iξ,t−l ) can be defined containing all lagged (endogenous) variables in the system for each i. Then, I define an N Mk · 1 vector γ k it,j = vec(g k it,j ) containing all columns, stacked into a vector 9 , of the matrices A(L), B(L),B(L), and C(L) for each pair of countries (i, j) for a given k, with g k it,j = (A i1,j , A i2,j , . . . , A iM,j , B i1,j , B i2,j , . . . , B iQ,j ,B i1,j¨ ,B i2,j¨ , . . . ,B iQ,j¨ , C i1,jφ , C i2,jφ , . . . , C iΞ,jφ ) , and γ t = (γ 1t , γ 2t , . . . , γ Nt ) denoting the time-varying coefficient vectors, stacked for i, for each country-variable pair. With these specifications, I can express the model (1) in a simultaneous-equation form: where Y t = (Y 1m,t , . . . , Y Nm,t ) and E t = (ε 1m,t , . . . , ε Nm,t ) are N M · 1 vector containing the observable variables of interest and the random disturbances of the model for each i for a given m, respectively, andX t = (I N M ⊗ X t ) contains all lagged time-varying variables within the system stacked in X t . Now, because the coefficient vectors in γ t vary in different time periods for each country-variable pair and there are more coefficients than data, it is impossible to eliminate γ t . Thus, to avoid the curse of dimensionality, I adapt the framework in Pacifico (2019b) using a hierarchical factor structure to effectively transform the overparametrized multicountry Structural VAR in (1) into a parsimonious SNLR model. More precisely, I take a flexible Bayesian viewpoint to weakly restrict the coefficients to depend on a low-dimensional vector of observable and time-varying factors (the loadings), which correspond to linear combinations of the right-hand variables of (1), where the forecast errors feature a timevariant heteroskedastic structure. These factors would capture-for example-coefficient variations that are specific to a variable (variable-specific effects), specific to a country given a set of observable and directly measured variables (country-specific effects), common among countries and sectors (common effects), and specific to additional time-varying effects among countries and sectors directly affecting the outcomes in Y t (misspecified effects due to endogeneity issues). Such a reparametrization has three appealing features. (i) First, it reduces the problem of estimating high dimensional (potential) combinations of time-varying coefficients into the problem of estimating a smaller number of loadings on some linear combinations of the right-hand variables of (1). (ii) Second, since the loadings of the SNLR model are observable and time-varying linear combinations of the right-hand variables of (1), an estimable hierarchical structure is feasible and suitable for policy purposes and strategies. (iii) Third, representing the main thrust of this study, the parsimonious SNRL is able to identify additional common or heterogeneous effects between different countries and sectors that vary over time and directly affect the variables of interest in Y t .
In this study, I assume γ t to have the following factor structure: where F N Mk and dim(β f t ) dim(γ t ) by construction, G f = [G 1 , G 2 , . . . , G F ] are N Mk · κ f matrices obtained by multiplying the matrix coefficients (g k it,j ), stacked in the vector γ t , by conformable matrices D f with elements equal to zero and one, with κ f being a numerical index that depends on the typology of the factorization, u t is an N Mk · 1 vector of unmodeled variations present in γ t , and E(u t u t ) = Σ u = Σ e ⊗ V, with Σ e denoting the covariance matrix of the vector E t that includes time-varying log-volatilities and V = (σ 2 I k ) as in Kadiyala and Karlsson (1997). In this framework, unobserved heterogeneity and functional forms of misspecification are absorbed in the κ f · 1 time-varying coefficient vectors β f t . They are observable smooth linear functions of the lagged variables and thus can be easily estimated with a gain in efficiency and accuracy. The idea is to shrink γ t to a much smaller dimensional vector β t , with β t = (β 1t , β 2t , . . . , β Ft ) , containing all the regression coefficients stacked into a vector. In this way, further investigations (e.g., policy regime shifts and interactions, international business cycles, and economic-institutional linkages) can be performed. Finally, the factorization of γ t becomes exact as long as σ 2 converges to zero. In Equation (4), all factors are permitted to be time-varying, and thus time-variant structures can be obtained via implementations of MCMC algorithms. Moreover, time variations in the variance of shocks u t to the factors β f t are also allowed so that Y t can capture (potential) structural changes and policy shifts among countries and variables. Running Equations (3) and (4) for Equation (1), the factorization is: Given the factorization in Equation (5), the reduced-form SPBVAR-MTV model in Equation (3) can be transformed into a Structural Normal Linear Regression 10 model with an error covariance matrix of an Inverse-Wishart (IW) distribution 11 . By Equations (3) and (4), the SNLR model can be written as where χ f t ≡X t G f is an N M · κ f matrix that stacks all coefficients and their possible interactions in the SPBVAR-MTV model in (1), with χ t = diag(χ 1t , χ 2t , . . . , χ Ft ), and η t ≡X t u t + E t ∼ N(0, σ t · Σ u ) has a particular heteroskedastic covariance matrix that needs to accounted for, with σ t = (I N + Σ ht ⊗X tX t ).
To complete the specification, I suppose the following state-space structure for the time-varying regression coefficients: where β t = (β 1t , β 2t , . . .) , P t = diag(P 1t ,P 2t , . . . ,P Ft ) is a block diagonal matrix, andP f t = (p f t · I k ), where p f t controls the tightness (stringent conditions) of the factorization ( f ) of the time-varying coefficient parameters (β t ) in order to make them estimable. Moreover, the computational costs involved in using this specification are moderate since the dimension of β t is considerably smaller than the dimensionality of Y t , and the block diagonality of P t guarantees the identifiability of the factors (β t ). Here, some considerations on the innovations are in order: (i) the errors E t , u t , and v t are mutually independent; (ii) the error terms η t andṽ t are allowed to be correlated between them; and (iii) v t and η t are correlated between them by construction. Finally, if the factorization in Equation (5) is exact, σ 2 → 0 and one has to act on three competing models:

•
Model I (M I ): A benchmark model with no change-points, denoting the 'General Case'.
• Model II (M I I ): A benchmark model with change-points in the only log-volatilities, denoting the 'Special Case'.
Here, η t ∼ N(0,Σ) and depends on the only disturbances contained in E t , withΣ = σ t . The M I I refers to the case of structural breaks because of (potential) unmodeled dynamics 12 in γ t , with h 0 = 0, h t evolving over time, Σ t = Σ, and uncorrelatedness between the innovations η t andṽ t . • Model III (M I I I ): A benchmark model with change-points in either time-varying parameters or log-volatilities, denoting the 'Full Case'.
Here, η t ∼ N(0,Σ) and depends on the disturbances contained in v t and u t , witḧ Σ = σ t · diag(Σ e , Σ e , . . . , Σ e ). The M I I I refers to the case of structural breaks because of both unmodeled dynamics and policy regime shifts, with h 0 = 0, Σ t = Σ, and h t evolving over time.
Finally, once the (conditional) marginal likelihood 13 is obtained for any model, the exact and final solution can be obtained via MCMC integrations, corresponding to the best 14 model solution with higher log Bayes factor (lBF): The methodology does not include studies focused on Markov-switching dynamics, modeling covariance matrices of country-specific Markov chains 16 , because of this paper aims to extend and improve recent works developed when studying macroeconomicfinancial linkages in multicountry dynamic panel setups. Nevertheless, to solve potential overfitting problems and cross-unit unobserved heterogeneity, Markov-switching models follow similar hierarchical prior specification strategy proposed in those works, where their empirical results met with positive feedback in the empirical analysis.

Model Features
To illustrate the conformation of the time-varying SPBVAR-MTV in (1), the exact form and meaning of the β f t 's and D f 's in (4), and the selection of the conformable matrices D f , I suppose there are M = 2 endogenous variables and Q =Q = Ξ = 2 additional observable and proxy variables varying over time for every N = 2 countries. For convenience, I suppose one lag and no intercept. Thus, the SPBVAR-MTV in (1) assumes the form:  11,11 B 11,21 B 11,12 B 11,22 B 21,11 B 21,21 B 21,12 B 21,22 B 12,11 B 12,21 B 12,12 B 12,22 B 22,11 B 22,21 B 22,12 11,11 C 11,21 C 11,12 C 11,22 C 21,11 C 21,21 C 21,12 C 21,22 C 12,11 C 12,21 C 12,12 C 12,22 , vec(C iξ,jφ ) be the 64 · 1 vector containing all columns (stacked) of the matrices A(L), B(L),B(L), and C(L) for each pair of countries (i, j), with ( , ,¨ ,φ) denoting the variables in each equation observed for j and independent of i, and let Z 22,t−1 ) be the 1 · 16 vector containing all lagged variables in the system for each i for a given k, with Concerning the more recent studies of international business cycles and spillover effects (see, e.g., Ciccarelli et al. (2018), Canova et al. (2012), Pacifico (2019aPacifico ( , 2020a, it is typical for authors to define a country-specific indicator for Y t for each i, a variable-specific indicator for Y t for each m, and a common indicator for Y t for each pair of countries and variables (i, m).
According to the empirical implementation addressed in this study, I assess three additional terms in the factorization assuming time-varying log-volatilities in ε im,t : (i) an indicator to account for the role of additional transmission channels (e.g., in W iq,t−1 ) affecting spillover effects, the impact of macroeconomic-institutional interdependencies (e.g., in Z iξ,t−1 ) in driving the transmission of country-specific shocks, and policy interactions (e.g., inẄ iq,t−1 ) when dealing with policy regime shifts; (ii) an indicator to investigate heterogeneity, interdependence and commonality among all lagged variables (Y m,t−1 , W q,t−1 ,Ẅq ,t−1 , Z ξ,t−1 ) affecting shock transmission; and (iii) an indicator to highlight different reactions and co-movements among countries and sectors (Y im,t−1 , W iq,t−1 , W iq,t−1 , Z iξ,t−1 ) due to an unexpected common shock.

Dynamic Analysis
Before specifying prior assumptions and posterior distributions, I recall the state-space structure of the reparametrized SPBVAR-MTV model in (6):

Hierarchical Prior Setups and Assumptions
Supposing exact factorization in (5), in order to complete the model, I need to define prior moments on (Σ −1 e , Σ −1 ht , p f 0 , h 0 , β 0 ). Thus, collecting them in a vector φ 0 , with φ 0 = (Σ −1 e , Σ −1 ht , p f 0 , h 0 , β 0 ) being prior densities, the conditional likelihood function can be derived from the sampling density p(Y|φ 0 ) by using a mixture hierarchical distribution. In other words, (i) a Normal distribution for factors β and log-volatilities h; (ii) a Wishart distribution for Σ −1 e ; and (iii) an Inverse-Gamma distribution for Σ ht and p f , where p f = vec(P t ). That is, withΣ depending on the benchmark model: is the sum of the squared errors, withβ = (Σ t χ t χ t ) −1 · (Σ t χ t Y t ) referring to the OLS estimate of β, and S h denotes the least squares estimate of σ h based on the (satured) model, with σ h = vec(Σ ht ). Finally, the Equation (18) corresponds to the proposal distribution obtained by MCMC integration (such as Metropolis-Hastings algorithm 17 ).
Since the above-mentioned hierarchical prior specification strategy is affected by common or subjective beliefs because of the marginal effect of economic-financial variables, I assume Independent Normal-Wishart Prior so as to allow cross-equation independence of the coefficients distributions and remove the dependence of β on Σ e .
Given the state-space structure in Equations (14) and (15), MCMC methods and implementations (such as Gibbs sampling, Kalman Filter algorithm, and Metropolis-Hastings algorithm) can be computed numerically and joint distributions characterised analytically. The first step is to suppose that data run from (t = 0) to (t = T) in order to obtain a training sample (−τ, 0) and then to estimate the features of the priors. When such a sample is unavailable, it is just sufficient to modify the expressions for the prior moments in Equations (16)- (20) as: Here, N() stands for a Normal distribution, iW() denotes an Inverse-Wishart distribution, IG() indicates an Inverse-Gamma distribution, and F −1 refers to the information available at time −1. The prior for β 0 in (26) and the law of motion for the factors imply that: whereβ t−1|t−1 and R t−1|t−1 denote mean and variance-covariance matrix of the conditional distribution ofβ t|t , respectively. All hyperparameters are known. More precisely, collecting them in a vector δ, where δ = (z 1 , β 1 ,ω 0 ,S 0 ,S 0 , α 0 , V 0 ,β 0 , R 0 ), they are treated as fixed and are either obtained from the data to tune the prior to the specific applications (such as z 1 ,ω 0 , α 0 , andβ 0 ) or selected a priori to produce relatively loose priors (such as β 1 ,S 0 ,S 0 , V 0 , and R 0 ). In this context, the only fully Bayesian approach that leads to analytical results requires the use of a natural conjugate prior. According to Equations (16), (17), and (19), the natural conjugate prior has the form 18 : whereβ t|t and R t|t are hyperparameters collected in the vector δ, and ω α and S V are parameters to be estimated. If P t = 0, allowing for time-variant factors and volatilities, draws of p f t and σ 2 hi can be taken from Normal-Inverse-Gamma distributions.
According to the natural conjugate prior (28), β t depends on Σ e and Σ ht . Thus, β t , Σ e , and Σ ht are not independent of one another. To allow different equations in the VAR to have different explanatory variables, previous specifications have to be modified. More precisely, given the SNLR model in (14), general priors that do not involve the restrictions inherent in the natural conjugate prior are the Independent Normal-Wishart (INW) and the Independent Inverted Gamma (IIG) distributions. The latter has different scale and shape parameters with respect to p f t and is obtained by maximum likelihood estimates. Thus, the natural conjugate priors (28)-(30) can be re-written as: Here, the hyperparametersω α andŜ V denotes scale and shape parameters, respectively, collected in the vector δ. Moreover, the prior (32), with P t = 0, allows for the prior covariance matrix R t|t to be anything the researcher chooses, rather than the restrictive (Σ e |Y T ⊗ R t|t + Σ ht ) form of the natural conjugate prior.

Conditional Likelihood and Kalman Filter Technique for Time-Varying Parameters
The posterior distributions for φ = (Σ −1 e , Σ −1 ht , p f t , h t , {β t } T t=1 ) are calculated by combining the prior with the (conditional) likelihood for the initial conditions of the data. The resulting function is then proportional to where ) refers to the unknowns whose joint distribution needs to be found, with φ −k standing the vector φ, excluding the parameter k.
Despite the dramatic parameter reduction obtained with Equation (14), the analytical computation of posterior distributions (φ|Y T ) is unfeasible. Thus, a variant of the Gibbs sampler approach-such as Kalman-Filter technique-can be used through MCMC integrations. More precisely, for the conditional posterior distribution of (β 1 , . . . , β T |Y T , φ −β t ), it gives the following forward recursions for posterior means and the covariance matrix, respectively:β where Starting from β T|T and R T|T , the marginal distributions of β t can be computed by averaging over draws in the nuisance dimensions, and the Kalman filter backward can be run to characterise posterior distributions for φ: The Equations (42) and (44) refer to the variance-covariance matrix of the conditional distribution ofβ t|t+1 and the Generalized Least Square (GLS) estimator, respectively. By rearranging the terms, Equation (41) can be rewritten as whereβ t|t+1 and R t|t+1 denote the smoothed one-period-ahead forecasts of β t and of the variance-covariance matrix of the forecast error, respectively. The above output of the Kalman filter is used to generate a random trajectory for {β t } by using the backward recursion starting with a draw of {β t } from N(β T|T , R T|T ) 19 . Given (40), the other posterior distributions can be defined as: Here, some considerations are in order.
In Equation (48),ã h,m = α 0 ·ω is obtained by Metropolis-Hastings algorithm 21 and V h,m = V 0 +Σ h * is computed by MCMC-based EM algorithm 22 , with α 0 and V 0 denoting the arbitrary degree of freedom and the arbitrary scale parameter, respectively, andΣ h * = (Σ 1T , . . . ,Σ NT ) referring to the estimated covariance matrix for each i in a regimem given t. In this analysis,ã h,m is constructed to be close to zero, withω ∼ = 0.1 · exp(N M), and h * is an arbitrary vector givenm regimes according to the state-transition Equations (2) and (15).
In Equation (49), , withS 0 and β f t denoting the arbitrary scale parameter and the f th subvector of β t , respectively. In this analysis,S 0 ∼ = 0.1, f refers to the factors described in Equation (5), and k denotes the number of all matrix coefficients in each equation of the SPBVAR-MTV model in (1).
Finally, the last two hyperparameters to be defined in the vector δ areβ 0 =β 0 , witĥ β 0 denoting the Ordinary Least Squares (OLS) estimates of Equation (14), and R 0 = I k .

Metropolis-Hastings Algorithm for h it
Supposem regimes, withm = 0, 1, . . . , s, and use Metropolis-Hastings algorithm to draw posteriors for h it from the proposal density distribution δ * (h it ), with probability αm equals: According to the SNLR described in Equation (14), let β * t denote the time-varying coefficient vectors whenm = 0 (some forms of break occur), the probability function takes the form: where In (52), α(β * t , β t ) is the probability to accept or reject a draw 23 . In addition, since the posterior distribution corresponds -by construction -to a multivariate normal distribution, the Optimal Acceptance Rates (OARs) 24 are:

Analytical Integration for Integrating out the Time-Varying Volatilities
Given the proposal density distribution δ * (β * t , β t ) in Section 3.2.2, with probability α(β * t , β t ) in (52), one needs to integrate out h it using importance sampling. More precisely, in this study, I approximate the log conditional marginal density log[p(h|Y, P t , Σ ht , β 0 , h 0 )] by using a Gaussian density, which is then used as the importance sampling density. Thus, the Expectation-Maximization (EM) algorithm is used to find the maximum of the log conditional marginal density and consists of two steps: the Expectation step (E-step) and the Maximization step (M-step) 25 .

Expectation Step (E-Step)
The E-step is implemented by computing the following conditional expectation: where the expectation is taken with respect to p(β|Y, h * , P t , Σ ht , β 0 , h 0 ) for an arbitrary vector h * as Kroese and Chan (2014). More precisely, given the SNLR described in (14) and letP β denote the first difference matrix, the state-transition Equation (15) can be rewritten as:P whereᾱ β = (β 0 , 0, . . . , 0) , Υ β = I T ⊗ P t , andP β is a lower triangular matrix of dimension k. Thus, the conditional marginal density for β would be distributed normally according to standard linear regression results: Here, the precision sampler of Chan and Jeliazkov (2009) can be used to sample from N(β, Φ −1 β ) efficiently. In other words, the mean vectorβ and the precision matrix Φ β are computed using h * .
The expectation in (54) can then be written in terms of an explicit expression: where tr() is the trace operator, c is a constant independent of h, 1 NT is a vector of ones, Π h is a lower triangular matrix of dimension N, α 0 ∼ = Π −1 h ·α 0 , withα 0 = (h 0 , 0, . . . , 0), and ε i = Y − (X t G)β is the error term.

Maximization Step (M-step)
The M-step consists of maximizing the function Ψ(h|h * ) with respect to h by using the Newton-Raphson method 26 . Thus, the gradient (g Ψ ) and the Hessian (H Ψ ) are, respectively: where ' ' refers to the entry-wise product andθ = (s 2 1T +ε 2 1T , . . . , Here, H Ψ is negative definite for all h, ensuring fast convergence of the N-R method, and Φ β , since it is a band matrix, guarantees that its Cholesky factor L β can be obtained without further effort. More precisely, given Φ β = L β · L β , with L β denoting a lower triangular matrix, G = L β · (X t G) can be obtained by solving the linear system L β G = (X t G) for G. Thus, the diagonal elements of (X t G) · Φ −1 β · (X t G) will be the row sums of the squares of G.
Repeat steps (A)-(D) until some convergence criterion is met at the OARs in (53).
Thus, terminate the iteration and set h ν = h (m,ν−1) , denoting that a certain changepoint among time-varying coefficient vectors and log-volatilities has been assessed correctly.

Data Description and Empirical Model
The SPBVAR-MTV model in (1) contains 17 country-specific models, including the United States, 8 Central-Eastern Europe (CEE) economies 27 and eight Western Europe (WE) economies 28 . The CEE and WE countries-except for SL-also refer to European emerging and advanced economies, respectively. Moreover, all the European countries are Eurozone members, with the exception of CZ, HU, and PO, and thus inter-sector and inter-country linkages can be investigated in depth.
The dataset contains the following collection of variables (Table A1). (i) Six endogenous variables are involved to describe real economy (real it,j ) and financial markets ( f in it,j ): three real variables to capture real business cycles (general government spending, gross fixed capital formation, GDP growth rate) and three financial variables to highlight the situation in the lending markets (bank leverage, flow of credit into economy, inflation rate). (ii) Bilateral flows of trade (rweights it,j ) and financial transactions ( f weights it,j ) are used to deal with endogeneity issues when studying international spillover effects among countries and variables. (iii) Three policy variables (policy it,j ) are used to investigate monetary and fiscal policy implications and interactions among countries and sectors (international interest policy rate, general government debt, and current account balance). (iv) Five (directly) observed variables are used as proxy variables (structures it,j ) to evaluate economic-institutional implications in driving the evolution of international spillovers and transmission of shocks over time among countries and variables: three indicators to deal with internal imbalances (financial consumption expenditure, private sector consumption, change in unemployment rate); one indicator to capture competitiveness developments and catching-up effects (nominal labour cost); and one indicator to monitor the probabilities of transitions between expansion/recession phases in business cycles and potential macroeconomic imbalances (house price indices). (v) The real GDP per capita in logarithmic form (productivity) is used to evaluate the size and the spreading of international spillover effects over time among countries and sectors given an unexpected shock. The weights it,j 29 , the policy it,j , and the structures it,j components are treated endogenously and used to jointly deal with endogeneity issues, structural changes, and policy regime shifts.
The series are expressed in standard deviations with respect to the same quarter of the previous year (q t /q t−1 ), and seasonally and calendar adjusted. All variables are used in year-on-year growth rates and all data comes from OECD data source.
The estimation sample covers the period from December 1994 to December 2018. It amounts, without restrictions, to 26, 384 regression parameters. More precisely, each equation of the time-varying SPBVAR-MTV in (1) has k = [17(6 + 2 + 3 + 5)] × 1 = 272 coefficients, and there are 97 equations in the system. Given the structural conformation of the model and a sufficiently large number of quarters describing economic-institutional and policy implications, it is able to capture: (i) endogeneity issues because of unobserved heterogeneity and misspecified dynamics across the sample; (ii) interdependency, commonality, and homogeneity because of potential international macroeconomic-financial linkages among countries and sectors; and (iii) relevant monetary and fiscal policy interactions, and contagion measures.
Given the factor structure in (4), I assume that the coefficient vector γ t depends on 10 factors in order to investigate sufficient macroeconomic-financial linkages and economic-institutional implications when performing conditional forecasting and policymaking. Thus, where, stacking for t, β f = (β 1 , β 2 , β 3 , β 4 , β 5 , β 6 , β 7 , β 8 , β 9 , β 10 ) contains all time-varying coefficient vectors to be estimated. Given the factorization in (63), the SNLR model in (14) can be written as: According to diagnostic tests (Table 1), the marginal (conditional) likelihood estimation confirms the exact γ t 's factorization in (63) and the estimates are asymptotically consistent given the absence of serial correlations across the residuals. Thus, the specified factors in (64) can be made clearer and estimated in terms of posterior means.
The indicators χ 1t β 1t and χ 2t β 2t are N M · 1 vectors of observable country-specific indicators for Y t , and account for the only real it,j and f in it,j components, respectively, in order to evaluate international spillover effects and transmission of shocks among countries in real economy and financial markets.
The indicators χ 3t β 3t and χ 4t β 4t are N M · 1 vectors of observable country-specific effects for Y t , and account for one additional component: (1) real it,j with policy it,j and (2) f in it,j with policy it,j . They are able to investigate monetary and fiscal policy implications and interactions among countries in the real and the financial dimensions, respectively.
The indicators χ 5t β 5t and χ 6t β 6t are N M · 1 vectors of observable country-specific effects for Y t , and account for two components further: (1) real it,j with rweights it,j and policy it,j , and (2) f in it,j with f weights it,j and policy it,j . They are able to jointly evaluate how interna-tional transmission channels and policy issues affect the size and spreading of spillover effects given an unexpected shock among countries in real economy and financial markets, respectively.
The indicators χ 7t β 7t and χ 8t β 8t are N M · 1 vectors of observable country-specific effects for Y t , and account for one component further: (1) real it,j with rweights it,j , policy it,j , and structures it,j , and (2) f in it,j with f weights it,j , policy it,j , and structures it,j . They are able to assess the role that international transmission channels, macroeconomic-institutional implications, potential macroeconomic imbalances, and policy implications play in allowing shocks to spill over among countries in real and financial sectors, respectively.
The indicator χ 9t β 9t is a N M · M v vector of observable cross-country variable-specific effects for Y t , where M v = (M v1 , M v2 , M v3 , M v4 ) denotes the number of variable groups: (i) M v1 = real it,j and policy it,j ; (ii) M v2 = f in it,j and policy it,j ; (iii) M v3 = real it,j , rweights it,j , policy it,j , and structures it,j ; and (iv) M v4 = f in it,j , f weights it,j , policy it,j , and structures it,j . The variable-specific factor is able to investigate endogeneity issues, policy regime shifts, and multivariate structural breaks among variables in real economy and financial markets.
Finally, the indicator χ 10t β 10t is a N M · M c vector of observable common effects for Y t , where M c = (M c1 , M c2 ) denotes the number of common groups: (i) M c1 = real it,j , f in it,j , and policy it,j and (ii) M c2 = real it,j , f in it,j , weights it,j , policy it,j , and structures it,j . The common factor is able to assess idiosyncratic spillover effects due to different reactions or co-movements among countries and variables for a given common unexpected shock in the real and financial dimensions. This latter (M c2 ) is also used to investigate and then quantify contagion measures during triggering events and policy regime shifts. Dynamic analyses have been conducted via accurate MCMC implementations. The total number of draws was 5000 + 1000 = 6000, which corresponds to the sum of the final number of draws to discard and draws to save, respectively. A total of 1000 draws has been used to conduct posterior inference at each t. The outcomes absorb the conditional forecasts computed for a time frame of 9 quarters (2 years and a quarter) in order to also address potential findings concerning the impact of an ongoing pandemic crisis on the global economy. The natural conjugate prior refers to four subsamples: (i) 1994q4-2008q3 and (ii) 2008q4-2015q4 in order to evaluate how monetary policy regimes affect the dynamics of the GDP growth; and (iii) 2006q1-2009q4 and (iv) 2010q1-2018q4 in order to highlight the impact of the most recent financial crisis and fiscal consolidation when investigating international spillover effects.
According to the log Bayes Factor in (8) and the exact factorization in (63), most of the time-varying estimated coefficient vectors (β f t ) embrace the 'Full Case' (M I I I ), where structural changes and policy regime shifts hold in either time-varying parameters or log-volatilities (see Table 2). It accounts for: two of the country-specific factors (χ 7tβ7t and χ 8tβ8t ); the cross-country variable-specific factor (χ 9tβ9t ) belonging to the variable groups M v3 and M v4 ; and the common factor (χ 10tβ10t ) belonging to the common group M c2 . All remaining empirical results embrace the 'Special Case' (M I I ), except for two factors concerning the 'General Case' (M I ). They correspond to the first two country-specific indicators (χ 1tβ1t and χ 2tβ2t ). These findings highlight the performance and then the potential of the SPBVAR-MTV model pointing out that: (i) change-points and policy regime shifts need to be taken into account when dealing with macroeconomic-financial linkages in multicountry dynamic panel setups; (ii) multiple structural changes in time-varying log-volatilities occur when evaluating international transmission channels and policy implications among countries and sectors in both the real and the financial dimensions; and (iii) change-points and policy regime shifts in either time-varying coefficients or logvolatilities occur when accounting for economic-institutional implications to investigate unobserved heterogeneity and misspecified dynamics among country-and variable-specific factors and common features.  9,1tβ9,1t , χ 9,2tβ9,2t 6 ≤ lBF ≤ 10 lBF >10 2 ≤ lBF ≤ 6 χ 9,3tβ9,3t , χ 9,4tβ9,4t 0 ≤ lBF ≤ 2 6 ≤ lBF ≤ 10 lBF >10 χ 10,1tβ10,1t 6 ≤ lBF ≤ 10 lBF >10 2 ≤ lBF ≤ 6 χ 10,2tβ10,2t 0 ≤ lBF ≤ 2 6 ≤ lBF ≤ 10 lBF >10

Time-Varying Factors
The first column denotes the time-varying factors and the other three columns refer to all three (potential) best benchmark models. The best model solution is highlithed in bold and corresponds to the highest log Bayes Factor according to the generalized version of the Kass and Raftery (1995)'s scale of evidence in (53).

Macroeconomic-Financial Linkages with Structural Changes and Policy Regime Shifts: A Counterfactual Assessment
The aim of this empirical analysis is to improve the existing literature on macroeconomicfinancial linkages in multicountry dynamic panel setups when dealing with either endogeneity or volatility issues. Thus, the SPBVAR-MTV model is appropriate to be used for investigating in depth how structural changes and policy regime shifts affect the intensity and the evolution (or dynamics) of international spillover effects among countries and sectors.
More precisely, intra-CEWE dynamics are assessed in four related contexts. (i) Firstly, international spillovers and policy issues are evaluated in an international and broader European setting (Section 5.1). (ii) Then, the empirical results are reevaluated accounting for additional time-varying factors to address endogeneity issues because of unobserved heterogeneity and misspecified dynamics (Section 5.2). (iii) A deepened investigation is further conducted accounting for multivariate change-points and policy regime shifts during different phases of financial cycles, where many emerging market economies have experienced a large surge of capital inflow, following the notably expansionary monetary policies of major advanced countries, and fiscal consolidation adjustments, playing a central role in the disinflation process (Section 5.3). (iv) Finally, policy implications and suggestions for decision makers are addressed according to all aforementioned findings (Section 5.4).

International Spillovers and Policy Issues among CEWE Economies
In Figure 1, where I consider the first two country-specific indicators (χ 1t , χ 2t ) for the overall sampled time-series, all CEE economies tend to be net receivers (or inward spillovers) in the real dimension and thus would be affected by the conditional impulse responses received from the European advanced countries (net senders). Overall, the size of the spillover effects is larger in the financial dimension because of highly strong cross-country interdependencies. These results find confirmation in previous related works such as Pacifico (2019bPacifico ( , 2019aPacifico ( , 2020a and Curcio et al. (2020).
However, contrary to them, the findings highlight that a consistent cross-country heterogeneity across the spillovers' dynamics would matter more in financial markets (Figure 1b), while a persistent degree of homogeneity and larger co-movements among countries tend to occur in the real dimension despite stronger inter-country linkages in the financial one (Figure 1a). The results confirm the presence of potential functional form of misspecifications that need to be investigated thoroughly when studying macroeconomicfinancial linkages.
(a) χ 1tβ1t Factor-Overall period (b) χ 2tβ2t Factor-Overall period In Figure 2, I account for the two country-specific indicators dealing with policy issues and their interactions (χ 3t , χ 4t ). In contrast to the previous results, most CEE economies tend to be net senders (outward international spillovers) in their real dimension (Figure 2a). Cross-country heterogeneity follows to be consistent and stronger in real economy and even more in financial markets (Figure 2b). In addition, larger commonality and homogeneity matter across the spreading and the intensity of spillover effects. The findings confirm the importance to account for either endogeneity and volatility issues.
From a global perspective, the same dynamic behaviour is observed in the transmission of US financial shocks, with outward spillover effects. The results are consistent and robust with the more recent literature on multicountry dynamic panel setups. More precisely, they confirm that US seem to be an important driver in allowing unexpected shocks to spill over and thus affecting European financial markets, mainly regarding CEE economies with inward spillovers. Then, intra-country shocks directly affect a country's own output growth in the real economy because of consistent cross-country interdependencies.
(a) χ 3tβ3t Factor-Overall period (b) χ 4tβ4t Factor-Overall period Established that structural changes and policy shifts affect macroeconomic-financial linkages among countries in an international and broader context, I consider the first two group-variable factors (χ 9,1t , χ 9,2t ) in order to examine in depth how monetary policy regimes and fiscal implications drive international shocks among real and financial sectors ( Figure 3). Here, the countries are grouped in three clusters: WE, CEE, and BLS 30 .
During the ISE Regime, most countries tend to be net receivers and net senders in the real and the financial dimensions, respectively (Figure 3a). Moreover, larger homogeneity in the spreading and the intensity of international spillovers would matter more among CEE economies given an unexpected financial shock. From a policy perspective, since in that period (1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008) the only country joined in with EU was Slovenia, the results highlight that the transmission of shocks among sectors are mainly affected by highly strong cross-country interdependencies rather than policy implications (e.g., because of high persistent inflation among emerging and then CEE countries).
During ZIRE Regime, emerging economies become net senders in real economy with larger spillover effects than advanced economies (Figure 3b). In financial markets, international spillovers show higher intensity than real economy due to a deeper state of severe fiscal reforms. In addition, higher co-movements among sectors in both the real and the financial dimensions matter than the ISE period because of persistent financial measures to foster the output stabilization.
From a modeling perspective, inward spillovers during the ISE Regime highlight that CEE countries are less competitive than WE economies, requiring appropriate emergency programs in order to be up against triggering events. From a policy perspective, outward spillovers during the ZIRE Regime in real economy point out that more stringent fiscal constraints would need to support developing economies in absorbing the effects of unexpected financial shocks (misspecified dynamics).

Unobserved Heterogeneity and Misspecified Dynamics Accounting for Additional Time-Variant Factors
Accounting for additional time-variant factors in order to investigate in depth policy regime shifts and structural breaks along with endogeneity issues, relevant empirical results and policy perspectives are derived ( Figure 4).
As regards weights it,j component (standing for omitted factors), most countries follow to show inward and outward spillovers in the real and the financial dimensions, respectively. Despite a consistent heterogeneity persists in their own output growth responses, larger co-movements matter accounting for additional shock transmission channels in real economy and even more in financial markets because of stronger cross-country financial linkages (Figure 4a,b). From a global perspective, outward spillovers in US confirm the importance about international spillover effects affecting European financial shocks (see, for instance, Pacifico (2020a) and Curcio et al. (2020)). From a modeling perspective, capital flows tend to matter more than trade flows in allowing shocks to spill over among countries (see, for instance, Pacifico (2019bPacifico ( , 2020a). However, higher intra-CEWE heterogeneity in the financial dimension, in terms of spillovers' intensity and spreading, emphasizes more consistent difference among financial markets due to tighter monetary policies.
Concerning structures it,j component (standing for unobserved heterogeneity), the intensity of spillover effects tends to increase confirming that economic-institutional linkages significantly affect countries' responses (Figure 4c,d). Cross-country commonality would be larger in real economy and thus if capital flows tend to matter more in driving shock transmission among financial markets, trade flows would matter more in affecting the spreading of spillover effects among countries. Moreover, output responses over time are larger in WE countries despite catching-up effects in CEE economies due to persistent and consistent cross-country heterogeneity. From a policy perspective, the results face a situation of trade-off. More precisely, if on one hand the adoption of sounder macroeconomic policies and economic-institutional changes-put in place to foster consolidated policy actions-have helped to bring inflation in emerging (and then CEE) economies back under control, on the other hand, in case of a noteworthy unexpected financial shock-without appropriate coordinated structural reforms in trade, product, and labour markets-outward government benefits will be not able for supporting the process of international financial integration among countries and boosting the output to potential growth.  Established that policy regime shifts and endogeneity issues (both omitted factors via additional transmission channels and unobserved heterogeneity via economic-institutional linkages) affect the spreading and the intensity of international spillovers, the same analysis is conducted by focusing on the last two cross-country variable-specific factors (χ 9,3t , χ 9,4t in Figure 5).
In this context, some main considerations are in order. First, during ZIRE Regime, CEE and BLS countries from net senders become net receivers in the financial dimension. It highlights that, even if substantial structural reforms in terms of radical fiscal adjustments were able to absorb unexpected financial shocks (outward spillovers), consistent crosscountry interdependencies among financial sectors-because of EA's common monetary policy-brought about 'pseudo-shock' in the short term to catch up with the economic growth of the other euro partecipants 31 (inward spillovers). Second, structural-institutional implications along with policy reforms affect the intensity (or volatility) of spillover effects in CEE and even more in BLS countries-because of larger current account deficits and lower real economic convergence-via international transmission channels, that allow in turn financial shocks to spill over. Third, persistent cross-country heterogeneity during monetary policy regimes emphasizes that the fairly well synchronized business cycles among emerging and advanced economies might be unlikely, mainly on account of triggering events in the long run. Thus, the increasing need of consistent reforms of the international financial system to accelerate well-suited financial integration in developing countries. These findings are against existing studies that support similarity across business cycles in CEWE economies because of dealing with too short periods. For instance, they consider up to seven years or less, implying that only a single business cycle would be covered by the available data.

Policy Interactions, Common Features, and Contagion Measures among Countries and Sectors
In the aftermath of the Great Recession and an ongoing postcrisis consolidation, the intensity of spillover effects becomes larger in real and even more in financial dimension because of stronger inter-country linkages among financial markets behind stringent fiscal adjustments (Figure 6a,b). More precisely, the spreading and the size of spillover effects tend to be higher among CEE and even more among BLS countries-because of extensive reforms-in real economy due to radical policy actions, and among WE countries in financial markets due to stronger interdependencies. Despite a consistent homogeneity holds among CEWE economies, different countries' responses matter during financial crisis and even more fiscal consolidation periods due to coordinated but not fairly flexible fiscal actions, mainly among emerging economies suffering from lower competitiveness.
Finally, according to all aforementioned findings, I compute the Total Contagion Index (TCI) on the only common indicator χ 10t,2β10t,2 , so as to investigate in depth commonality among sectors and countries in real economy and financial markets. To do it, the cumulative impulse responses are restricted in the interval [0, 1] and the (individual) spillover effects are restricted in the interval [−1, +1] so that the index will be bound between 0 and 100 (or between −100 and 0 if negative effects occur). Thus, the TCI is so obtained: where, IR y i →y j denotes individual (out) spillover effects and N −ṽ refers to the degrees of freedom depending on the needs of the investigation, withṽ accounting for the terms chosen in the factorization (63).
Here, some considerations are in order (Figure 6c). During crisis period, emerging and advanced economies show inward and outward spillovers in the financial dimension, respectively. Contrary to postcrisis consolidation periods, where CEE and BLS countries become net senders and WE countries are net receivers. It highlights the presence of consistent policy interactions: an unexpected shock in financial markets (e.g., because of inflationary pressures, unsustainable credit boom, stiffening of banking supervision) affects the real economy through fiscal adjustments (e.g., public expenditure cuts, lowers increasing tax). Then, stringent economic-institutional linkages cause a 'pseudo-shock' among CEE and BLS economies because of larger fiscal adjustments-mainly in the last two decade-to catch up with the economic growth of the other advanced EA economies (from net receivers to net senders).
Finally, to better highlight additional findings can be drawn because of the hierarchical structural framework used in this analysis, I address a cross-country spillover analysis supposing an unexpected volatility change in financial markets (Figure 7a). Let the estimation sample cover the period from December 1994 to December 2018, and let the current coronavirus outbreak have been officially declared a global pandemic crisis in 2020 32 , the aim of this analysis is to emphasize the performance of the estimating procedure of forecasting potential monetary policy strategies-along with fiscal tools-during (unexpected) dramatic structural breaks. Thus, I use conditional projections to include forecasts from 2019q1 to 2021q1, and focus on the cross-country indicators dealing with macroeconomic-financial linkages (χ 2tβ2t ), policy shifts and international transmission channels (χ 6tβ6t ), and endogeneity issues and policy shifts (χ 8tβ8t ).
The findings highlight three different scenarios ( Figure 7). (i) First, without accounting for unobserved effects and policy implications, high-income countries (advanced economies, WE) have better a credit rating than low-income countries (emerging economices, CEE and BLS). Let spillover effects among countries be positive (outward spillovers) and heterogeneous per intensity and spreading (high-income countries' responses larger than the low-income ones), the results highlight that a consistent crosscountry heterogeneity across the spillovers' dynamics persists in financial markets despite stronger interdependencies (Figure 7a). They confirm the presence of potential functional form of misspecifications (due to endogeneity and/or volatility issues) that need to be investigated thoroughly when studying macroeconomic-financial linkages. (ii) Second, when accounting for fiscal policy tools and not directly observed (omitted) factors, the results find confirmation with the recent literature 33 (Figure 7b). More precisely, advanced economies have better a credit rating than emerging ones. That is, a country's credit rating affects its ability to pursue expansionary fiscal policies (outward spillovers in credit) and countercyclical fiscal policies are not common in countries with higher credit risk (inward spillovers in cpi). Compared to Romer and Romer (2018), countries with lower debt-to-GDP (larger debit sustainability) tend to use fiscal policy more aggressively during crisis (inward spillovers in mkt per WE). (iii) Third, jointly dealing with endogeneity and policy issues, the results contrast the previous findings (Figure 7c). During dramatic structural breaks, either advanced or emerging economies would address unconventional monetary policy measures increasing monetary policy space to help central banks meet their output and inflation goals (inward spillovers in cpi), mitigating limitations to monetary transmission that may hamper the provision of credit where it is most needed (outward spillovers in credit), and supporting liquidity in financial markets or expanding fiscal space (inward spillovers in mkt). Even if low-income countries with more developed capital markets and effective transmission via interest rates should be more likely to benefit more from unconventional monetary policy, the results show smaller benefits than in advanced economies (lower spillover effects). Thus, in emerging economies, countercyclical fiscal policy has been conducted but with delay. The limit to resort to fiscal policy during recession in low-income countries is due to their limited ability in using traditional monetary tools. These findings reflect the recent reports on the conduct of monetary policy during the coronavirus pandemic crisis (e.g., European Central Bank, 19 October 2020 andInternational Monetary Fund, 23 September 2020). In the next section, based on the full estimation sample, I emphasize the results conducting a further analysis on the economic outlook amid the COVID-19 pandemic shock.

Lessons and Matters for Future Policy Efforts
In summary, all aforementioned results lead to four important chain-effect findings: (i) given an unexpected shock in financial markets, countries' responses show higher heterogeneity among real sectors, pointing out non-homogeneous real economic convergence among countries (endogeneity issues); (ii) the related spillover effects show larger intensity among developing economies due to sever overheating periods, mainly during the recent financial crisis due to radical structural fiscal adjustments (policy-regime shifts); (iii) at the same time, higher intensity in the spreading of spillovers bring about larger volatilities among either sectors or countries over time (structural changes); and (iv) these volatility issues highlight an unlikely international business cycle synchronization among emerging economies and thus a solid but not properly achieved integration within EU, increasing the cost of partecipation in the European and Monetary Union (EMU).
Since developing countries tend to bear the brunt of triggering events due to their relatively low economic weight (in terms of international trade exposures), 'quasi-flexible' policies should be conducted in order to ensure in a not-too-distant future: (i) the restora-tion of the confidence in financial systems, still recovering from the recent financial crisis; (ii) higher homogeneity across countries' responses in real economy given an unexpected financial shock so as to safeguard the inter-country real convergence; and (iii) stronger cross-correlations among CEWE economies when facing international shocks transmission.
In this context, 'quasi-flexible' policies stand for coordinated structural policy actions among foreign and domestic sectors along with more pointed fiscal adjustments according to country-specific requirements. Furthermore, the analysis highlights that, in case of a noteworthy unexpected shock in real economy, outward government benefits would be really beneficial for supporting the European integration and boosting the output to potential growth. Thus, the need of examining international spillovers accounting for both model misspecification problems and implied volatility changes.
In Figure 8, I display the generalized Entropy Index from 1994q4 to 2021q1. It corresponds to the Theil's Entropy, calculated by weighing the GDP with the population in terms of proportions with respect to the total, and can be used to measure the degree of divergence and economic inequality among countries. Here, forecasts from 2019q1 to 2021q1 correspond to conditional projections of each variable drawn in the SPBVAR-MTV in (1) and thus are able to point out the impact of an ongoing pandemic crisis on the global economy.
The coronavirus (or COVID-19) pandemic is a major global crisis negatively affecting sustainable development, economic growth, and stability and security across the globe. It constitutes an unprecedented challenge with very severe socio-economic consequences and highly strong deterioration of already existing humanitarian crises. In this study, the findings confirm the radical decrease of the economy in the last two quarters of the current year pursuant to the pandemic of coronavirus disease. However, a hint of the economic recovery shows up among countries in the next quarter. Thus, some considerations can be addressed. (i) First, coordinated and radical policy actions are necessary to deal with health emergency needs, support inter-country economic activity, and face the ground for the recovery. (ii) These adjustments should be implemented combining short, medium and long-term initiatives, but taking into account the dynamics of international spillovers and the cross-country economic-financial linkages so as to preserve confidence, stability, and financial integration (where highly strong heterogeneity and volatility matter). (iii) Moreover, even if several measures have already been taken at the national and EU levels, temporary and targeted discretionary fiscal stimulus have to keep on being adopted in a coordinated manner. More precisely, public resources and structural reforms in trade, product, and labour markets have to be directed to strengthen the healthcare sector and support affected economic-financial sectors. (iv) As regards monetary policy, closed resolute actions have to be taken by the European Central Bank to support liquidity and finance conditions to households and banks in order to preserve the smooth provision of credit to the economy. (v) Finally, to overcome the financing pressures faced by banks and households, all these policy adjustments need to be implemented by closely monitoring the evolution of the situation in each country and coordinating country-specific European and national measures. However, if an increasing degree of divergence should overlook among countries, further and different actions, including legislative measures, will have to be taken-where appropriate-to mitigate the impact of COVID-19. It corresponds to the Theil's Entropy and is computed by weighing the GDP with the population in terms of proportions with respect to the total. The conditional projections of each variable drawn in the SPBVAR-MTV in (1) have been used to perform forecasting from 2019q1 to 2021q1.

Estimating Procedure through Monte Carlo Simulations
In this analysis, I address a simulated example to highlight the performance of the estimating procedure of SPBVAR-MTV model in (1) using Monte Carlo simulations. It is performed accounting for two of the three described competing models in order to compare the results to related works: M I ('General Case'), with no structural breaks and volatility changes (hereafter 'unobserved effects'); and M I I I ('Full Case'), with unobserved effects in either time-varying parameters or log-volatilities. The idea is to highlight the thrust of the estimation method of performing better conditional forecasts when studying macroeconomic-financial linkages among countries and sectors with a high dimensional estimation sample covering-for example-triggering events (unobserved effects). Then, I compare the estimation method to the following related works: (i) multicountry Bayesian VAR (BVAR) as in Canova and Ciccarelli (2009); (ii) multicountry Panel Bayesian VAR (PBVAR) as in Ciccarelli et al. (2018); (iii) Structural Panel Bayesian VAR (SPBVAR) as in Pacifico (2019b); and (iv) Large Bayesian VAR with Stochastic Volatility (LBVAR-SV) as in Carriero et al. (2019) 34 .
According to model features in Section 2.1, I use a simplified version of (1)-supposing one lag and no intercept-and construct two distinct sets: a training set by simulating 20 variables as independent standard normal vectors (M s = 2, Q s = 6,Q s = 6, and Ξ s = 6), with s stands for 'simulated', for N = 15 and T = 100 time-series data; and a prediction set by generating 100 additional observations in the same manner. Then, I suppose two time-varying (theoretical) coefficient factors (β s f t ), with f = 1, 2, so defined. (i) β s 1t = (β 1t,1 , β 1t,2 , β 1t,3 , β 1t,4 , β 1t,5 , β 1t,6 ) s denoting variable-specific effects for m s = 1 and m s = 2 stacked in six ( (1) and SNLR in (6) assume the form: Without restrictions, the simulated sample amounts to 30,000 (theoretical) regression parameters (close to the estimation sample in the empirical analysis). More precisely, each equation of the supposed SPBVAR-MTV in (66) has k s = [15 × 20] × 1 = 300 coefficients, and there are 100 equations in the system 35 .
Let the supposed factorization in (67) be exact and the hyperparameters in δ be all known just as in the empirical application, (theoretical) posterior distributions are computed according to Equation (40) for β s f t , (46)-(49) for moment distributions in the state-space factorization structure, and (57)-(59) for ε s t . To draw projections for β s f t as close as possible to the empirical ones, I simulate random coefficient vectors by using-as initial value of the random number seed-the standard deviations (on average) of the (observed) variables 36 in (64). I use 1000 until 5000 draws and find that convergence is obtained at about 3000 draws 37 . The total number of draws has been 2000 + 3000 = 5000, which corresponds to the sum of the final number of draws to discard and save, respectively. A total of 5000 retained replications have been used to conduct posterior inference at each t. Table 3 displays the (theoretical) posterior estimates obtained simulating all five supposed Bayesian approaches: BVAR, PBVAR, and SPBVAR for M I ; and LBVAR-SV and SPBVAR-MTV for M I I I . The best model solution accounts for the only (theoretical) factors highlighted in bold and corresponds to the highest log Bayes Factor (lBF > 10) according to the generalized version of the Kass and Raftery (1995)'s scale of evidence in (53). Here, some considerations are in order. (i) According to absence of unobserved effects (M I ), SPBVAR model performs better (higher significance) because of its structural framework.
(ii) Posterior estimates in LBVAR-SV are close to the ones in SPBVAR-MTV, and would perform better (higher PMPs although lower significance) when dealing with the only (theoretical) macroeconomic-financial linkages (β s 1t,1 andβ s 1t,2 ). It would be because of its more general Bayesian procedure in which the estimation algorithm can be easily inserted into any pre-existing algorithm for sampling the posterior of the VAR's coefficients and then constructing coincident and leading indicators. (iii) SPBVAR-MTV would perform better when dealing with common and additional factors, and policy shifts because of its hierarchical (structural) prior specification strategy, multivariate time-varying volatility, and multidimensional data (panel data analysis). (iv) When jointly dealing with endogenity issues and unobserved effects, volatility changes need to be modeled through appropriate estimation methods (e.g., MCMC integrations). (v) When investigating unobserved effects, modeling time-varying log-volatilities as alternative to the stochastic volatility specification, posterior estimates seem to display either higher significance or consistency fitting the data (larger PMPs).
In Figure 9, I draw the corresponding (theoretical) density forecast combinations of the supposed SNLR model in (67)  s 2t,3 ). They correspond to the (theoretical) projections of every supposed variables drawn in the simulated sample. The yellow and red curves denote the 95% confidence bands, and the blue and purple curves denote the (theoretical) conditional 38 and unconditional 39 projection of outcomes Y s t for each N in the simulated sample, respectively. Here, the training sample is used to run hierarchical conjugate priors and then obtain posterior distributions, and the test sample to run forecasting accuracy.  (53). The significance codes stand for: (*) significance at 10%; (**) significance at 5%; and (***) significance at 1%.
The results emphasize the findings obtained in Table 3. (i) (Theoretical) conditional projections lie in the confidence interval; conversely, (theoretical) unconditional projections tend to diverge over T. Thus, when investigating macroeconomic-financial linkages, cross-unit lagged interdependencies and structural time variations-along with dynamic feedback and interactions-have to be accounted for. (ii) Density forecasts in (theoretical) total economy (third plot) tend to show a similar evolution to the ones according to (theoretical) financial economy (second plot) because of highly stringent cross-country interdependencies. (iii) Density forecasts in (theoretical) real economy tend to diverge (larger projections) because of stronger cross-country heterogeneity (first plot). (iv) Overall, the persistent divergence (different dynamics) observed in (theoretical) real and financial dimensions and even more in (theoretical) total economy-along with wide confidence intervals-highlights the relevant impact of consistent endogeneity and volatility issues. Thus, the need for forecasters and policymakers to account for additional factors and policy shifts with time-varying volatilities when formulating policies and forecasting in a large set of multidimensional data. s 1t,6 ('theoretical' financial markets); and χ s 2t,3β s 2t,3 ('theoretical' total economy). They correspond to (theoretical) conditional (blue line) and unconditional (purple line) projections of each supposed variable drawn in (67).

Concluding Remarks
This paper provides new empirical insights in order to give a relevant contribution to the more recent literature on international macroeconomic-financial linkages when jointly modeling and quantifying multicountry data using the information contained in a large set of endogenous and economic-financial variables. A multicountry SPBVAR with Multivariate Time-varying Volatility is developed to jointly deal with issues of endogeneity, because of omitted factors and unobserved heterogeneity, and volatility, because of policy regime shifts and structural changes. The two main differences with respect to a standard SPBVAR lie in an additional component to investigate fiscal and monetary policy implications and interactions, and in the variance-covariance matrix allowed to be time-variant. The latter is an useful way of modeling time-varying conditional second moments to provide an alternative to the stochastic volatility specification. The computational costs involved in using that specification are moderate since the high dimensionality is avoided via Bayesian inference and Monte Carlo Markov Chain (MCMC) implementations.
An empirical application is developed by accounting for the Central, Eastern, and Western European countries, with particular emphasis to the Great recession and successive post-crisis periods. The United States are included in the analysis to assess international spillover effects and possible contagion measures among financial markets. In this study, I focus on the latest two alternative monetary policy regimes that have been in place since the 1990: (1) the Inflation Stabilization Era from 1994 to 2008 and (2) the Zero Interest Rate Era from 2008 to 2015. Two more additional periods are also considered: (1) 2006q1 − 2009q4 to investigate possible commonality between financial markets and real economy during the Great recession and (2) 2010q1 − 2018q4 to evaluate fiscal implications and policy perspectives during post-crisis consolidation.
A simulated experiment -compared to related existing approaches -is also addressed to highlight the performance of the estimating procedure developed in this study using some Monte Carlo simulations. The findings prove that the hierarchical structural framework with time-varying log-volatilities perform better conditional forecasts when studying macroeconomic-financial linkages with structural breaks and volatility changes.
From a global perspective, the same dynamic behaviour is observed in the transmission of US financial shocks, with outward spillover effects. The findings are consistent and robust with the more recent literature on multicountry dynamic panel setups. More precisely, they confirm that US seem to be an important driver in allowing unexpected shocks to spill over and thus affecting European financial markets, mainly concerning CEE economies with inward spillovers. Then, intra-country shocks directly affect a country's own output growth in the real economy because of consistent cross-country interdependencies.
From a modeling perspective, the presence of highly strong intra-CEWE heterogeneity, in terms of intensity and Figure 9. The plot draws (theoretical) density forecast combinations for outcomes Y s t given N for the indicators that account for additional factors and policy shifts with unobserved effects. They are: χ s 1t,5β s 1t,5 ('theoretical' real economy); χ s 1t,6β s 1t,6 ('theoretical' financial markets); and χ s 2t,3β s 2t,3 ('theoretical' total economy). They correspond to (theoretical) conditional (blue line) and unconditional (purple line) projections of each supposed variable drawn in (67).

Concluding Remarks
This paper provides new empirical insights in order to give a relevant contribution to the more recent literature on international macroeconomic-financial linkages when jointly modeling and quantifying multicountry data using the information contained in a large set of endogenous and economic-financial variables. A multicountry SPBVAR with Multivariate Time-varying Volatility is developed to jointly deal with issues of endogeneity, because of omitted factors and unobserved heterogeneity, and volatility, because of policy regime shifts and structural changes. The two main differences with respect to a standard SPBVAR lie in an additional component to investigate fiscal and monetary policy implications and interactions, and in the variance-covariance matrix allowed to be time-variant. The latter is an useful way of modeling time-varying conditional second moments to provide an alternative to the stochastic volatility specification. The computational costs involved in using that specification are moderate since the high dimensionality is avoided via Bayesian inference and Monte Carlo Markov Chain (MCMC) implementations.
An empirical application is developed by accounting for the Central, Eastern, and Western European countries, with particular emphasis to the Great Recession and successive post-crisis periods. The United States are included in the analysis to assess international spillover effects and possible contagion measures among financial markets. In this study, I focus on the latest two alternative monetary policy regimes that have been in place since the 1990: (1) the Inflation Stabilization Era from 1994 to 2008 and (2) the Zero Interest Rate Era from 2008 to 2015. Two more additional periods are also considered: (1) 2006q1-2009q4 to investigate possible commonality between financial markets and real economy during the Great Recession and (2) 2010q1-2018q4 to evaluate fiscal implications and policy perspectives during post-crisis consolidation.
A simulated experiment-compared to related existing approaches-is also addressed to highlight the performance of the estimating procedure developed in this study using some Monte Carlo simulations. The findings prove that the hierarchical structural framework with time-varying log-volatilities perform better conditional forecasts when studying macroeconomic-financial linkages with structural breaks and volatility changes.
From a global perspective, the same dynamic behaviour is observed in the transmission of US financial shocks, with outward spillover effects. The findings are consistent and robust with the more recent literature on multicountry dynamic panel setups. More precisely, they confirm that US seem to be an important driver in allowing unexpected shocks to spill over and thus affecting European financial markets, mainly concerning CEE economies with inward spillovers. Then, intra-country shocks directly affect a country's own output growth in the real economy because of consistent cross-country interdependencies.
From a modeling perspective, the presence of highly strong intra-CEWE heterogeneity, in terms of intensity and spreading of spillover effects, emphasize more consistent difference among financial markets due to tighter monetary policies. In the aftermath of the Great Recession and an ongoing post-crisis consolidation, despite a consistent homogeneity holds among CEWE economies, different countries' responses tend to matter due to coordinated but not fairly flexible fiscal actions, mainly among emerging economies suffering from lower competitiveness. The findings confirm the need of examining international spillovers accounting for both misspecification problems and implied volatility changes.
From a policy perspective, the empirical results face a situation of trade-off. More precisely, if on one hand the adoption of sounder macroeconomic policies and economicinstitutional changes-put in place to foster consolidated policy actions-have helped to bring inflation in emerging economies back under control, on the other hand, in case of a noteworthy unexpected financial shock-without appropriate coordinated structural reforms in trade, product, and labour markets-outward government benefits will be not able for supporting the process of international financial integration among countries and boosting the output to potential growth.
Funding: The APC was funded by the Knowledge Unlatched initiative.

Data Availability Statement:
The datasets analysed during the current study are available from the corresponding author on reasonable request.