Capitulo4 (RBC) 4

Capitulo 4
Modelo RBC con trabajo constante*

Hamilton Galindo
Arizona State University (ASU)
hamilton.galindo@asu.edu
Alexis Montecinos
Massachusetts Institute of Technology (MIT)
alexis.montecinos@sloan.mit.edu
Borrador: 8 de julio de 2017
Indice
1. Introduccion 3
2. Construccion del modelo 3

2.1. Familias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.2. Empresas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.3. Equilibrio de mercado y definicion del choque . . . . . . . . . . . . . . . . . 8
2.4. Ecuaciones principales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3. Calibracion 8
4. Estado estacionario 10
5. Log-linealizacion 12
5.1. Efecto sustitucion y efecto ingreso de la tasa de interes . . . . . . . . . . . . 16
6. Solucion del sistema lineal 24

6.1. Metodo de coeficientes indeterminados . . . . . . . . . . . . . . . . . . . . . 24
6.2. Analisis de elasticidades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
7. Representacion de series de tiempo 31

7.1. Serie de tiempo del capital . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
7.2. Serie de tiempo del producto . . . . . . . . . . . . . . . . . . . . . . . . . . 32
7.3. Serie de tiempo del consumo . . . . . . . . . . . . . . . . . . . . . . . . . . 32
7.4. Serie de tiempo de la tasa de interes real bruta . . . . . . . . . . . . . . . . 33
7.5. Serie de tiempo de la inversion . . . . . . . . . . . . . . . . . . . . . . . . . 33
Todos los derechos reservados.

*
Disclaimer: cualquier error u omision es responsabilidad de los autores.
1
Macrodinamica (Tomo I) Modelos de ciclos economicos reales
8. Funciones impulso-respuesta 33
9. Simulacion de las variables endogenas 41
10.Componente cclico de las variables simuladas 43
11.Calculo de los momentos teoricos 43
12.Comparacion modelo teorico con los datos empiricos 47
13.Codigos 48
2
1. Introduccion
El objetivo de este captulo es comprender detalladamente el proceso de construccion
y solucion de un modelo de ciclos economicos reales. Ademas, entender como se construye
la simulacion de las variables y como se obtiene la funcion impulo-respuesta. Para ello, en
este captulo, se analiza en detalle uno de los modelos propuestos por Campbell (1994).
El modelo base propuesto por Campbell (1994) es un modelo estacionario (sin ten-
dencia), pero con crecimiento diferente de cero en el estado estacionario. Este modelo es
un extension del modelo de crecimiento estocastico, el cual permite rastrear los efectos
dinamicos de cualquier evento aleatorio (choque).
No obstante, la solucion del modelo estocastico es difcil principalmente por las no-
linealidades que emergen del mismo modelo, las cuales se derivan de la interaccion entre ele-
mentos multiplicativos (funcion de produccion Cobb-Douglas) y elementos aditivos (ley de
movimiento del capital). Un caso especial es el modelo propuesto por Long y Plosser (1983),
descrito en el captulo 3. En ese modelo las no-linealidades desaparecen debido al supuesto
no realista que la depreciacion es total; es decir, la tasa de depreciacion es igual a uno
( = 1), y que ademas, la funcion de utilidad es logartmica (u(ct , ht ) = lnct + ln(1 ht )).
En este caso el modelo llega a ser lineal y puede ser resuelto analticamente; en los demas
casos una solucion aproximada es requerida.
En lnea con lo anterior, Campbell (1994) menciona que un paper tpico en la literatura
RBC plasma el modelo y luego se mueve directamente a la discusion de las propiedades de
solucion, sin especificar como se llego a dicha solucion. Lo anterior no permite que el lector
entienda el proceso para obtener dichas propiedades de solucion, ni la solucion en s misma.
Ante ello, el autor propone un enfoque analtico simple del modelo de crecimiento
estocastico, cuya version log-lineal puede ser resuelto analticamente para mostrar el me-
canimo de solucion lo mas transparente posible. Con el fin de ilustrar el metodo de solucion,
Campbell (1994) lo aplica a cuatro modelos: [1] modelo con oferta ed trabajo fija, [2] mo-
delo con oferta de trabajo variable y con funcion de utilidad aditivamente separable, [3]
modelo con oferta de trabajo variable y con funcion de utilidad no aditivamente separable,
y [4] el segundo modelo extendido con un choque de gasto publico.
Este captulo se centra en el primer modelo (oferta de trabajo constante), dejando para
el siguiente captulo el modelo con oferta de trabajo variable.
2. Construccion del modelo

Este modelo esta compuesto por familias y empresas en un entorno de economa ce-
rrada, en la cual existe un unico bien. Por un lado las familias tienen trabajo fijo; es decir,
todas las familias estan empleadas. Por otro lado las familias son duenas del capital y por
tanto demandan bienes para invertir lo cual a su vez crea una oferta de capital. Asimismo,
las familias demandan bienes de consumo.
3
De otro lado, las empresas tienen un tecnologa para producir el unico bien en la
economa en funcion del capital. Por ello las empresas demandan capital. En la figura [1]
se esquematiza el modelo.
Figura 1: Esquema del modelo de oferta de trabajo constante
Economa cerrada
Familias
Mercado de capital
(competencia perfecta)
Mercado de bienes
(competencia perfecta)

Empresas
2.1. Familias
En este modelo se asume que la economa esta poblada por un conjunto de familias
identicas que tienen vida infinita. La familia representativa busca maximizar su funcion
de utilidad descontada:

X
Max E0 t u(ct ) (1)
{ct ,kt+1 }t=0
t=0
Donde ct es el consumo del periodo t y es el factor de descuento. Ademas, la funcion

de utilidad instantanea esta descrita por la siguiente forma funcional:
ct1
u(ct ) = (2)
1
4
Propiedades de la funcion de utilidad
La funcion de utilidad previa tiene un coeficiente de aversion al riesgo igual a y una

elasticidad de sustitucion intertemporal (del consumo) = 1/
c
Calculo de ESIt+1,t ():
uct = c
t

c uct ct
T M gSIt+1,t = Et = Et
uct+1 ct+1
c
ln( ct+1
ct )
c
T M gSIt+1,t 1 1
ESIt+1,t = c = ct+1 c
T M gSIt+1,t
=
ln(T M gSIt+1,t ) ct
ct+1

ct
La elasticidad de sustitucion intertemporal del consumo () se entiende como la disposi-

cion de la familia de sustituir consumo hoy ( ct ) por consumo de manana ( ct+1 ). Cuando
se dice que dicha elasticidad es fuerte ( es grande), se entiende que el consumidor esta
dispuesto a reducir su consumo hoy en mayor cantidad.
De otro lado, se asume que la familia es duena del capital fsico (kt ), cuya dinamica de
acumulacion esta representada por la ley de movimiento del capital:
kt+1 = (1 )kt + it (3)
Dicho capital (kt ) es alquilado a las empresas a una tasa de interes real rt . Este flujo
(rt kt ) positivo representa los ingresos de la familia, los cuales son distribuidos entre el
consumo (ct ) y la inversion (it ). Esta equivalencia de flujos, para cada periodo de tiempo,
esta representada en la restriccion presupuestaria.
ct + it = rt kt (4)
Problema de optimizacion
El problema de optimizacion de la familia representativa es el siguiente:

X c1
Max E0 t t
{ct ,kt+1 }t=0 1
t=0
sujeto a la restriccion presupuestaria:
ct + kt+1 (1 )kt = rt kt
Donde la inversion (it ) se ha reemplazado por su expresion derivada de la ley de movimiento

del capital (ecuacion (3)). Ademas, cabe mencionar que las variables de control, en este
problema de optimizacion, son: ct y kt+1 .
El problema de optimizacion de las familias puede ser escrito como una funcion de La-
grange:

X
t u(ct ) + t rt kt (ct + kt+1 (1 )kt )

L = E0
t=0
5
Donde, de manera similar al capitulo 3 (modelo Lomg y Plosser (1983)), la version

extendida de la funcion de Lagrange se puede expresar de la siguiente manera:

L = E0 0 u(c0 ) + 0 r0 k0 (c0 + k1 (1 )k0 ) +

1 u(c1 ) + 1 r1 k1 (c1 + k2 (1 )k1 ) +

2 u(c2 ) + 2 r2 k2 (c2 + k3 (1 )k2 ) +

3 u(c3 ) + 3 r3 k3 (c3 + k4 (1 )k3 ) +

4 u(c4 ) + 4 r4 k4 (c4 + k5 (1 )k4 ) +

... +
t u(ct ) + t rt kt (ct + kt+1 (1 )kt ) +

t+1 u(ct+1 ) + t+1 rt+1 kt+1 (ct+1 + kt+2 (1 )kt+1 ) +

...+
...
Las condiciones de primer orden, en el periodo t, son:

L t

= 0 = E0 uct + t (1) = 0
ct
uct = t (5)

L t
t+1

= 0 = E0 t (1)] + t+1 (rt+1 + (1 )) = 0
kt+1
t = Et t+1 (rt+1 + (1 )) (6)

Reemplazando le ecuacion (5) en la ecuacion (6) se obtiene la ecuacion de Euler:
uct = Et uct+1 (rt+1 + (1 ))

c
t = Et c
t+1 (rt+1 + (1 )) (7)
En lnea con Campbell (1994), se define la variable Rt como la tasa bruta de interes
de la inversion en capital de un periodo, el cual es igual a la tasa de interes real neta (rt )
mas el capital no depreciado (1 ). En el periodo t+1 esta relacion se expresa de la
siguiente manera.
Rt+1 = rt+1 + (1 ) (8)

Considerando la expresion anterior, la ecuacion de Euler tendra la forma siguiente:
c
t = Et ct+1 Rt+1 (9)
La ecuacion de Euler expresa una comparacion beneficio/costo marginal de consumir
una unidad del bien. Por un lado se tiene el costo marginal de dejar de consumir una
6
unidad adicional del bien, el cual es expresado por la utilidad marginal uct . Por otro lado
se tiene el beneficio marginal de no consumir dicha unidad del bien en t, la cual en el
periodo siguiente t + 1 se convierte en 1(1 + rt+1 ) unidades de bien. Esto se debe a
que existe una tasa de interes y una tasa de depreciacion. La utilidad marginal que brinda
esta unidad adicional en t+1 es uct+1 Rt+1 . Sin embargo, para compararlo con el costo
marginal en t es necesario traerlo a valor presente por medio del factor de descuento
. Por tanto, el beneficio marginal, en t, es igual a uct+1 (Rt+1 ). Esto se observa en
la siguiente ecuacion.
uct = Et uct+1 (rt+1 + (1 ))

