EconBase
← Back to paper

The Origin and the Resolution of Nonuniqueness in Linear Rational Expectations

Extracted main text — title through conclusion, appendix excluded. This is what our citation measures are computed over, published so the extraction can be checked by eye.

99,712 characters · 17 sections · 49 citation commands

Rendered from LaTeX for readability, not typeset faithfully. Citation keys are highlighted; maths is left as source; figures, tables and equation environments are summarised rather than reproduced; unrecognised commands are greyed out so nothing is silently dropped. Email addresses are removed.

The Origin and the Resolution of Nonuniqueness in Linear Rational Expectations

abstractThe nonuniqueness of rational expectations is explained: in the stochastic, discrete-time, linear, constant-coefficients case, the associated free parameters are coefficients that determine the public's most immediate reactions to shocks. The requirement of model-consistency may leave these parameters completely free, yet when their values are appropriately specified, a unique solution is determined. In a broad class of models, the requirement of least-square forecast errors determines the parameter values, and therefore defines a unique solution. This approach is independent of dynamical stability, and generally does not suppress model dynamics. Application to a standard New Keynesian example shows that the traditional solution suppresses precisely those dynamics that arise from rational expectations. The uncovering of those dynamics reveals their incompatibility with the new I-S equation and the expectational Phillips curve. \begin{flushleft} {\bf JEL Classification:} E17, C53 \\ {\bf Keywords:} rational expectations; model-consistency; nonuniqueness; indeterminacy, New Keynesian models. \end{flushleft}

Introduction

The rational expectations hypothesis plays a central role in macroeconomics, but it has long been recognized that it need not determine expectations uniquely, and the resulting indeterminacy has raised objections on philosophical, mathematical, and practical grounds. This paper traces nonuniqueness to its economic origins, in the dynamics of expectation formation.

It has in effect been widely supposed that the heart of the problem of nonuniqueness might lie in an infinite regression (if it is possible to speak of “regression” into the future): if the values of endogenous variables depend on contemporary forecasts of their future values, then those forecasts of future values depend on future forecasts, which in turn depend on still later forecasts; and so forth. This apparent regression has suggested a need for a terminal boundary condition, and raised the hope that such a condition might help in determining a unique solution. But this paper shows that the origin of indeterminacy lies not in the infinitely remote future, but in the immediate present: in the stochastic, linear, discrete-time, constant-coefficients case, the associated free parameters determine the immediate response of expectations to shocks. The requirement of rationality -- or model-consistency -- generally leaves those parameters unconstrained; yet when their values are appropriately specified, a unique solution is determined.

In particular, if the requirement of model-consistency is strengthened with the additional requirement that forecasting errors be minimized in the least-squares sense, then, for a broad class of models, the free parameters are determined, and solutions are unique.

RePEc:oup:oxford:v:34:y:2018:i:1-2:p:43-54. has expressed the view that rational expectations are “insufficiently inertial”: it is shown here that the traditional approach to rational expectations, in producing a unique solution to a standard New Keynesian model, suppresses exactly those eigenvalues that arise from rational expectations. In contrast, the methods of this paper do not presuppose any stability or instability properties, nor do they generally suppress any model dynamics. The uncovering of the dynamics of rational expectations in the New Keynesian model reveals incompatibility with the new I-S equation and the expectational Phillips curve.

Taylor's example

The nature of the free parameters that underlie nonuniqueness is illustrated by an early example of RePEc:ecm:emetrp:v:45:y:1977:i:6:p:1377-85. Consider the following equation,

align*[align* omitted — 70 chars of source]

where $\/\delta_1\/$ is a nonzero real constant, $\/u_t\/$ is a sequence of independent, identically distributed, zero-mean random variables with finite variance, $\/p_t\/$ is an endogenous variable (specifically, proportional to the logarithm of the price of output)\footnote{This version of Taylor's equation has been linearized through a change of coordinates: here, $\/p_t\/$ denotes $\/p_t - \delta_0/\delta_1\/$ in Taylor's coordinates.}, and, for $\/i = 1,2\/$, $\/\hat{p}_{i,t}\/$ denotes a forecast of $\/p_{t+i}\/$, formulated at time $\/t\/$ (that is, formulated in terms of $\/u_\tau\/$, for $\/\tau \leq t\/$).

Seeking solutions of the form

align*[align* omitted — 54 chars of source]

Taylor imposes rational expectations by setting

align*[align* omitted — 141 chars of source]

Substituting these expressions into the model equation, he finds:

align*[align* omitted — 106 chars of source]

The coefficient $\/\pi_0\/$ is determined,\footnote{The analysis of appendix (ref) shows that $\/\pi_0\/$ is determined because the model equation contains no unlagged forecasts.} but $\/\pi_1\/$ is free, and its value determines all other coefficients. Taylor then turns to the imposition of a finite-variance condition (dynamical stability) and a minimum-variance condition, as means of limiting the possible values of $\/\pi_1\/$.

But what is the economic significance of the quantity $\/\pi_1\/$? According to the parameterization of $\/\hat{p}_{1,t-1}\/$,

align*[align* omitted — 75 chars of source]

so the coefficient $\/\pi_1\/$ determines the effect of $\/u_{t-1}\/$ on this forecast. In other words, this coefficient determines the immediate effect of shocks; it models the integration of new information into the forecast. This is an undeniably important parameter of the “expectations mechanism.”

Yet, $\/\pi_1\/$ is completely free under the assumption of rational expectations. As strong an assumption as rational expectations is, it does not determine, either in whole or in part, the integration of new information into the forecast. The results of this paper generalize this finding; they show that the free variables that account for the nonuniqueness of rational expectations are precisely the coefficients that govern the immediate effects of shocks on forecasts. Consequently, nonuniqueness is resolved if and only if those immediate effects are modeled unambiguously -- by other means.

The prevailing approach is to constrain such effects indirectly, principally by imposing a terminal condition, which requires dynamical stability. The use of the stability criterion is arbitrary, because it bears no particular relationship to the cause of nonuniqueness. It is also unrealistic, because it depends on infinite-precision cancellation of unstable dynamics. Taylor shows that, if $\/\delta_1\/$ is positive, then his model is stable if and only if $\/\pi_1\/$ is exactly zero. He derives the following recurrence:

align*[align* omitted — 119 chars of source]

The moving averages in $\/p_t\/$ and $\/u_{t}\/$ are of the same form if $\/\pi_1 = 0\/$, leading to a so-called `pole-zero cancellation' that suppresses the unstable eigenvalue at $\/1 + \delta_1\/$. The cancellation is displayed more explicitly by bringing in a `left-shift' operator $\/z\/$, so that $\/z x_t\/$ stands for $\/x_{t+1}\/$, and $\/z^{-1} x_t\/$ for $\/x_{t-1}\/$ (so $\/z^{-1}\/$ is a `right-shift' or `lag' operator). Then, assuming that normal algebraic operations can be applied to the operator,\footnote{This operator notation is formalized with the use of the z-transform for the purposes of the rest of the paper.} the above equation can be rewritten as,

align*[align* omitted — 111 chars of source]

Hence, if $\/\pi_1 = 0\/$,

align*[align* omitted — 112 chars of source]

The zero of the rational function in $\/z\/$ coincides with (and cancels) the unstable pole (the root of the denominator polynomial, which coincides with the characteristic polynomial) if and only if $\/\pi_1 = 0\/$. In that case, the resulting model has reduced-order dynamics -- and therefore, no dynamics at all.

This method of ensuring a unique solution depends on the model's possessing an appropriately unstable eigenvalue. It also requires the immediate response of the public's aggregate forecast $\/\hat{p}_{1,t}\/$ to the input $\/u_t\/$ to equal zero, exactly.

Another way of specifying the same unique value for $\/\pi_1\/$ in Taylor's example is to require in addition that the variance of the forecast errors be minimized. This paper shows that, in a broad class of models, that requirement determines a unique solution, regardless of stability properties, and generally without entailing pole-zero cancellation.

Background

The source of nonuniqueness of rational-expectations models lies in the reason for modeling expectations in the first place: forecasts have a bearing on the behavior of economic variables; moreover, they may affect the very quantities being forecast. The study of this self-referential phenomenon dates at least to Tinbergen:horizon, who examined the effect of forecast horizons on the movement of commodity prices. The work of RePEc:ucp:jpolec:v:62:y:1954:p:465, on the public prediction of social events, makes the circularity of the problem explicit, treating model-based expectations as fixed points. This property was later summarized by RePEc:eee:moneco:v:4:y:1978:i:1:p:1-44 as that of an “expectations mechanism which `reproduces itself' in [the] model.” The fixed-point characterization immediately raises the possibilities of nonexistence and nonuniqueness.

Muth:rational applied similar ideas to market dynamics, hypothesising that there were no systematic discrepancies between the forecasts of market participants and the predictions of economic theory. He called such forecasts {\em rational expectations\/}.\footnote{RePEc:ecj:econjl:v:101:y:1991:i:408:p:1245-53 compares Muth's contribution to Tinbergen's.} RePEc:eee:jetheo:v:4:y:1972:i:2:p:103-124,RePEc:eee:crcspp:v:1:y:1976:i::p:19-46 applied the rational expectations hypothesis to macroeconomic models, to show how changes in anticipated policy may give rise to changes in the behavior of economic agents. The advent of rational expectations is often described as a revolution in macroeconomics, and the approach has been developed by many economists, of different schools of thought. But as soon as it was applied to models with significant dynamics, and corresponding forecasts of the future outcomes of those dynamics, it was found that expectations were not defined uniquely RePEc:ecm:emetrp:v:45:y:1977:i:6:p:1377-85,RePEc:eee:moneco:v:4:y:1978:i:1:p:1-44;\footnote{RePEc:eee:jetheo:v:7:y:1974:i:1:p:53-65 focused on nonuniqueness, but under the relatively strong requirement that the initial conditions consistent with continual economic equilibrium, and a bounded rate of inflation, be unique.} the potential for nonuniqueness inherent in Grunberg and Modigliani's fixed-point characterization was realized.

