Finance and Economics Discussion Series: 2007-14 Screen Reader version ^{♣}

Keywords: bifurcation, perturbation, relative price distortion, optimal monetary policy

Abstract:

In perturbation analysis of nonlinear dynamic systems, the presence of a bifurcation implies that the first-order behavior of the economy cannot be characterized solely in terms of the first-order derivatives of the model equations. In this paper, we use two simple examples to illustrate how to
detect the existence of a bifurcation. Following the general approach of Judd (1998), we then show how to apply l'Hospital's rule to characterize the solution of each model in terms of its higher-order derivatives. We also show that in some cases the bifurcation can be eliminated through
renormalization of model variables; furthermore, renormalization may yield a more accurate first-order solution than applying l'Hospital's rule to the original formulation.

JEL Classification: C63; C61; E52.

In recent analysis of nonlinear dynamic macroeconomic models, the characterization of their first-order dynamics has been an important step for understanding theoretical implications and evaluating empirical success. However, the presence of a bifurcation in perturbation analysis of nonlinear dynamic systems implies that the first-order behavior of the economy cannot be characterized solely in terms of the first-order derivatives of the model equations.

In this paper, we use two simple macroeconomic models to address several issues regarding bifurcations. In particular, the bifurcation problem would emerge in conjunction with the price dispersion generated by staggered price setting in the part of firms. We then show how to apply l'Hospital's rule to characterize the solution of each model in terms of its higher-order derivatives. We also show that in some cases the bifurcation can be eliminated through renormalization of model variables; furthermore, renormalization may yield a more accurate first-order solution than applying l'Hospital's rule to the original formulation.

Before presenting our results, it is noteworthy that our definition of bifurcation is distinct from the one analyzed in Benhabib and Nishimura (1979). In particular, their analysis on bifurcation is associated with time evolution of dynamic systems. However, our concern with bifurcation arises in the process of approximating nonlinear equations, as discussed in Judd (1998).

We proceed as follows. Section 2 describes the two examples and illustrates how to detect the existence of a bifurcation problem. Section 3 follows the general approach of Judd (1998) and applies l'Hospital's rule to characterize the first-order behavior of each model. Section 4 shows how the bifurcation can be eliminated through renormalization of model variables. Section 5 concludes.

This section discusses how we can detect the existence of bifurcation in two simple economies. In both models, Calvo-style price setting behavior of firms can be summarized by the following law of motion for the relative price distortion:

where the distortion index is defined as The parameters and represent the percentage of firms that cannot change their price in each period and the elasticity of substitution across goods , respectively. The variable is the gross inflation rate of the price index aggregated over firms.

To discuss the issue of bifurcation, we have to close the model with another equation. In the first example, we simply assume that inflation follows an exogenous stochastic process,

where the logarithm of follows a mean zero process. We can rationalize this process in terms of monetary policy by a version of strict inflation targeting around the exogenous process or a version of strict output-gap targeting in a model with cost-push shocks.By combining the two equations, we now have a single-equation model:

Since this equation is backward looking, this exact nonlinear form can be used for any dynamic analysis. However, we suppose that we have to rely on approximation methods to analyze this model as would be the case when there are forward-looking equations.

Woodford (2003) and Benigno and Woodford (2005) pointed out that, when deviations of the (net) inflation rate from its zero steady state are of first order in terms of exogenous variations, deviations of the distortion index from one is of second order. Based on this observation, one can naturally approximate the system with respect to the square root of the logarithm of the relative price distortion index. Note that this distortion index is unity at the steady state with zero inflation rate. We follow the convention of using lower cases for log deviations,

Specifically, corresponds to the approximation in Woodford (2003) and Benigno and Woodford (2005). It will be also shown in Section 4 that can be used as the basis of an alternative approximation.

Under the choice of as the approximation variable, (2.2) can be rewritten as follows:

Now let's see what happens if we try a Taylor approximation of this system with respect to and . It is easy to see that the derivative with respect to the endogenous variable would be zero at the steady state. Based on this zero derivative, we can diagnose the bifurcation problem in this case. Put in an alternative way, the implicit function theorem cannot be applied when the derivative with respect to the endogenous variable is zero.

It is noteworthy to find out what would happen if we feed this case into computer codes commonly available for dynamic macroeconomic analysis. The Dynare package (version 3.05) produces an message saying ``Warning: Matrix is singular to working precision`', and AIM (developed by Gary
Anderson and George Moore, and widely used at the Federal Reserve Board) returns a code indicating ``Aim: too many exact shiftrights`'. The routine developed by Christopher Sims (`gensys.m`) ends without any output or error message.