|{z} | {z }
costo marginal beneficio marginal
Por tanto, la ecuacion de Euler indica que la familia esta dispuesta a sacrificar consumo
hoy hasta que el costo marginal de dejar de consumir una unidad del bien hoy sea igual
al beneficio marginal de dicha unidad del bien traido a valor presente.
2.2. Empresas
Se asume que las empresas se desarrollan en un contexto de competencia perfecta tanto
en el mercado de bienes como en el mercado de factores de produccion. En este escenario,
la empresa representativa maximiza su funcion de beneficios sujeta a su tecnologa (funcion
de produccion). Dicho problema de optimizacion esta descrito de la siguiente manera:
Max t = yt rt kt
{kt }
t=0
Sujeto a la funcion de produccion:
yt = at kt1 (10)
La funcion de produccion solo depende de la productividad at y del capital kt debido
a que se asume que el trabajo ht es constante (fijo). Ademas, debido a que la empresa no
toma decisiones intertemporales, su problema de optimizacion se realiza para cada uno de
los periodos. Por tanto, el problema de optimizacion se puede realizar en t y extender el
resultado para los siguientes periodos.
Introduciendo la funcion de produccion en la funcion objetivo y derivando esta ultima

con respecto a la unica variable de control (kt ), se obtiene la siguiente expresion.

(at kt1 rt kt ) at
= 0 = = 0 = (1 ) rt = 0
kt kt kt
De esta condicion de primer orden se obtiene la demanda de capital:

at
rt = (1 ) (11)
kt
7
2.3. Equilibrio de mercado y definicion del choque

Para completar el modelo antes descrito es necesario especificar dos ecuaciones adi-
cionales. La primera describe el equilibrio en el mercado de bienes; es decir, todo lo que
se produce en la economa debe encontrar su contraparte en los diferentes componentes
del gasto agregado. La segunda especifica el comportamiento de la productividad. Con
respecto a esta ultima, usualmente se supone que es estacionaria en media y que tiene una
varianza constante. La forma estandar de representarla es asumiendo que la productividad
sigue un proceso autorregresivo de orden uno.
En este modelo en particular, se asume que no existe gasto de gobierno (gt = 0) y

que la economa es pequena y cerrada. Por tanto, toda la produccion tendra dos posibles
destinos: el consumo (ct ) y la inversion (it ). En ese sentido, la condicion de equilibrio esta
descrita por la siguiente ecuacion:
yt = ct + it (12)
De otro lado, la productividad sigue un comportamiento estacionario AR(1), en la cual
el choque esta representado por el ruido blanco t , que tiene una funcion de distribucion
normal con media cero y varianza constante [N (0, 2 )]. En estado estacionario, se asu-
me que dicho ruido blanco toma el valor de su media. Asimismo, cuando se dice que la
economa ha sufrido un choque en t = 0 significa que en dicho periodo el ruido blanco
(t ) ha dejado de ser cero y ha tomado, solo en ese periodo, algun valor proporcional a
su desviacion estandar (n ). Usualmente, se considera que n es igual a uno. La ecuacion
(13) describe el comportamiento de la productividad.
lnat = lnat1 + t (13)

Cabe subrayar que el logartmo de la productividad se comporta como un AR(1) y
no la productividad en s misma. Esto es importante porque permite que, en el estado
estacionario, la productividad sea igual a uno, lo cual evita cualquier division entre cero.
2.4. Ecuaciones principales

Las ecuaciones principales del modelo se resumen en cuadro [1]:
Este conjunto de ecuaciones representan un sistema de ecuaciones en diferencias no
lineales y estocasticas. Para resolver dicho sistema, como es usual en la literatura, se trans-
forma en un sistema de ecuaciones lineales. Esto es debido a que las tecnicas matematicas
de solucion de sistemas lineales son ampliamente conocidas en la literatura. La solucion del
sistema lineal sera una aproximacion de la solucion del sistema no-lineal. Cabe mencionar
que un paso previo a la linealizacion de sistema de ecuaciones es la asignacion de valores
a los parametros (calibracion) y el calculo del estado estacionario.
3. Calibracion
Calibracion es una metodologa emprica, la cual consiste en asignar un valor a los
parametros del modelo de equilibrio general basado en una diversidad de fuentes. Segun
Heer y Mauner (2009), las fuentes mas comunes son las siguientes:
8
Cuadro 1: Sistema de ecuaciones no lineal del modelo

Ecuaciones Descripcion
c
t = Et ct+1 Rt+1 Ecuacion de Euler
yt = at kt1 Funcion de produccion

rt = (1 ) aktt Demanda del capital
Rt = rt + (1 ) Rt es la tasa de interes real (bruta)
rt es la tasa de interes real (neta)
que considera la depreciacion
yt = ct + it Equilibrio mercado de bienes
kt+1 = (1 )kt + it Ley de movimiento del capital
lnat = lnat1 + t Choque de productividad
Nota: Estas 7 ecuaciones se pueden escribir directamente en un mod
en Dynare para obtener la solucion del modelo y los IRFs.
1. El uso del promedio del nivel de variables economicas de series de tiempo o el pro-
medio de los ratios de dichas variables.
2. La estimacion econometrica de una ecuacion.
3. Referencia a estudios econometricos basados en datos microeconomicos o macro-

economicos.
4. Ajustar los parametros para que el modelo replique ciertos hechos empricos como
segundo momentos de los datos o impulso-respuesta de un VAR estructural.
La forma de evaluar el poder del modelo para capturar la realidad es por medio de la
comparacion de los valores de los segundos momentos y de las funciones impulso-respuesta
con los valores obtenidos empricamente. En el cuadro [2] se indica los valores de los
parametros del modelo, los cuales estan basados en Campbell (1994).
Cuadro 2: Calibracion (valores base)

Parametro Nombre Sustento anual
= 0.667 (1) es la participacion del capital
en el producto
= 0.025 Tasa de depreciacion 10 % anual
ln(Rss ) = 0.015, lleva a Rss = Tasa de interes real bruta de estado 6.184 % anual:
1.015 y por tanto: = 0.9852 estacionario 4
(1 + 0.015) 1
= 0.2 Tlasticidad de sustitucion intertem-
poral del consumo
= 0.95 Persistencia del choque
= 1 Desviacion estandar del choque
9
4. Estado estacionario
Para el calculo del estado estacionario se considera que la variable xt se mantiene
constante. Entonces, en el estado estacionario se tiene que xt = xt+1 = xss . Esta ultima
condicion se aplica a todas las variables endogenas. Ademas, en el estado estacionario el
choque ss toma su valor promedio, que es igual a cero.
Para la ecuacion de Euler se tiene lo siguiente:
c
t = Et c
t+1 Rt+1
c
ss = c
ss Rss
1 = Rss
1
Rss = (14)

Para la funcion de produccion:
yt = at kt1
yss = ass kss
1
(15)
Para la demanda de capital:
at

rt = (1 )
kt
ass

rss = (1 ) (16)
kss
De la ecuacion de la tasa de interes bruta:
Rt = rt + (1 )
Rss = rss + (1 )
por la ecuacion (14):
1
= rss + (1 )

1
rss = (1 ) (17)

Para la ecuacion de equilibrio en el mercado de bienes:
yt = ct + it
yss = css + iss (18)
10
De la misma manera para la ley de movimiento del capital:
kt = (1 )kt + it
kss = (1 )kss + iss
iss = kss (19)
Finalmente para la ecuacion de comportamiento de la productividad:
lnat = lnat1 + t
lnass = lnass + ss
|{z}
=0(valor de su media)
lnass = lnass
ln(ass ) = ln(ass )
ass = ass (20)
Al igual que en el modelo de Long y Plosser (1983), dos valores de ass podran resolver
esta ultima ecuacion (20): ass = 1 o ass = 0. Sin embargo, solo cuando ass = 1, el lnass
existe. Por tanto, la solucion correcta es ass = 1. La ventaja de considerar la ecuacion del
choque de productividad en logaritmos es que evita que la productividad en estado esta-
cionario pueda ser cero. Esto es importante porque evita que en las ecuaciones de estado
estacionario y en las ecuaciones log-lineales se encuentre algun numero o variable divida
por cero.
Hasta aqu se ha encontrado el valor de estado estacionario de la tasa de interes bruta

Rss , de la tasa de interes neta rss y de la productividad ass ; sin embargo para encontrar
el estado estacionario para las demas variables se tiene que hacer algunas operaciones
algebraicas adicionales. De la ecuacion (16) se tiene:
ass

rss = (1 )
kss
Como ya se conoce el valor de rss por la ecuacion (17) y de ass , entonces se puede
conocer el valor del capital kss .
ass

rss = (1 )
kss
1
rss
kss = ass (21)
(1 )
Debido a que ya se conoce kss , entonces se puede hallar el valor del producto yss , de
la inversion iss y del consumo css :
yss = ass kss

1
, de la ecuacion (15) (22)
iss = kss , de la ecuacion (19) (23)
css = yss iss , de la ecuacion (18) (24)
11
En el cuadro [3] se resume la expresion del estado estacionario de cada variable del
modelo.
Cuadro 3: Estado estacionario

Estado estacionario (forma recursiva) Estado estacionario (forma parametrica)
Rss = 1 = 1
rss = Rss (1 ) = 1 (1 )
ass = 1 =1
1 1 1
(1)
rss
kss = ass (1) = 1
1 (1)
(1)
yss = ass kss
1 =
1
1 1
(1)

iss = kss = 1
1 1 1
+1 (1)

css = yss iss = 1 1
Nota: El calculo de los estados estacionarios se encuentran en Campbell Lfijo.m (seccion 1).
5. Log-linealizacion
El sistema de ecuaciones que describe el modelo de Campbell (1994) es no lineal. Es-
ta caracterstica del modelo dificulta la forma de encontrar la solucion de dicho sistema.
Una forma estandar de abordar esta dificultad es log-linealizar cada ecuacion; es decir,
convertir una ecuacion no lineal en una ecuacion lineal en terminos de log desviacion de
la variable con respecto a su estado estacionario. Ademas, para pequenas desviaciones del
estado estacionario, la log desviacion de una variable tiene una interpretacion economica
importante: ella es aproximadamente igual a la desviacion, en porcentaje, del estado esta-
cionario (Uhlig, 1995).
La ventaja de aplicar log-linealizacion es que convierte el sistema no lineal en lineal, al

cual se le puede aplicar los metodos matematicos estandar para resolver dichos sistemas
(Blanchard y Kahn, 1981).
En primer lugar, se define la variable en log-desviaciones:
bt = lnxt lnxss
x (25)
En segundo lugar, despejando la variable xt de la ecuacion [25] se tiene:
xt = xss exbt (26)

En tercer lugar, se hace una aproximacion de Taylor de primer orden de exbt con respecto
al estado estacionario, en el cual x
bt = 0; es decir, xt = xss :
12
exbt xbt =0
= exbt =0 + exbt =0 (b

xt 0)
exbt

x
bt =0 = 1+x bt
ex
bt
= 1+x
bt (27)
Esta ultima ecuacion se reemplaza en la ecuacion (26):
xt = xss exbt
= xss (1 + x
bt ) (28)
De la ecuacion (28) se despeja x

bt :
xt xss
bt
x = (29)
xss
Por tanto, la variable en log-desviaciones es aproximadamente igual a la desviacion, en
porcentaje, del estado estacionario. De un punto de vista practico, se puede reemplazar
cada variable por su expresion log-lineal y luego se aplica la aproximacion de primer orden
segun la ecuacion (27).
Log-linealizando la ecuacion de Euler se tiene:
c
t = Et c
t+1 Rt+1
ct