In the absence of an explanation of the source of this nonuniqueness, the predominant response has been to impose a terminal boundary condition, limiting the growth rates of trajectories, in the hope that it is satisfied by a unique solution. Indeed, the term “solution” has almost come to mean “stable solution” Funovits201747. If the underlying model has appropriate unstable eigenstructure, its instabilities may be sufficient to narrow the space of stable solutions to a singleton.\footnote{See, for example, RePEc:ecm:emetrp:v:41:y:1973:i:6:p:1043-48,RePEc:eee:moneco:v:4:y:1978:i:1:p:1-44,MINFORD1979117,RePEc:ecm:emetrp:v:48:y:1980:i:5:p:1305-11,RePEc:cup:etheor:v:13:y:1997:i:06:p:877-888_00,RePEc:ier:iecrev:v:39:y:1998:i:4:p:1015-26,RePEc:kap:compec:v:20:y:2002:i:1-2:p:1-20,RePEc:aea:aecrev:v:94:y:2004:i:1:p:190-217.} This approach entails the cancellation of unstable dynamics in one fashion or another, which effectively requires the public to act, in aggregate, with infinite precision. The effect of the cancellation is to obliterate dynamical features of the model. The result of these practices has been that proponents and critics alike have confined their studies of rational expectations to pathologically special instances. Moreover, these ad hoc methods may still be insufficient to resolve nonuniqueness: for such cases, a variety of other ideas, such as minimum-variance solutions, or minimal state-variable realizations, have also been proposed RePEc:ecm:emetrp:v:45:y:1977:i:6:p:1377-85,BASAR1989591,RePEc:kap:itaxpf:v:6:y:1999:i:4:p:621-639,EH:learning. However, like dynamical stability, such criteria fail to get to the heart of the matter.

Overview

The main point of the article is to explain nonuniqueness by identifying the associated free parameters, which turn out to be parameters of the expectations mechanism itself. For the sake of generality, the analysis proceeds exclusively from simple, minimal assumptions, that are satisfied by traditional approaches to rational expectations. Beyond a standard, technical assumption, it is supposed only that forecasts depend linearly on initial conditions and on the model's driving variables, their dependence on the driving variables being representable by means of a linear difference equation with constant coefficients. This mild assumption echoes assumption 3 of Muth:rational; it is satisfied by conventional rational-expectations solutions; and it is necessary for the preservation of the linear, constant-coefficient structure of the model equations. In this paper, it is formalized in the form of equations called {\em forecasting mechanisms\/}.

The problem of deriving a general rational-expectations solution consists in solving for a forecasting mechanism subject to the constraint of model-consistency -- that is, subject to the constraint that the forecasting mechanism “reproduce itself in the model.” That problem in turn reduces to the solution of some generally singular, deterministic matrix difference equations that describe the interrelationships of model parameters. Because these equations do not lend themselves to a simple recursion either forward or backward in time, it is useful to solve them in the frequency domain, by means of the z-transform (and without loss of generality). Uniqueness demands only the appropriate specification of a parameter of the forecasting mechanism that governs the immediate response of forecasts to shocks; it is that parameter that distinguishes one fixed-point forecasting mechanism from another. In other words, it is that parameter that distinguishes different rational-expectations solutions.

The determination of the necessary form of any model-consistent forecasting mechanism (in sections (ref) and (ref)) leads to a necessary and sufficient condition for existence, and to a characterization of the general solution in terms of the aforementioned parameter (section (ref)). This result explains the nature of the nonuniqueness of rational-expectations models, with reference to the dynamics of expectation formation. Because of the minimal set of assumptions on which it is based, all of this analysis is general: it encompasses any reasonable approach to rational expectations under linear, constant-coefficient models.

This general picture can be tidied up through the assumption of a simple structural condition, called well-posedness. Well-posedness ensures existence of a model-consistent forecasting mechanism, and permits the realization of that forecasting mechanism in the relatively robust form of a predictor that incorporates feedback (section (ref)).

Conventional stability criteria, as ad hoc means of consuming the undesired degrees of freedom in the general solution, are antithetical to the spirit of this paper, but it is nevertheless shown for purposes of comparison (section (ref)) that the framework of this paper can reproduce well known observations on the “determinacy” of a small New Keynesian model. RePEc:oup:oxford:v:34:y:2018:i:1-2:p:43-54. has expressed the view that rational expectations are insufficiently inertial: it turns out that the pole-zero cancellation required for stability suppresses precisely those dynamics of the New Keynesian model that arise from expectations (section (ref)).

For a broad class of models, nonuniqueness can be eliminated simply by taking the assumption of rational expectations -- or unbiased forecasts -- one step further, and assuming that forecast errors are not only zero-mean, but are minimized, in the least-squares sense (section (ref)). The strengthened assumption is arguably milder than that of unstable pole-zero cancellations, and does not require any particular stability or instability properties of the model. Nor does it generally result in the cancellation of dynamics.

Because the methods of the paper are independent of dynamical stability considerations, they allow for the study of stabilization via policy. This point is illustrated in section (ref), where it emerges that the standard formulation of the New Keynesian model is incompatible with the dynamics of rational expectations.

A brief review is given in section (ref) of the vast related literature, and the paper concludes with some general suggestions for research. Appendices show how key results of the paper extend to more general models, give details of some proofs omitted from the main body of the paper, and outline relevant mathematical background, for consultation as necessary.

Problem formulation

The exposition is based on a simple, abstract model of CHO2015160, though the principles generalize (see appendix (ref)).

For all $\/t \in \mathbb{Z}\/$,

align[align omitted — 133 chars of source]

The matrices $\/A,\hat{A} \in \mathbb{R}^{n \times n}\/$ ($\/\hat{A} \neq 0\/$), $\/B \in \mathbb{R}^{n \times m}\/$ are constants.

The variables include the independent variable $\/t \in \mathbb{Z}\/$, representing discrete time instants, a vector of endogenous variables $\/x_t \in \mathbb{R}^{n \times 1}\/$, the vector $\/\hat{x}_{1,t} \in \mathbb{R}^{n \times 1}\/$ representing a one-time-step “forecast” of the value of $\/x_t\/$, and a vector of exogenous inputs $\/u_t \in \mathbb{R}^{m \times 1}\/$, driven by a sequence $\/w_t \in \mathbb{R}^{m \times 1}\/$ of real-valued, independent, zero-mean random variables, with finite variance,\footnote{Specifically, all entries of the covariance (or variance) matrix are finite.} defined on a common probability space. The forecast $\/\hat{x}_{1,t}\/$ may depend on the initial conditions $\/x_{-1}\/$, $\/\hat{x}_{1,-1}\/$ and $\/u_{-1}\/$, which are formally treated as constants, and on the sequence of random variables $\/u_0\/$, $\/u_1\/$,\ \ldots\ , $\/u_t\/$ -- or equivalently, the sequence $\/w_0, w_1,\ \ldots\ , w_t\/$.

It will be assumed that the polynomial matrix $\/[z^2 \hat{A} - z I + z]\/$ is {\em regular\/} -- meaning that its determinant does not vanish for all $\/z \in \mathbb{C}\/$, in which case the model too will be called regular. This is a common assumption, which serves in this instance to rule out a source of nonuniqueness that is unrelated to expectations.\footnote{The model is intended by Cho and McCallum to capture a local, linear approximation around the steady state of a dynamic, stochastic, general-equilibrium (DSGE) model. For a concrete, numerical instance, see section (ref). Under the conventional rational-expectations paradigm, the Cho-McCallum model has been used to capture equations with arbitrary finite numbers of expectation terms, having arbitrary expectational leads and lags broze_gouri�roux_szafarz_1995,RePEc:cup:etheor:v:13:y:1997:i:06:p:877-888_00,RePEc:eee:dyncon:v:31:y:2007:i:4:p:1376-1391; see Appendix (ref) for a direct extension of the results of the paper to such general models. }

Of course, the random variables $\/\hat{x}_{1,t}\/$ and $\/x_t\/$ are undefined, in the absence of additional equations. For the sake of generality, it will be assumed only that the forecasts $\/\hat{x}_{1,t}\/$ depend linearly on the initial conditions and the driving variables $\/w_t\/$ -- and moreover, that their dependence on the driving variables can be represented by a linear, constant-coefficient, stochastic difference equation. This is in essence a version of assumption 3 of Muth:rational, and it is satisfied by other approaches to rational expectations. The theory of linear, constant-coefficient difference equations then tells us that expectations must then obey an equation of the following form:

align*[align* omitted — 167 chars of source]

where $\/\tilde{F}_t \in \mathbbm{R}^{n \times m}\/$ vanishes for negative $\/t\/$, and $\/\overline{x}_{t+1}\/$ depends linearly on the initial conditions, but does not depend on the $\/w_t\/$. Such an equation will be called a {\em forecasting mechanism\/}. The convolution kernel determines the matrix of impulse responses of the forecasting mechanism, and the limits of the convolution sum ensure that forecasts are based on the appropriate information set.

The variables $\/w_t\/$ are in essence merely a convenient means of defining the stochastic structure of the $\/u_t\/$, so it is also of interest to consider forecasting mechanisms driven by the economic, exogenous variables $\/u_t\/$:

align*[align* omitted — 181 chars of source]

For every $\/t \geq 0\/$, $\/F_t \in \mathbbm{R}^{n \times m}\/$. The sequence $\/\tilde{u}_t := \sum_{\tau = 0}^t R^{t-\tau} w_\tau = u_t - R^{t+1} u_{-1}\/$ denotes the component of the sequence of endogenous variables $\/u_t\/$ that depends only on the $\/w_t\/$ and not on the initial conditions.

All of the analysis carried out in the next section follows from the above mild assumption on the form of forecasts, together with the aforementioned regularity property. It therefore represents a general analysis of rational expectations, within the context of linear, constant-coefficient, stochastic difference equations.

Of specific interest are forecasting mechanisms that are unbiased. Let $\/y\/$ be any square-integrable random variable defined on the common probability space of the $\/w_t\/$. For $\/t \geq -1\/$, $\/E_t(y)\/$ denotes the expected value of $\/y\/$, conditioned on the driving variables $\/w_0\/$, $\/w_1\/$,\ \ldots\ , $\/w_t\/$ -- and subject to the full model, including the forecasting mechanism, and its initial conditions. Given consistent initial conditions, a forecasting mechanism ((ref)) (respectively, ((ref))) is {\em model-consistent\/} if the full model ((ref)--(ref)) (respectively, ((ref),(ref),(ref))) satisfies

align*[align* omitted — 128 chars of source]

Such a forecasting mechanism embodies the “strong” rational-expectations hypothesis, whereby economic agents behave, in aggregate, as if they have access to all relevant information about the economy, and, on that basis, form expectations that do not incorporate any systematic errors.

It should be emphasized that the appropriate conditional expectations are subject to the forecasting mechanism, even if that forecasting mechanism is initially unknown to the modeler or analyst. This is a crucial point, which the author considers to be logically implied by, and entirely in the spirit of, the strong rational-expectations hypothesis. In economic terms, it represents the assumption that economic actors (in their aggregate) behave as if they know not only how the economy responds to shocks and to actors' expectations, but also how actors' own expectations are formed, and revised in response to shocks.