The second example is a case with multiple equations. Our example is a prototypical Calvo-style sticky-price model, and the optimal policy problem is to maximize the household welfare subject to the following four constraints: the law of motion for relative price distortions, the social resource constraint, the firms' profit maximization condition, and the present-value budget constraint of the household. However, it is shown in Yun (2005) that the optimal policy problem can be reduced to minimizing the index for relative price distortion (2.1). At the optimum, we have the following relationship:

Therefore, the solution to the optimal policy problem can be represented with the following bivariate nonlinear system:

As in the single-equation case, we start with a normalization according to which and are endogenous variables and is exogenous:

As explained in Judd (1998) and Judd and Guu (2001), the bifurcation problem can be resolved by using l'Hospital's rule.

To understand the approximated behavior of in the singe-equation example, we need to compute and where is defined as an implicit function as follows:

In cases for which regular perturbation analysis could be applied, the first-order approximation of would come from the implicit function theorem as follows: The number in the parenthesis indicates the order of approximation. However, the assumption of the implicit function theorem does not hold in our case since . We need to adopt an advanced asymptotic method--the bifurcation method in this case.Noting that the derivatives in the numerators are also zero at the steady state, we apply l'Hospital's rule to the two ratios in the form of and obtain the following first-order
approximation:^{1}

This is an example of the transcritical bifurcation.

In this single-equation model, it is easy to avoid the bifurcation problem when we consider the following equation that is equivalent to (2.3),

To illustrate how we can invoke the bifurcation method in the multi-dimensional example, we substitute the second equation in (2.5) into the first to obtain:

0 | ||

Were the assumptions of the bifurcation theorem to hold, then differentiation of the implicit expression with respect to would produce the equation

However, since both derivatives on the right-hand side are zero at the steady state, we need to apply l'Hospital's rule to compute . The first-order solution for is

and the second-order accurate expression for inflation isNote that the dependence of on is purely quadratic (i.e. the zero coefficient for the linear term) around the steady state with zero inflation rate.

The presence of bifurcations is not only related to the economic model in hand, but also to the choice of the variable with respect to which the Taylor approximation is applied. This section shows that the bifurcation can be eliminated through renormalization of model variables; furthermore, renormalization may yield a more accurate first-order solution than applying l'Hospital's rule to the original formulation.

In the single-equation setting, if we can approximate the model with respect to
and
instead of
and
, then the bifurcation problem would not emerge.^{3} To see this, rewrite (2.2) as follows:

Under this renormalization, the expression for the relative price distortion is richer--and more accurate--than (3.6) derived using l'Hospital's rule. Another renormalization that produces a solution similar to (3.6) is to approximate with respect to (instead of ). This alternative way is based on the interpretation that the initial relative price distortion is of second order. Specifically, we rewrite the model as

In the multi-dimensional case, the two ways of renormalization would correspond to

According to the first renormalization, the second-order approximation of (2.1) is

and the logarithmic transformation of (2.4) is Therefore, the second-order solution of this problem would be

It is noteworthy to point out that, according to this renormalization, the first-order relationship between inflation and relative price distortions ( ) replicates the exact nonlinear relationship (2.4).

The alternative renormalization consistent with the timeless perspective is to adopt (instead of ) as an exogenous variable. Based on this choice of an expansion parameter, Woodford (2003) concluded that the optimal inflation rate is zero to the first order in the absence of cost-push shocks. Under this normalization, the two model equations are approximated as follows:

The second-order solution to this system of equations would be purely quadratic

The first-order approximation of this solution is consistent with the optimality of zero inflation, as derived in the linear-quadratic approximation by Woodford (2003), Benigno and Woodford (2005), and Levine, Pearlman and Pierse (2006). Furthermore, the second-order solution for inflation is equivalent to the one via the bifurcation method, (3.7).

After presenting two different renormalizations, it is natural to compare approximation errors for these two methods.^{6} For this purpose, we use as a
reference point the closed-form solution to the optimal policy problem (2.5). Specifically, as shown in Yun (2005), the exact nonlinear solution for the optimal inflation rate is

It is noteworthy that this closed-form solution is feasible only when the relative price distortion is the only distortion--due to the assumption that there is an optimal subsidy and there are no cost-push shocks. The optimal rate of inflation is less than zero as long as there are initial price distortions .

