Macroeconometric Analysis: Chetan Dave David N. Dejong
Macroeconometric Analysis: Chetan Dave David N. Dejong
Macroeconometric Analysis: Chetan Dave David N. Dejong
analysis. The …rst step of the preparation stage is the construction of a linear approximation
The purpose of this chapter is to demonstrate the completion of this …rst step for three
prototypical model environments that will serve as examples throughout the remainder of
the text. This will set the stage for Part II, which outlines and demonstrates alternative
The …rst environment is an example of a simple real business cycle (RBC) framework,
patterned after that of Kydland and Prescott (1982). The foundation of models in the RBC
tradition is a neoclassical growth environment, augmented with two key features: a labor-
leisure trade-o¤ that confronts decision makers; and uncertainty regarding the evolution
of technological progress. The empirical question Kydland and Prescott (1982) sought to
address was the extent to which such a model, bereft of market imperfections and featuring
fully ‡exible prices, could account for observed patterns of business cycle activity while
capturing salient features of economic growth. This question continues to serve as a central
focus of this active literature; an overview is available in the collection of papers presented
in Cooley (1995).
Viewed through the lens of an RBC model, business cycle activity is interpretable as re-
Such interpretations are not without controversy. Alternative interpretations cite the exis-
1
tence of market imperfections, costs associated with the adjustment of prices, and other
nominal and real frictions as potentially playing important roles in in‡uencing business cy-
cle behavior, and giving rise to additional sources of business cycle ‡uctuations. Initial
skepticism of this nature was voiced by Summers (1986); and the collection of papers con-
tained in Mankiw and Romer (1991) provide an overview of DSGE models that highlight the
role of, e.g., market imperfections in in‡uencing aggregate economic behavior. As a com-
plement to the RBC environment, the second environment presented here (that of Ireland,
2004) provides an example of a model within this neo-Keynesian tradition. Its empirical
purpose is to simultaneously evaluate the role of cost, demand and productivity shocks in
driving business cycle ‡uctuations. Textbook references for models within this tradition are
The realm of empirical applications pursued through the use of DSGE models extends
well beyond the study of business cycles. The third environment serves as an example of
this point: it is a model of asset-pricing behavior adopted from Lucas (1978). The model
consumption in the face of exogenous stochastic movements in income and dividends earned
from asset holdings. Viewed through the lens of this model, two particular features of asset-
pricing behavior have proven exceptionally di¢ cult to explain. First, Shiller (1981) used a
version of the model to underscore the puzzling volatility of prices associated with broad
indexes of assets (such as the Standard & Poor’s 500), highlighting what has come to be
known as the ‘volatility puzzle’. Second, Mehra and Prescott (1985) used a version of the
model to highlight the puzzling dual phenomenon of a large gap observed between aggregate
returns on risky and risk-free assets. When coupled with exceptionally low returns yielded
2
by risk-free assets, this feature came to be known as the ‘equity premium’puzzle. The texts
of Shiller (1989) and Cochrane (2001) provide overviews of literatures devoted to analyses
of these puzzles.
1.1 Environment
is to maximize U , the expected discounted ‡ow of utility arising from chosen streams of
The household is equipped with a production technology that can be used to produce a
yt = zt f (kt ; nt ); (2)
where kt and nt denote quantities of physical capital and labor assigned by the household
to the production process, and zt denotes a random disturbance to the productivity of these
3
inputs to production (that is, a productivity or technology shock).
Within a period, the household has one unit of time available for division between labor
1 = nt + lt : (3)
In addition, output generated at time t can either be consumed or used to augment the stock
of physical capital available for use in the production process in period t + 1. That is, output
y t = c t + it ; (4)
where it denotes the quantity of investment. Finally, the stock of physical capital evolves
according to
where denotes the depreciation rate. The household’s problem is to maximize (1) subject
Implicit in the speci…cation of the household’s problem are two sets of trade-o¤s. One is a
consumption/savings trade-o¤: from (4), higher consumption today implies lower investment
(savings), and thus from (5), less capital available for production tomorrow. The other is
a labor/leisure trade-o¤: from (3), higher leisure today implies lower labor today and thus
explicit functional forms for u( ) and f ( ), and to characterize the stochastic behavior of
4
the productivity shock zt . We pause before doing so to make some general comments. As
noted, an explicit goal of the RBC literature is to begin with a model speci…ed to capture
important characteristics of economic growth, and then to judge the ability of the model to
capture key components of business cycle activity. From the model builder’s perspective,
the former requirement serves as a constraint on choices regarding the speci…cations for
u( ); f ( ) and the stochastic process of zt . Three key aspects of economic growth serve as
constraints in this context: over long time horizons the growth rates of fct ; it ; yt ; kt g are
roughly equal (balanced growth), the marginal productivity of capital and labor (re‡ected
by relative factor payments) are roughly constant over time, and flt ; nt g show no tendencies
Beyond satisfying this constraint, functional forms chosen for u( ) are typically strictly
increasing in both arguments, twice continuously di¤erentiable, strictly concave and satisfy
@u(ct ; lt ) @u(ct ; lt )
lim = lim = 1: (6)
c!0 @ct l!0 @lt
Functional forms chosen for f ( ) typically feature constant returns to scale and satisfy similar
limit conditions.
Finally, we note that the inclusion of a single source of uncertainty in this framework,
via the productivity shock zt , implies that the model carries nontrivial implications for
the stochastic behavior of a single corresponding observable variable. For the purposes of
this chapter, this limitation is not important; however, it will motivate the introduction of
5
1.1.1 Functional Forms
The functional forms presented here enjoy prominent roles in the macroeconomics lit-
erature. Instantaneous utility is of the constant relative risk aversion (CRRA) form:
1
c't lt1 '
u(ct ; lt ) = ; (7)
1
where > 0 determines two attributes: it is the coe¢ cient of relative risk aversion, and also
1
determines the intertemporal elasticity of substitution, given by (for textbook discussions,
see e.g., Blanchard and Fischer, 1998; or Romer, 2001).1 Note that the larger is , the
Also, ' 2 (0; 1) indicates the importance of consumption relative to leisure in determining
instantaneous utility.
yt = zt kt n1t ; (8)
where 2 (0; 1) represents capital’s share of output. Finally, the log of the technology shock
2
"t N ID(0; ); 2 ( 1; 1): (10)
1
When = 1; u(:) = log(:):
6
The solution to the household’s problem may be obtained via standard application of the
theory of dynamic programming (e.g., as described in Stokey and Lucas, 1989). Necessary
conditions associated with the household’s problem expressed in general terms are given by
The intratemporal optimality condition (11) equates the marginal bene…t of an additional
unit of leisure time with its opportunity cost: the marginal value of the foregone output
resulting from the corresponding reduction in labor time. The intertemporal optimality con-
dition (12) equates the marginal bene…t of an additional unit of consumption today with its
opportunity cost: the discounted expected value of the additional utility tomorrow that the
corresponding reduction in savings would have generated (higher output plus undepreciated
capital).
Consider the qualitative implications of (11) and (12) for the impact of a positive produc-
tivity shock on the household’s labor/leisure and consumption/savings decisions. From (11),
higher labor productivity implies a higher opportunity cost of leisure, prompting a reduction
in leisure time in favor of labor time. From (12), the curvature in the household’s utility
serves to increase output, thus a¤ording an increase in consumption; but since the marginal
utility of consumption is decreasing in consumption, this drives down the opportunity cost of
savings. The greater is the curvature of u(:); the more intense is the consumption-smoothing
objective, and thus the greater will be the intertemporal reallocation of resources in the face
7
of a productivity shock.
Dividing (11) by the expression for the marginal utility of consumption, and employing
1 ' ct kt
= (1 )zt (13)
' lt nt
( " #)
1
'(1 ) 1 (1 ')(1 ) '(1 ) 1 (1 ')(1 ) nt+1
ct lt = Et ct+1 lt+1 zt+1 + (1 ) : (14)
kt+1
1 ' ct kt
= (1 )zt (15)
' lt nt
( " #)
1
nt+1
ct lt = Et ct+1 lt+1 zt+1 + (1 ) (16)
kt+1
yt = zt kt nt1 (17)
y t = c t + it (18)
1 = nt + lt (20)
may be computed analytically from this system. These are derived by holding zt to its steady
8
state value z, which we set to 1:
y c i 1 k
= ; = ; = ; n= ; l=1 n; = ;
n n n 1+ 1 1 '
1 1 n
1 '
(22)
where
1
1
=
1= 1+
= :
Note that in steady state the variables fyt ; ct ; it ; kt g do not grow over time. Implicitly,
these variables are represented in the model in terms of deviations from trend, and steady
state values indicate the relative heights of trend lines. To incorporate growth explicitly,
!t = !t 1 + "t : (24)
Note that, absent shocks, the growth rate of zt is given by g; and that removal of the
trend component (1 + g)t from zt yields the speci…cation for log zt given by (21). Further,
the reader is invited to verify that under this speci…cation for zt , fct ; it ; yt ; kt g will have a
g
common growth rate given by 1
:Thus the model is consistent with the balanced-growth
requirement, and as speci…ed, all variables are interpreted as being measured in terms of
9
There is one subtlety associated with the issue of trend removal that arises in dealing with
the dynamic equations of the system. Consider the law of motion for capital (19). Trend
g t
removal here involves division of both sides by 1 + 1
; however, the trend component
g t+1
associated with kt+1 is 1 + 1
; so the speci…cation in terms of detrended variables is
g
1+ kt+1 = it + (1 )kt : (25)
1
Likewise, there will be a residual trend factor associated with ct+1 in the intertemporal
optimality condition (16). Since ct+1 is raised to the power = '(1 ) 1; the residual
g
factor is given by 1 + 1
:
( " #)
1
g nt+1
ct lt = Et 1+ ct+1 lt+1 zt+1 + (1 ) : (26)
1 kt+1
1
With negative (insured by < 1; i.e., an inelastic intertemporal elasticity of substitution
towards t:
Exercise 1 Rederive the steady state expressions (22) by replacing (19) with (25), and (16)
with (26). Interpret the intuition behind the impact of g on the expressions you derive.
1.3 Linearization
The linearization step involves taking a log-linear approximation of the model at steady
state values. In this case, the objective is to map (15)-(21) into the linearized system
Axt+1 = Bxt + C t+1 +D t+1 for eventual empirical evaluation. Regarding D, dropping
10
Et from the Euler equation (16) introduces an expectations error in the model’s second
Regarding A and B, using the solution methodology discussed in Chapter 2, these can
1 '
0 = log( ) + log c0 log l0 log(1 ) log z 0 log k + log n0 (27)
'
exp [(1 ) log n0 ]
0= log c + log l log log c0 log l0 log exp(log z 0 ) + (1 )
exp [(1 ) log k 0 ]
(28)
0 = log y 0 log fexp [log (c0 )] + exp [log (i0 )]g (30)
The mapping from (15)-(21) to (27)-(33) involves four steps. First, logs of both sides of each
equation are taken; second, all variables not converted into logs in the …rst step are converted
using the fact, e.g., that y = exp(log(y)); third, all terms are collected on the right-hand side
of each equation; fourth, all equations are multiplied by 1: Derivatives taken with respect
to log y 0 ; etc. evaluated at steady state values yield A; and derivatives taken with respect
to log y; etc. yield B: Note that capital installed at time t is not productive until period
11
t + 1; thus, e.g., k rather than k 0 appears in (29).
Having obtained A; B; C and D, the system can be solved using any of the solution
methods outlined in Chapter 2 to obtain a system of the form xt+1 = F ( )xt + et+1 . This
system can then be evaluated empirically using any of the methods described in Part II.
0
yt ct it nt lt kt zt
xt = log ; log ; log ; log ; log ; log ; log
y c i n l k z
and
show that the steady state values of the model are y = 0:9; c = 0:7; i = 0:2; n =
0:47; l = 0:53 and k = 3:5 (and take as granted z = 1): Next, use a numerical gradient
procedure to derive
2 3
6 0 1 0 0:33 1 0 1 7
6 7
6 7
6 0 1:5 0 0:12 0:5 0 0:17 7
6 7
6 7
6 7
6 1 0 0 0:67 0 0 1 7
6 7
6 7
6 7
A=6
6 1 0:77 0:23 0 0 0 0 7 7
6 7
6 7
6 7
6 0 0 0:18 0 0 1 0 7
6 7
6 7
6 7
6 0 0 0 0:47 0:53 0 0 7
6 7
4 5
0 0 0 0 0 0 1
12
2 3
6 0 0 0 0 0 0:33
0 7
6 7
6 7
6 0 1:5 0 0 0:5 0:9 0 7
6 7
6 7
6 7
6 0 0 0 0 0 0:33 0 7
6 7
6 7
6 7
B=6
6 0 0 0 0 0 0 0 77:
6 7
6 7
6 0 0 0 0 0 0:77 0 7
6 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 7
6 7
4 5
0 0 0 0 0 0 0:9
Exercise 3 Rederive the matrices A and B given the explicit incorporation of growth in the
model. That is, derive A and B using the steady state expressions obtained in Exercise 1,
Policy
This section outlines a model of imperfect competition featuring ‘sticky’ prices. The
model includes three sources of aggregate uncertainty: shocks to demand, technology and the
competitive structure of the economy. The model is due to Ireland (2004), who designed it to
determine how the apparent role of technology shocks in driving business-cycle ‡uctuations
From a pedological perspective, the model di¤ers in two interesting ways relative to the
RBC model outlined above. While the linearized RBC model is a …rst-order system of
di¤erence equations, the linearized version of this model is a second-order system. However,
13
as we demonstrate, it is possible to represent a system of arbitrary order into the …rst-
order form taken by Axt+1 = Bxt + C t+1 +D t+1 , given appropriate speci…cation of the
into the behavior of non-stationary data is revisited from an alternative perspective than that
adopted in the discussion of the RBC model. Rather than assuming the actual data follow
stationary deviations around deterministic trends, here the data are modelled as following
drifting random walks; stationarity is induced via di¤erencing rather than detrending.2
2.1 Environment
The economy once again consists of a continuum of identical households. Here, there
are two distinct production sectors: an intermediate-goods sector and a …nal-goods sector.
di¤erentiated products which serve as factors of production in the …nal-goods sector. While
…rms in this sector have the ability to set prices, they face a friction in doing so. Finally,
2
Nelson and Plosser (1982) initiated an intense debate regarding the nature of non-stationarity that
characterizes macroeconomic variables. The issue has proven di¢ cult to resolve, thus the use of alterna-
tive detrending methods remains common practice. For an overview of the debate, see, e.g., DeJong and
Whiteman (1993) and Stock (1994).
14
2.1.1 Households
The representative household maximizes lifetime utility de…ned over consumption, money
( )
X
1
mt nt
t
max U = E0 at log ct + log (34)
ct ;mt ;nt
0
pt
bt
s:t: pt c t + + m t = m t 1 + bt 1 + t + wt nt + dt ; (35)
rt
where 2 (0; 1) and 1: According to the budget constraint (35), the household divides
its wealth between holdings of bonds bt and money mt ; bonds mature at the gross nominal
rate rt between time periods. The household also receives transfers t from the monetary
authority and works nt hours in order to earn wages wt to …nance its expenditures. Finally,
the household owns an intermediate-goods …rm, from which it receives a dividend payment
dt . Note from (34) that the household is subject to an exogenous demand shock at that
Recognizing that the instantaneous marginal utility derived from consumption is given by
at
ct
; the …rst-order conditions associated with the household’s choices of labor, bond holdings
wt at 1
= nt (36)
pt ct
1 at+1 1 at
Et = (37)
pt+1 ct+1 rt p t ct
1
mt 1 at+1 1 at
+ Et = : (38)
pt pt+1 ct+1 pt ct
15
Exercise 4 Interpret how (36)-(38) represent the optimal balancing of trade-o¤s associated
2.1.2 Firms
There are two types of …rms, one that produces a …nal consumption good yt ; which sells
at price pt , and a continuum of intermediate-goods …rms that supply inputs to the …nal-good
…rm. The output of the ith intermediate good is given by yit ; which sells at price pit : The
intermediate goods combine to produce the …nal good via a constant elasticity of substitution
(CES) production function. The …nal-good …rm operates in a competitive environment and
Z1
F
max t = pt yt pit yit di (39)
yit
0
81 9 t
<Z t 1 = t 1
s:t: yt = yit t
di : (40)
: ;
0
The solution to this problem yields a standard demand for intermediate inputs and a price
aggregator:
t
pit
yit = yt (41)
pt
81 91 1
<Z = t
pt = p1it t
di : (42)
: ;
0
Notice that t is the markup of price above marginal cost; randomness in t provides the
16
Intermediate-goods …rms are monopolistically competitive. Since the output of each …rm
enters the …nal-good production function symmetrically, the focus is on a representative …rm.
The …rm is owned by the representative household, thus its objectives are aligned with the
household’s. It manipulates the sales price of its good in pursuit of these objectives, subject
X
1
at dt
I t
max it = E0 ; (43)
pit
0
ct pt
s:t: yit = zt nit (44)
t
pit
yit = yt (45)
pt
2
pit
c(pit ; pit 1 ) = 1 yt ; > 0; (46)
2 pit 1
where is the gross in‡ation rate targeted by the monetary authority (described below),
dt pit yit wt nt
= c(pit ; pit 1 ) : (47)
pt pt
t t 1
pit yt pit wt yt 1
( t 1) = t
pt pt pt pt zt pt
pit yt at+1 ct pit+1 yt+1 pit+1
1 Et 1 : (48)
pit 1 pit 1 at ct+1 pit p2it
The left-hand side of (48) re‡ects the marginal revenue to the …rm generated by an increase
in price; the right-hand side re‡ects associated marginal costs. Under perfect price ‡exibility
17
( = 0) there is no dynamic component to the …rm’s problem; the price-setting rule reduces
wt wt
to pit = t
t
1 zt
; which is a standard markup of price over marginal cost zt
. Under ‘sticky
prices’( > 0) the marginal cost of an increase in price has two additional components: the
direct cost of a price adjustment, and an expected discounted cost of a price change adjusted
by the marginal utility to the households of conducting such a change. Empirically, the
The monetary authority chooses the nominal interest rate according to a Taylor Rule.
With all variables expressed in terms of logged deviations from steady state values, the rule
is given by
2
ret = et 1
rr + et + et
gg + et
oo + "rt ; "rt IIDN (0; r ); (49)
where et is the gross in‡ation rate, get is the gross growth rate of output, and oet is the output
gap (de…ned below). The i parameters denote elasticities. The inclusion of ret 1 as an input
into the Taylor Rule allows for the gradual adjustment of policy to demand and technology
The output gap is the logarithm of the ratio of actual output yt to capacity output ybt .
Capacity output is de…ned to be the ‘e¢ cient’ level of output, which is equivalent to the
18
level of output chosen by a benevolent social planner who solves:
8 0 19
X
1 >
< Z1 >
=
1@
max US = E0 t
at log ybt nit diA (50)
ybt ;nit
0
>
: >
;
0
0 1 t
Z1 t 1
t 1
1
ybt = at zt : (52)
In addition to the monetary policy shock "rt introduced in (49), the model features a
demand shock at , a technology shock zt , and a cost-push shock t. The former is IID; the
drifting random walk. This induces similar behavior in the model’s endogenous variables,
and necessitates the use of an alternative to the detrending method discussed above in the
context of the RBC model. Here, stationarity is induced by normalizing model variables
19
by zt . For the corresponding observable variables, stationarity is induced by di¤erencing
rather than detrending: the observables are measured as deviations of growth rates (logged
di¤erences of levels) from sample averages. Details are provided in the linearization step
discussed below.
The model is closed through two additional steps. The …rst is the imposition of symmetry
The second is the requirement that the money and bond markets clear:
mt = mt 1 + t (57)
bt = bt 1 = 0: (58)
In its current form, the model consists of twelve equations: the household’s …rst-order
conditions and budget constraint; the aggregate production function; the aggregate real div-
idends paid to the household by its intermediate-goods …rm; the intermediate-goods …rm’s
…rst-order condition; the stochastic speci…cations for the structural shocks; and the expres-
sion for capacity output. Ireland’s (2004) empirical implementation focused on a linearized
Rule (speci…ed in linearized form in (49)); the three exogenous shock speci…cations; and
de…nitions for the growth rate of output and the output gap.
20
The reduced system is recast in terms of the following normalized variables:
:: yt :: ct :: ybt pt
yt = ; ct = ; ybt = ; t =;
zt zt zt pt 1
:: (dt =pt ) :: (wt =pt ) :: (mt =pt ) :: zt
dt = ; wt = ; mt = ; zt = :
zt zt zt zt 1
Using the expression for real dividends given by (47), the household’s budget constraint is
rewritten as
:: :: 2 ::
t
y t = ct + 1 yt: (59)
2
at at+1 1 1
:: = rt Et :: :: : (60)
ct ct+1 z t+1 t+1
Next, the household’s remaining …rst-order conditions, the expression for the real dividend
payment (47) it receives, and the aggregate production function can be combined to elim-
inate wages, money, labor, dividends and capacity output from the system. This serves to
introduce the expression for the output gap into the system:
::
yt y
ot = 1t : (61)
ybt
at
Finally, normalizing the …rst-order condition of the intermediate-goods …rm and the stochas-
21
tic speci…cations leads to the following non-linear system:
:: :: 2 ::
t
y t = ct + 1 yt (62)
2
at at+1 1 1
:: = rt Et :: :: (63)
ct ct+1 z t+1 t+1
:: :: ::
ct :: 1 t t ct at+1 t+1 t+1 y t+1
0 = 1 t+ t y 1 + Et :: 1 :: (64)
at t ct+1 at yt
:: ::
zty
gt = :: t (65)
yt 1
::
yt yt
ot = = 1 (66)
ybt
at
log(at ) = (1 a ) log(a) + a log(at 1 ) + "at (67)
::
log(z t ) = log(z) + "zt (69)
2.3 Linearization
Log-linearization proceeds with the calculation of steady state values of the endogenous
variables:
1 1
z 1 1
r= ; c=y= a ; o= ; (70)
(62)-(69) are then log-linearized around these values. As with Model I, this can be accom-
plished through the use of a numerical gradient procedure. However, as an alternative to this
approach, here we follow Ireland (2004) and demonstrate the use of a more analytically ori-
ented procedure. In the process, it helps to be mindful of the re-con…guration Ireland worked
22
with: an IS curve; a Phillips curve; the Taylor Rule; the shock processes; and de…nitions of
As a …rst step, the variables appearing in (62)-(69) are written in logged form. Log-
::
yt
linearization of (62) then yields yet log y
=e
ct , since the partial derivative of yet with
respect to et (evaluated at steady state) is zero.3 Hence upon linearization, this equation is
Relating output and the output gap via the log-linearization of (66),
1
yet = eat + oet ; (72)
the term Et yet+1 yet may be substituted out of (71), yielding the IS curve:
1
oet = Et oet+1 (rt Et et+1 ) + 1 (1 a )e
at : (73)
Similarly, log-linearizing (64) and eliminating yet using (72) yields the Phillips curve:
( 1) 1e
where = and eet = t. This latter equality is a normalization of the cost-push
3
Recall our notational convention: tildes denote logged deviations of variables from steady state values.
23
shock; like the cost-push shock itself, the normalized shock follows an AR(1) process with
1
persistence parameter = e; and innovation standard deviation e = .
The resulting IS and Phillips curves are forward looking: they include the one-step-ahead
equations to include lagged variables of the output gap and in‡ation in order to enhance the
empirical coherence of the model. This …nal step yields the system he analyzed. Dropping
oe = e
oo + (1 e0
o )Et o (r Et e0 ) + 1 1
(1 a )e
a (75)
e = e + (1 )Et e0 + oe ee (76)
1 0
oe0 = ye0 e
a (78)
re0 = e+
rr e0 + e0
gg + e0
oo + "0r (79)
a0 =
e ae
a + "0a (80)
ee0 = e+
ee "0e (81)
where the structural shocks t = f"rt ; "at ; "et ; "zt g are IIDN with diagonal covariance matrix
. The additional parameters introduced are o 2 [0; 1] and 2 [0; 1]; setting o = =0
The augmentation of the IS and Phillips curves with lagged values of the output gap and
in‡ation converts the model from a …rst- to a second-order system. Thus a …nal step is re-
24
quired in mapping this system into the …rst-order speci…cation Axt+1 = Bxt +C t+1 +D t+1 :
This is accomplished by augmenting the vector xt to include not only contemporaneous ob-
servations of the variables of the system, but also to include lagged values of the output gap
and in‡ation:
xt [e
ot oet 1 et et 1 yet ret get e
at eet zet ]0 :
This also requires the introduction of two additional equations into the system: e0 = e0 and
oe0 = oe0 . Specifying these as the …nal two equations of the system, the corresponding matrices
2 3
6 (1 0) 1 1 0 0 0 0 0 0 0 7
6 7
6 7
6 0 (1 ) 1 0 0 0 0 0 0 7
6 7
6 7
6 7
6 0 0 0 0 1 0 1 0 0 1 7
6 7
6 7
6 7
6 1 0 0 0 1 0 0 1
0 0 7
6 7
6 7
6 7
6 0 0 0 1 0 0 0 7
6 0 g 7
A=6
6
7
7 (83)
6 0 0 0 0 0 0 0 1 0 0 7
6 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 0 1 0 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 0 0 1 7
6 7
6 7
6 7
6 0 0 0 1 0 0 0 0 0 0 7
6 7
4 5
0 1 0 0 0 0 0 0 0 0
25
2 3
1
6 0 o 0 0 0 1 0 (1 )(1 a) a 0 0 7
6 7
6 7
6 0 0 0 0 0 0 0 1 0 7
6 7
6 7
6 7
6 0 0 0 0 1 0 0 0 0 0 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 0 0 0 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 0 0 7
6 r 7
B=6
6
7:
7 (84)
6 0 0 0 0 0 0 0 0 0 7
6 a 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 0 e 0 7
6 7
6 7
6 7
6 0 0 0 0 0 0 0 0 0 1 7
6 7
6 7
6 7
6 0 0 1 0 0 0 0 0 0 0 7
6 7
4 5
1 0 0 0 0 0 0 0 0 0
0
Further, de…ning t =[ 1t 2t ] ; where 1t+1 = Et oet+1 oet+1 and 2t+1 = Et et+1 et+1 ;
2 3
6 04x4 7
6 7
6 7
6 1 0 0 0 7 2 3
6 7
6 7
6 7 6 1 0 1 7
6 0 1 0 0 7 6 7
6 7 6 7
C=6
6
7;
7 D=6
6 0 (1 ) 7
7 (85)
6 0 0 1 0 7 6 7
6 7 4 5
6 7 02x8
6 7
6 7
6 0 0 0 1 7
6 7
4 5
02x4
The …nal step needed for empirical implementation is to identify the observable variables
of the system. For Ireland, these are the gross growth rate of output gt ; the gross in‡ation
rate t; and the nominal interest rate rt (all measured as logged ratios of sample averages).
26
Under the assumption that output and aggregate prices follow drifting random walks, gt and
t are stationary; the additional assumption of stationarity for rt is all that is necessary to
Construct vectors (xt ; et ) and matrices ( ; ) so that the model may be re-cast in …rst-order
form as
xt+1 = xt + et+1 :
Exercise 6 Solve the linearized system (75)-(82) using any of the methods outlined in Chap-
ter 2. Note that the vector of deep parameters is now given by:
0
= [z ! x r g x a a z r] :
Exercise 7 Consider the following CRRA form for the instantaneous utility function for
Model II:
mt c mt nt
u(ct ; ; nt ) = at t + log :
pt pt
27
3 Model III: Asset Pricing
ior. Alternative versions of the model have played a prominent role in two important strands
of the empirical …nance literature. The …rst, launched by Shiller (1981) in the context of a
single-asset version of the model, concerns the puzzling degree of volatility exhibited by prices
associated with aggregate stock indexes. The second, launched by Mehra and Prescott (1985)
in the context of a multi-asset version of the model, concerns the puzzling coincidence of a
large gap observed between the returns of risky and risk-free assets, and a low average risk-
free return. Resolutions to both puzzles have been investigated using alternate preference
speci…cations. After outlining single- and multi-asset versions of the model given a generic
role of preferences in the equity-premium literature are provided by Kocherlakota (1996) and
Cochrane (2001); and in the stock-price volatility literature by DeJong and Ripoll (2004).
The model features a continuum of identical households and a single risky asset. Shares
held during period (t 1), st 1 ; yield a dividend payment dt at time t; time-t share prices
an exogenous stochastic dividend stream, proceeds from sales of shares, and an exogenous
is given by
X
1
t
max U = E0 u(ct ); (86)
ct
t=0
28
where 2 (0; 1) again denotes the discount rate, and optimization is subject to
ct + pt (st st 1 ) = dt st 1 + qt : (87)
Since households are identical, equilibrium requires st = st 1 for all t, and thus ct = dt st +qt =
u0 (dt+1 + qt+1 )
pt = Et (dt+1 + pt+1 ): (88)
u0 (dt + qt )
From (88), following a shock to either dt or qt , the response of pt depends in part upon
the variation of the marginal rate of substitution between t and t + 1. This in turn depends
upon the instantaneous utility function u( ). The puzzle identi…ed by Shiller (1981) is that
pt is far more volatile than what (88) would imply, given the observed volatility of dt :
The model is closed by specifying stochastic processes for (dt ; qt ). These are given by
29
3.2 Multi-Asset Environment
An n-asset extension of the environment leaves the household’s objective function intact,
but modi…es its budget constraint to incorporate the potential for holding n assets. As a
special case, Mehra and Prescott (1985) studied a two-asset speci…cation, including a risk-
free asset (ownership of government bonds) and risky asset (ownership of equity). In this
where pet denotes the price of the risky asset, set represents the number of shares held in the
asset during period t 1; and pft and sft are analogous for the risk-free asset. The risk-free
asset pays one unit of the consumption good at time t if held at time t 1 (hence the loading
First-order conditions associated with the choice of the assets are analogous to the pricing
u0 (ct+1 ) pet+1 + dt
Et = 1 (93)
u0 (ct ) pet
u0 (ct+1 ) 1
Et 0 = 1: (94)
u (ct ) pft
e pet+1 + dt
rt+1 =
pet
f 1
rt+1 = f;
pt
30
Mehra and Prescott’s identi…cation of the equity premium puzzle centers on
u0 (ct+1 ) f
Et r = 1 (95)
u0 (ct ) t+1
u0 (ct+1 ) h e f
i
Et 0 rt+1 rt+1 = 0; (96)
u (ct )
where (96) is derived by subtracting (94) from (93). The equity premium puzzle has two
components. First, taking fct g as given, the average value of re rf is quite large: given
CRRA preferences, implausibly large values of the risk-aversion parameter are needed to
account for the average di¤erence observed in returns. Second, given a speci…cation of u(c)
that accounts for (96), and again taking prices as given, the average value observed for rf is
far too low to reconcile with (95). This second component is the risk-free rate puzzle.
As noted, alternative preference speci…cations have been considered for their potential in
resolving both puzzles. Here, in the context of the single-asset environment, three forms for
the instantaneous utility function are presented in anticipation of the empirical applications
self control preferences. The presentation follows that of DeJong and Ripoll (2004), who
sought to evaluate empirically the ability of these preference speci…cations to make headway
31
3.3.1 CRRA
c1t
u(ct ) = ; (97)
1
thus > 0 measures the degree of relative risk aversion, and 1= the intertemporal elasticity
(dt+1 + qt+1 )
p t = Et (dt+1 + pt+1 ): (98)
(dt + qt )
Notice that, ceteris paribus, a relatively large value of will increase the volatility of price
responses to exogenous shocks, at the cost of decreasing the correlation between pt and dt
(due to the heightened role assigned to qt in driving price ‡uctuations). Since fdt g and fqt g
are exogenous, their steady states d and q are simply parameters. Normalizing d to 1 and
q
de…ning = d
, so that = q; the steady state value of consumption (derived from the
p= d= : (99)
1 1
Letting = 1=(1+r), where r denotes the household’s discount rate, (99) implies p=d = 1=r.
Thus as the household’s discount rate increases, its asset demand decreases, driving down
the steady state price level. Empirically, the average price/dividend ratio observed in the
32
Exercise 8 Linearize the pricing equation (98) around the model’s steady state values.
3.3.2 Habit/Durability
Following Ferson and Constantinides (1991) and Heaton (1995), an alternative speci…-
cation of preferences that introduces habit and durability into the speci…cation of preferences
is parameterized as
ht1
u(ht ) = ; (100)
1
with
where 2 (0; 1), hdt is the household’s durability stock, and hht its habit stock. The stocks
X1 X1 X
1
j d j i
hht = (1 ) ht 1 j = (1 ) ct 1 i (103)
j=0 j=0 i=0
where 2 (0; 1) and 2 (0; 1). Thus the durability stock represents the ‡ow of services
from past consumption, which depreciates at rate . This parameter also represents the
as a weighted average of the durability stock, where the weights sum to one. Notice that
more recent durability stocks, or more recent ‡ows of consumption, are weighted relatively
The variable ht represents the current level of durable services net of the average of past
services; the parameter measures the fraction of the average of past services that is netted
33
out. Notice that if = 0, there would only be habit persistence, while if = 0 only durability
survives. Finally, when = 0, the habit stock includes only one lag. Thus estimates of these
" j 1
#
X
1 X X
1
j i j i 1
ht = ct + (1 ) ct j j ct j ; (104)
j=1 i=0 j=0
P
1
j P
1
j i ct+1+j i
j=0 i=0
p t = Et (dt+1 + pt+1 ) ; (105)
P
1
j P
1
j i ct+j i
j=0 i=0
To see how the presence of habit and durability can potentially in‡uence the volatility of
When there is a positive shock to say qt , ct increases by the amount of the shock, say
2
q. Given (89)-(90), ct+1 would increase by q q, ct+2 would increase by q q, etc. Now,
examine the …rst term in parenthesis both in the numerator and the denominator. First,
q + 1 q 7 q. Thus, whether the share price pt increases by more than in the standard
CRRA case depends ultimately on whether q + 1 7 1. Notice that if j = 0 for j > 0, the
34
equation above reduces to the standard CRRA utility case. If we had only habit and not
durability, then 1 < 0, and thus the response of prices would be greater than in the CRRA
which strengthens the smoothing motive relative to the time-separable CRRA case.
Alternatively, if there was only durability and not habit, then 0 < 1 < 1, but one still
would not know whether + 1 7 1. Thus with only durability, we cannot judge how the
volatility of pt would be a¤ected: this will depend upon the sizes of and 1. Finally, we
also face indeterminacy under a combination of both durability and habit: if is large and
is small enough to make + 1 < 1, then we would get increased price volatility. Thus
this issue is fundamentally quantitative. Finally, with respect to the steady state price, note
Exercise 9 Given that the pricing equation under Habit/Durability involves an in…nite num-
ber of lags, truncate the lags to 3 and linearize the pricing equation (106) around its steady
state.
Consider next a household that every period faces a temptation to consume all of its
wealth. Resisting this temptation imposes a self-control utility cost. To model these prefer-
ences we follow Gul and Pesendorfer (2004), who identi…ed a class of dynamic self-control
preferences. In this case, the problem of the household can be formulated recursively as
W (s; P ) = max
0
fu(c) + v(c) + EW (s0 ; P 0 )g max
0
v(e
c); (107)
s se
35
where P = (p; d; e); u(:) and v(:) are Von Neuman-Morgenstern utility functions; 2 (0; 1);
c represents temptation consumption; and s0 denotes share holdings next period. While u(:)
e
is the momentary utility function, v(:) represents temptation. The problem is subject to the
c = ds + e p(s0 s) (108)
e
c = ds + e s0
p(e s): (109)
given that the agent has chosen c. With v(c) speci…ed as strictly increasing, the solution for
max
0
c) is simply to drive e
v(e c to the maximum allowed by the constraint e
c = ds + e s0
p(e s),
se
W (s; P ) = max
0
fu(c) + v(c) + EW (s0 ; P 0 )g v(ds + e + ps) (110)
s
subject to
and since
36
the optimality condition becomes
Notice that when v(:) = 0, there is no temptation, and the pricing equation reduces to
the standard case. Otherwise, the term u0 (d0 + e0 ) + v 0 (d0 + e0 ) v 0 (d0 + e0 + p0 ) represents
tomorrow’s utility bene…t from saving today. This corresponds to the standard marginal
utility of wealth tomorrow u0 (d0 + e0 ), plus the term v 0 (d0 + e0 ) v 0 (d0 + e0 + p0 ) which
represents the derivative of the utility cost of self-control with respect to wealth.
DeJong and Ripoll (2004) assume the following functional forms for the momentary and
c
v(c) = ; (117)
this preference speci…cation for the stock-price volatility issue. To understand why, rewrite
37
(118) as
2 3
(d0 +e0 )
(d+e)
+ (d + e) (d0 + e0 ) 1
(d0 + e0 + p0 ) 1
p = E [d0 + p0 ] 4 1+
5: (119)
1 + (d + e)
Suppose > 1; so that v(:) is convex, and consider the impact on p of a positive endowment
(d + e) (d0 + e0 ) 1
(d0 + e0 + p0 ) 1
in the numerator. Both e¤ects imply that relative to the CRRA case, in which = 0, this
speci…cation reduces price volatility in the face of an endowment shock, which is precisely
the opposite of what one would like to achieve in seeking to resolve the stock-price volatility
puzzle.
The mechanism behind this reduction in price volatility is as follows: a positive shock to d
or e increases the household’s wealth today, which has three e¤ects. The …rst (“smoothing”)
captures the standard intertemporal motive: the household would like to increase saving,
which drives up the share price. Second, there is a “temptation” e¤ect: with more wealth
today, the feasible budget set for the household increases, which represents more temptation
to consume, and less willingness to save. This e¤ect works opposite to the …rst, and reduces
price volatility with respect to the standard case. Third, there is the “self-control”e¤ect: due
to the assumed convexity of v(:), marginal self-control costs also increase, which reinforces
the second e¤ect. As shown above, the last two e¤ects dominate the …rst, and thus under
38
In contrast, price volatility would not necessarily be reduced if v(:) is concave, and thus
(d + e) (d0 + e0 ) 1
(d0 + e0 + p0 ) 1
increases. On the other hand, if 1+ > 0, i.e., if the risk-aversion parameter > 1,
the denominator also increases. If the increase in the numerator dominates that in the
denominator, then higher price volatility can be observed than in the CRRA case.
To understand this e¤ect, note that the derivative of the utility cost of self-control with
respect to wealth is positive if v(:) is concave: v 0 (d0 + e0 ) v 0 (d0 + e0 + p0 ) > 0. This means
that as agents get wealthier, self-control costs become lower. This explains why it might
be possible to get higher price volatility in this case. The mechanism behind this result
still involves the three e¤ects discussed above: smoothing, temptation, and self-control. The
di¤erence is on the latter e¤ect: under concavity, self-control costs are decreasing in wealth.
This gives the agent an incentive to save more rather than less. If this self-control e¤ect
dominates the temptation e¤ect, then these preferences will produce higher price volatility.
Notice that when v(:) is concave, conditions need to be imposed to guarantee that W (:)
is strictly concave, so that the solution corresponds to a maximum (e.g., see Stokey and
1 2 2
(d + e) + ( 1) (d + e) (d + e + p) <0 (120)
which holds for any d, e, and p > 0, and for > 0, > 0, and 0 < < 1. The empirical
39
implementation in Part II of the text proceeds under this set of parameter restrictions.
Finally, from the optimality conditions under self-control preferences, steady-state temp-
tation consumption is e
c = 1 + + p. From (118), the steady-state price in this case is given
by
" #
1 1
(1 + ) + (1 + ) (1 + + p)
p= (1 + p) 1
: (121)
(1 + ) + (1 + )
Regarding (121), the left-hand-side is a 45-degree line. The right-hand side is strictly con-
cave in p, has a positive intercept, and a positive slope that is less than one at the intercept.
Thus (121) yields a unique positive solution for p for any admissible parameterization of
the model. (In practice, (121) can be solved numerically, e.g., using GAUSS’s quasi-Newton
algorithm NLSYS; see Judd, 1998, for a presentation of alternative solution algorithms.) An
increase in causes the function of p on the right-hand-side of (121) to shift down and ‡at-
ten, thus p is decreasing in . The intuition for this is again straightforward: an increase in
This drives down its demand for asset shares, and thus p. Note the parallel between this
Linearize the asset-pricing equation (119) using the steady state values for p; d; q implied
40
References
[1] Benassy, J. 2002. The Macroeconomics of Imperfect Competition and Nonclearing Mar-
Press.
[3] Clarida, R., J. Gali and M. Gertler. 2000. Monetary policy rules and macroeconomic
stability: evidence and some theory. Quarterly Journal of Economics 115: 147-180.
[5] Cooley, T. F. 1995. Frontiers of Business Cycle Research. Princeton: Princeton Univer-
sity Press.
[6] DeJong, D. N. and M. Ripoll. 2004. Self-control preferences and the volatility of stock
[7] DeJong, D.N. and C.H. Whiteman. 2003. Unit roots in U.S. macroeconomic time series:
a survey of classical and Bayesian methods. In D. Brillinger et al., Eds. New Directions
[8] Ferson, W. E. and G. M. Constantinides. 1991. Habit formation and durability in ag-
[9] Gul, F. and W. Pesendorfer. 2004. Self-control, revealed preference, and consumption
41
[10] Heaton, J. 1985. An empirical investigation of asset pricing with temporally dependent
[11] Ireland, P. 2004. Technology shocks in the new Keynesian model. Review of Economics
[13] Kocherlakota, N. R. 1996. The equity premium: it’s still a puzzle. Journal of Economic
[14] Kydland, F. E. and E. C. Prescott. 1982. Time to build and aggregate ‡uctuations.
[15] Lucas, R. J. 1978. Asset prices in an exchange economy. Econometrica. 46: 1429-1445.
[16] Mankiw, N. G. and D. Romer. 1991. New Keynesian Economics. Cambridge: MIT Press
[17] Mehra, R. and E. C. Prescott. 1985. The equity premium: a puzzle. Journal of Monetary
[19] Rotemberg, J. 1982. Sticky prices in the United States. Journal of Political Economy.
90: 1187-1211.
[20] Shiller, R. J. 1981. Do stock prices move too much to be justi…ed by subsequent changes
42
[22] Stock, J. H. 2004. Unit roots, structural breaks, and trends. In R. F. Engle and D. L.
[23] Stokey, N. L. and R. E. Lucas. 1989. Recursive Methods in Economic Dynamics. Cam-
[24] Summers L. 1986. Some skeptical observations on real business cycle theory. Federal
[25] Woodford, M. 2003. Interest and Prices: Foundations of a Theory of Monetary Policy.
43