The main results of the paper include a general existence-and-uniqueness result for model-consistent forecasting mechanisms (Theorem (ref), page (ref)), and identification of a simple structural condition that ensures existence, as well as realizability in the form of combined feedforward/feedback implementations (Theorem (ref), page (ref)). Finally, it is shown that uniqueness can be ensured, for a broad class of models, by requiring not only that forecasting errors be zero-mean, but that they be minimized, in the least-squares sense (Corollary (ref), page (ref)).

The general model-consistent solution

To derive a general representation of model-consistent forecasting mechanisms, we consider separately the case where the initial conditions are zero-valued, and that where the driving variables are zero-valued; by linearity, we then superimpose the resulting solutions.

It is in the first case that the reason for nonuniqueness becomes apparent.

Zero-state response

Suppose that the initial conditions $\/x_{-1},\/$ $\/\hat{x}_{1,-1}\/$, and $\/u_{-1}\/$ are zero. In the terminology of linear, time-invariant systems, the resultant response is called the {\em zero-state response\/}.

By definition, the zero-state response of the forecasting mechanism must have the form of a convolution. It is convenient in the first instance to find a description in terms of the driving variables $\/w_t\/$:

align[align omitted — 158 chars of source]

Because the model ((ref),(ref)) is then effectively a system of linear, constant-coefficient equations, the sequence of endogenous-variable vectors must have the form of a similar convolution:

align*[align* omitted — 149 chars of source]

Naturally, $\/\tilde{G}_t\/$ will depend on $\/\tilde{F}_t\/$, and vice versa. The model and the model-consistency condition will furnish two equations describing their relationship.

Substitute the convolution sums ((ref),(ref)) into the equation ((ref)), and use ((ref)) to eliminate $\/u_t\/$:

align*[align* omitted — 444 chars of source]

Applying the conditional expectation operator $\/E_0\/$,

align*[align* omitted — 122 chars of source]

But here, $\/w_0 \in \mathbb{R}^{m \times 1}\/$ is arbitrary, so it must be that

align*[align* omitted — 167 chars of source]

For a second equation relating the two impulse responses, bring in the model-consistency condition for $\/t \geq 0\/$:

align*[align* omitted — 263 chars of source]

Once again taking expected values conditioned on $\/w_0\/$, and factoring out $\/w_0\/$ on the grounds that the resulting equation must hold for arbitrary $\/w_0\/$,

align*[align* omitted — 130 chars of source]

The calculation of the zero-state response amounts to solving the system ((ref),(ref)), and thus effectively reduces to the solution of the equation obtained by substituting $\/\tilde{G}_{t+1}\/$ for $\/\tilde{F}_t\/$ in ((ref)):

align*[align* omitted — 152 chars of source]

For this, it is assumed that the matrix polynomial $\/[z^2 \hat{A} - z I + A]\/$ is {\em regular\/} -- that its determinant is not identically zero. This is a standard assumption, for which there is ample justification. See, for example, RePEc:ier:iecrev:v:39:y:1998:i:4:p:1015-26,RePEc:eee:ecolet:v:61:y:1998:i:2:p:143-147. In particular, it is necessary for the uniqueness of solutions, irrespective of expectations. Under such a regularity assumption, any solution that exists is unique. Moreover, it is of exponential order, and therefore possesses a unilateral z-transform: consider that ((ref)) can be rewritten as a first-order matrix difference equation, via {\em linearization\/} of the above matrix polynomial Lanc:lin; and see, for example, the results of Bruell2009.

Applying the unilateral z-transform, and its left- and right-shift properties, to ((ref),(ref)), and keeping in mind that $\/\tilde{G}_t\/$ must vanish for negative $\/t\/$, one finds

align*[align* omitted — 267 chars of source]

The appearance of $\/\tilde{G}_0\/$ in the equations reflects the need for an additional boundary condition; nonuniqueness results when, in effect, only the condition $\/\tilde{G}_{-1} = 0\/$ is applied. But, given the form of the model ((ref),(ref)), such nonuniqueness can only arise from that of model-consistent forecasting mechanisms. The additional boundary condition should therefore be expressed in terms of $\/\tilde{F}_t\/$: setting $\/t = 0\/$ in ((ref)) and ((ref)),

align*[align* omitted — 144 chars of source]

Substituting $\/\hat{A}\tilde{F}_0 + B\/$ for $\/\tilde{G}_0\/$ then, solve for the transforms $\/\tilde{F}[z]\/$ and $\/\tilde{G}[z]\/$ in terms of $\/\tilde{F}_0\/$:

align*[align* omitted — 321 chars of source]

This establishes that, if suitable solutions of ((ref),(ref)) exist, they must have transforms of the above forms, and therefore must be unique. But the form of $\/\tilde{F}[z]\/$ also determines existence:

propositionSuppose that $\/[z^2 \hat{A} - zI + A]\/$ is regular. Then for any given value of the product $\/\hat{A} \tilde{F}_0 \in \mathbbm{R}^{n \times m}\/$, there exists a solution of the system ((ref),(ref)), such that $\/\tilde{F}_t\/$ and $\/\tilde{G}_t\/$ vanish for negative $\/t\/$, and $\/\hat{A}\tilde{F}_0\/$ has the specified value, if and only if the rational matrix \begin{align*} \tilde{F}[z] = [z^2 \hat{A} - z I + A]^{-1} \left [[zI - A](\hat{A}\tilde{F}_0 + B)[zI - R] - z^2 B \right ] [I - R z^{-1}]^{-1} \end{align*} is proper. In that case, the inverse z-transforms \begin{align*} \tilde{F}_t & :={\cal Z}^{-1}\left \{ [z^2 \hat{A} {-} z I {+} A]^{-1} \left [[zI {-} A](\hat{A}\tilde{F}_0 + B)[zI {-} R] - z^2 B \right ] [I {-} R z^{-1}]^{-1} \right \}\\ \tilde{G}_t & := {\cal Z}^{-1}\left \{ [z^2 \hat{A} {-} z I {+} A]^{-1} \left [ z \hat{A} (\hat{A}\tilde{F}_0 + B) [zI {-} R] - z B \right ] [I {-} R z^{-1}]^{-1} \right \} \end{align*} comprise the unique such solution.
proofSuppose that an appropriate solution of the system exists. Then $\/\tilde{F}_t\/$ and $\/\tilde{G}_t\/$ vanish for negative $\/t\/$. Their respective z-transforms are then proper rational matrices; and, by the preceding discussion, $\/\tilde{F}[z]\/$ and $\/\tilde{G}[z]\/$ have the form given in the statement of the proposition, where $\/\tilde{F}_0\/$ is the initial value of $\/\tilde{F}_t\/$. This establishes uniqueness (by z-transform inversion), and the necessary condition for existence. Conversely, for any specified value of $\/\hat{A} \tilde{F}_0\/$, if the given matrix $\/\tilde{F}[z]\/$ is proper, then so is \begin{align*} \tilde{G}[z] & = [z^2 \hat{A} - z I + A]^{-1} [z^2 \hat{A} (\hat{A} \tilde{F}_0 + B) - z B [I - R z^{-1}]^{-1}] \\ & = (\hat{A} \tilde{F}_0 + B) + [z^2 \hat{A} - z I + A]^{-1}[[zI {-} A] (\hat{A} \tilde{F}_0 + B) - z B [I {-} R z^{-1}]^{-1}] \\ & = (\hat{A}\tilde{F}_0 + B) + z^{-1} \tilde{F}[z]\ . \end{align*} Both $\/\tilde{F}[z]\/$ and $\/\tilde{G}[z]\/$ are therefore unilateral z-transforms. Because $\/z^{-1} \tilde{F}[z]\/$ is then strictly proper, the initial value of the inverse transform of $\/\tilde{G}[z]\/$ must equal $\/\hat{A} \tilde{F}_0 + B\/$. It follows that ((ref)) and ((ref)) are satisfied. Transforming back to the time domain then shows that ((ref)) and ((ref)), and consequently ((ref)), are satisfied. Applying ((ref)) at $\/t=0\/$ gives $\/\tilde{G}_0 = \hat{A} \tilde{G}_{1} + B\/$, so $\/\hat{A} \tilde{G_1}\/$ equals the specified value of $\/\hat{A} \tilde{F}_0\/$. But by ((ref)), $\/\tilde{F}_t = \tilde{G}_{t+1}\/$ for all $\/t \geq 0\/$. So the product $\/\hat{A}\tilde{F}_t\/$ indeed has the specified value at $\/t=0\/$. This proves the sufficient condition for existence. \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

The above result leads to a necessary condition on the form of model-consistent forecasting mechanisms ((ref)):

corollaryLet the model ((ref),(ref)) be regular. For any given value of the product $\/\hat{A} \tilde{F}_0 \in \mathbbm{R}^{n \times m}\/$, any model-consistent forecast mechanism ((ref)) must satisfy \begin{align*} \tilde{F}_t & ={\cal Z}^{-1}\left \{[z^2 \hat{A} {-} z I {+} A]^{-1} \left [[zI {-} A](\hat{A}\tilde{F}_0 + B)[zI {-} R] - z^2 B \right ] [I {-} R z^{-1}]^{-1} \right \}\ . \end{align*} with $\/\tilde{F}[z]\/$ proper. If $\/x_{-1}\/$, $\/\hat{x}_{1,-1}\/$, and $\/u_{-1}\/$ are zero, and $\/\hat{x}_{1,t} = \sum_{\tau = 0}^t \tilde{F}_{t - \tau} w_\tau\/$, then the model ((ref),(ref)) satisfies \begin{align*} x_t = \sum_{\tau = 0}^t \tilde{G}_{t - \tau} w_\tau\ , \end{align*} where \begin{align*} \tilde{G}_t & := {\cal Z}^{-1}\left \{ [z^2 \hat{A} - z I + A]^{-1} \left [ z \hat{A} (\hat{A}\tilde{F}_0 + B) [zI - R] - z B \right ] [I - R z^{-1}]^{-1} \right \} \ . \end{align*} The forecast error realized at time $\/t \geq 0\/$ is $\/x_t - \hat{x}_{1,t-1} = (\hat{A} \tilde{F}_0 + B) w_{t}\/$.
proofThe necessary conditions on the form of model-consistent forecast mechanisms, and on the solutions of the resulting models, are direct consequences of Proposition (ref), by the foregoing discussion. For the satisfaction of the model equations, we have, again by Proposition (ref), that equation ((ref)) is satisfied. By linearity, and the fact that $\/\tilde{G}_t\/$ vanishes for negative $\/t\/$, so then is ((ref)). It follows that the $\/\hat{x}_{1,t}\/$ and the $\/x_t\/$ in the statement of the corollary solve the model equations ((ref),(ref)). For the forecast dated at $\/t=-1\/$, then, \begin{align*} x_0 - \hat{x}_{1,-1} = x_0 = \tilde{G}_0 w_0 = (\hat{A} \tilde{F}_0 + B) w_0\ ; \end{align*} and by equation ((ref)), for all $\/t \geq 0\/$, \begin{align*} x_{t+1} - \hat{x}_{1,t} = \sum_{\tau = 0}^{t+1} \tilde{G}_{t+1 - \tau} w_\tau - \sum_{\tau = 0}^t \tilde{F}_{t - \tau} w_\tau = \tilde{G}_0 w_{t+1} = (\hat{A} \tilde{F}_0 + B) w_{t+1}\ . \addtocounter{equation}{1}\tag{\theequation} \end{align*} \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

