A note on the stability of the Swedish Phillips curve

We use Bayesian techniques to estimate bivariate VAR models for Swedish unemployment rate and inflation. Employing quarterly data from 1995Q1 to 2018Q3 and new tools for model selection, we compare models with time-varying parameters and/or stochastic volatility to specifications with constant parameters and/or covariance matrix. The evidence in favour of a stable dynamic relationship between the unemployment rate and inflation is mixed. Model selection based on marginal likelihood calculations indicates that the relation is time varying, whereas the use of the deviance information criterion suggests that it is constant over time; we do, however, note consistent evidence in favour of stochastic volatility. An out-of-sample forecast exercise is also conducted, but similarly provides mixed evidence regarding which model to favour. Importantly though, even if time-varying parameters are allowed for, our results do not suggest that the Phillips curve has been flatter in more recent years. This finding thereby questions the explanation that a flatter Phillips curve is the cause of the low inflation that Sweden has experienced in recent year.


Introduction
Inflation in Sweden has been stubbornly low over the last years. This development is similar to that in several other inflation-targeting countries where inflation has been moderate, and increasing slowly, despite historically low policy-interest rates and developments in the real economy which many argue should have generated a stronger inflationary pressure; see, for example, Jansson (2017) and Yellen (2017). One explanation for the low inflation which has been put forward is that the Phillips curve has become flatter, for example, due to globalisation and digitalisation. 1 Other explanations have, however, also been suggested. For instance, the oil price, aggregate demand, monetary policy and demographic trends have-just to mention a few variables-been put forward as having some importance; see, for example, Conti et al. (2017), Juselius and Takáts (2018) and Williamson (2018).
In this paper we contribute to the discussion regarding the properties of the Phillips curve by providing evidence based on Swedish data. Employing a Bayesian VAR (BVAR) framework, we estimate bivariate models using quarterly data on unemployment rate and inflation. The models are estimated under different assumptions concerning the dynamics and covariance matrix. As noted by a growing literature, time variation in both dynamics and volatility appears to be important features of macroeconomic relationships; see, for example, Cogley and Sargent (2005), Koop et al. (2009), Franta et al. (2014, Chan et al. (2016), Knotek and Zaman (2017) and Akram and Mumtaz (2019). We therefore estimate models with time-varying parameters and/or stochastic volatility and compare these models to specifications with constant parameters and/or covariance matrix. Four combinations are considered: (i) constant parameters and covariance matrix, (ii) time-varying parameters and constant covariance matrix, (iii) constant parameters and stochastic volatility and (iv) timevarying parameters and stochastic volatility. Relying on new tools for model selection developed by Chan and Eisenstat (2018), we formally assess which model is preferred by the data. This constitutes a step forward relative to the vast majority of previous related research. Since model selection is a non-trivial issue when models with timevarying parameters and stochastic volatility are involved, it has often simply been assumed that it is reasonable to employ a model with such features. 2 In this paper we instead evaluate this assumption using marginal likelihoods and the deviance information criterion (DIC) and can provide statistical evidence on the stability of the Swedish Phillips curve.
The focus in our analysis is put on bivariate models estimated using price inflation and the unemployment rate. While the original analysis conducted by Phillips (1958) was based on wage inflation, the more recent debate regarding the flattening of the Phillips curve has typically been focused on price inflation and we believe that this focus is reasonable given the many central banks-including Sveriges Riksbank-that target price inflation. However, we do not wish to go too far from Phillips' original idea and the intuitive bivariate setting; after all, we believe that there is value in assessing a reasonably clean and simple relation between inflation and the unemployment rate. That said, we also conduct analysis based on wage inflation and use different measures of resource utilization in order to assess the sensitivity of our findings from the benchmark specification (with price inflation and the unemployment rate). As part of this sensitivity analysis, the model is also augmented beyond the bivariate case with some key variables. Finally, we conduct an out-of-sample forecast exercise in order to further evaluate the relative merits of the different assumptions regarding the dynamics and covariance matrix.
Our main results suggest that the evidence in favour of a stable relation between the unemployment rate and inflation is mixed. Model selection based on marginal likelihoods indicates that the relation is time varying; if we instead rely on the DIC, the relation is judged to be constant over time. Importantly though, we do not-regardless of model specification-find any evidence suggesting that the Phillips curve has been flatter than usual during the last years of our sample and accordingly conclude that the low inflation in the last few years has not been caused by a changing slope of the Phillips curve. Our findings instead indicate that low trend inflation may have contributed to this development. Regarding the results from the sensitivity analysis, we argue that our main conclusion-that is, that there is no evidence of the Swedish Phillips curve being flatter during the last few years of the sample-is very robust. The results from the out-of-sample forecast exercise are similar to those from the withinsample analysis-that is, mixed. No model shows superior forecasting performance at all investigated horizons.
The rest of this paper is organised as follows: In Sect. 2, we describe the Bayesian VAR models which we employ. The results from our main empirical analysis relying on price inflation and the unemployment rate are presented and discussed in Sect. 3. In Sect. 4, we show results from various alternative specifications in order to assess how sensitive our main results are. An out-of-sample forecast exercise is conducted in Sect. 5. We provide an overall discussion of our results in Sect. 6. Finally, Sect. 7 concludes.

The Bayesian VAR models
We rely on BVARs for our analysis since the inflation equation of the BVAR can be seen as a "dynamic generalization of the Phillips curve" (King and Watson 1994, p. 172). While we pay special attention to the inflation equation, our analysis is mainly based on a bivariate system since important aspects of the dynamic relation between the variables otherwise could be lost.
Defining the vector of dependent variables as y t u t π t , where u t is the unemployment rate and π t is inflation, we specify the BVAR in its most general form-that is, with time-varying parameters and stochastic volatility ("TVP-SV")-in Eq. (1): B 0t y t γ t + B 1t y t−1 + . . . + B pt y t− p + ε t (1) where B 0t is a 2 × 2 lower triangular matrix with ones on the diagonal, γ t contains the time-varying intercepts and the matrices B 1t , . . . , B pt describe the dynamics. The vector of disturbances, ε t , follows ε t ∼ N (0, Σ t ), where Σ t diag(exp(h 1t ), exp(h 2t )); the fact that the covariance matrix is diagonal and that the matrix B 0t premultiplies y t means that the model is specified in a "structural" form. 3 The lag length is set to p 4. Collecting the free parameters of γ t and B it in the 19 × 1 parameter vector θ t , we specify the processes for the time-varying parameters and log-volatilities as random walks: Having presented the most general version of the VAR model above, restrictions are then imposed when estimating the three alternative specifications, that is, (i) constant parameters and covariance matrix ("constant"), (ii) time-varying parameters and constant covariance matrix ("TVP") and (iii) constant parameters and stochastic volatility ("SV"). When estimating a model with constant parameters, Σ θ is restricted to be zero. Similarly, Σ h is restricted to be zero when a model with constant covariance matrix is estimated.
When it comes to the issue of model selection-that is, assessing which model is preferred by the data-we rely on marginal likelihoods as well as the DIC. In a Bayesian setting, the marginal likelihood is the appropriate measure of how well the model and prior agree with the data; the model with the highest marginal likelihood is the one preferred by the data. While the marginal likelihood is the prior expected value of the likelihood, the DIC (Spiegelhalter et al. 2002) combines a measure of fit, the posterior expectation of the log-likelihood, with a penalty term for the number of parameters and can thus be argued to be less sensitive to the choice of prior distributions.
As we rely on both marginal likelihoods and the DIC for model choice, we take care to ensure that the prior information is comparable across models. For the regression parameters, θ , we use an uninformative prior for the overall location together with a relatively tight prior on the evolution over time. Specifically, for models with constant parameters, the prior is θ ∼ N (0, V θ ) with V θ 10 I. For models with time-varying parameters, the same normal prior is used for the initial state, θ 0 , and the evolution is governed by the variance of the increments. The diagonal elements of Σ θ are given inverse Gamma priors, σ 2 θi ∼ i G(v θi , S θi ), with v θi 5 and prior means of 0.01 for the intercepts and 0.0001 for the other regression parameters. In expectation this yields a prior variance for the final state close to 11 for intercepts and close to 10 for the other parameters.
The specification of the prior for the error variances or volatilities requires more care as conditional conjugacy suggests inverse Gamma priors for the diagonal elements of Σ in the constant variance models and a normal prior for the initial state, h 0 , or a log-normal prior for the initial variance in the stochastic volatility models. As the first step, we tailor the prior to match the scale and variation of the data. Let s 2 i be the residual variance from a univariate AR model. Employing a normal prior, h 0 ∼ N μ h , c I , we set μ h,i ln s 2 i − c 2 to match the prior mean of exp(h it ) with the residual variance. Similarly, for the constant variance case and inverse gamma priors, Finally, c and v 0i are selected so that the shapes of the prior distributions match while keeping them relatively uninformative, resulting in c 0.25 and v 0i 5. Finally, the prior for the diagonal elements of Σ h is inverse Gamma, σ 2 hi ∼ i G(v hi , S hi ), with v hi 5 and S hi 0.04 for a prior mean of the variance of 0.01.
While the models with time-varying parameters and/or stochastic volatility are inherently more flexible than the VAR with constant parameters and covariance matrix, the properties of the data implied by the priors are quite comparable across models. The prior means of the regression parameters are the same for the different models, and the increase of the prior variance with t is moderate for the models with time-varying parameters. While the prior expectation of the log-volatilities is constant over time with only a moderate increase in the variance, the log-normal form does imply that the volatility is increasing with t a priori. We have E(exp(h it )) exp μ h,i + tσ 2 hi /2 , and at the prior mean of σ 2 hi , this implies that the prior expectation of the volatility at the end of the sample is 60% higher than the beginning of the sample. This is, however, more an issue of model specification than prior specification where most of the tractable model specifications that result in constant prior expectation of the volatilities imply decreasing prior expectation of the log-volatilities. For simplicity we have opted for the current, quite standard, specification of the time-varying volatilities.
For posterior inference, we rely on the Markov Chain Monte Carlo sampler employed by Chan and Eisenstat (2018). The sampler proceeds by sampling from the full conditional posteriors The full conditional posterior of θ θ 1 , . . . , θ T is normal both in models with constant and time-varying parameters. In models with stochastic volatility, h is sam-pled using the auxiliary mixture sampler of Kim et al. (1998). The diagonal elements of Σ θ and Σ h are conditionally independent with inverse Gamma full conditional posteriors. Finally, θ 0 and h 0 are conditionally independent with normal full conditional posteriors. For the models with constant variance, the last three steps are replaced by an update of the diagonal elements of Σ which are conditionally independent with inverse Gamma full conditional posteriors.
The marginal likelihood and DIC calculations are carried out using the methods developed in Chan and Eisenstat (2018). Write the marginal likelihood for the TVP-SV model as where ξ collects the parameters θ 0 , h 0 , θ and h of the state Eqs. (2) and (3). The integral with respect to θ is analytic, and importance sampling is used to integrate out h. Having integrated out the latent time-varying parameters and volatilities, this yields what Chan and Eisenstat term the integrated likelihood The remaining integration with respect to the parameters in ξ is then carried out in an outer importance sampling loop to obtain m( y) p( y|ξ ) p(ξ )dξ . For the TVP model the integrated likelihood calculation skips the integration with respect to h and in the outer loop ξ collects the parameters θ 0 , θ and diag( ), while the SV model omits the integration with respect to θ in the integrated likelihood and ξ collects the parameters θ , h 0 and h . That is, the constant regression parameters are integrated out in the outer loop. For the constant parameter VAR, there are no latent parameters, the integrated likelihood coincides with the marginal likelihood, and the integration with respect to θ and diag( ) in ξ is analytic.
The DIC is based on the deviance which, for the purpose of model comparison, can be defined as D(ξ ) −2 ln p( y|ξ ). Model fit is measured by the posterior mean deviance That is, in contrast to the marginal likelihood, the integration is over the posterior distribution of ξ and not the prior. Model complexity is measured by the effective number of parameters, whereξ is an estimate of ξ , here the posterior mean. Combining these yields the DIC as The calculation of the integrated likelihood in (9) is done in the same way as for the marginal likelihood, and the integration with respect to ξ in (10) is carried out in an outer importance sampling loop.

Empirical findings: unemployment rate and price inflation
In this section, we use data on seasonally adjusted unemployment rate and CPIF 4 inflation ranging from 1995Q1 to 2018Q3. The unemployment rate refers to the age group 16-64 years. CPIF inflation is calculated as π p t 100(P t /P t−4 − 1), where P t is the CPIF index at time t. We choose 1995Q1 as the starting point since it formally constitutes the start of Sweden's inflation-targeting regime and this gives the model with constant parameters and covariance matrix a fair chance to be considered a relevant modelling choice. 5 Data are shown in Fig. 1.
We estimate the four different models presented in Sect. 2. Evaluation measures-in terms of log marginal likelihoods and DICs-are shown in Table 1. Judging by the log marginal likelihoods, the model with time-varying parameters and stochastic volatility is the preferred one. Based on this evidence, we would conclude that the Phillips curve has not been stable. However, Table 1 also shows that the DIC indicates that the Based on bivariate BVARs with y t u t π p t , where u t is the unemployment rate and π p t is price inflation. The table gives the natural logarithm of the marginal likelihood of the models. "DIC" is the deviance information criterion. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters best model is that with constant parameters and stochastic volatility; perhaps somewhat surprisingly, it also suggests that the model with time-varying parameters and stochastic volatility is the worst. The method used for model selection obviously matters, and we conclude that this exercise has left us with contradictory results regarding whether there is time variation in the Swedish Phillips curve. It can be noted though that regardless of which selection tool that is used, a model with stochastic volatility comes out on top.
If the Phillips curve is stable over time-as the DIC suggests-the dynamic aspects that we are primarily interested in can quite easily be summarised. Looking at the preferred model-that is, the model with constant parameters and stochastic volatility-the first aspect of interest is the impulse-response function which describes how inflation reacts to a shock in the unemployment rate. This is shown in Fig. 2. As can be seen, a one-standard-deviation shock to the unemployment rate lowers-in line with our expectations-inflation for a number of quarters. The largest effect is found at the two-quarter horizon. It should be noted that the minor differences in the impulseresponse function that can be found between different dates are only due to the fact that the size of the shock is time varying. The size of the shock-that is, its standard deviation-is illustrated in the impulse-response function which describes the effect that a shock to the unemployment rate has on the unemployment rate itself; see Fig. 11 in "Appendix". As can be seen, there is limited variation in this standard deviation-it ranges from 0.20 to 0.25. The evidence for stochastic volatility in the model appears to be driven by the volatility of the shocks to the inflation equation which is shown in Fig. 12. 6 The second aspect of the Phillips curve that is of primary interest is its slope. In line with other related research-see, for example, Knotek and Zaman (2017) and Karlsson and Österholm (2018)-we report the sum of the coefficients on the (1) is given by where δ t B −1 0t γ t , A it B −1 0t B it and e t B −1 0t ε t . It can be noted that for the model with constant parameters and stochastic volatility, there is no time variation in γ t , B 0t or B it ; the slope of the Phillips curve is accordingly constant during the investigated sample, and it is estimated to be −0.06. This is a flatter Phillips curve than what Karlsson and Österholm (2018) found for most periods for the USA (comparing point estimates) but very close to the slope that Knotek and Zaman (2017) found for the US Phillips curve in the two time points that they focused on, namely 1999Q3 and 2017Q3.
As shown in Table 1 though, the marginal likelihoods indicated that the model with time-varying parameters and stochastic volatility was the preferred one. We therefore next turn our attention to the properties of the Phillips curve based on this model. The impulse-response function describing the effect that a shock to the unemployment rate has on inflation is given in Fig. 3 7 ; the other impulse-response functions are shown in Figs. 14, 15 and 16 in "Appendix". Also in this model, an unexpectedly high unemployment rate tends to decrease inflation at short horizons (apart from the three and four quarter horizons, the latter only in the beginning of the sample). The impulse-response function is reasonably stable over time despite the fact that the  Fig. 17 in "Appendix" for the unemployment rate equation. 8 It can also be noted that the time variation in the equation for inflation is similar in magnitude to that in the equation for the unemployment rate. We thus conclude that the time variation suggested by the marginal likelihoods is not due to time variation in the equation for the unemployment rate alone.
Turning to the slope of the Phillips curve, this is plotted in Fig. 5. Taking the 68% credible interval into account, one might claim that the slope of the Phillips curve has not changed dramatically during the sample. However, looking at the point estimate, the story is somewhat different as it ranges from − 0.10 to − 0.35. Interestingly though, the slope has not been unusually close to zero between 2011 and 2016 indicating that Sweden's low inflation in this period cannot be explained by a flat Phillips curve.
The final aspect of the evolving properties of the model with time-varying parameters and stochastic volatility which we consider is a concept called "trend inflation". 8 Looking at Fig. 4, it is also worth noting that the sum of the coefficients on lagged inflation is close to zero. Other things equal, a sum closer to zero (rather than closer to unity) means that inflation's deviation from what we refer to as "trend inflation" below will be more short-lived. One explanation for this finding is the moderate persistence that Swedish CPIF inflation shows in general; the first-order autocorrelation is 0.77. This can be contrasted with US data where inflation is often modelled as having a unit root; see, for example, Stock and Watson (2009). Given that the model with time-varying parameters also can remove a persistent component from the series through the time-varying intercept, we do not consider our finding of a sum close to zero on the coefficients on lagged inflation to be very surprising. It is also in line with the findings of modest inflation gap persistence in the USA found by Cogley and Sbordone (2009) for the period 1984-2003-a period during which monetary policy generally is considered to have been focused on inflation stabilization. In contrast, Cogley and Sbordone show that between 1960 and 1983, inflation gap persistence was substantially higher. Further discussions concerning inflation gap persistence in the USA can be found in, for example, Benati and Surico (2008) and . In line with the terminology of Faust and Wright (2013) and Clark and Doh (2014), we define trend inflation as the value to which the inflation forecasts from the BVAR model converge. 9 In order to define the concept algebraically, we first rewrite the reduced form VAR in Eq. (13) in companion form as andẽ t e t , 0 , . . . , 0 Trend inflation is given as the second element of φ t (I − A t ) −1δ t . Estimated median trend inflation at each point in time is shown in Fig. 6 and is interesting to discuss from a monetary policy perspective.
Around 2010-2014, the Riksbank was "leaning against the wind" in order to dissuade households from taking on housing-related debt. This contributed to the low inflation outcomes during this period and sparked an intense debate regarding whether the Riksbank's policy was expansive enough; see, for example, Svensson (2014). As can be seen, trend inflation fell while the Riksbank was leaning against the wind. By 2014 trend inflation was between 0.6 and 1%. This coincides with when the Riksbank abandoned the policy of leaning against the wind and declared that its focus was on achieving the inflation target-partly because they were concerned about deanchoring of inflation expectations. 10 The model results hence give some support for the Riksbank's decision to change its policy in 2014 since it suggests that the low inflation outcomes had become problematic as its long-run inflation forecasts were not consistent with the inflation target. In order to fight low inflation and low inflation expectations, the Riksbank lowered the repo rate to zero in October 2014 and kept it at − 0.5% between February 2016 and December 2018. 11 After the change in the Riksbank's policy, inflation has slowly but surely increased and trend inflation with it; this development has occurred despite the weak development in oil prices that took place during part of this period (see Fig. 25 in "Appendix"), suggesting that a monetary policy focused on achieving the inflation target has been effective in moving both inflation outcomes and expectations towards the target.
Summarising the results so far, we note that there is mixed evidence on whether the Phillips curve is stable or not. However, even if we allow for the possibility that it has been unstable, it does not show signs of having been flatter than usual in the last few years of the analysed sample.

Empirical findings: sensitivity analysis
In order to assess the robustness of our findings in the previous section, we next conduct additional analysis that covers some relevant aspects given our empirical setting. First, we vary some of the prior settings of the models. Second, data on wage inflation are used instead of price inflation. Third, we replace the unemployment rate with alternative measures of resource utilization. Finally, we move beyond the bivariate framework and estimate models which include additional variables that could be of particular interest when studying Swedish inflation.

Priors
We believe the choice of priors described in Sect. 2-which were used when estimating the models underlying the results presented in Sect. 3-to be uncontroversial. The prior settings are, for example, close to the ones used by Chan and Eisenstat (2018). Still, given the stated purpose of facilitating model comparison, we have to some extent deviated from what can be seen as "common practice" with different model classes. While Minnesota style priors are common with constant parameter VARs, we do not use this type of prior as the increased prior shrinkage towards zero with increasing lag , where u t is the unemployment rate and π p t is price inflation. Table gives the natural logarithm of the marginal likelihood of the models. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters length to some extent is in conflict with the nature of time-varying parameter models. Similarly, while training sample priors are common in the literature using VARs with time-varying parameters, we do not use this type of prior since the monetary policy regime was quite different in the period before our sample and there are clear level breaks in both inflation and unemployment just prior to the start of our sample. Given this, a training sample prior would put constant parameter models at a disadvantage, which we want to avoid. It is nevertheless of interest to see how sensitive our results are with respect to the priors. More specifically, we vary the prior settings by halving and doubling key prior parameters one at a time. Specifically, we vary the prior variance of the regression parameters (or their initial state), V θ , and the scale parameters S θi and S hi governing the prior mean of the innovation variances in the state Eqs. (2) and (3).
In Tables 2 and 3, we present log marginal likelihoods and DICs for the four specifications under consideration using six different variations in the priors; the benchmark specification (V θ 10 I, S θi 0.04 for constants and 0.0004 for the other regression parameters and S hi 0.04) used in Sect. 3 is given in the result column furthest to the left. As can be seen from the two tables, the results with respect to model selection are extremely robust: Based on marginal likelihoods, the model with time-varying parameters and stochastic volatility is always the preferred one. Similarly, when the DIC is employed, the model with constant parameters and stochastic volatility is always judged as being superior to the others.
Based on the results above, we conclude that our results regarding model selection are not particularly sensitive with respect to our choice of prior parameters.

Wage inflation
We next turn our attention to the properties of the Phillips curve when it is estimated using data on wage inflation instead of price inflation. This has the benefit of taking our analysis in the direction of Phillips' (1958) original observation which-as stated above-was based on wage inflation. But perhaps more importantly, it also means that we can scrutinise what may be described as a key relationship through which central , where u t is the unemployment rate and π p t is price inflation. Table  gives the DIC-that is, the deviance information criterion-of the models. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters banks assume that their goals are achieved: Expansive monetary policy is assumed to stimulate the real economy, thereby lowering the unemployment rate; this tighter labour market pushes up wage growth (wage inflation) which in turn increases price inflation. Assessing the stability of the relation between the unemployment rate and wage inflation is thus of considerable interest.
Price inflation in the bivariate BVARs is accordingly replaced with wage inflation; this is calculated as π w t 100(W t /W t−4 − 1), where W t is the hourly wage at time t. 12 Our vector of dependent variables thus becomes y t u t π w t . As before, our sample ranges from 1995Q1 to 2018Q3. Data are shown in Fig. 7. Based on bivariate BVARs with y t u t π w t , where u t is the unemployment rate and π w t is wage inflation. Table gives the natural logarithm of the marginal likelihood of the models. "DIC" is the deviance information criterion. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters As was the case when the models were estimated using price-inflation data, we find that the model with time-varying parameters and stochastic volatility is the preferred one when looking at marginal likelihoods. However, the DICs once again conclude that the model with constant parameters and stochastic volatility is most suitable ( Table 4).
The model with constant parameters and stochastic volatility obviously implies that the Phillips curve has been stable during the studied period and it is therefore not particularly interesting to study it in detail. In Fig. 18 in "Appendix", we nevertheless show the impulse-response function describing the effect that a shock to the unemployment rate has on wage inflation. 13 It can be noted though that while the effect in the short run is the expected-that is, wage inflation decreases-the slope of the Phillips curve for this model is 0.00. So while the model suggests that the Phillips curve has been stable, it can be questioned whether there actually exists a Phillips curve in this case.
If we instead turn to the preferred model according to marginal likelihoods, Fig. 8 shows the effect that a shock to the unemployment rate has on wage inflation. 14 In this model with time-varying parameters and stochastic volatility, a shock to the unemployment rate tends to decrease wage inflation; the maximum effect of approximately −0.06 is reached after five quarters. And similarly to the case when using price inflation, the impulse-response function is reasonably stable over time.
In Fig. 9, the slope of the Phillips curve is presented. As can be seen, this shows a tendency to have become steeper over time, though one should of course recall that there is a fair bit of uncertainty associated with the estimates as indicated by the 68% credible interval. Importantly though, there is no evidence of the Phillips curve having flattened during the last few years of the sample.
Finally, in Fig. 10 we show the model's estimate of trend inflation. This indicates that there have been noticeable differences in trend inflation over time. While it hovered between 4 and 5% from 1996 to 2002, it has varied around much more modest levels 13 The shock size is one standard deviation, which varies between 0.19 and 0.25. 14 The shock size is one standard deviation, which varies between 0.10 and 0.16. where u t is the unemployment rate and π w t is wage inflation. Size of impulse is one standard deviation. Effect in percentage points on vertical axis. Horizon in quarters and dates on horizontal axes Fig. 9 Estimated "slope" of the Phillips curve. Wage inflation and unemployment rate data. Model with time-varying parameters and stochastic volatility. Note: Based on bivariate BVAR with y t u t π w t , where u t is the unemployment rate and π w t is wage inflation. Coloured band is 68% equal-tail credible interval after the global financial crisis; in fact, trend inflation has not been more than 3% since 2009. Given that actual wage inflation has been low since 2010, this is not particularly surprising though. , where u t is the unemployment rate and π w t is wage inflation. Coloured band is 68% equal-tail credible interval Overall, the results based on wage inflation are largely in line with what we found using price inflation: The same two types of models were preferred by the data, and the evidence regarding the stability of the Phillips curve is mixed. Also, there is no support to a claim of a Phillips curve that has been flatter than usual in the last few years. However, when using wage inflation data, the model with constant parameters and stochastic volatility finds little support for a Phillips curve as the slope is estimated to be zero. 15 Given Phillips' (1958) original observation, the unemployment rate is the most natural choice of variable to relate to inflation when conducting the kind of analysis that we do in this paper. However, during the process of becoming popular, the meaning of the term "Phillips curve" has become somewhat blurred. It now comes in many different forms and using a range of different variables. 16 This fact is reflected in a description of the Phillips curve provided by McCracken (2006, p. 1127), who mean that it broadly can be defined as "… a model relating inflation to the unemployment rate, output gap, or capacity utilization". 15 On the other hand, when using price inflation data in Sect. 3, the 68% credible interval for the slope was all in the negative range, from −0.09 to −0.02. 16 One of the most important versions in modern macroeconomics is the so-called New Keynesian Phillips curve; see, for example, Galí and Gertler (1999) and Rudd and Whelan (2005). Recent empirical literature taking a less structural appoach and relying on various techniques from time series economectrics includes Svensson (2015), Busetti and Caivano (2016), Chan et al. (2016) and Knotek and Zaman (2017).

Other measures of resource utilization
Seeing that different measures of capacity utilization are considered in the empirical literature related to the Phillips curve, we next investigate the sensitivity of our results with respect to this aspect. 17 We do this by estimating bivariate models with price inflation but using either the unemployment gap or the output gap instead of the unemployment rate. That is, we define the vector of dependent variables in Eq. (1) as where u gap t is the unemployment gap and g gap t is the output gap.
The unemployment gap is defined as the actual unemployment rate minus the equilibrium (or "natural") unemployment rate. If the equilibrium rate is constant, then it is irrelevant whether the unemployment gap or the unemployment rate is used. However, there are both theoretical and empirical arguments suggesting that the equilibrium rate is time varying; see, for example, Hall (2003). We account for this fact by using the unemployment gap as calculated by the National Institute of Economic Research (NIER). 18 This gap is based on the NIER's best estimate of a time-varying equilibrium unemployment rate in Sweden. 19 Data on both the unemployment gap and the equilibrium unemployment rate are given in Fig. 19 in "Appendix". 20 As can be seen, the equilibrium rate is assumed to move very moderately over time.
Turning to the output gap, this is also based on data from the NIER. It is calculated as the percentage deviation between actual GDP and the NIER's best estimate of potential GDP. Data on the output gap can also be found in Fig. 19 in "Appendix".
Results from the model selection exercise are given in Table 5. As can be seen, we find exactly the same results as before, regardless of capacity utilization measure: According to marginal likelihoods, the model with time-varying parameters and stochastic volatility is preferred whereas the DIC chooses the model with constant parameters and stochastic volatility.
Focusing on the slope of the Phillips curve in the models with constant parameters and stochastic volatility, we can note that the model with the unemployment gap has a slope of −0.07, which is very close to what we found in our main analysis in Sect. 3 (that is, −0.06). The model estimated with output gap data has a positive slope as expected, namely 0.03. Turning to the models with time-varying parameters and stochastic volatility, the slope of the Phillips curve when the model was estimated using the unemployment gap is shown in Fig. 20 in "Appendix". This suggests that 17 Examples include broad unemployment, underemployment, the unemployment gap, the output gap, a labour market conditions indicator gap and labour underutilization; see Riggi and Venditti (2015), Bulligan et al. (2017), Bell and Blanchflower (2018) and Conti and Gigante (2018) for some recent examples. To some extent the use of alternative measures has been driven by a failure to find support for a Phillips curve when using specifications based on a "traditional" unemployment rate, unlike what we do in this paper. It is unclear why the unemployment rate appears to carry informational value for inflation in Sweden but not in several other countries. Given the complexity of the issue though, we refrain from suggesting explanations rather than providing what would have to be speculations. 18 The NIER is a Swedish government authority under the Ministry of Finance which does analysis and forecasting. https://www.konj.se/. 19 See, for example, National Institute of Economic Research (2017) for a discussion. is the output gap and π p t is price inflation. Table gives the natural logarithm of the marginal likelihood of the models. "DIC" is the deviance information criterion. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters our previous findings are robust; there are no signs that the Phillips curve has been flatter than usual the last few years. The same conclusion is drawn also when using output gap data-see Fig. 21 in "Appendix". It can be noted though that the Phillips curve has the wrong slope in the beginning of the sample. However, given the dramatic developments in the Swedish economy after the financial crisis of the early 1990s, it is perhaps not completely surprising that some counterintuitive results emerge at that point in time.
Summing up the results, we note that using the unemployment gap or output gap as measures of resource utilization yields results that are quite similar to when the unemployment rate is used. We find no evidence of a flattening Phillips curve during the last years of the sample.

Additional factors
The models that we have relied upon so far have exclusively been bivariate. While we believe that this has many benefits, it is obviously the case that other factors than the unemployment rate might affect Swedish inflation. We therefore want to investigate the sensitivity of our findings to the use of somewhat larger models, keeping in mind that the purpose of this paper clearly is not to find the best analysis or forecasting model for Swedish inflation. In this subsection, we consider two aspects that appear particularly relevant for inflation, namely foreign variables and Swedish interest rates.
We first augment the model with the repo rate, i r t , which is the policy rate of the Riksbank. Data are shown in Fig. 22 in "Appendix". This makes the vector of dependent variables y t u t π p t i r t . Trivariate models based on the unemployment rate, inflation and an interest rate are commonly employed tools in the related empirical literature using Bayesian VARs-see, for example, Cogley and Sargent (2005) and Primiceri (2005)-and accordingly seem like a natural extension in this case.
Results from the model selection exercise are shown in Table 6, and as can be seen, the model with time-varying parameters and stochastic volatility is chosen when we Based on trivariate BVARs with y t u t π p t i r t and y t u t π p t i s t , where u t is the unemployment rate, π p t is price inflation, i r t is the repo rate and i s t is the shadow rate of De Rezende and Ristiniemi (2018) based on a three-factor model. Table gives the natural logarithm of the marginal likelihood of the models. "DIC" is the deviance information criterion. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters rely on marginal likelihoods; the DIC, however, once again selects the model with constant parameters and stochastic volatility. The slope of the Phillips curve-where we still measure this as the sum of the coefficients on lagged unemployment in the inflation equation-for the latter model is −0.11. For the model with time-varying parameters and stochastic volatility, the slope is shown in Fig. 23 in "Appendix". Comparing it to the slope from the benchmark model (in Fig. 5), it can be seen that adding the repo rate to the system does not substantially affect the estimate.
However, in the sample employed here, traditional monetary policy has been restricted by the "zero lower bound" which lead the Riksbank to also conduct unconventional monetary policy, for example, by buying government bonds in large quantities. 21 In order to take both conventional and unconventional monetary policy into account, we also estimate a model with a shadow rate (instead of the repo rate) 22 ; we accordingly define y t u t π p t i s t , where i s t is the shadow rate. 23 Results from the estimation are shown in Table 6.
As can be seen, we find results that echo our previous findings: Marginal likelihoods favour the model with time-varying parameters and stochastic volatility, and the DIC selects the model with constant parameters and stochastic volatility. The slope of the model with time-varying parameters and stochastic volatility is shown in Fig. 24 in "Appendix" and shows a striking resemblance to the results based on the model using the repo rate (in Fig. 23). The model with constant parameters and stochastic volatility has an estimated slope of −0.13, that is, only a minor difference relative to what was found above when employing the repo rate. 21 It should be noted that the concept "zero lower bound" is something of a misnomer. In Sweden, the repo rate was lowered to − 0.5% and other central banks have also had negative policy rates. 22 We use the shadow rate of De Rezende and Ristiniemi (2018) based on a three-factor model. For other recent examples of analysis using shadow rates, see Wu and Xia (2016) and Conti (2017). 23 Note that the data on the shadow rate range from 1996Q1 to 2018Q3. The model using these data is accordingly estimated on a marginally shorter sample than previous models. and y t π oil t u t π p t , where g gapUS t is the output gap in the USA, π oil t the year-on-year percentage change in the oil price, u t is the unemployment rate and π p t is price inflation. Table gives the natural logarithm of the marginal likelihood of the models. "DIC" is the deviance information criterion. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters We next turn to the foreign variables. Through the inclusion of these, we aim to-at least to some extent-take into account the fact that Sweden is a small open economy and that it has been established that global factors have been shown to be empirically relevant determinants of inflation in several economies; see, for example, Conti et al. (2015) and Bobeica and Jarocinski (2017). We estimate two trivariate models-one where the US output gap is added to our benchmark bivariate system and the other where we instead add the percentage change (year-on-year) in the oil price. Data are shown in Fig. 25 in "Appendix". 24 Results regarding model selection are shown in Table 7. Turning to the model including the US output gap first-that is, where we have set y t g gapUS t u t π p t -it can be seen that the model with time-varying parameters and stochastic volatility is once again chosen when looking at marginal likelihoods. 25 The DIC on the other hand suggests that the model with constant parameters and covariance matrix is the preferred one. The difference to the "usual winner" when using the DIC-that is, the model with constant parameters and stochastic volatility-is extremely small though. Looking at the slope of the model with time-varying parameters and stochastic volatility in Fig. 26 in "Appendix", we can note that it is very similar to what we found in our benchmark specification (shown in Fig. 5). The model with constant parameters and covariance matrix has a very moderate slope −0.02, which can be compared to −0.06 which we found for the model with constant parameters and stochastic volatility in Sect. 3. 24 The US output gap data range from 1995Q1 to 2018Q1. Similar to the the model using the shadow rate, it is accordingly estimated on a marginally shorter sample than the other models. 25 We have chosen to place the US output gap first in the system. This reflects the fact that Sweden is a small open economy and accordingly should not affect a US variable. Obviously, the choice of ordering in the system does not alone make the US output gap exogenous with respect to the Swedish variables; Swedish variables can here affect the US output gap through lags. This could be changed, for example, through the use of additional hyperparameters which can forcefully shrink the coefficients of the Swedish variables in the equation for the US output gap to zero; see, for example, Villani and Warne (2003).
Finally, we note that for the model including oil price changes-that is, where y t π oil t u t π p t -model selection based on marginal likelihoods suggests that the model with constant parameters and stochastic volatility is the most appropriate. A model with constant parameters and covariance matrix is supported by the DIC. The estimated slope for these two models is −0.04 and −0.03, respectively. The model selection hence points in a somewhat different direction than what we have seen for other specifications. Overall, it does not, however, change our main conclusion from the analysis conducted up until here: Even though the evidence in favour of a stable Phillips curve in Sweden is mixed, there is no evidence that the Swedish Phillips curve has become flatter in the last few years of the sample.

Out-of-sample forecasts
Above we have shown-perhaps somewhat ironically-that the mixed evidence in favour of a stable Phillips curve is a robust finding over many different specifications. The conclusion ultimately depends on which model selection tool one prefers. As a final empirical exercise, we will try to shed some light on the model selection issue by conducting an out-of-sample forecast exercise. There are admittedly pros and cons with this approach. In practice, we aim to establish whether there is Granger causality and Ashley et al. (1980) mean that out-of-sample forecasting is the natural way to do this. Diebold (2015), on the other hand, points out that any out-of-sample exercise makes inefficient use of the available information; while not completely dismissing out-of-sample forecasting performance as a tool for model selection, he concludes (p. 8) that "In general, full-sample model comparison procedures appear preferable".
While we believe that Diebold has a valid point, we also believe that there might be additional information to be gained from an out-of-sample forecast exercise. First, there is a general risk that more flexible models are more prone to over-fitting the data and, as is widely known, over-fitting tends to hurt forecast performance. 26 In this case, our primary concern is that the model with both time-varying parameters and stochastic volatility-which is the most flexible specification-is most subject to this problem. Assessing this issue seems relevant. Second, and from a more general perspective, the number of out-of-sample forecast exercises conducted using BVARs with drifting parameters and/or stochastic volatility is fairly limited. 27 It therefore seems useful to also add to this literature. 26 See, for example, Clark (2004) or Rossi and Sekhposyan (2011) for related discussions. It can also be noted that the traditional Minnesota prior employed in many Bayesian VARs-see, for example, Litterman (1986)-was developed to address the problem with over-fitting in VARs. 27 We briefly mention a few important contributions here. Clark (2011) assessed the performance of a BVAR with stochastic volatility (and constant parameters and employing a steady-state prior) using US data on GDP growth, unemployment rate, inflation and the Federal funds rate. He concluded that there were modest gains to be made when it came to the accuracy of point forecasts but that relevant improvements could be made when it came to density forecasts. Using a BVAR with time-varying parameters and stochastic volatility, and employing US data on inflation, the unemployment rate and a short-term interest rate, D' Agostino et al. (2013) concluded that on average the model's forecasting performance was better than that of the alternative models. The benefits from using the model with time-varying parameters and stochastic volatility were particularly large when forecasting inflation. Barnett et al. (2014) considered a fairly wide range of models when forecasting GDP growth, inflation and the 3-month treasury bill rate in the UK. One of their findings The out-of-sample forecasting exercise is conducted using all four VAR models described earlier, that is, (i) constant covariance matrix and parameters, (ii) constant covariance matrix and time-varying parameters, (iii) stochastic volatility and constant parameters and (iv) stochastic volatility and time-varying parameters. In addition, we also introduce a univariate Bayesian AR model with constant parameters and covariance matrix for inflation-which is the variable we will focus upon-in order to get a simple benchmark model. The lag length is selected as five by the Akaike (1974) information criterion except for one forecast origin, 2005Q1, where the chosen lag length is one. 28 Regarding the variables used, it can be noted that we return to our main specification from Sect. 3 with the unemployment rate (age group 16-64 years) and CPIF inflation.
Out-of-sample forecasts are generated based on an expanding sample. We first estimate all models using data from 1995Q1 to 2004Q4 and then use the estimated models to generate forecasts up to eight quarters ahead with the median of the predictive distribution as the forecast. The sample is then extended one period, the models are re-estimated, and new forecasts up to eight quarters ahead generated. This continues until we generate the last forecasts based on the sample 1995Q1 to 2018Q2. We thus have a different number of forecasts to evaluate, from 55 forecasts at the one-quarter horizon to 48 at the eight-quarter horizon.
The results from the out-of-sample forecast exercise-in terms of the root mean square forecast errors (RMSFEs)-are given in Table 8. As can be seen, there is no single model that performs best at all horizons. The model with time-varying parameters and stochastic volatility is the best VAR model at most horizons (one to six quarters). It can be noted though that at the longest horizon, this model actually has the highest RMSFE of all models. At the two longest horizons (seven and eight quarters), it is instead the model with constant parameters and covariance matrix which performs best of the four VAR models. Differences are in general not particularly large though. For example, comparing the model with time-varying parameters and stochastic volatility to that with constant parameters and covariance matrix, the difference is never larger than 0.08 percentage points. Conducting a Diebold and Mariano (1995) test between these two models-under the assumption of quadratic loss-we cannot find statistical significance at any horizon. 29 Finding that the model with time-varying parameters and stochastic volatility has difficulties outperforming the model with constant parameters Footnote 27 continued was that a model with drifting parameters and stochastic volatility could generate improvements in forecast accuracy over a simple autoregressive benchmark (with constant parameters and variance) but that it was difficult to beat a VAR (with constant parameters and covariance matrix) estimated on a rolling sample. Aastveit et al. (2017) estimated models of various dimensions using US data but focused on the evaluation of the forecast accuracy with respect to GDP growth, unemployment rate and inflation. They found that the model with time-varying parameters and stochastic volatility had difficulties outperforming a model with constant parameters and covariance matrix. 28 As an alternative, we have also used the Bayesian information criterion (Schwarz 1978) to select the lag length, with lag length one chosen for the first 42 forecasts (up to forecast origin 2015Q2) and five lags for the remaining forecast origins. This has only a marginal effect on the RMSEs of the forecasts, and the conclusions about the relative performance of the models are the same. Results are not reported but are available upon request. 29 The absolute value of the test statistic is never larger than 0.98. Results are not reported but are available upon request. It can be noted that we use a "standard" Diebold and Mariano (1995) test even though the validity of this test in the present setting can be questioned. However, Diebold (2015) has pointed out that where u t is the unemployment rate and π p t is price inflation. "Constant" is the BVAR model with constant covariance matrix and parameters. "TVP" is the BVAR model with constant covariance matrix and time-varying parameters. "SV" is the BVAR model with stochastic volatility and constant parameters. "TVP-SV" is the BVAR model with stochastic volatility and time-varying parameters. "AR" is a univariate AR model for inflation. Out-of-sample forecasts generated using an expanding sample; the first sample is 1995Q1 to 2004Q4 and the last 1995Q1 to 2018Q2 and covariance matrix in out-of-sample forecasting, our results bear some resemblance to those of Aastveit et al. (2017).
It should, however, also be noted that the univariate AR model has the lowest RMSFE of all models at horizons one, six, seven and eight quarters. This finding is in line with what has been commonly pointed out in the literature on inflation forecasting, namely that simple univariate models tend to do quite well; see, for example, Stock and Watson (2009).
Given that the differences between the models are small, they should not be overinterpreted. What seems clear though is that the strong statements regarding the superiority of the model with time-varying parameters and stochastic volatility-for example, when it is compared to the model with constant parameters and covariance matrix-that the marginal likelihood comparisons in Table 1 show do not carry over to a clear improvement in out-of-sample forecast performance. One can also note regarding the VAR models that the model with constant parameters and stochastic volatility-which in previous sections received a lot of support from the DIC-never is the preferred specification according to the out-of-sample forecast analysis.
We conclude this section by noting that while marginal likelihood calculations tend to clearly favour the model with time-varying parameters and stochastic volatility, out-of-sample forecasts are less supportive of that model (even if it does perform best at several horizons). The fact that the differences in forecasting performance between the various models are small means that it is difficult to assess what the differences are due to. However, the finding that both the univariate AR model and the VAR with constant parameters and covariance matrix perform well in out-of-sample forecasting Footnote 29 continued if one uses Diebold and Mariano-type tests for model selection, one might as well use the "standard" test and associated critical values.
gives some support to old concern that in-sample measures might favour more complex models at the expense of forecasting performance.

Discussion
The analysis presented in Sects. 3 and 4 paints a fairly consistent picture. While the evidence is mixed regarding the stability of the Swedish Phillips curve, we consistently find that the Phillips curve has not been unusually flat during the last few years of the sample when inflation-despite efforts by the Riksbank to make monetary policy highly expansive-has been low. A flattened Phillips curve is, according to our analysis, hence not the cause for the low inflation in recent years.
Obviously, it is useful to be able to dismiss a potential explanation for the low inflation. But we would also of course like to know what caused it. The analysis presented in this paper suggests at least one explanation based on the model with time-varying parameters (and stochastic volatility). According to the model, the low inflation during the last few years of the sample can to some extent be explained by the low trend inflation that we found regardless of whether price or wage inflation was used. This might perhaps be perceived as a somewhat circular line of reasoning since the low outcomes for actual inflation are one reason for why the estimated trend inflation is low. But the findings of the model can be given an economic interpretation along the following lines: In the aftermath of the global financial crisis, the Riksbank has had substantial problems with making inflation reach the target rate. As can be seen from Fig. 1, CPIF inflation did not reach 2% for 6 years in a row, 2011-2016. In light of this development, economic agents in Sweden began to lose confidence in the inflation target which manifested itself in falling inflation expectations among other things. As pointed out above, this was a main reason for the Riksbank to abandon its leaning against the wind policy in 2014. The model's low trend inflation can be seen as reflecting the low inflation expectations that we saw during this period. There is evidence that the low inflation outcomes affected the wage formation process; as employers and employees no longer expected inflation around the target, wage agreements with low nominal increases were signed. Seeing that wage agreements tend to run for several years, this means that such effects based on expectations will be fairly long-lasting. As price inflation has increased, so has trend inflation in the model relying on prices. As new wage agreements will be negotiated in a period with higher price inflation, it is reasonable to assume that wage increase claims will be higher and that both wage inflation and the accompanying measure of trend inflation will increase. 30 30 That higher wage inflation will drive up price inflation is in line with our theoretical expectations as discussed above. It also finds support in modelling work that we have conducted. Estimating yet another bivariate BVAR with time-varying parameters and stochastic volatility, now with wage and price inflation as the variables in the system-that is, we set y t (π w t π p t ) -we find that a shock to wage inflation increases price inflation; see Fig. 27 in "Appendix". But while this empirical evidence points in the right direction, it can also be noted that the variance decomposition from this model with wage and price inflation indicates that shocks to wage inflation explain little of the forecast error variance of price inflation; see Fig. 28 in "Appendix". These results are in line with the findings of Peneva and Rudd (2017) who established that changes in labour costs seemed to have little effect on price inflation in the USA.
Summing up, we argue that part of the problem with low inflation that Sweden has experienced recently may have been self-inflicted and a consequence of the Riksbank's policy of leaning against the wind. That said, we also want to point out that other factors of course have affected inflation during this period. It can be noted that the broad findings put forward in this paper are in line with those of Karlsson and Österholm (2018) regarding the US Phillips curve; together with the fact that inflation developments over the last few years have been similar in several countries, this indicates that there are international determinants of inflation at work as well, for instance the commonly suggested globalization and digitalization.

Conclusions
In this paper we have assessed the stability of the Swedish Phillips curve. Conducting model selection using new methods, we found mixed evidence in favour a stable relationship between inflation and the unemployment rate. Based on this finding, we conclude that it is unclear whether the Phillips curve has been stable over time. However, while support for time-varying parameters (and stochastic volatility) can be found when looking at the models' marginal likelihood, our results also clearly indicate that the low inflation in Sweden in recent years is not due to a flatter Phillips curve.
In terms of explaining these low inflation outcomes in recent years, the models with time-varying parameters suggest that low trend inflation-reflecting low inflation expectations which likely affected the wage formation process-has contributed. This is in line with the findings of Clark (2014) and Watson (2014) who both point to the importance of a well-anchored trend level of inflation when it comes to explaining inflation outcomes in the USA. Other factors, such as globalisation and digitalisation, are likely to be important parts of the explanation as well given that many other countries also have seen low inflation outcomes during the same period. The results presented in this paper thereby add to our knowledge both concerning the properties of the Swedish Phillips curve but also international developments, although there is admittedly further work needed. From a more methodological viewpoint, it is worth pointing out that such further work should-to a larger extent than what is being done currently-consider the possibility of relying on frameworks that allow for time variation in the relationships between macroeconomic variables as well as the variance of the shocks that hit the economy. Price inflation data, model with constant parameters and stochastic volatility. Note: Based on bivariate BVAR with y t u t π p t , where u t is the unemployment rate and π p t is price inflation. Size of impulse is one standard deviation. Effect in percentage points on vertical axis. Horizon in quarters and dates on horizontal axes Fig. 13 Impulse-response function: Effect of shocks to inflation on the unemployment rate. Price inflation data, model with constant parameters and stochastic volatility. Note: Based on bivariate BVAR with y t u t π p t , where u t is the unemployment rate and π p t is price inflation. Size of impulse is one standard deviation. Effect in percentage points on vertical axis. Horizon in quarters and dates on horizontal axes Price inflation data, model with time-varying parameters and stochastic volatility. Note: Based on bivariate BVAR with y t u t π p t , where u t is the unemployment rate and π p t is price inflation. Size of impulse is one standard deviation. Effect in percentage points on vertical axis. Horizon in quarters and dates on horizontal axes

Fig. 27
Impulse-response function: Effect of shocks to wage inflation on price inflation. Model with constant parameters and stochastic volatility. Note: Based on bivariate BVAR with y t π w t π p t , where π w t is wage inflation and π p t is price inflation. Size of impulse is one standard deviation. Effect in percentage points on vertical axis. Horizon in quarters and dates on horizontal axes