The difference between the two methods is that the expansion parameter of the first renormalization is
, while that of the second is
. Figure 1 compares the accuracy of the two normalizations based on the first-order solution under each normalization.^{7} The black solid line represents the exact closed-form solution for annualized inflation (
) in terms of initial relative distortion (
). The blue line with crosses is the linear approximation of this nonlinear solution. This corresponds to the first-order approximation of when the expansion parameter is
--that is,
. It is evident that this approximation is more accurate than
: the first-order approximation with
as the expansion parameter, depicted by the red circles.

We can provide an intuitive understanding about the improved accuracy of the approximation with respect to as follows. Since is the square of , the first-order approximation with respect to is equivalent to the second-order approximation with respect to :

Note that the equality holds because no linear terms are included in with zero steady-state inflation rate.We have illustrated how to detect the existence of a bifurcation and demonstrated how to apply l'Hospital's rule to characterize the solution. We have also shown that the bifurcation can be eliminated through renormalization of model variables; furthermore, renormalization may yield a more accurate first-order solution than applying l'Hospital's rule to the original formulation. This paper has focused on the consequences of renormalization on the treatment of bifurcations. However, the renormalization is also associated with the welfare evaluation of different policies as in Benigno and Woodford (2005).

"The hopf bifurcation and the existence and stability of closed orbits in multisector models of optimal economic growth." Journal of Economic Theory, December, 1979, Vol. 21(3), pp. 421-444.

" Inflation Stabilization and Welfare: The Case of A Distorted Steady State." *Journal of the European Economic Association*, December, 2005, Vol. 3(6), pp. 1185-1236.

"Staggered Prices in a Utility Maximizing Framework." *Journal of Monetary Economics*, 1983, 12 (3), pp. 383-398.

"Asymptotic Methods for Asset Market Equilibrium Analysis." *Economic Theory*, 2001, 1 (18), pp. 127-157.

"Monetary Policy under Uncertainty in Micro-Founded Macroeconometric Models" , in M. Gertler and K. Rogoff, eds., NBER Macroeconomics Annual, 2005. Cambridge, MA: MIT Press, 2006.

"Linear-Quadratic Approximation, External Habit and Targeting Rules." unpublished manuscript, University of Surrey, October, 2006.

"Optimal Simple and Implementable Monetary and Fiscal Rules." forthcoming in *Journal of Monetary Economics*,
2006.

"Optimal Monetary Policy with Relative Price Distortions." *American Economomic Review*, March, 2005, 95 (1), pp. 89 - 109.

* We have benefited from discussion with Gary Anderson, Jean Boivin, Chris Sims, and participants at the Canadian Macroeconomics Study Group meeting in Montreal. The views in this paper are solely the responsibility of the authors and should not be interpreted
as reflecting the views of the Board of Governors of the Federal Reserve System or any other person associated with the Federal Reserve System. Return to Text

† Corresponding Author: Division of Monetary Affairs, Mailstop 71, Federal Reserve Board, Washington, DC 20551. Tel: (202) 452-2981, E-mail: jinill.kim@frb.gov. Return to Text

‡ Same mailing address, E-mail: andrew.levin@frb.gov. Return to Text

§ Same mailing address, E-mail: tack.yun@frb.gov. Return to Text

1. A detailed derivation of this first-order approximation is available upon request. Return to Text

2. See Judd (1998, Ch. 15) and Judd and Guu (2001) for details on the bifurcation method. Return to Text

3. Examples include Schmitt-Grohe and Uribe (2006) and Levin et al. (2006). Return to Text

4. In perturbation analysis, the timeless perspective has implications for the choice of normalization variables. This point can be equivalently addressed in a linear-quadratic analysis. The timeless perspective is also different from the Ramsey optimal
policy in terms of how to deal with the lagged Lagrange multiplier for the behavior of models and their welfare implications. Return to Text

5. To use this solution recursively, we have to make the set of equations complete in the dynamic context by adding an equation (
) that relates the contemporaneous value of
and
. There are two ways to express the solution of this complete system. One way is to replace
with
. Alternatively, we can derive the linear dynamics of
as
. Return to Text

6. Benigno and Woodford (2005) also comment on the difference between the two methods in Footnotes 30 and 52. Return to Text

7. Other ways to measure accuracy include checking Euler equation errors and calculating welfare costs caused by approximation errors. Return to Text