These forecast errors are zero-mean, confirming model-consistency.

The solution can be expressed in terms of the exogenous inputs $\/u_t\/$ rather than the driving variables $\/w_t\/$, by simply right-multiplying $\/\tilde{F}[z]\/$ and $\/\tilde{G}[z]\/$ by $\/[I - R z^{-1}]\/$:

align*[align* omitted — 206 chars of source]

The time-domain counterpart of right-multiplying by $\/[I - R z^{-1}]^{-1}\/$ to recover $\/\tilde{F}[z]\/$ and $\/\tilde{G}[z]\/$ is convolution with $\/R^t\/$. It follows that $\/F_0 = \tilde{F}_0\/$. It also follows that convolution of $\/F_t\/$ (respectively, $\/G_t\/$) with $\/u_t\/$ is equivalent to convolution of $\/\tilde{F}_t\/$ (resp., $\/\tilde{G}_t\/$) with $\/w_t\/$ (by the associativity of convolution). Because $\/[I - R z^{-1}]\/$ and its rational-matrix inverse are both proper, and because the initial value of each of their inverse z-transforms is the identity matrix, multiplication by either of them does nothing to alter properness.

The following counterpart of Corollary (ref) is immediate:

corollarySuppose that the model ((ref),(ref)) is regular. Then, for any specified value of $\/ \hat{A} F_0 \in \mathbbm{R}^{n \times m}\/$, any model-consistent forecasting mechanism ((ref)) has \begin{align*} F_t = {\cal Z}^{-1}\{F[z]\} = {\cal Z}^{-1}\{[z^2 \hat{A} - z I + A]^{-1} [[zI {-} A](\hat{A}F_0 + B)[zI {-} R] - z^2 B]\}\ , \end{align*} where $\/F[z]\/$ is proper. If $\/x_{-1}\/$, $\/\hat{x}_{1,t}\/$, and $\/u_{-1}\/$ are all zero, and $\/\hat{x}_{1,t} {=} \sum_{\tau = 0}^t {F}_{t - \tau} u_\tau\/$, then the resulting full model ((ref)-(ref)) satisfies \begin{align*} x_t = \sum_{\tau = 0}^t {G}_{t - \tau} u_\tau\ , \end{align*} where \begin{align*} G_t & = {\cal Z}^{-1}\{G[z]\} = {\cal Z}^{-1}\{[z^2 \hat{A} - z I + A]^{-1} z [ \hat{A} (\hat{A}F_0 + B)[zI - R] - B]\} \ . \end{align*} The forecast error realized at time $\/t \geq 0\/$ is $\/x_{t} - \hat{x}_{1,t-1} = (\hat{A} {F}_0 + B) w_{t}\/$. \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

If $\/{F}[z]\/$ and $\/{G}[z]\/$ are proper, they constitute what are commonly called the {\em transfer matrices\/}, or matrices of {\em transfer functions\/}, of a forecasting mechanism and of the corresponding full model, respectively. When transfer matrices are proper, they have state-space realizations. Such realizations are not unique, and their choice may depend on the details of the structure of the rational matrices. But realizations can be found easily with the aid of numerical routines such as the command {\tt tf2ss\/} of the MATLAB Control Systems Toolbox, or of Scilab scilab. See section (ref) for concrete, numerical examples.

The especially alert reader will have noticed that the formulas for $\/\tilde{G}[z]\/$ and $\/G[z]\/$ implicitly assume the exact cancellation of the term $\/[I - A z^{-1}]^{-1}\/$, arising from the model equation ((ref)) via equation ((ref)), by the term $\/[I - Az^{-1}]\/$ arising in $\/\tilde{F}[z]\/$. (To see the cancellation, substitute the solution for $\/\tilde{F}[z]\/$ into equation ((ref)).) However, the implementation of $\/\tilde{F}[z]\/$ given in section (ref) ensures that the latter term also arises from the model equation itself, via feedback, in which case the cancellation is a matter of algebra, and does not demand infinite precision in the aggregate actions of the public.