= Et css ebct+1 Rss eRt+1

css eb b
ebct = Et ebct+1 eRt+1

b
ebct = Et ebct+1 +Rt+1

b

1 b
ct = Et 1 b ct+1 + Rbt+1
1b
ct
b = Et bct+1 R t+1 (30)

Haciendo lo mismo para la funcion de produccion:
yt = at kt1
b 1
yss eybt = ass ebat kss ekt

yss eybt = ass ebat kss

1 (1)kt
e
b
eybt = ebat +(1)kt

b
1 + ybt = 1 + bat + (1 )b
kt
ybt = bat + (1 )kt
b (31)
Con respecto a la demanda de capital:
13

at
rt = (1 )
kt
ass ebat

rss erbt = (1 )
kss ebkt
ass ebat

rss erbt = (1 )
kss kt
e
b

ass
rss erbt = (1 ) (e(bat kt ) )
b
kss
erbt = e(bat kt )
b
at b
1 + rbt = 1 + (b kt )
at kt )
rbt = (b b (32)
En el caso de la tasa bruta de interes, su forma log-lineal se obtiene de la siguiente

manera:
Rt = rt + (1 )
Rss eRt = rss erbt
b
Rss (1 + Rbt ) = rss (1 + rbt )

bt = rss rbt
R (33)
Rss
En el equilibrio de mercado de bienes:
yt = ct + it
yss eybt = css ebct + iss eit
b
yss (1 + ybt ) = css (1 + b

ct ) + iss (1 + bit )
yss + yss ybt = css + css bct + iss + issbit
yss ybt = css b
ct + issbit
css iss b
ybt = ct +
b it (34)
yss yss
La ley de moviento de capital en su forma log-lineal quedara:
kt+1 = (1 )kt + it
kss ekt+1 = (1 )kss ekt + iss eit
b b b
kt+1 ) = (1 )kss (1 + b
kss (1 + b kt ) + iss (1 + bit )
kss + kss bkt+1 = (1 )kss + (1 )kss b kt + iss + issbit
kt+1 = (1 )kss b
kss b kt + issbit
iss b
kt+1 = (1 )b
b kt + it (35)
kss
14
Finalmente, la ecuacion de la productividad:
lnat = lnat1 + t
lnass ebat = lnass ebat1 + t
lnass + b
at = lnass + b
at1 + t
at = b
b at1 + t (36)
El cuadro [4] resume las ecuaciones log-lineal del modelo:
Cuadro 4: Ecuaciones log-lineal

Ecuaciones log-lineal Descripcion
ct+1 1 R

[1] ct = Et b
b bt+1 Ecuacion de Euler
[2] ybt = b at + (1 )b kt Funcion de produccion
[3] rbt = [bat b kt ] Demanda de capital
[4] Rbt = rss rbt Tasa de interes bruta
Rss
[5] ybt = ycss
ss
ct + yiss
b ss
bit Equilibrio en el mercado de bienes
[6] kt+1 = (1 )b
b kt + kiss
ss
bit Ley de movimiento del capital
[7] at = b
b at1 + t Choque de productividad
Nota: Para obtener directamente la solucion del modelo con Dynare se puede
utilizar el mod Campbell Lfijo Dynare.mod
El numero de ecuaciones del cuadro [4] se puede resumir en cinco, para ello se introduce
la ecuacion de equilibrio del mercado de bienes (ecuacion 5) en la ecuacion del movimiento
del capital (ecuacion 6). La variable que relaciona ambas ecuaciones es la inversion. En
primer lugar se despeja la inversion de la ecuacion 5:

bit = ybt css b
ct
yss
yss iss
En segundo lugar, se introduce esta ecuacion en la ley de movimiento de capital:

iss css yss
kt+1 = (1 )kt +
b b ybt ct
b
kss yss iss
Ademas, se introduce la ecuacion de la funcion de produccion (yt ):

iss css yss
kt+1 = (1 )kt +
b b bat + (1 )kt
b ct
b
kss yss iss
Ordenando los terminos algrebraicos se tiene:

yss b yss css
kt+1 = (1 ) + (1 )
b kt + at b
b ct (37)
iss i iss
| {z } | {zss}
1 2
De los coeficientes de la ecuacion (37) se demuestra que:

css
= 1 1 2
iss
15
Por tanto, la ecuacion final es:
kt+1 = 1 b
b at + (1 1 2 )b
kt + 2 b ct (38)
De otro lado, la ecuacion [3] (demanda de capital) se introduce en la ecuacion [4] (tasa
de interes bruta):
bt = rss rbt
R
Rss
rss
Rt =
b at b
[b kt ]
Rss
R at b
bt = 3 [b kt ] (39)
Donde en la ecuacion previa se ha definido el coeficiente 3 :

rss
3 =
Rss
El cuadro [5] resume las cinco principales ecuaciones log-lineal del modelo de trabajo
fijo de Campbell (1994).
Cuadro 5: Ecuaciones log-lineal (sistema reducido)

Ecuaciones log-lineal
ct+1 1 R

[1] ct = Et b
b bt+1
[2] at + (1 )b
ybt = b kt
[3] Rt = 3 [b
b at kt ]
b
[4] kt+1 = 1 b
b kt + 2 bat + (1 1 2 )b
ct
[5] at = b
b at1 + t
5.1. Efecto sustitucion y efecto ingreso de la tasa de interes

Antes de resolver el sistema log-lineal es importante analizar el impacto de la tasa de
interes real sobre el consumo. Para abordar este analisis es muy util utilizar las ecuaciones
log-lineales.
La teora del consumidor sugiere que cuando el precio (pt ) de un bien (qt ) cambia hay
dos efectos sobre el consumidor: primero, el precio de qt relativo a otros productos cambia.
Segundo, debido al cambio en pt , el ingreso real del consumidor tambien cambia. El cambio
del consumo optimo como resultado de un cambio en el precio contiene ambos efectos.
El efecto sustitucion es el efecto obtenido solo por el cambio de precios relativos, man-
teniendo constante el ingreso real. Mientras que el efecto ingreso es el efecto obtenido solo
por el cambio en el ingreso real.
La tasa de interes representa el precio relativo de la canasta en el periodo t+1

(ct+1 ) con respecto a hoy (ct ). Por tanto, un cambio en la tasa de interes producira dos
16
efectos: sustitucion e ingreso.
Efecto sustitucion (ES): un incremento en la tasa de interes real hace que el consumo
de manana ct+1 sea relativamente menos costoso comparado con el consumo de hoy ct .
Esto se debe a que el ahorro es mas rentable para alcanzar el mismo monto de consumo
manana; es decir, el consumidor necesita sacrificar menos consumo hoy. Por tanto, el efecto
sustitucion se resume en:
Efecto Sustitucion
Rt ct y ct+1
Cabe mencionar que la ecuacion de Euler refleja el efecto sustitucion del consumo.
Ademas, es la elasticidad de sustitucion intertemporal del consumo.
1b
b ct+1 R
ct = Et b t+1

La magnitud del efecto sustitucion es controlado por , mientras mas grande sea
mayor sera el efecto sustitucion; es decir:
Efecto Sustitucion
Rt ct y ct+1
Efecto ingreso (EI): un incremento de la tasa de interes produce un efecto ingreso. Si
el consumidor tiene activos (bonos o ahorro), un incremento de la tasa de interes pro-
duce mayores ganancias por esos activos y por tanto mayor ingreso. Este efecto tiende a
incrementar el consumo en todos los periodos.
Efecto Ingreso
Rt ct y ct+1
Cabe mencionar que la restriccion presupuestaria refleja el efecto ingreso:
ct + it = rt kt
Un incremento de la tasa de interes produce dos efectos:

ES ct y ct+1 (Ecuacion de Euler)
EI ct y ct+1 (Restriccion presupuestaria)

ET Depende de ESI y ct+1
Efecto total (ET): para observar el efecto final de la tasa de interes sobre el consumo
nos basaremos en la restriccon presupuestaria y la ecuacion de Euler (de las variables en
niveles).
17
ct + it = rk kt (40)
pero se sabe :
kt+1 = (1 )kt + it
despejando it :
it = kt+1 (1 )kt (41)
(41) en (40) :
ct + kt+1 (1 )kt = rk kt
ct + kt+1 = (rk + (1 ))kt
| {z }
Rt
ct + kt+1 = Rt kt (42)
Como se sabe el ingreso de la familia representativa en t es Rt kt , la cual se resumira

en At . De igual forma para el ingreso en t+1: Rt+1 kt+1 = At+1 . Reescribiendo la ecuacion
(42) en terminos de ingreso se tiene:
ct + kt+1 = Rt kt
Rt+1 kt+1
ct + = Rt kt
Rt+1
At+1
ct + = At (43)
Rt+1
La ecuacion (43) es una ecuancion en diferencias, la cual se puede resolver iterando

hacia adelante. Por induccion matematica hacemos lo siguiente:
At+1
At = ct + (44)
Rt+1
At+2
At+1 = ct+1 + (45)
Rt+2
At+3
At+2 = ct+2 + (46)
Rt+3
Luego la ecuacion (46) se reemplaza en (45):
At+2
At+1 = ct+1 +
Rt+2
1 At+3
At+1 = ct+1 + (ct+2 + )
Rt+2 Rt+3
ct+2 At+3
At+1 = ct+1 + + (47)
Rt+2 Rt+2 Rt+3
La ecuacion (47) se reemplaza en (44):
18
At+1
At = ct +
Rt+1
1 ct+2 At+3
At = ct + (ct+1 + + )
Rt+1 Rt+2 Rt+2 Rt+3
ct+1 ct+2 At+3
At = ct + + + + (48)
Rt+1 Rt+1 Rt+2 Rt+1 Rt+2 Rt+3
Dividiendo toda la ecuacion (48) por Rt para hacer una generalizacion (en sumatoria)
mas sencilla:
At ct ct+1 ct+2 At+3

= + + +
Rt Rt Rt Rt+1 Rt Rt+1 Rt+2 Rt Rt+1 Rt+2 Rt+3
resumiendo : en una sumatoria...
2
At c At+3
Qs t+s
X
= + Q3 (49)
Rt j=0 Rt+j j=0 Rt+j
s=0
generalizando para n :
n
At X ct+s At+(n+1)
= Qs + Qn+1
s=0
aplicando Limite cuando :n

At c At+(n+1)
Qs t+s
X
= + Limn Qn+1
s=0
| {z }
=0(por transversalidad)

At X ct+s
= Qs (50)
Rt j=0 Rt+j
s=0
Para encontrar la relacion de la tasa de interes con el consumo de hoy es necesario

encontrar la relacion del ct+s con el consumo actual ct , para ello se usa la ecuacion de
Euler (abstrayendo el operador expectativa) para t, t + 1 y t + 2:
c
t = c
t+1 Rt+1
c
t+1 = ct+2 Rt+2
c
t+2 = ct+3 Rt+3
Multiplicando estas ecuaciones se tiene:
19
c 3
t ct+1 ct+2 = ct+1 Rt+1 ct+2 Rt+2 ct+3 Rt+3
Rt
c
t = 3 c
t+3 Rt+1 Rt+2 Rt+3
Rt
Y 3
3 ct+3
c
t = Rt+j
Rt
j=0
generalizando para s :
c
s
c t+s
Y
t = s Rt+j
Rt
j=0

ct+s R
= Qs t
ct s j=0 Rt+j
despejando ct+s :
1
Rt
ct+s = Qs ct (51)
s j=0 Rt+j
Introduciendo la ecuacion (51) en la ecuacion (50):
1

s
QsRt ct
At X j=0 Rt+j
= Qs
Rt j=0 Rt+j
s=0
s
At X s Y 1 1 1/
= Rt+j
ct Rt
Rt
s=0 j=0

X s
At s Y 1 1 1/
= ct Rt+j
Rt (52)
Rt
s=0 j=0
Caso simplificado: para analizar el efecto de la tasa de interes sobre el consumo de hoy
ct se supone que la tasa de interes es la misma en todos los periodos; es decir, Rt = Rt+1 =
Rt+2 = ... = Rt+j = R. Introduciendo este supuesto en la productoria de la ecuacion (52)
se tiene:
s
Y
Rt+j = Rs+1
j=0
Reemplazando la expresion anterior en la ecuacion (52) se tiene:
20

X
At s s+1
1 1 1/
= ct R
R
R
s=0

X
At s
(s+1) 1 1
= ct R R1/
R
s=0

X
At s
(s( 1 1)+ 1 1
= ct R R1/
R
s=0

X
At s
(s( 1 1)
= ct R R1
R
s=0

X
s
(s( 1 1)
At = ct R
s=0

X s
1 1
1
At = ct R
(53)
s=0
s
1 1
P 1
Por progresion geometrica de s=0 R se tiene que:
s 2 3
X 1 1
1 1 1
1 1 1
1 1 1
1 1
R
= 1+ R
+ R
+ R
... = 1 1 (54)
1
s=0 1 R
Reemplazando la expresion (54) en la ecuacion (53) se tiene:

1
At = ct 1 1
1
1 R
1 1
1
ct = At 1 R (55)
Aplicando logartmo a la ecuacion (55) se tiene:

1 1
1
ln(At ) = ln(ct ) + ln 1 R (56)
1 1
1
Aplicando la aproximacion de Taylor de primer orden a ln 1 R se tiene que:
1 1 1 1
1 1
ln 1 R R (57)
Reemplazando (57) en (56):
1 1
1
ln(At ) = ln(ct ) + R (58)
Tomando diferencial a la ecuacion (58) y considerando que At no cambia, y ademas,
1
= (ESI), entonces:
ln(ct ) = ( 1) R R
ln(ct )
= ( 1) R (59)
R
21
La ecuacion (59) refleja el efecto final sobre el consumo de hoy un movimiento de

la tasa de interes real. Una conclusion importante es que el efecto final depende de la
elasticidad de sustitucion intertemporal del consumo (). La expresion siguiente muestra
el efecto final sobre el consumo dependiendo del valor de la ESI:
ln(ct )
< 1 > 0 ct
R
ln(ct )
= 1 = 0 R no afecta el consumo
R
ln(ct )
> 1 < 0 ct
R
Caso general: considerando la ecuacion (52) y desarrollandola se tiene:

X s
At s Y 1 1 1/
= ct Rt+j
Rt
Rt
s=0 j=0
X s
At s Y 1 1
11/
= ct Rt+j
Rt s=0 j=0
siendo explcito en la sumatoria:

At 1 1
1 2 1
1 3 1
1
11/
= ct 1 + (Rt Rt+1 )
+ (Rt Rt+1 Rt+2 )
+ (Rt Rt+1 Rt+2 Rt+3 )
...
Rt | {z }
Nt
11/
At = ct Rt [1 + Nt ]
11/ 11/
At = ct Rt + ct Rt Nt (60)
Diferenciando la ecuacion (60) con respecto a Rt+1 y considerando que Rj (j 6= 1) no

depende de Rt+1 :
At 11/ ct ct 11/ 11/ Nt

= Rt + Rt Nt + ct Rt (61)
Rt+1 Rt+1 Rt+1 Rt+1
Nt
Desarrollando el diferencial: Rt+1 ,
22
Nt 1 1 1
2 1 2 1
2
= 1 (Rt Rt+1 ) Rt + 1 (Rt Rt+1 Rt+2 ) Rt Rt+2 +
Rt+1
1 3 1
2
1 (Rt Rt+1 Rt+2 Rt+3 ) Rt Rt+2 Rt+3 + ...

multiplicando y dividiendo por Rt+1

1 1 1 1
2 1 2 1
2
= 1 (Rt Rt+1 ) Rt Rt+1 + 1 (Rt Rt+1 Rt+2 ) Rt Rt+1 Rt+2 +
Rt+1

1 3 1
2
1 (Rt Rt+1 Rt+2 Rt+3 ) Rt Rt+1 Rt+2 Rt+3 + ...

1 1 1 1
1 2 1
1 3 1
1
= 1 (Rt Rt+1 )
+ (Rt Rt+1 Rt+2 )
+ (Rt Rt+1 Rt+2 Rt+3 )
+ ...
Rt+1
1 1
= 1 Nt
Rt+1
1 Nt
= 1 (62)
Rt+1
Introduciendo la ecuacion (62) en la ecuacion (61):
At ct
11/ ct 11/ 11/ Nt
= Rt + R Nt + ct Rt (63)
Rt+1 Rt+1 Rt+1 t Rt+1
At 11/ ct ct 11/ 11/ 1 Nt
= Rt + Rt Nt + ct Rt 1
Rt+1 Rt+1 Rt+1 Rt+1
11/ 11/ 11/ 1 Nt
At = Rt ct + ct Rt Nt + Rt+1 ct Rt 1
Rt+1
Se sabe que At = 0, entonces:
11/ 11/ 11/ 1 Nt

0 = Rt ct + ct Rt Nt + Rt+1 ct Rt 1
Rt+1
11/ 1 Rt+1
0 = Rt ct + ct Nt + 1 ct Nt
Rt+1
1 Rt+1
0 = ct + ct Nt + 1 ct Nt
Rt+1
1 Rt+1
0 = ct [1 + Nt ] + 1 ct Nt
Rt+1
Ordenando algebraicamente los terminos, se tiene:
1 Rt+1
ct [1 + Nt ] = 1 ct Nt
Rt+1
ct [1 + Nt ] 1 Rt+1
= 1
ct Nt Rt+1
ct 1 Nt Rt+1
= 1 (64)
ct 1 + Nt Rt+1
23
De la ecuacion (64) se puede concluir que el impacto de la tasa de interes del perio-
do siguiente sobre el consumo de hoy es gobernada por la elasticidad de sustitucion del
consumo ( 1 = ), tal como se observo en el caso simplificado.
6. Solucion del sistema lineal

En el captulo 1 y 3 se senalo que en la literatura existen varios metodos para solucionar
sistemas de ecuaciones lineales. En el captulo 3 se ilustro el metodo de Blanchard y
Kahn (1981) y dada la naturaleza del modelo de Longy Plosser (1983) se pudo obtener
la solucion analticamente tambien. En este captulo se utilizara el metodo de coeficientes
indeterminados de Uhlig (1999) con el fin de tener un panorama de los distintos metodos
de solucion.
6.1. Metodo de coeficientes indeterminados

El metodo de coeficientes indeterminados busca que las variables de control esten en
funcion de las variables de estado (b
kt ) y de la variable exogena (b
at ). Es decir, de la misma
forma que el metodo de Blanchard y Kahn, este metodo busca la funcion de poltica y la
funcion de estado.
Al analizar si cada ecuacion log-lineal se encuentra en funcion del capital (b kt ) y de la

productividad (b at ) se observa, en el cuadro [5], que la ecuacion [2] (funcion de produccion)
y la ecuacion [3] (demanda de capital que considera la tasa de interes bruta) dependen de
dichas variables. Ademas la ecuacion [5] describe la productividad.
Al introducir la demanda de capital en la ecuacion de Euler, dicha ecuacion estara en

funcion del capital y de la productividad:
b ct+1 3 (b
ct = Et (b at+1 b
kt+1 )) (65)
De otro lado, la ley de movimiento del capital contiene a la variable de estado y al choque:
kt+1 = 1 b
b at + (1 1 2 )b
kt + 2 b ct (66)
Por tanto, si encontramos el bct y b

kt+1 en funcion de (b
kt , b
at ), el sistema estara solucionado.
Para ello, bajo el metodo de coeficientes indeterminados, se propone la siguiente solucion:
ct = ck b
b kt + ca b
at (67)
b kt + ka b
kt+1 = kk b at (68)
En este contexto, el problema radica en encontrar los valores de los coeficientes: ck ,

ca , kk , ka . Con este fin, el analisis se realizara en cinco pasos:
[1] Ecuacion de Euler: al reemplazar la solucion propuesta en la ecuacion de Euler (65)
se obtiene una expresion para los coeficientes ca y ck :
ka (3 + ck ) 3
ca = ca = f (ka , ck ) (69)
1
kk 3
ck = ck = f (kk ) (70)
1 kk
24
[2] Ecuacion del capital: al reemplazar la solucion propuesta en la ecuacion de movi-

miento del capital (66) se obtiene una expresion para los coeficientes kk y ka :
kk = 1 + (1 1 2 )ck kk = f (ck ) (71)

ka = 2 + (1 1 2 )ca ka = f (ca ) (72)
[3] Primer coeficiente: para hallar ck elegimos (70) y (71):
ck = f (kk ) :
kk 3
ck = (73)
1 kk
kk = f (ck ) :
kk = 1 + (1 1 2 )ck (74)
[4] Encontrando ck : la ecuacion (74) se reemplaza en (73), de la cual se obtiene:
2
Q2 ck + Q1 ck + Q0 = 0 (75)
Donde, en primer lugar las dos raices de esta ecuacion representan los dos valores que
puede tomar ck . En segundo lugar, el valor de este coeficiente permite obtener el valor de
los tres restantes, y finalmente, los valores de Qi son:
Q2 = 1 1 2
Q1 = 1 1 + 3 (1 1 2 )
Q0 = 1 3
Al resolver la ecuacion (75) se obtiene los dos valores de ck :

p
Q1 + Q21 4Q2 Q0
ck1 =
2Q
p 2
Q1 Q21 4Q2 Q0
ck2 =
2Q2
El signo de ck que se debe de elegir es positivo porque esto permite que kk sea menor
a uno, lo cual indica que la ecuacion del capital es estable (no explosiva). Para ello, se
evalua el signo de cada Qi :
Q2 < 0 (porque 1 > 1 y 2 > 0)
Q0 > 0
Q1 > 0 (Q1 = 1 1 + Q2 Q0 /1 )
De lo anterior, se demuestra que ck2 tiene signo positivo, por tanto se elige esta raz. Esto
permite obtener los dos coeficientes ck y kk :
p
Q1 Q21 4Q2 Q0
ck = (76)
2Q2
kk = 1 + (1 1 2 )ck (77)
25
[5] Coefientes restantes: para hallar los dos coeficientes restantes ca y ka se elige la
ecuacion (69) y (72):
ka = 2 + (1 1 2 )ca ka = f (ca )
ka (3 + ck ) 3
ca = ca = f (ka , ck )
1
ka y ca :
ck 2 + 3 ( 2 )
ca =
1 + (1 1 2 )(ck + 3 )
ka = 2 + (1 1 2 )ca
Con los parametros calibrados para el modelo base se obtiene que: ck = 0.3253,
ca = 0.2643, kk = 0.9841 y ka = 0.0551. Finalmente, la solucion del modelo para cada
una de las variables endogenas son:
Solucion para el consumo:

ct = ck b
b kt + ca b
at (78)
Solucion para el capital:
kt+1 = kk b
b kt + ka b
at (79)
Solucion para el producto:
ybt = (1 )b
kt + b
at (80)
Solucion para la inversion:
css iss b
ybt = ct +
b it
yss yss
yss css
bit = ybt ct
b
iss yss
Reemplazando (78) y (80):
yss css yss css
bit = (1 ck )b
kt + ( ca )b
at (81)
iss yss iss yss
Solucion (tasa de interes neta):
at b
rbt = (b kt ) (82)
Solucion (tasa de interes bruta):
Rbt = rss (b
at b
kt ) (83)
Rss
6.2. Analisis de elasticidades

Los coeficientes de la solucion de cada una de las variables representan elasticidades.
Esto se debe a que las variables estan expresandas en logaritmos. Por ejemplo para el caso
del consumo se tiene:
ct = ck b
b kt + ca b
at
26
ct = ln( ccsst ) y de manera similar para las demas variables se tiene:

Dado que b
ct kt at
ln( ) = ck ln( ) + ca ln( )
css kss ass
ln(ct ) ln(css ) = ck (ln(kt ) ln(kss )) + ca (ln(at ) ln(ass ))
ln(ct ) = [ln(css ) + ln(kss ) + ln(ass )] + ck ln(kt ) + ca ln(at )
Tomando diferencial con respecto al capital (kt ) se tiene:
ln(ct ) = ck ln(kt )
ct kt
= ck
ct kt
ct
ct
kt
= ck
kt
Elasticidadct ,kt = ck (84)
Como se puede observar la ecuacion (84), ck refleja la elasticidad del consumo ante
un cambio del capital. En particular, ck mide el efecto del capital (kt ) sobre el con-
sumo actual (ct ), manteniendo constante la productividad (at ); es decir, si el capital
aumenta 1 %, el consumo aumenta en ck %. De esta forma se lee todos los coeficientes de
la solucion del sistema log-lineal. El cuadro [6] resume las elasticidades.
Cuadro 6: Coeficientes (elasticidades) de la solucion del modelo lineal

Elasticidad Expresion Valor

Q1 Q21 4Q2 Q0
Elasticidad del consumo al capital: ck = 2Q2 0.3253
ck
ck 2 +3 (2 )
Elasticidad del consumo a la pro- ca = 1+(11 2 )(ck +3 ) 0.2643
ductividad: ca
Elasticidad del capital de manana al kk = 1 + (1 1 2 )ck 0.9841
capital de hoy: kk
Elasticidad del capital de manana a ka = 2 + (1 1 2 )ca 0.0551
la productividad: ka
Nota: La expresion de las elasticidades y sus valores estan en Campbell Lfijo.m (seccion 2).
En el analisis de elasticidades dos parametros son importantes: la elasticidad de sus-

titucion intertemporal del consumo y la persistencia del choque . Para ver como estos
parametros influyen sobre las elasticidades vamos a revisar cada una de las elasticidades.
Revisando 1 , 2 y 3 :
27
yss
1 = (1 ) + (1 )
i
ss
1
= (1 ) + (1 ) kss

rss
= (1 ) + (1 )
1
1
= (1 ) + ( (1 ))

1
=

1 = F () (85)
yss
2 =
i
ss
1
= kss

rss
=
1
1
= ( (1 ))
1
2 = F (, , ) (86)
rss
3 =
Rss
1
(1 )
= 1

= (1 (1 )
3 = F (, , ) (87)
Revisando Q0 , Q1 y Q2 :
Q2 = 1 1 2
1 1
= 1( ) ( (1 ))
1
1
+ 1

=
1
Q2 = F (, , ) (88)
Q1 = 1 1 + 3 (1 1 2 )
= 1 1 + 3 (1 1 2 )
|{z} | {z }
F () F (,,)
(+)
Q1 = F ( , , ) (89)
28
Q0 = 1 3
= 1 3
|{z} |{z}
F () F (,,)
(+)
Q0 = F ( , , ) (90)
De que parametros dependen las elasticidades (ck y kk )?
Debido a que Q2 es negativo, el componente dentro del radical es positivo. En ese caso
, que afecta positivamente a Q0 y Q1 , tiene un impacto positivo sobre ck . De otro lado,
Q1 que se encuentra fuera del radical tambien traslada el efecto positivo de sobre ck .
Cabe mencionar que ck no depende de la persistencia del choque ().
p
Q1 Q21 4Q2 Q0
ck = = F ( (+) , , ) (91)
2Q2
De lo anterior se concluye la siguiente observacion:
Observacion 1: ck se incrementa a medida que se incrementa la elasticidad de sustitu-
cion del consumo ().
De otro lado, al analizar el coeficiente kk se obtiene lo siguiente:
kk = 1 + (1 1 2 )ck
= 1 + (1 1 2 ) ck
|{z} | {z } |{z}
F () = ci ss F ( (+) ,,)
ss
css
= 1 ck
|{z} iss |{z}
F () F ( (+) ,,)
css
= 1 ck
|{z} iss |{z}
F () |{z} F ( (+) ,,)
F (,,)
kk = F ( () , , ) (92)
De la ecuacion (92) se concluye las siguientes observaciones:

Observacion 2: ck y kk no dependen de .
Observacion 3: kk se reduce a medida que se incrementa la elasticidad de sustitucion
intertemporal del consumo ().
De que parametros dependen las elasticidades (ca y ka )?
ck 2 + 3 ( 2 )
ca = = F (, , , , ) (93)
1 + (1 1 2 )(ck + 3 )
ka = 2 + (1 1 2 )ca = F (, , , , ) (94)
29
De la expresion (94) se puede ver que ca tiene una relacion no lineal con y . De
manera similar para ka . De lo anterior se concluye las siguientes observaciones:
Observacion 4: ca se incrementa a medida que aumenta para valores bajos de

( 6 1), pero se reduce para valores altos ( > 1).
Observacion 5: kk y ka se reducen a medida que se incrementa la elasticidad de susti-
tucion del consumo ().
Figura 2: Elasticidades (coeficientes de la solucion)
ck kk
1.8 1
1.6
0.98
1.4
0.96
1.2
1 0.94
0.8 0.92
0.6
0.9
0.4
0.88
0.2
0 0.86
0 2 4 6 8 10 0 2 4 6 8 10
Elasticidad de sustitucin intertemporal del consumo () Elasticidad de sustitucin intertemporal del consumo ()
ca ka
0.15 0.1
0.1 0.095
0.05 0.09
0 0.085
0.05 0.08
0.1 0.075
0.15 0.07
=0 =0.4 =0.8 =0 =0.4 =0.8
0.2 0.065
0 2 4 6 8 10 0 2 4 6 8 10
Elasticidad de sustitucin intertemporal del consumo () Elasticidad de sustitucin intertemporal del consumo ()
Nota: Cabe mencionar que estos graficos se obtienen del codigo Campbell Lfijo Sim Parametros.m
En el cuadro [7] se mencionan tres casos de especial interes.
30
Cuadro 7: Casos especiales

Funcion de
Caso Valor de Elasticidad Serie de tiempo
utilidad
No existe efecto
Caso 1 =0 sustitucion kk = 1 ln(ct ) es un random walk, y
intertemporal ln(kt ) y ln(kt ) cointegran con
el ln(ct )
Funcion de
utilidad
Caso 2 =1 El efecto sustitucion y el efec-
logaritmica:
to ingreso se anulan.
u(ct ) = ln(ct )
Funcion de
Caso 3 = utilidad lineal: kk = 0, ka = kt se comporta como un
u(ct ) = ct AR(1), mientras ct y yt
se comportan como un AR-
MA(1,1)
7. Representacion de series de tiempo

Debido a que se tiene la solucion del modelo; es decir, cada variable endogena en funcion
de la variable de estado (capital) y de la variable exogena (productividad), considerando
ademas que la productividad se comporta como un proceso AR(1), entonces se puede
hallar la representacion de series de tiempo ARMA (p,q) de cada variable.
7.1. Serie de tiempo del capital

De la solucion del modelo, en particular de la ecuacion que describe el comportamiento
del capital en t + 1 en funcion del capital en t y de la productividad se tiene:
kt+1 = kk b
b kt + ka b
at
Donde los coeficientes kk y ka han sido hallados previamente. De esta ecuacion se puede
encontrar la forma autorregresiva del capital (b
kt+1 ):
(1 kk L)b
kt+1 = ka b
at
ka
kt+1 =
b at (95)
1 kk L
b
Ademas, at = b
b at1 + t
Considerando que b
at se puede expresar como:
t
at = (96)
1 L
Entonces se tiene que:
ka t
kt+1 =
b (97)
(1 kk L) (1 L)
La expresion anterior demuestra que el capital se comporta como un AR(2): dos raices
reales ( y kk ) y menores a 1 (kt+1 es estable). La expresion AR(2) del capital es:
31
ka t
kt+1 =
b
(1 kk L) (1 L)
(1 kk L)(1 L)b kt+1 = ka t
2 b
(1 kk L L + kk L )kt+1 = ka t
kt+1 kk b
b kt b
kt + kk b
kt1 = ka t
b kt kk b
kt+1 = ( + kk )b kt1 + ka t (98)
7.2. Serie de tiempo del producto

De igual manera que en el caso del capital, para encontrar la expresion de series de
tiempo del producto se parte de la solucion del modelo:
at + (1 )b
ybt = b kt (99)
Para encontrar el modelo de series de tiempo del producto (yt ) se reemplaza en la

ecuacion previa (98) la expresion de la productividad (en funcion del error) y la expresion
del capital (en funcion de la productividad). Esta ultima corresponde a la ecuacion (96).
ybt = bat + (1 )b
kt
et ka
ybt = + (1 ) at1
1 L (1 kk L)
b
et ka et1
ybt = + (1 )
1 L (1 kk L) (1 L)
et ka L et
ybt = + (1 ) (100)
1 L (1 kk L) (1 L)
La ecuacion (100) sugiere que el producto se comporta como un ARMA(2,1):

+ [(1 )ka kk ]L
ybt = t (101)
(1 kk L)(1 L)

(1 kk L)(1 L)b
yt = + [(1 )ka kk ]L t
(1 kk L L + kk L2 )b
yt = t + [(1 )ka kk ]t1
ybt kk ybt1 b yt2 = t + [(1 )ka kk ]t1
yt1 + kk b
yt1 kk b
ybt = (kk + )b yt2 + t + [(1 )ka kk ]t1
| {z } | {z }
AR(2) M A(1)
7.3. Serie de tiempo del consumo

De la solucion del modelo:
ct = ck b
b kt + ca b
at
El consumo se comporta como un ARMA(2,1)

ca + (ck ka ca kk )L
ct = t (102)
(1 kk L)(1 L)
b
32
7.4. Serie de tiempo de la tasa de interes real bruta

De la solucion del modelo:
R at+1 b
bt+1 = 3 (b kt+1 )
La tasa de interes se comporta como un ARMA(2,1)

Rbt+1 = 3 (1 ka kk L) t (103)
(1 kk L)(1 L)
7.5. Serie de tiempo de la inversion

De la solucion para la inversion (ecuacion (81)):
bit = (1 css ck ) b
kt + (
css
ca ) b
at
y y
| {z ss } | {zss }
ik ia
bit = ik b
kt + ia b
at
bit = ka t1 t
ik + ia
(1 kk L) (1 L) 1 L

bit = ka L t
ik + ia
(1 kk L) 1 L

bit = ik ka L + ia (1 kk L) t
(1 kk L) 1 L

bit = ia + (ik ka ia kk )L t
(1 kk L) 1 L
Por tanto, la inversion se comporta como un ARMA(2,1):

bit = ia + (ik ka ia kk )L t (104)
(1 kk L)(1 L)
8. Funciones impulso-respuesta
La construccion de la funcion impulso-respuesta de las variables endogenas requiere de
dos etapas. En la primera etapa se transforma la forma autorregresiva del capital AR(2)
a su version de media moviles MA(). En la segunda etapa se cuantifica el impacto, en
cada periodo, de un choque temporal (de un solo periodo) en cada variable endogena.
Primera etapa: se obtiene la forma MA() del capital.
k + kk b
kt+1 = ( + kk )b k + ka t
| {z } t | {z } t1
b
1 2
kt+1 = 1 b
b kt1 + 2 b
kt + ka t
2 b
(1 1 L 2 L )kt+1 = ka t
33
Calculando las raices del AR(2):
1 1 L 2 L2 = 0
En factores:
(L y1 )(L y2 ) = 0
Factorizando y1 del primer factor e y2 del segundo:

1 1
y1 L 1 y2 L1 =0
y1 y2
La expresion se reduce a:

1 1
L1 L1 =0
y1 y2
|{z} |{z}
1 2
Multiplicando por (-) a ambos terminos:
(1 1 L)(1 2 L) = 0
Por tanto: equivalencia de raices
(L y1 )(L y2 ) = (1 1 L)(1 2 L) = 0
Donde:
1
1 = y1
1
2 = y2
Utilizando la equivalencia de raices del AR(2):
(L y1 )(L y2 )b
kt+1 = ka t
(1 1 L)(1 2 L)kt+1 = ka t
b
1
kt+1 =
b ka t
(1 1 L)(1 2 L)
| {z }
(L)
Donde:
(L) = 1 + 1 L + 2 L2 + 3 L3 + ... + k Lk + ...

k
1j 2kj
X
k = =
j=0
Version MA() del capital:
kt+1 = (1 + 1 L + 2 L2 + 3 L3 + ...)ka t
b (105)
34
Con esta expresion calculamos la funcion impulso-respuesta. La version extendida de la

ecuacion (105) es:
kt+1 = (1 + 1 L + 2 L2 + 3 L3 + ...)ka t
b
kt+1 = ka t + (1 ka )t1 + (2 ka )t2 + (3 ka )t3 + ...
b (106)
Segunda etapa: en esta etapa se calcula la funcion impulso-respuesta del capital ante un
choque de productividad. En este caso para el calculo de la funcion impulso respuesta del
capital se considera que el impulso o choque t se realiza en un solo periodo (el periodo
uno) y que toma el valor de una desviacion estandar , el cual se asume que es igual
a uno; es decir, en t = 1, 1 = = 1. El error (t ) tomar el valor de cero durante los
periodos antes del choque y despues del choque. El cuadro [8] muestra la construccion de
la funcion impulso-respuesta del capital.
Cuadro 8: Construccion de la funcion impulso-respuesta del capital

t t Version MA() de b
kt+1 IFR de bkt+1
0 0 = 0 k1 = ka 0 +(1 ka ) 1 +...
b k1 = ka 0
b
|{z} |{z}
=0 =0
1 1 = 1 k2 = ka 1 +(1 ka ) 0 +...
b k2 = ka 1
b
|{z} |{z}
=1 =0
2 2 = 0 k3 = ka 2 +(1 ka ) 1 +(2 ka ) 0 +...
b k3 = 1 ka 1
b
|{z} |{z} |{z}
=0 =1 =0
3 3 = 0 k4 = ka 3 +(1 ka ) 2 +(2 ka ) 1 +(3 ka ) 0 +...
b k4 = 2 ka 1
b
|{z} |{z} |{z} |{z}
=0 =0 =1 =0
4 4 = 0 ... k5 = 3 ka 1
b
En t = 0 todas las variables se encuentran en su estado estacionario. El capital en

t = 1, el cual se determina en t = 0, tambien se encuentra en estado estacionario. Tal es
as que se cumple la ley de movimiento del capital: k1 = (1)k0 +i0 , donde k1 = k0 = kss .
El choque de productividad se realiza en el periodo t = 1 produciendo los siguientes efectos:
1er Efecto (sobre las empresas): un incremento de la productividad produce un in-

cremento en la funcion de produccion para cada nivel de capital. El capital se hace mas
productivo en t = 1; es decir, con el mismo capital se puede producir mas. Por tanto, la
demanda de capital aumenta.
2do Efecto (sobre las empresas): el aumento de la demanda de capital permite que
la tasa de interes en t = 1 se incremente: rt (r0 r1 ), r1 > r0 . Esto se debe a que la
oferta de capital en t = 1 se mantiene constante ya que no se ve afectada por el choque
de productividad.
3er Efecto (sobre las familias): el incremento de la tasa de interes real produce un
efecto ingreso sobre el consumo:
rt (r1 > r0 ) r1 k1 > r0 k0 c1
35
Figura 3: Efecto sobre la funcion de produccion

Y
Y1
Y0
B
Y1
Y0 A
Kss = K0 = K1 K
Figura 4: Efecto sobre la demanda de capital

Ok0 = Ok1
r
r1 B
r0 A
Dk1
Dk0
Kss = K0 = K1 K
4to Efecto (sobre las familias): el incremento de la tasa de interes incentiva el ahorro,
el cual en economa cerrada es igual a la inversion. Entonces la inversion pasa de i0 a i1
(i1 > i0 ). El impacto de una mayor inversion se observa en el incremento de la oferta de
capital en el siguiente periodo (t = 2)
k2 = (1 )k1 + i1
Por tanto:
i1 > i0 k2 > k1
5to Efecto (sobre las empresas y las familias): como el impacto del choque de pro-
ductividad tiene persistencia; es decir, sus efectos son positivos aunque cada vez menores
en el tiempo. En t = 2 la funcion de produccion se incrementa generando que la demanda
de capital tambien se incremente, pero en menor magnitud que lo observado en t = 1. Esto
produce que la tasa de interes real en t = 2 sea menor que en t = 1 (r2 < r1 ); sin embargo,
sigue siendo mayor que el valor en t = 0. Entonces dado que el individuo compara su
36
Figura 5: Efecto sobre la oferta y demanda de capital

Ok1 = Ok0 Ok2
r
B
r2 r1 C
r0 A
Dk2
Dk1
Dk0
Kss = K0 K2 K
= K1
Cuadro 9: Valores de la funcion impulso-respuesta (variales log-lineales)
t ybt kt+1
b ct
b bit R
bt at
b
0 0 0 0 0 0 0
1 0.667 0.05512 0.26429 2.20468 0.02636 1
2 0.652 0.10660 0.26901 2.11441 0.02359 0.95
3 0.63747 0.15464 0.27321 2.02834 0.02098 0.9025
4 0.62337 0.19943 0.27691 1.94626 0.01852 0.85738
. . . . . . .
Nota: Debido a que el choque se realiza en el primer periodo (t = 1),
el valor de las variables en t = 0 es cero. Cabe mencionar que estos
valores se obtienen del codigo Campbell Lfijo.m (seccion 4)
situacion en cada periodo con respecto a t = 0 (estado estacionario), entonces esta mayor
tasa de interes (r2 > r0 ) produce dos efectos sobre el consumo:
r2 > r0 : Efecto sustitucion c1 c2
r2 > r0 : Efecto ingreso r2 k2 > r1 k1 c2

Por tanto el efecto final de la tasa de interes sobre el consumo, para pequeno, es:
EI > |ES| c1 c2
El cuadro [9] muestra los valores de la funcion impulso-respuesta de las variables
endogenas del modelo. Para leer correctamente estos valores se debe de recordar que estas
funciones corresponden a las variables log-lineales, las cuales por ejemplo para el produc-
to esta expresada de la siguiente manera: ybt = ln(yt ) ln(yss ) o en su forma reducida
ybt = ln yysst .

37
Figura 6: Funcion impulso-respuesta de las variables macroeconomicas log-lineales
Producto Capital Consumo

0.7 0.7 0.35
0.6 0.6 0.3
0.5 0.5 0.25
0.4 0.4 0.2
0.3 0.3 0.15
0.2 0.2 0.1
0.1 0.1 0.05
0 0 0
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200
Inversin Tasa de inters real bruta Productividad

2.5 0.03 1
0.025
2 0.02 0.8
0.015
1.5 0.6
0.01
0.005
1 0.4
0
0.5 0.005 0.2

0.01
0 0.015 0
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200
Nota: Estas funciones impulso-respuesta corresponde a las variables log-lineales; es decir a ybt , b
kt , b
ct , bit ,
at . Cabe mencionar que estos graficos se obtienen del codigo Campbell Lfijo.m (seccion 4).
rbt y b
En lnea con lo anterior, segun el cuadro [9] el valor del producto (log-lineal) en t = 0
es igual a cero. Es decir, yb0 = ln yyss0 = 0. La unica solucion para esta expresion es que
y0
yss = 1, lo cual conlleva a que y0 = yss . Esto quiere decir que cuando la variable log-lineal
ybt se encuentra en el valor cero, esto significa que la variable en niveles yt se encuentra en
su estado estacionario.
De otro lado, en t =1 el valor del producto (log-lineal) es igual a 0.667, en el cual se
cumple: yb1 = 0.667 = ln yyss1 . Resolviendo la segunda igualdad se tiene que yyss1 = e0.667
1 + 0.667. Por tanto, yyss1 = 1 + 0.667, lo cual conlleva finalmente a y1 = (1 + 0.667)yss .
En t = 1 7 yb1 = 0.667 7 y1 = (1 + 0.667)yss

| {z } | {z }
variable log-lineal variable en niveles
Por tanto, el valor (0.667) de la funcion impulso respuesta en t = 1 significa que la

variable producto en niveles (y1 ) esta 66.7 % por encima de su nivel de estado estacionario
(yss ).
38
En la figura [6] y el cuadro [9] se puede observar lo siguiente:
1. En t = 0 (antes de choque) todas las variables permanecen en su estado estacionario.

xss = ln xxss

Por tanto, las variables log-lineales en t = 0 son iguales a cero (b ss
=
ln(1) = 0).
2. En el periodo del choque (t = 1), 1 toma el valor de su desviacion estandar, en este

caso igual a 1.
3. El primer efecto del choque de productividad es un incremento en la funcion de

produccion, la cual incrementa de productividad marginal del capital P M gkt ; es
decir, la demanda del capital en t (Dk ).
4. El incremento de la demanda de capital aumenta la tasa de interes de hoy (R bt ). Esto

se debe a que la oferta del capital es perfectamente inelastica (vertical) porque es
fijada en el periodo anterior b
kt .
5. R
bt produce un Efecto Ingreso (EI): (R kt )
bt b
6. El efecto ingreso incrementa el ct y it
7. it expande kt+1 (oferta del capital de t+1).
8. Lo anterior produce una cada de la tasa de interes en t+1 ( rt+1 ), pero aun esta
por encima de su estado estacionario; es decir, es mas alta que la tasa de interes
antes del choque R b0 , lo cual incentiva a la familia trasladar consumo de hoy t
hacia manana t + 1. Es decir, existe un efecto sustitucion que es gobernado por
la elasticidad de sustitucion del consumo. Para poder ver esta relacion revisemos la
ecuacion de Euler log-lineal:
1b
ct = Et b
b ct+1 Rt+1

Aqu se puede observar que si la tasa de interes de t+1 se incrementa en 1 % entonces
el consumo hoy t se reduce en 1 (elasticidad de sustitucion del consumo). Todo
ello es el efecto sustitucion que produce la tasa de interes.
9. R
bt+1 (pero por encima del estado estado estacionario) produce dos efectos: Efecto
Sustitucion (ES) y Efecto Ingreso (EI).
10. Efecto sustitucion (de la tasa de interes):
R bt b
bt+1 > R ct
11. Efecto ingreso (de la tasa de interes):
R bt b
bt+1 > R ct+1
De la figura [7] se puede concluir algunas ideas. La primera conclusion es que el capital
es mas grande en unidades que cualquier otra variable. Por ejemplo, el valor de estado
estacionario del capital es de 23.88 unidades, la cual es mayor en gran medida con respecto
a las demas variables (el valor de estado estacionario del producto es de 2.87 unidades).
39
Figura 7: Funcion impulso-respuesta de las variables macroeconomicas en niveles

6 50 3.1
3
5.5 45
2.9
5
40 2.8
4.5 2.7
35
4 2.6
30 2.5
3.5
2.4
3 25
2.3
2.5 20 2.2
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200

6 1.045 2.8
1.04 2.6
5
1.035 2.4
4 1.03 2.2
1.025 2
3
1.02 1.8
2 1.015 1.6
1.01 1.4
1
1.005 1.2
0 1 1
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200
Nota: Estas funciones impulso-respuesta corresponde a las variables en niveles; es decir a yt , kt , ct , it , Rt

y at . Cabe mencionar que esta grafica se obtienen del codigo Campbell Lfijo.m (seccion 5).
Para entender por que el stock de capital en estado estacionario es grande se debe de
revisar los parametros de los cuales depende:
1 1
(1 )
kss =
1
Aplicando el signo del exponente se tiene:
1
1
kss = 1
(1 )
Entonces kss es una funcion de , y . En primer lugar, el exponente 1 es mayor a
uno porque es menor a uno (= 0.667), mientras mas pequeno sea mas grande sera el
exponente y mayor sera el numerador incrementando el kss . En segundo lugar, un incre-
mento de la tasa de depreciacion reduce kss , lo cual hace sentido debido a que el capital
se consume a una mayor depreciacion. Por ejemplo, si el capital se deprecia totalmente
( = 1), entonces kss = 0.1880. Finalmente, una mayor tasa de descuento incrementa kss .
40
La segunda conclusion es que la inversion es muy pequena en comparacion con el

capital. Esto se debe a que en estado estacionario la inversion iss es igual a una proporcion
del capital kss , es mas es igual a 2.5 %; es decir, la inversion en estado estacionario (=
0.597) es igual al 2.5 % del capital. Una tercera conlusion es que dado que los valores de la
funcion impulso-respuesta de las variables log-lineales cumplen con la siguiente expresion:
css iss b
ybt = ct +
b it
yss yss
Entonces se puede obtener una relacion entre los niveles de las variables (en la funcion
impulso-respuesta de las variables en niveles):
css iss b
ybt = ct +
b it
yss y
ss
yt css ct iss it
ln = ln + ln
yss yss css yss iss
css iss css iss
ln(yt ) = (ln(yt ) ln(css ) ln(iss )) + ln(ct ) + ln(it ) (107)
yss yss yss yss
De la figura [8] (grafica de la derecha) se desprende una conclusion importante: ante un

choque de productividad, la inversion reacciona fuertemente superando al producto y al
consumo. Es mas, la inversion se incrementa un poco mas del 200 % del valor de su estado
estacionario. Ademas, las variables demoran mas de 100 periodos (trimestres) en volver a
su estado estacionario debido a que el choque tiene una alta persistencia ( = 0.95).
9. Simulacion de las variables endogenas

Para la simulacion del capital usaremos su representacion autorregresiva AR(2):
kt+1 = 1 b
b kt + 2 b
kt1 + ka t
Asumiremos que la variable inicia de su estado estacionario: b

k0 = 0. Ademas, se asume
que la variable en periodos previos se ha mantenido en steady state, entonces: b
k1 = 0.
Cuadro 10: Simluacion del capital log-lineal

t t Representacion AR(2) del kt+1
0 0 = 0 (steady state) k 1 = 1 b
b k0 + 2 b
k1 + ka 0
1 1 = valor aleatorio k2 = 1 k1 + 2 k0 + ka 1
b b b
2 2 = valor aleatorio k 3 = 1 b
b k2 + 2 b
k1 + ka 2
b k3 + 2 b
k2 + ka 3
b k4 + 2 b
k3 + ka 4
Para la simulacion de las variables macroeconomicas como el producto, consumo e

inversion, primero se necesita la serie simulada de la productividad b
at y del capital b
kt , las
cuales se muestran en el cuadro 11. Para esto ultimo se utiliza la solucion del sistema de
ecuaciones log-lineal:
41
Figura 8: Funcion impulso-respuesta (comparacion de las variables log-lineal vs en niveles)
Variables en niveles Variables loglineal

6 2.5
Producto Producto
Consumo Consumo
Inversin Inversin
SSProducto
SSConsumo
5 SSInversin
2
1.5
0.5
1
0 0
0 50 100 150 200 0 50 100 150 200
Nota: Cabe mencionar que esta grafica se obtiene del codigo Campbell Lfijo.m (seccion 5).
at = b
b at1 +
kt+1 = 1 b
b k t + 2 b
kt1 + ka t
Para la simulacion de las demas variables macroeconomicas (b

yt , b
ct , bit y R
bt ) se utiliza
la solucion:
ybt = b at + (1 b kt )
ct = ck kt + ca b
b b at
yss css yss css
bit = (1 ck )b
kt + ( ca )b
at
iss yss iss yss
R
bt = 3 (b at b
kt )
Donde: ik = yiss
ss
(1 ycss ) y ia = yiss
ss ck
ss
( ycss
ss
ca )
La simulacion de las variables pueden tomar valores negativos debido
a que se encuen-
tran expresadas en log-desviaciones de su estado estacionario (ln xxsst ). El valor negativo
42
Cuadro 11: Simulacion de la productividad y del capital (log-lineal)
t t at
b kt+1
b
0 0 = 0 a0
b =0 k1
b =0
0 1 = valor aleatorio de una N(0,1) a1
b = 0.1832 k2
b = 0.0101
b = 0.8557 k3
b = 0.0372
b = 0.1363 k4
b = 0.0291
b = 0.4366 k5
b = 0.0046

Cuadro 12: Simulacion de las variables macroeconomicas log-lineales

t ybt ct
b bit Rbt
0 y0 = a0 + (1 k0 ) c0 = ck k0 + ca a0 i0 = ik k0 + ia a0 R0 = 3 (a0 k0 )
. . . . .
. . . . .
de la variable log-lineal significa que la variable en niveles se encuentra por debajo de

su estado estacionario. Ademas, era de esperar que la variable simulada log-lineal tenga
valores negativos porque su media es igual a cero. Sin embargo, la variable simulada en
niveles tiene solo valores positivos.
10. Componente cclico de las variables simuladas

Para hallar el componente cclico de las variables log-lineales se aplica el filtro Hodrick-
Prescott (Filtro HP). Este filtro permite separar la serie en dos componentes: el compo-
nente tendencial y el componente cclico. La figura [10] muestra el componente cclico y
tendencial para cada variable simulada.
Lo que nos interesa evaluar del modelo son los momentos del componente cclico de
cada variable simulada. Es decir, la varianza, la autocorrelacion y la correlacion con otras
variables.
11. Calculo de los momentos teoricos

Al realizar la simulacion del modelo en Matlab, lo que se obtiene son los momentos
teoricos para las variables log-lineales. Sin embargo, lo que se necesita es volver a las
variables en niveles debido a que los momentos empricos corresponden a las variables
en niveles. Para calcular los momentos teoricos de las variables en niveles se utiliza la
relacion entre la variable log-lineal y la variable en niveles:
43
Figura 9: Simulacion de las variables macroeconomicas log-lineales

1 1 0.5
0 0 0
1 0.5
1
2 1
2
3 1.5
3
4 2
4
5 2.5
5
6 3
6 7 3.5
7 8 4
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200

10 0.25 4
0.2
5 2
0.15
0 0
0.1
5 0.05 2
0
10 4
0.05
15 6
0.1
20 0.15 8
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200
Nota: El capital se comporta como AR(2), mientras que el producto, el consumo, la inversion y la tasa
de interes real se comportan como ARMA(2,1). Ademas, la productividad sigue un proceso AR(1). Este
grafico se obtiene del codigo Campbell Lfijo.m (seccion 6).

xt
x
bt = ln (108)
xss
Media: la media de la variable log-lineal es cero. Luego de los artificios algrebraicos se
concluye la relacion (109), la cual indica que la media del logaritmo de la variable x es
igual al logaritmo de la variable en estado estacionario.
44
Figura 10: Aplicacion del filtro HP a las variables simuladas

2 1 1
1 0
0 1 0
1 2
1
2 3
3 4
2
4 5
5 6 3
6 7
y Tendencia Ciclo k Tendencia Ciclo c Tendencia Ciclo
7 8 4
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200

10 0.25 4
0.2
5 2
0.15
0 0
0.1
5 0.05 2
0
10 4
0.05
15 6
0.1
i Tendencia Ciclo R Tendencia Ciclo a Tendencia Ciclo
20 0.15 8
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200
Nota: Debido a que los parametros calibrados corresponden a datos trimestrales, entonces cada periodo
en esta figura se entiende como un trimestre. Esto sugiere que se debe de usar el parametro de suaviza-
miento del filtro HP correspondiente a datos trimestrales ( = 1600). Este grafico se obtiene del codigo
Campbell Lfijo.m (seccion 7).

xt
x
bt = ln
xss

xt
Eb
xt = Eln
xss
xbt = E(ln(xt ) ln(xss ))
xbt = E(ln(xt )) ln(xss )
xbt = ln(xt ) ln(xss )
0 = ln(xt ) ln(xss )
ln(xt ) = ln(xss ) (109)
Varianza: la varianza de la variable log-lineal es igual a la varianza del logartmo de la

variable en niveles (ver la ecuacion (110)). Por su defecto la desviacion estandar de ambas
45
variables es la misma.

xt
V ar(b
xt ) = V ar ln
xss
Por propiedad V ar(ax + b) = a2 V ar(x)
:

xt ) = V ar ln(xt ) ln(xss )
V ar(b

V ar(b
xt ) = V ar ln(xt )
xb2t 2
= ln(x t)
(110)
Correlacion: la correlacion entre dos variables log-lineales es igual a la correlacion entre

el logaritmo de dichas variables en niveles.
PN
i=1 (b
xi xbi )(b
yi ybi )
corr(b
x, yb) =
xbyb
Pero : xbi = 0 y xbi = 0
PN
i=1 (b
xi )(b
yi )
corr(b
x, yb) =
xbyb
PN xi yi
i=1 (ln xss )(ln yss )
=
xbyb
PN
i=1 (ln(xi ) lnxss )(ln(yi ) lnyss )
=
xbyb
PN
i=1 (ln(xi ) ln(xt ) )(ln(yi ) ln(yt ) )
=
ln(xt ) ln(yt )
corr(b
x, yb) = corr(ln(xt ), ln(yt )) (111)
De lo anterior se concluye que los segundos momentos de las variables log-lineales son
exactamente iguales a los segundos momentos del logartmo natural de las variables en
niveles. Por tanto, podemos concluir que los segundos momentos del componente cclico
de las variables log-lineales (simuladas) son identicas a los segundos momentos del com-
ponente cclico del logartmo de las variables en niveles. Entonces, basta comparar los
momentos (del componente cclico) obtenidos por el modelo (simulacion) con los corres-
pondientes obtenidos de los datos (variables en logartmos).
Para realizar la comparacion se hace lo siguiente. En primer lugar se fija el periodo

de estudio de las variables macroeconomicas (por ejemplo: 1980.I hasta 2016.III). Este
periodo esta compuesto por 147 trimestres. En segundo lugar se realizara 500 simulaciones
considerando que cada variable tenga un periodo de 147 trimestres. En tercer lugar, para
cada simulacion se aplica el filtro HP y se abstrae el componente cclico. En cuarto lugar, se
calcula la media, varianza, autocorrelacion de primer orden y correlacion del componente
cclico para cada simulacion. En quinto lugar, se construye una distribucion para la media,
la varianza, autocorrelacion de primer orden y la correlacion. Se elige el valor medio de cada
distribucion, y ese valor representa el valor de cada segundo momento teorico. Finalmente,
dichos valores se comparan con los datos.
46
Figura 11: Distribuciones de la desviacion estandar del modelo teorico

18 25 18
16 16
14 20 14
12 12
15
10 10
8 8
10
6 6
4 5 4
2 2
0 0 0
0.8 1 1.2 1.4 0.1 0.2 0.3 0.4 0.5 0.25 0.3 0.35 0.4 0.45
Inversin Tasa de inters real Productividad

18 20 18
16 16
14 14
15
12 12
10 10
10
8 8
6 6
5
4 4
2 2
0 0 0
2 2.5 3 3.5 4 0.025 0.03 0.035 0.04 0.045 0.5 1 1.5 2
Nota: Estas distribuciones se obtiene de simular 100 veces las variables considerando un periodo de 150
trimestres. El valor que se toma de cada distribucion es el valor promedio. Este grafico se obtiene del codigo
Campbell Lfijo Sim Variables.m (seccion 7).
12. Comparacion modelo teorico con los datos empiricos

Del cuadro [13] se desprende algunas conclusiones. La primera es que los datos sugieren
que la volatilidad de la inversion es mayor que el producto y el consumo, tal es as que la
desviacion estandar de la inversion es 4.8 veces la volatilidad del producto y 6.4 veces la
del consumo. Sin embargo, el modelo teorico no captura todas estas relaciones encontra-
da en los datos. En primer lugar, aunque mantiene que la volatilidad de la inversion es
mayor que la del consumo, sobreestima la magnitud (la volatilidad de la inversion es 8.2
veces la del consumo). En segundo lugar, la volatilidad para las tres variables (producto,
consumo e inversion) esta muy por debajo de lo observado en los datos. De otro lado,
una bondad del modelo es que captura la direccion correcta en cuanto a volatilidad ():
inversion > producto > consumo .
Una segunda conclusion es que la correlacion temporal del producto con el consumo
y la inversion esta sobrevaluada por el modelo. No obstante, el modelo se acerca mucho
mas a la correlacion producto-inversion que la correpondiente con el consumo.
47
Cuadro 13: Comparacion del comportamiento cclico del modelo teorico con los datos
empiricos
Datos empricos USA Modelo teorico
Variable DesEst ( %) Corr. con producto (t) DesEst ( %) Corr. con producto (t)
Producto 1.72 1 0.85 1
(0.093)
Consumo 1.27 0.83 0.34 0.9874
(0.0382)
Inversion 8.24 0.91 2.82 0.9972
(0.3065)
Nota: Los valores empricos han sido tomado de Cooley y Prescott (1995), los cuales han sido
calculados bajo el periodo muestral de 1954.I hasta 1991.III. Mientras los valores teoricos se han
obtenido de una simulacion de 100 veces considerando un periodo de 150 trimestres. Los valores
mostrados en el modelo teorico son los valores promedios de cada distribucion. Estos valores se
obtiene del codigo Campbell Lfijo Sim Variables.m (seccion 8).
13. Codigos
La solucion del modelo como las funciones impulso-respuesta y la simulacion de las
variables se ha desarrollado directamente en Matlab (mediante la construccion de varios
m-file) y tambien por medio de Dynare (mediante la construccion de un mod-file). El
resultado de ambos caminos es el mismo, la ventaja de construir directamente un m-file
es que se puede hacer explcito muchos detalles en la solucion y simulacion del modelo, los
cuales ya estan programados en Dynare.
48
Cuadro 14: Codigos en Matlab y Dynare

Codigos Descripcion
Matlab
Campbell Lfijo.m Este m-file calcula el estado estacionario,
calcula los coeficientes de la solucion de-
rivados del metodo de coeficientes inde-
terminados. Ademas, calcula la funcion
impulso-respuesta y realiza una simula-
cion de las variables. Finalmente, aplica
el filtro HP para obtener la tendencia y el
ciclo de cada variable simulada
Campbell Lfijo Sim Parametros.m Este m-file simula los valores de las elas-
ticidades (coeficientes de la solucion) en
funcion de la persistencia () y de la ESI
del consumo ().
Campbell Lfijo Sim Variables.m Este m-file simula 100 veces las variables
(forma ARMA(p,q)) para 150 periodos.
Cabe mencionar que ambos parametros
pueden ser cambiados por el usuario pa-
ra realizar distintas simulaciones.
Dynare
Campbell Lfijo Dynare nolineal log.mod Este archivo mod contiene el modelo no-
lineal con variables en logartmo. Este ar-
chivo es el que usa Dynare para resolver
el modelo.
Campbell Lfijo Dynare lineal log.mod Este archivo mod contiene el modelo lineal
con variables en logartmo. Este archivo es
el que usa Dynare para resolver el modelo.
Cabe mencionar que ambos archivos mod
brindan el mismo resultado en cuanto a
las funciones de poltica y de estado.
49

Capitulo4 (RBC) 4

Cargado por

Copyright:

Formatos disponibles

Capitulo4 (RBC) 4

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Capitulo4 (RBC) 4

Cargado por

Copyright:

Formatos disponibles

Capitulo 4

Modelo RBC con trabajo constante*

Borrador: 8 de julio de 2017

2. Construccion del modelo 3

6. Solucion del sistema lineal 24

7. Representacion de series de tiempo 31

Todos los derechos reservados.

9. Simulacion de las variables endogenas 41

10.Componente cclico de las variables simuladas 43

11.Calculo de los momentos teoricos 43

12.Comparacion modelo teorico con los datos empiricos 47

2. Construccion del modelo

Figura 1: Esquema del modelo de oferta de trabajo constante

Donde ct es el consumo del periodo t y es el factor de descuento. Ademas, la funcion

Propiedades de la funcion de utilidad

La funcion de utilidad previa tiene un coeficiente de aversion al riesgo igual a y una

La elasticidad de sustitucion intertemporal del consumo () se entiende como la disposi-

El problema de optimizacion de la familia representativa es el siguiente:

sujeto a la restriccion presupuestaria:

Donde la inversion (it ) se ha reemplazado por su expresion derivada de la ley de movimiento

Donde, de manera similar al capitulo 3 (modelo Lomg y Plosser (1983)), la version

1 u(c1 ) + 1 r1 k1 (c1 + k2 (1 )k1 ) +

2 u(c2 ) + 2 r2 k2 (c2 + k3 (1 )k2 ) +

3 u(c3 ) + 3 r3 k3 (c3 + k4 (1 )k3 ) +

4 u(c4 ) + 4 r4 k4 (c4 + k5 (1 )k4 ) +

t+1 u(ct+1 ) + t+1 rt+1 kt+1 (ct+1 + kt+2 (1 )kt+1 ) +

Las condiciones de primer orden, en el periodo t, son:

t = Et t+1 (rt+1 + (1 )) (6)

uct = Et uct+1 (rt+1 + (1 ))

Rt+1 = rt+1 + (1 ) (8)

uct = Et uct+1 (rt+1 + (1 ))

Sujeto a la funcion de produccion:

Introduciendo la funcion de produccion en la funcion objetivo y derivando esta ultima

2.3. Equilibrio de mercado y definicion del choque

En este modelo en particular, se asume que no existe gasto de gobierno (gt = 0) y

lnat = lnat1 + t (13)

2.4. Ecuaciones principales

Cuadro 1: Sistema de ecuaciones no lineal del modelo

2. La estimacion econometrica de una ecuacion.

3. Referencia a estudios econometricos basados en datos microeconomicos o macro-

Cuadro 2: Calibracion (valores base)

Para la ecuacion de Euler se tiene lo siguiente:

Para la demanda de capital:

De la ecuacion de la tasa de interes bruta:

De la misma manera para la ley de movimiento del capital:

Finalmente para la ecuacion de comportamiento de la productividad:

Hasta aqu se ha encontrado el valor de estado estacionario de la tasa de interes bruta

yss = ass kss

Cuadro 3: Estado estacionario

La ventaja de aplicar log-linealizacion es que convierte el sistema no lineal en lineal, al

En primer lugar, se define la variable en log-desviaciones:

xt = xss exbt (26)

Esta ultima ecuacion se reemplaza en la ecuacion (26):

De la ecuacion (28) se despeja x

Log-linealizando la ecuacion de Euler se tiene:

ebct = Et ebct+1 eRt+1

ebct = Et ebct+1 +Rt+1

yss eybt = ass ebat kss

eybt = ebat +(1)kt

lnat = lnat1 + t (13)