A related concern is that the overall model may be very sensitive to the value of $\/\hat{A}\/$: mathematically, variations in that value could give rise to singular perturbations of the nominal matrix polynomial $\/[z^2 \hat{A} - zI + A]\/$. A thorough treatment of the latter robustness issue is beyond the scope of this paper; but an element of a resolution might be to treat $\/\hat{A}\/$ solely as a parameter of the forecasting mechanism, whose output would then be the product $\/\hat{A} \hat{x}_{1,t}\/$ -- the `economic effect' of the forecast $\/\hat{x}_{1,t}\/$, as opposed to the forecast itself.

Zero-input response

To characterize model-consistent forecasting mechanisms fully, the effect of nonzero initial considerations must be considered. That of the driving terms having been analyzed, it suffices now (by linearity) to consider the case where the initial conditions may have nonzero values, but all of the driving variables $\/w_t\/$ vanish. In system-theoretic terms, the corresponding solution of the model is called its “zero-input response.”

By the definition of forecasting mechanisms, in this case, $\/\hat{x}_{1,t} = \overline{x}_{t+1}\/$, where $\/\overline{x}_{t+1}\/$ depends linearly on the initial conditions, but not at all on the $\/w_t\/$. For any $\/t \geq 0\/$, model-consistency therefore requires that

align*[align* omitted — 61 chars of source]

In other words, $\/\hat{x}_{1,t}\/$ must be an exact forecast of $\/x_{t+1}\/$.

The sequence $\/x_t\/$ must therefore be the solution of the following “perfect-foresight” model, which captures the case where the driving variables $\/w_t\/$ are zero-valued for all $\/t \geq 0\/$.

align*[align* omitted — 231 chars of source]

Substituting for $\/\hat{x}_{1,t}\/$ in the first equation,

align*[align* omitted — 151 chars of source]

Taking unilateral z-transforms and solving for $\/X[z]\/$,

align*[align* omitted — 207 chars of source]

Here, the second equation of the perfect-foresight model has been used to replace $\/x_0\/$ with the initial condition $\/\hat{x}_{1,-1}\/$.

propositionSuppose that the model ((ref),(ref)) is regular. Let \begin{align*} \overline{x}_t & := {\cal Z}^{-1}\{\overline{X}[z]\} \\ & = {\cal Z}^{-1}\left \{[z^2 \hat{A} - z I + A]^{-1} [z^2 \hat{A} \hat{x}_{1,-1} - z A x_{-1} - z B R [I - R z^{-1}]^{-1} u_{-1}] \right \}\ . \end{align*} Then, a perfect-foresight solution exists if and only if $\/\overline{X}[z] - \hat{x}_{1,-1}\/$ is strictly proper. In that case, the unique such solution has $\/x_t = \overline{x}_t\/$ and $\/\hat{x}_{1,t-1} = \overline{x}_{t+1}\/$.
proofBy the above discussion, any solution $\/x_t\/$ of ((ref)) must have a unilateral z-transform of the form of $\/\overline{X}[z]\/$. As a unilateral z-transform, $\/\overline{X}[z]\/$ must be proper. Moreover, because $\/x_0\/$ must equal $\/\hat{x}_{1,-1}\/$, the matrix $\/\overline{X}[z] - \hat{x}_{1,-1}\/$ must be strictly proper. This establishes the necessary condition for existence of a perfect-foresight solution. Now, suppose that $\/\overline{X}[z] - \hat{x}_{1,-1}\/$ is strictly proper. That implies that $\/\overline{X}[z]\/$ is proper, and that $\/\hat{x}_{1,-1}\/$ is the value of the inverse transform $\/\overline{x}_t\/$ at $\/t=0\/$. Rearrange the expression for $\/\overline{X}[z]\/$: \begin{align*} \overline{X}[z] & = A z^{-1}[\overline{X}[z] + z x_{-1}] + \hat{A} z [\overline{X}[z] - \hat{x}_{1,-1}] + B R [I - R z^{-1}]^{-1} u_{-1}\ . \end{align*} Transforming to the time domain, that inverse transform $\/\overline{x}_t\/$ is seen to satisfy ((ref)). The uniqueness of this solution follows from that of the transform $\/\overline{X}[z]\/$. If, in addition, $\/\hat{x}_{1,t-1} = \overline{x}_{t}\/$, for all $\/t \geq 0\/$, then ((ref),(ref)) is satisfied. The full model ((ref)--(ref)) then satisfies $\/x_t = \overline{x}_t\/$, if $\/w_t = 0\/$, $\/\forall t \geq 0\/$. \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

As a consequence of Proposition (ref), the initial conditions will be said to be {\em consistent\/} if $\/\overline{X}[z] - \hat{x}_{1,-1}\/$ is strictly proper; and in that case, a model-consistent forecasting mechanism must have $\/\overline{x}_t = {\cal Z}^{-1}\{\overline{X}[z]\}\/$.\footnote{Like the state-space realization of proper rational matrices $\/F[z]\/$ and $\/G[z]\/$, the inversion of the proper matrix $\/\overline{X}[z]\/$ can be carried out with the use of standard software tools; it amounts to computing the impulse response of a system with a given proper transfer matrix.}

By the last part of the proposition, $\/\overline{x}_t\/$ can be said, in system-theoretic terminology, to represent the {\em zero-input response\/} of the full model ((ref)-(ref)), under model-consistent forecasts. It shows how the system responds to nonzero initial conditions, when any driving terms vanish for $\/t \geq 0\/$.

However, the nature of the derivation should be borne in mind in potential applications that might otherwise exceed the limitations of the results. If the model has been evolving through negative time instants, and the model parameters that hold for $\/t \geq 0\/$ differ from those in effect for $\/t = -1\/$, it is plausible to suppose that realized values of $\/x_{-1}\/$ and $\/u_{-1}\/$ are valid initial conditions for the zero-input response for $\/t \geq 0\/$. But it is less clear what are the implications of realizations of the forecast $\/\hat{x}_{1,-1}\/$. “Perfect foresight” need not apply under an unforeseen change in model parameters.

Like the formulas derived in the previous section, the above solution for $\/\overline{X}[z]\/$ implicitly assumes that a term $\/[I - A z^{-1}]^{-1}\/$ arising from the model is exactly cancelled by a term $\/[I - A z^{-1}]^{-1}\/$ derived from the forecast mechanism. But this concern can be addressed by realizing the forecast mechanism using feedback from the model equation, as in section (ref).

Total response

The total response of the full system is the sum of its zero-state and zero-input responses. Indeed, fix any model-consistent forecasting mechanism; then $\/x_t\/$ and $\/\hat{x}_{1,t}\/$ must be linear functions of the random variables $\/w_\tau\/$, $\/0 \leq \tau \leq t\/$, and of the initial conditions $\/x_{-1}\/$, $\/\hat{x}_{1,-1}\/$, and $\/u_{-1}\/$. It follows by linearity that they must then be obtained by summing the separate respective responses to the driving variables and to the initial conditions -- namely, the zero-state and zero-input responses.

theoremSuppose that the model ((ref),(ref)) is regular. For any $\/\hat{A} F_0 \in \mathbbm{R}^{n \times m}\/$, define \begin{align*} {F}_t & = {\cal Z}^{-1}\{F[z]\} = {\cal Z}^{-1}\{[z^2 \hat{A} - z I + A]^{-1} [[zI {-} A](\hat{A}F_0 + B)[zI {-} R] - z^2 B]\}\ , \\ {G}_t & = {\cal Z}^{-1}\{ [z^2 \hat{A} - z I + A]^{-1} z[ \hat{A} (\hat{A}F_0 + B)[zI {-} R] - B]\} \ ,\ and \\ \overline{x}_t & = {\cal Z}^{-1}\{\overline{X}[z]\} \\ & = {\cal Z}^{-1}\left \{[z^2 \hat{A} - z I + A]^{-1} [z^2 \hat{A} \hat{x}_{1,-1} - z A x_{-1} - z B R [I {-} R z^{-1}]^{-1} u_{-1}] \right \}\ . \end{align*} Suppose that the initial conditions are consistent ($\/\overline{X}[z] - \hat{x}_{1,-1}\/$ is strictly proper). Then there exists a model-consistent forecasting mechanism ((ref)) for ((ref),(ref)) if and only if $\/F[z]\/$ is proper. In that case, in terms of the above inverse transforms, the unique model-consistent forecasting mechanism ((ref)) is \begin{align*} \hat{x}_{1,t} = \sum_{\tau = 0}^t {F}_{t-\tau} \tilde{u}_\tau + \overline{x}_{t+1} \ ,\ \forall t \geq 0\ ; \addtocounter{equation}{1}\tag{\theequation} \end{align*} and the resulting full model ((ref)--(ref)) satisfies \begin{align*} x_t = \sum_{\tau = 0}^t {G}_{t-\tau} \tilde{u}_\tau + \overline{x}_t \ ,\ \forall t \geq 0. \addtocounter{equation}{1}\tag{\theequation} \end{align*} The forecast error realized at time $\/t \geq 0\/$ is \begin{align*} (\hat{A} F_0 + B) w_t = (\hat{A} F_0 + B)(u_t - R u_{t-1})\ . \end{align*}
proofIf the initial conditions are consistent, then the necessary form of a model-consistent forecasting mechanism, and that of the unique solution of the full model resulting from a forecasting mechanism of that form, follow from Corollary (ref) and Proposition (ref), by the linearity of the model and of conditional expectations. The terms derived from the zero-input solution have no effect on the forecasting errors, so (again by linearity) the overall forecasting error realized at time $\/t\/$ is $\/(\hat{A} F_0 + B) w_t\/$, as derived in section (ref). It follows that the specified forecasting mechanism ((ref)) is indeed model-consistent. \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

It should be emphasized that, assuming that a suitable value of $\/\hat{A} F_0\/$, representing the immediate economic effects of shocks on endogenous variables via expectations, can be specified, the results of this section resolve the nonuniqueness of rational expectations: as long as a model-consistent forecasting mechanism exists, $\/\hat{A} F_0\/$ determines $\/G_0\/$, and “rationality” then determines $\/F_t\/$ and $\/G_t\/$ for all $\/t > 0\/$.

In the next section, a simple assumption is introduced that ensures existence for any value of $\/\hat{A} F_0\/$.

Well-posedness, existence, and feedback

A regular model ((ref),(ref)) will be called {\em well-posed\/} if the inverse of the “characteristic matrix” $\/[-z \hat{A} + I - A z^{-1}]\/$ is proper -- or equivalently, if $\/[z^2 \hat{A} - z I + A]^{-1}\/$ is strictly proper. For example, this is so whenever $\/\hat{A}\/$ is nonsingular, but not when $\/\hat{A}\/$ is nilpotent.

Well-posedness admits a simple sufficient condition for consistency of the initial conditions. The initial conditions will be called {\em weakly consistent\/} if

align*[align* omitted — 79 chars of source]

Indeed, this condition is plainly necessary if a solution of ((ref),(ref)) is to exist when $\/w_t = 0\/$ and the prediction $\/\hat{x}_{1,-1}\/$ is exact. Together, well-posedness and weak consistency imply the existence of a model-consistent forecasting mechanism for any possible value of $\/\hat{A} F_0\/$. In economic terms, well-posedness obviates any assumption that economic actors use their presumed aggregate knowledge of the model to “choose” parameter values $\/\hat{A} F_0\/$ for which solutions exist.

The theoretical definition of forecasting mechanisms is an unrealistic one, in the sense that, incorporating no feedback, a direct implementation would be completely lacking in robustness. But it turns out that well-posedness also allows for the realization of model-consistent forecasting mechanisms in the form of feedforward/feedback interconnections with the rest of the model, as represented by ((ref),(ref)). (The condition implies the “well-posedness” of that interconnection, in the specific sense in which that term is applied to feedback systems.) The use of feedback resolves key robustness issues with respect to the model parameter $\/A\/$, but sensitivity to the parameter $\/\hat{A}\/$ remains an issue (discussed at the end of section (ref)).

For brevity, proofs for this section are relegated to appendix (ref). The results themselves are summarized in the following:

theoremIf the model ((ref),(ref)) is regular and well-posed, and the initial conditions are weakly consistent, then for every possible value of the product $\/\hat{A} F_0 \in \mathbbm{R}^{n \times m}\/$, there exists a unique model-consistent forecasting mechanism. That forecasting mechanism can be realized by the following feedforward/feedback law: \begin{align*} \hat{x}_{1,t} & = \sum_{\tau = 0}^{t} \varPhi_{t-\tau} [A x_\tau + B R u_\tau] - \sum_{\tau=0}^{t} \Psi_{t-\tau} \hat{A} F_0 w_\tau - \Psi_t (\hat{x}_{1,-1} - A x_{-1} - B R u_{-1})\/ . \end{align*} Here, $\/\varPhi_t := {\cal Z}^{-1}\{ [I - z \hat{A}]^{-1}\}\/$ and $\/\Psi_t := {\cal Z}^{-1}\{ [I - z \hat{A}]^{-1} z \hat{A} \hat{A}^g \}\/$, where $\/\hat{A}^g\/$ is a generalized inverse of $\/\hat{A}\/$ (s.t. $\/\hat{A} \hat{A}^g \hat{A} = \hat{A}\/$). \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

Conventional determinacy of a New Keynesian model

In this section it is shown -- purely for purposes of comparison -- how the general solution of the previous section lends itself to the reproduction of conventional results.

Consider for example the loglinearized New Keynesian DSGE model of RePEc:aea:aecrev:v:94:y:2004:i:1:p:190-217:

align*[align* omitted — 611 chars of source]

The scalar variables $\/y_t\/$, $\/\pi_t\/$, and $\/r_t\/$ respectively represent output, inflation, and the nominal interest rate, expressed as percentage deviations from a trend path or a steady state; $\/g_t\/$ and $\/z_t\/$ represent the effects of exogenous shifts on the first two equations. In accordance with our usual notation, $\/\hat{y}_{1,t}\/$ and $\/\hat{\pi}_{1,t}\/$ represent forecasts of $\/y_{t+1}\/$ and $\/\pi_{t+1}\/$, dated at time $\/t\/$. We shall write $\/w_t = [\epsilon_{g,t} \ \epsilon_{z,t} \ \epsilon_{r,t}]^\prime\/$, and consider distinct values of the resulting vector-valued sequence to be independent, identically distributed, and zero-mean.

The scalar coefficients are as follows: $\/\tau\/$ represents intertemporal substitution elasticity, $\/\beta\/$ is the households' discount factor, $\/\kappa\/$ is the slope of the expectational Phillips curve; the third equation describes the monetary authority's behavior, $\/\psi_1\/$ and $\/\psi_2\/$ being `Taylor-rule' coefficients.

Letting $\/x_t = [y_t \ \pi_t \ r_t]^\prime\/$, $\/u_t = [g_t \ z_t \ \epsilon_{r,t}]^\prime\/$, and $\/w_t = [\epsilon_{g,t} \ \epsilon_{z,t} \ \epsilon_{r,t} ]^\prime\/$, it is an easy matter to put the equations into the form ((ref),(ref)). Setting the coefficients equal to the mean values given in Table 1 of RePEc:aea:aecrev:v:94:y:2004:i:1:p:190-217, with $\/\beta = 0.99\/$, according to RePEc:eee:dyncon:v:28:y:2003:i:2:p:273-285, one finds

align*[align* omitted — 801 chars of source]

(the respective matrix coefficients being the values of $\/A\/$, $\/\hat{A}\/$, $\/B\/$, and $\/R\/$). This model satisfies well-posedness, because all nine entries of $\/[z^2 \hat{A} - zI + A]^{-1}\/$ are strictly proper, so a model-consistent forecasting mechanism exists for every possible value of $\/\hat{A} F_0\/$.

According to the chosen values of the coefficients, $\/\psi_1\/$, the coefficient of inflation in the interest-rate policy `Taylor rule,' has the value $\/1.10\/$. If $\/\psi_1 > 1\/$, then when inflation rises, the monetary authority policy raises the interest rate by a greater percentage (all other things being equal): such a policy is said to be {\em active\/}; if $\/\psi_1 < 1\/$, policy is called {\em passive\/}. With $\/\psi_1 = 0.90\/$, for example, the general solution has only one unstable eigenvalue, while for $\/\psi_1 = 1.10\/$, it has two.\footnote{Throughout the paper, computations were performed in Scilab (version 5.5.2, for Mac OS X).} In the conventional approach to rational expectations, the model is found to be “indeterminate” in the passive case, but “determinate” in the active case: only when there are two unstable eigenvalues does their suppression consume sufficient degrees of freedom to yield a unique solution. Specifically, when $\/\psi_1 = 1.10\/$, the general solution has unstable eigenvalues at $\/z = 1.4461829\/$ and $\/z = 1.0446352\/$. The corresponding left eigenvectors of the denominator polynomial are respectively,\footnote{The results of these and related calculations are displayed with seven or eight significant digits, as a reminder that the calculations really require infinite precision.}

align*[align* omitted — 152 chars of source]

On the other hand, if $\/\psi_1 = 0.9\/$, and policy is “passive,” then the model has only one unstable eigenvalue.

Again, the imposition of stability is carried out here strictly for the purpose of comparison. Following RePEc:aea:aecrev:v:94:y:2004:i:1:p:190-217, assume that all initial conditions are zero-valued, and therefore focus on the zero-state response. Because the matrix $\/R\/$ is stable, it suffices to consider the matrix polynomial $\/G[z]\/$. The free parameter $\/\hat{A} F_0\/$ appears only in the numerator of the matrix-fraction description of $\/G[z]\/$, so the only way its value can be chosen so as to stabilize an otherwise unstable model is by arranging for unstable `poles' of $\/G[z]\/$ to be canceled by `zeros.' In the multivariable case, this means that the numerator matrix polynomial of $\/G[z]\/$ must have the same unstable eigenvalues, with the same respective left eigenvectors, as the denominator matrix polynomial.

Note that $\/\lambda_i \in \mathbb{C}\/$ is an eigenvalue of the numerator matrix polynomial of $\/G[z]\/$, with left eigenvector $\/c_i\/$, for both $\/i = 1\/$ and $\/i = 2\/$, if and only if

align*[align* omitted — 175 chars of source]

Each row of the above equation represents three equations, each in a distinct pair of unknowns (the entries from the first two rows of a distinct column of $\/\hat{A} F_0 + B\/$). So a single pole-zero cancellation does not determine a unique solution, but the two simultaneous cancellations do. The first two columns of

align*[align* omitted — 62 chars of source]

are linearly independent, yielding a unique solution for the first two rows of $\/\hat{A} F_0 + B\/$; approximately,

align*[align* omitted — 126 chars of source]

This in turn yields the first two rows of $\/F_0\/$,

align*[align* omitted — 143 chars of source]

The values of the first two rows of $\/F_0\/$ determine a unique value for $\/\hat{A} F_0\/$ -- approximately,

align*[align* omitted — 199 chars of source]

So the requirement of dynamical stability of the model determines a unique value of $\/\hat{A} F_0\/$, and therefore, by Theorem (ref), a unique model-consistent forecasting mechanism.

The resulting matrix $\/G[z]\/$ is the following:

align*[align* omitted — 276 chars of source]

It can be realized in the form of the following state-space model:\footnote{In this instance, by means of the Scilab command {\tt tf2ss}.}

align*[align* omitted — 423 chars of source]

Whereas a minimal state-space realization of $\/G[z]\/$ is normally third-order (see section (ref) for an example), at this isolated point in the parameter space of $\/\hat{A}F_0\/$ it is first-order. This confirms the two pole-zero cancellations, and reproduces the results of the conventional approach to determinacy.

In order for the two unstable eigenvalues to be suppressed in this manner, the initial values of the impulse responses of the forecasts of output and inflation to the three shocks must equal, exactly, the values given (approximately) by the respective entries of the first two rows of $\/F_0\/$ given above: unless the aggregate forecasts of the public respond to the respective shocks in this manner, with infinite precision, the model will be unstable.\footnote{In double-precision floating-point, it takes only a multiplicative perturbation of $\/F_0\/$ on the order of $\/1 \pm10^{-9}\/$ to spoil the pole-zero cancellation.} The complete lack of robustness of this cancellation makes it untenable that it models any real-world phenomenon.

Moreover, the next section shows that the effect of this suppression is precisely to eliminate the eigenvalues (and the dynamics) that arise from expectations.

Rational expectations and stability

The relationship between expectations and stability is a longstanding concern AN:exp_stability. This section presents an asymptotic analysis of the effects of expectation terms on dynamics.

If a `small' scalar multiplicative weight $\/\epsilon\/$ is attached to the matrix coefficient $\/\hat{A}\/$ of the forecast term, then unless that matrix is nilpotent, the denominator matrix $\/[z^2 \epsilon \hat{A} - zI + A]\/$ of the overall model is a singular perturbation of that of the lower-order model with $\/\epsilon = 0\/$: the degree of the denominator polynomial is a fixed integer greater than $\/n\/$ for all $\/\epsilon > 0\/$; but for $\/\epsilon = 0\/$ it is equal to $\/n\/$.

It turns out that the overall model is always unstable when that weight is positive but sufficiently small. Indeed, if $\/z\/$ remains bounded as $\/\epsilon\/$ tends to zero, then the denominator polynomial tends to $\/- [zI - A]\/$, so $\/n\/$ of the eigenvalues tend toward those of the matrix $\/A\/$. But if $\/\hat{A}\/$ has nonzero eigenvalues, then the denominator polynomial has more than $\/n\/$ finite eigenvalues: in order for the degree to drop at $\/\epsilon = 0\/$, some of those finite eigenvalues must `escape' to infinity.

Indeed, suppose that $\/\hat{A}\/$ has $\/m > 0\/$ nonzero eigenvalues. To capture the behavior of eigenvalues that vary like $\/1/\epsilon\/$ as $\/\epsilon\/$ tends to zero, perform the change of variable $\/z = \lambda/\epsilon\/$. The matrix polynomial $\/[z^2 \epsilon \hat{A} - z I + A]\/$ becomes

align*[align* omitted — 149 chars of source]

As $\/\epsilon\/$ tends to zero, $\/m\/$ eigenvalues of the matrix polynomial in $\/\lambda\/$ approach the reciprocals of the nonzero eigenvalues of $\/\hat{A}\/$. Therefore, the moduli of $\/m\/$ eigenvalues of the matrix polynomial in $\/z = \lambda/\epsilon\/$ tend to infinity.

On the other hand, if $\/\hat{A}\/$ is nonsingular,\footnote{Eigenvalues at infinity complicate the establishment of upper bounds for the finite eigenvalues HIGHAM20035.} and the weight applied to expectations is sufficiently {\em large\/}, then the eigenvalues of the matrix polynomial in $\/z\/$ will be stable:

corollary\\ {\bf Small expectation gain:} In equation ((ref)) above, replace the coefficient matrix $\/\hat{A}\/$ with $\/\epsilon \hat{A}\/$, where $\/\epsilon\/$ is a real, positive scalar. Suppose that $\/\hat{A}\/$ has some nonzero eigenvalue. Then, for sufficiently small $\/\epsilon > 0\/$, the denominator matrix polynomial $\/[z^2 \epsilon \hat{A} - z I + A]\/$ has unstable eigenvalues; consequently, barring pole-zero cancellations, the full model ((ref)-(ref)) is dynamically unstable under model-consistent expectations. \\ {\bf Large expectation gain:} On the other hand, suppose that $\/\hat{A}\/$ is nonsingular. Then the modulus of any eigenvalue of $\/[z^2 \hat{A} - z I + A]\/$ is at most \begin{align*} \frac{1 + \sqrt{1 + 4 \lVert \hat{A}^{-1}\rVert^{-1} \lVert A \rVert}}{2 \lVert \hat{A}^{-1} \rVert^{-1}}\ . \end{align*} where $\/\lVert \ldots \rVert\/$ denotes any subordinate matrix norm.
proofThe “small-gain” result follows from Corollary 1 of ABG2004, incorporated into a comprehensive theory in AGM2014a. The upper bound for the case of nonsingular $\/\hat{A}\/$ is from Lemma 3.1 of HIGHAM20035. \nobreak \ifvmode \relax \else \ifdim\lastskip<1.5em \hskip-\lastskip \hskip1.5em plus0em minus0.5em \fi \nobreak \vrule height0.75em width0.5em depth0.25em\fi

The above analysis shows that the unstable pole-zero cancellation required to eliminate unstable eigenvalues of the New Keynesian model (of section (ref)) under “active” policy would have the effect of obliterating the dynamical features particularly associated with expectations terms. Indeed, the eigenvalues of $\/[z^2 \epsilon \hat{A} - z I + A]\/$ vary continuously with the parameter $\/\epsilon\/$, and letting $\/\epsilon\/$ range from one to zero, numerical computation shows that it is precisely the two eigenvalues that are unstable under active policy that tend to infinity as $\/\epsilon\/$ tends to zero (moving to the right along the real axis), while the other eigenvalues tend toward those of $\/A\/$. Now consider reversing the process, starting with $\/\epsilon = 0\/$ and therefore no expectation terms in the model, and gradually increasing $\/\epsilon\/$ to $\/1\/$ to restore the expectations. The eigenvalues of $\/A\/$ are $\/0\/$, $\/0\/$, and $\/0.417\/$; as $\/\epsilon\/$ increases from $\/0\/$ to one, the two eigenvalues at the origin remain fixed, while the third shifts from $\/0.417\/$ to $\/0.334\/$. One effect of expectations is therefore to produce a modest shift in the nonzero eigenvalue of $\/A\/$. The more important dynamical effect of expectations, then, is to bring into being two additional modes corresponding to the finite eigenvalues $\/1.045\/$ and $\/1.04\/$. However, these are precisely the modes that are suppressed in order to select a unique model-consistent forecasting mechanism under the conventional approach to rational expectations.\footnote{The suppression also cancels “zero dynamics” of the model, which can explain behavior such as “price puzzles” TM:nofreelunch.} That traditional approach to rational expectations therefore destroys the very dynamical features that arise from expectations.

The next section presents one method of determining a unique model-consistent forecasting mechanism without resorting to pole-zero cancellation. Indeed, it yields a unique solution, regardless of dynamical stability or instability.

The least-square-error solution

This section proposes a unique model-consistent solution that, at least superficially, seems a natural extension of rational expectations. It applies in all cases (regardless of dynamical stability or instability), and does not generally entail pole-zero cancellation. In addition to the requirement that forecast errors be zero-mean, it calls for all forecast errors be minimized, in the least-squares sense: the key to the approach is that a unique value of $\/\hat{A}F_0\/$ will achieve such minimization.

Specifically, let $\/e_t\/$ denote $\/(\hat{A} F_0 + B) w_t\/$, the forecast error realized at time $\/t\/$, under model-consistent expectations. Then it is assumed that $\/\hat{A} F_0\/$ is such that, given the sequence of the $\/w_t\/$, every squared-error term $\/e_t^\prime e_t\/$ is minimized.

This criterion is not proposed as a panacea. On the contrary, it has features that may be undesirable, at least in some contexts. Because $\/G_0 = \hat{A} F_0 + B\/$, error minimization amounts to a minimization (in the least-squares sense) of the immediate response of the endogenous variables in $\/x_t\/$ to shocks. It implies that $\/\hat{A}F_0\/$ is `selected' in such a way that the immediate response of the model to shocks via expectations $\/\hat{x}_{1,t}\/$ at least partially cancels its immediate response to shocks via the exogenous variables $\/u_t\/$. Partial `blocking' of the effect of shocks by expectation terms is a fixture of the pole-zero cancellation employed in the traditional solution of rational-expectations models, but that does not imply that it is realistic. If the column span of $\/\hat{A}\/$ contains that of $\/B\/$, then $\/\hat{A}F_0\/$ can be chosen so that $\/\hat{A} F_0 + B\/$ is zero, and the immediate effects of shocks are fully blocked: the model will then exhibit `perfect foresight,' or `self-fulfilling expectations,' and the minimum squared-error terms will of course all be zero. In particular, the above will apply whenever $\/\hat{A}\/$ is nonsingular -- as it must be in the scalar case, for instance. Another consequence of the vanishing of $\/\hat{A}F_0 + B\/$ is that $\/G_0\/$, the initial value of the impulse response of the overall model vanishes, and $\/G[z]\/$ is strictly proper.

But while this new assumption clearly represents a significant strengthening of the usual rational-expectations hypothesis, it is in a similar spirit; and in comparison with the infinite-precision, unstable pole-zero cancellations that are usually assumed in conjunction with rational expectations, it is arguably relatively mild. It preserves the key feature that expectations are policy-dependent, if, as in the example of section (ref), the policy parameters are reflected in the matrix coefficients of the model. But it does not require specific stability properties of the model, and it generally avoids the feature that the unique solution has reduced-order dynamics, owing to pole-zero cancellation.

To define the unique solution, note that each column $\/b\/$ of the matrix $\/B\/$ can be uniquely decomposed into $\/b_\parallel + b_\perp\/$, where $\/b_\parallel\/$ is its projection onto the column span of $\/\hat{A}\/$, and $\/b_\perp\/$ is orthogonal to that vector space. Let $\/B_\parallel\/$ consist of the projections $\/b_\parallel\/$ of the respective columns of $\/B\/$, and $\/B_\perp\/$ of the respective orthogonal components $\/b_\perp\/$. Then, for any $\/t \geq 0\/$,

align*[align* omitted — 377 chars of source]

Because this is a sum of nonnegative quantities, its minimum value over all possible $\/\hat{A}F_0\/$ must be at least $\/w_t^\prime B_\perp^\prime B_\perp w_t\/$.

Now, for any column $\/b\/$ of $\/B\/$, the corresponding column $\/f_0\/$ of $\/F_0\/$ can be chosen (not necessarily uniquely) so that $\/\hat{A} f_0 = - b_\parallel\/$. Let $\/F_0\/$ be composed of such columns. Then $\/\hat{A} F_0 = - B_\parallel\/$, and; for any $\/t \geq 0\/$,

align*[align* omitted — 75 chars of source]

This unique value $\/\hat{A} F_0 = - B_\parallel\/$ therefore achieves the minimum squared error, for any given $\/w_t\/$.

The foregoing discussion establishes the following

corollarySuppose that the model ((ref),(ref)) is regular. Let $\/\hat{A} F_0 = - B_\parallel\/$. Suppose that the resulting $\/F[z]\/$ is proper and the initial conditions are consistent. Then the forecasting mechanism given by Theorem (ref) is the unique model-consistent forecasting mechanism that gives rise to least-square forecasting errors. By Theorem (ref), a weaker sufficient condition is that the model be well-posed and the initial conditions weakly consistent.

To be sure, the approach outlined above requires the public, in its aggregate, to calibrate its immediate reaction to shocks so as to meet the requirement of least-square forecasting errors. But here, the unique solution is not as ill-conditioned or `brittle' as under the method of pole-zero cancellation: a small error in the value of $\/\hat{A}F_0\/$ generally means only that forecast errors will be slightly larger than necessary; whereas in the case of pole-zero cancellation, the slightest error means that the model will be dynamically unstable.\footnote{In fact, in instances where the least-square-error criterion calls for pole-zero cancellation, a practical means of producing a less pathological solution might be simply to apply a `small,' random perturbation to the value of $\/\hat{A}F_0\/$.}

New Keynesian example

For the example of section (ref), we find that,

align*[align* omitted — 192 chars of source]

Setting $\/\hat{A}F_0 = - B_\parallel\/$, and finding a minimal state-space realization for the resulting matrix $\/F[z]\/$, we arrive at a representation of the unique, least-squared-error, model-consistent forecasting mechanism:

align*[align* omitted — 903 chars of source]

The second matrix coefficient of the second equation being the initial value of the system's impulse response, it equals $\/F_0\/$; multiplying it by $\/\hat{A}\/$ yields the above $\/-B_\parallel\/$. The model is well-posed, so this forecasting mechanism could be alternatively be represented as a feedback/feedforward predictor, according to Theorem (ref).

Coupling the forecasting mechanism to the model equation ((ref)) yields a system that can be represented in the form of the following state-space model (obtained as a realization of the corresponding rational matrix $\/G[z]\/$):

align*[align* omitted — 870 chars of source]

The coefficient $\/\hat{A}F_0 + B\/$ that determines the forecast errors is the same as $\/G_0\/$, and therefore equal to the second matrix coefficient of equation ((ref)).

In contrast with the traditional solution, this one does not have a reduced dynamical order, and the preservation of model dynamics turns out to reveal difficulties with the formulation of the New Keynesian model. Indeed, the above state-space equations are at odds with the standard interpretation of key model equations. For instance, the New I-S equation ((ref)) is typically considered to assert that the higher the real interest rate $\/r_t - \hat{\pi}_{1,t}\/$, the lower output $\/y_t\/$; and by the same token, the greater $\/g_t\/$, the higher the output $\/y_t\/$. But the state-space equations disagree. Note that $\/g_t\/$ is the first component of the vector $\/u_t\/$, and $\/y_t\/$ is the first component of $\/x_t\/$. According to ((ref)), $\/g_t\/$ has no immediate effect upon any of the components of $\/x_t\/$; and by ((ref)) its only effect on the forecasts of interest is a negative effect on $\/\hat{y}_{1,t}\/$.

But it is not only the least-squared-error solution that disagrees with the usual interpretation of the New Keynesian model. It is shown in the next section that, barring cancellations of eigenvalues of $\/[z^2 \hat{A} - z I + A]\/$, an analysis of the model eigenvalues under any model-consistent forecasting system conflicts with the usual interpretations of both the I-S equation and the expectational Phillips curve.

Stabilization of the New Keynesian model

By decoupling rational expectations from dynamical stability, the methods of this paper give a clearer picture of the dynamics of macroeconomic models. An illustration is provided by the formulation of a stabilizing interest-rate policy for the New Keynesian model.

If the $\/R\/$ matrix is stable, then in the absence of pole-zero cancellations, the stability of the overall model depends only on the eigenvalues of the matrix polynomial $\/[z^2 \hat{A} - z I + A]\/$. Consequently, one can play at a toy version of central banking by altering the Taylor-rule parameters in such a way that all of that matrix polynomial's eigenvalues become stable. For instance, with $\/\psi_1 = 1.10 \/$ and $\/\psi_2 = -1.50\/$, the eigenvalues are $\/0.81 \pm 0.045 i\/$ and $\/0.76\/$ (while their moduli can be decreased further by decreasing the value of $\/\rho_r\/$). These indicate a stable model, with a response that is only slightly oscillatory. Under the conventional approach to rational expectations, these Taylor-rule parameters would not produce a `determinate' model, dynamical stability being incompatible with determinacy. But in general, there is no such difficulty: for instance, this stabilization procedure could be carried out in conjunction with the use of the above least-square-error criterion.

Note that the value of the above Taylor-rule coefficient $\/\psi_2\/$ applied to output is negative, meaning that the higher the level of output, the lower is the policy interest rate. The reason for this surprising feature lies in the peculiarity of the new I-S equation. The equation is derived from microfoundations, and specifically from the solution of optimal planning problems for members of different sectors of the economy. Households, for example, are assumed to schedule their consumption over all time, according to expected values of nonlinear functions of interest rates and other variables. Under that scenario, the higher the expected real interest rate at time $\/t\/$, the lower consumption will indeed be at the same time -- moreover, the lower consumption will be at earlier time instants as well -- the better to profit from a higher rate.

Linear model equations are then obtained by “loglinearizing” the first-order conditions associated with optimality -- a process of informal approximation that typically assumes, for instance, that expected values of nonlinear functions can be approximated as expected values of approximations of those functions. The new I-S equation is derived by identifying output with consumption in the resulting system of equations. Unlike the solution of the consumption planning problem, this equation of course does not imply that the higher the present real interest rate, the lower earlier levels of consumption. In fact, contrary to the usual interpretation, it is difficult to argue that it implies a negative relationship even with the current level of output: for any alteration in the current real interest rate there will almost surely be a cognate change in the forecast $\/\hat{y}_{1,t}\/$ of the next period's output, so it is not obvious how $\/y_t\/$ will be affected. Better to rearrange the equation as follows,

align*[align* omitted — 142 chars of source]

and see that, with $\/\tau > 0\/$, it means simply that, the larger the real interest rate, the larger the expected increase in output over the next period.

But more can be said in the context of the overall model: note that the denominator matrix polynomial $\/[z^2 \hat{A} - z I + A]\/$ is exactly the same in the case of the zero-input response (section (ref)) as for the zero-state response ((ref)). In order to study model eigenvalues it therefore suffices (in the absence of pole-zero cancellations) to consider the zero-input response. But if the sequences $\/\epsilon_{r,t}\/$, $\/\epsilon_{g,t}\/$, $\/\epsilon_{z,t}\/$ are identically zero, then, according to the analysis of section (ref), the I-S equation and the expectational Phillips curve effectively become,

align*[align* omitted — 261 chars of source]

The result is a perfect-foresight version of the New Keynesian model, that gives rise to the same eigenvalues as the general version. For any given initial conditions, it has a unique solution, according to which the values of $\/\pi_{t+1}\/$ and $\/y_{t+1}\/$ are determined by their previous values and by $\/r_{t}\/$, $\/g_{t}\/$ and $\/z_{t}\/$. It follows that the larger $\/r_t\/$, the {\em larger\/} $\/y_{t+1}\/$. So the greater the interest rate, the {\em greater\/} will be output at the next instant. Similarly, with $\/\beta\/$ essentially unity and $\/\kappa > 0\/$, the expectational Phillips curve ((ref)) effectively implies that the higher the level of output, the {\em lower\/} the expected rate of increase of inflation over the next period. For use with model-consistent expectations, both equations are therefore flawed.

In determining inflation in the full model, the two inversions essentially cancel each other, because the only transmission channel from the interest rate to inflation goes by way of output. That explains why the sign of $\/\psi_1\/$ in the above stabilizing policy rule conforms to economic intuition, while that of $\/\psi_2\/$ does not. On the other hand, if one crudely `corrects' the two key equations by simply changing the signs of $\/\tau\/$ and $\/\kappa\/$, one finds -- as expected -- that in stabilizing policy rules, both $\/\psi_1\/$ and $\/\psi_2\/$ are positive. Indeed, for stability, the policy reactions to both inflation and the output gap should be `active': with $\/\psi_2 = 1.5\/$, the model is stable provided $\/1.03 \leq \psi_1 \leq 1.49\/$; and with $\/\psi_1 = 1.03\/$, the model is stable provided $\/1.04 \leq \psi_2 \leq 1.50\/$. In the toy game of central banking, one might opt for the policy gains $\/\psi_1 = 1.10\/$ and $\/\psi_2 = 1.50\/$: that would give eigenvalues at $\/0.812 \pm 0.0453 i\/$ and $\/0.763\/$; these are about as far from the unit circle as can be arranged with those two policy parameters, and they produce a response that is only slightly oscillatory.\footnote{The moduli of the eigenvalues could be reduced further by reducing the value of $\/\rho_r\/$ that serves to smooth changes in the policy interest rate.}

In a recent, constructive critique of DSGE models, RePEc:oup:oxford:v:34:y:2018:i:1-2:p:43-54. asserts that the new I-S equation ((ref)) and the expectational Phillips curve ((ref)) are “deeply flawed.” The results of this section lend support to that statement, but not for all of the reasons that Blanchard cites. He describes rational expectations as insufficiently “inertial.” The fundamental reason for that assertion is likely the suppression of the associated dynamics that is part and parcel of the conventional application of rational expectations. However, that suppression is not a general feature of rational expectations, but only of the very particular approach that has held sway for the last several decades.

Related work

In seeking a way forward for DSGE models, RePEc:oup:oxford:v:34:y:2018:i:1-2:p:43-54. asks, “how can we deviate from rational expectations, while keeping the notion that people and firms care about the future?” This paper suggests that the problem is not that rational expectations have been tried and found wanting, but rather that they have never been tried (with apologies to G.K. Chesterton).

The main reason is that the cause of nonuniqueness was never before understood, chiefly because earlier work did not include any explicit parameterization of the forecasting mechanism. RePEc:ecm:emetrp:v:45:y:1977:i:6:p:1377-85 adapted the early method of Muth:rational to dynamic macroeconomic models by means of a parameterization of the overall model (see subsection (ref)) in the form of a Wold decomposition. While Taylor's formulation resembles the present one in some respects, it assumes dynamical stability, an unnecessary limitation. More to the point, Taylor did not explicitly parameterize the forecasting mechanism. RePEc:eee:moneco:v:4:y:1978:i:1:p:1-44 recognized that the key to the solution of rational-expectations models was to solve for the forecasting mechanism, but neither did he include such an explicit parameterization. The same is true of all other previous approaches; for recent examples, see TAN201595 and al-sadoon_2017. Yet the free parameters that underlie the nonuniqueness of solutions are essentially parameters of the forecasting mechanism itself.

Without insight into the true reason for nonuniqueness, and perhaps because of a generally perceived need for a terminal boundary condition, stability quickly became entrenched as the primary means of attempting to select a unique solution. RePEc:ecm:emetrp:v:45:y:1977:i:6:p:1377-85 and RePEc:aea:aecrev:v:69:y:1979:i:2:p:114-18 considered alternatives, but -- invoking the stability constraint himself only shortly afterward -- RePEc:aea:aecrev:v:71:y:1981:i:1:p:132-43 referred to its application as having already become “a standard if not entirely convincing practice.” It has remained so for decades, even though far from universally applicable, and in spite of the recognition that it obliterates model dynamics and alters stability properties RePEc:aea:aecrev:v:94:y:2004:i:1:p:190-217.

Consequently, standard methods of computing rational-expectations solutions are needlessly bound up with the question of stability. It is generally assumed that an equation of the following form holds:

align*[align* omitted — 72 chars of source]

The vector $\/y_t\/$ is made up of endogenous variables, including forecasts, $\/v_t\/$ is an exogenous random driving variable, and $\/\eta_t\/$ is a vector of zero-mean forecasting errors, that is determined as part of the solution process.\footnote{See, for example, RePEc:ecm:emetrp:v:48:y:1980:i:5:p:1305-11,RePEc:cup:etheor:v:13:y:1997:i:06:p:877-888_00,RePEc:eee:dyncon:v:24:y:2000:i:10:p:1405-1423,RePEc:kap:compec:v:20:y:2002:i:1-2:p:1-20,RePEc:aea:aecrev:v:94:y:2004:i:1:p:190-217.} A matrix decomposition is applied, to allow separate consideration of the respective eigenvalues, and where unstable eigenvalues are concerned, it is arranged (if possible) for the corresponding terms in $\/v_t\/$ and $\/\eta_t\/$ to cancel each other, giving rise to a zero that cancels the unstable eigenvalue. If the corresponding calculations determine the forecasting errors uniquely, then there is a unique rational-expectations solution -- among those that obey the stability constraint. As this article proves, this is merely an arbitrary and indirect means of unknowingly modeling the public's immediate responses to shocks. But the approach has become so closely identified with the solution of rational-expectations models that it should be stressed that the link with stability is completely unnecessary -- as is the unrealistic, exact cancellation of unstable eigenvalues, and the consequent suppression of dynamics.

RePEc:ecm:emetrp:v:45:y:1977:i:6:p:1377-85 considered, as an alternative means of ensuring uniqueness, the requirement that the variance of the endogenous variables be minimized. But BASAR1989591 appears to have been the first to study explicitly the minimization of the variance of forecast errors. Ba\c{s}ar considered a univariate model with a single forecast term $\/\hat{x}_{2,t-1}\/$, driven by an independent, zero-mean sequence with finite variance, and assumed that the information set underlying forecasts consisted either of $\/x_{t-1}\/$, or of a noisy measurement of that scalar variable, with all random variables being Gaussian in the latter case. He showed that a minimum-variance rational-expectations solution could be computed as the limiting case of solutions to finite-horizon problems -- under an assumption that amounts to the realness of the two eigenvalues of the overall model. Specifically, the finite-horizon problems called for minimization of a discounted sum of squared forecast errors, over a finite time interval, and subject to a terminal boundary condition. Apart from the restriction on information sets, Ba\c{s}ar made no assumptions as to the form of forecasting mechanisms, but his solutions satisfied linear recurrences with constant coefficients. General solutions forTaylor's and Ba\c{s}ar's models are given by the results of appendix (ref), together with formulas for the forecast errors. In both cases, the space of model-consistent solutions can be characterized by two parameters, $\/\tilde{G}_0\/$ and $\/\tilde{G}_1\/$, of which the first is determined (because the models contain no unlagged expectation terms) and the second is completely free (because the models satisfy the appropriate generalization of the well-posedness condition). Forecast-error variances are positive, and are minimized by setting $\/\tilde{G}_1 = 0\/$ -- this amounts to nulling the initial value of the impulse response of an unlagged one-step forecast.

Conclusion

This paper has shown that the free parameters that underlie the nonuniqueness of solutions of linear rational-expectations models are coefficients that capture the immediate reactions of economic actors to shocks. It has presented a general solution that has a unique instance for every appropriate specification of those free parameters. If the model is {\em well-posed\/}, and the initial conditions {\em weakly consistent\/}, then there exists a (unique) solution for every possible specification of the parameters.

If, in addition, the requirement of model-consistency is augmented with that of least-square forecast errors, then -- for a broad class of models -- unique values of the free parameters are determined, and therefore so is a unique solution. This result is offered as one concrete means of ensuring uniqueness. It avoids the most serious drawbacks of the prevalent approach, while preserving its advantages, but its suitability for a given application should be considered critically in the light of the paper's main results, and the corresponding form of the solution that it yields.

Indeed, the fundamental objective of the paper has been to explain the nature of the nonuniqueness of rational expectations, so that arbitrary, ad hoc solutions might be avoided. In that connection, the main finding is that model consistency is not as strong an assumption as may have been supposed: it generally implies nothing about the economic effects of the most immediate responses of forecasts to shocks. To put it another way, there is no fundamental reason for any model of the most short-term reactions of economic agents to entail systematic forecasting errors. A rough paraphrase suggests itself, in the terminology of Kahneman:thinking: rational expectations exemplify analytical, “slow” thinking, and as such do not capture the “fast” thinking that underlies the most immediate responses to economic shocks -- but clearly, any deliberate and reasoned process of expectation formation must necessarily take account of the economic effects of more instantaneous, reflexive actions. Like most loose paraphrases, the analogy probably should not be pushed too far, but it should be underlined that “fast thinking” does not necessarily imply on-the-fly improvisation, and could include the automatic application of a preconditioned response, such as might indeed be represented by the free parameters identified here.

By dispensing with the cancellation of unstable dynamics, the results of the paper remove substantial obstacles to the application and the study of rational expectations. They greatly expand the range of models to which rational expectations can be applied, and they avoid the suppression of model dynamics. Consequently, they allow the examination of fundamental questions of stability and of stabilization. They should also simplify problems that are central to the pertinence of the rational expectations hypothesis, such as those of model estimation and the “learning” of models by economic agents RePEc:eee:moneco:v:4:y:1978:i:1:p:1-44,RePEc:red:issued:13-148.

More generally, these new results provide a richer picture of the dynamics of macroeconomic models, and allow for more reasoned and realistic means “to attribute to individuals some view of the behavior of the future values of variables of concern to them” RePEc:eee:crcspp:v:1:y:1976:i::p:19-46.