Kenji Sato
2016-12-21
The master equation:
\[ \dot k = s f(k) - (\delta + g + n) k \]
Important variables in the steady state \( k^* \):
\[ \frac{K}{L} = A k^* \]
\[ \frac{Y}{L} = Af(k^*) \]
both grow at the rate of \( g \).
Let \( k = k^* \).
The common growth rate for
\[ K/L,\quad Y/L,\quad C/L \] is \[ g = \dot A / A. \]
The common growth rate for \[ K,\quad Y,\quad C \] is \[ g + n = \dot A/A + \dot L/L. \]
The situation in which important variables share growth rate is called balanced growth.
When \( k=k^* \), the economy is on the balanced growth path.
Comparative Statics/Dynamics is a common excercise of macroeconomics.
It is important to understand what happens after a (small) parameter change.
What happens after an increase of the saving rate?
\( s \) is an important policy variable for the government.
may have impact on \( s \).
library(ggplot2)
library(ggthemes)
s0 = 0.3
s1 = 0.4
alpha = 0.3
delta = 0.05
g = 0.02
n = 0.01
f = function(k) {
return(k^alpha)
}
k = seq(0.0, 25.0, by=0.01)
df = data.frame(k=k, f=f(k), s0f=s0*f(k), s1f=s1*f(k))
head(df)
k f s0f s1f
1 0.00 0.0000000 0.00000000 0.0000000
2 0.01 0.2511886 0.07535659 0.1004755
3 0.02 0.3092495 0.09277485 0.1236998
4 0.03 0.3492500 0.10477499 0.1397000
5 0.04 0.3807308 0.11421924 0.1522923
6 0.05 0.4070905 0.12212716 0.1628362
fig = ggplot(df) +
geom_line(aes(x=k, y=f)) + # Production Function
geom_line(aes(x=k, y=s0f), color='blue', size=1.5) + # For s0
geom_line(aes(x=k, y=s1f), color='red', size=1.5) + # For s1
geom_line(aes(x=k, y=(delta+g+n)*k)) # Break-Even
\( k^* \) is larger when \( s \) gets larger.
What about growth rate?
k0 = (s0 / (g + n + delta))^(1 / (1 - alpha)) # steady state
t0 = 10 # Change of policy
solow_update = function(t, k){
if (t < t0){
return(k + dt * (s0 * f(k) - (delta + g + n) * k))
} else {
return(k + dt * (s1 * f(k) - (delta + g + n) * k))
}
}
dt = 0.01 # controls precision of approximation
t_max = 100 # simulation for t_max years
t = seq(from=0, to=t_max, by=dt)
simulation = as.data.frame(t)
simulation[1, "k"] = k0
for (i in 2:nrow(simulation)){
simulation[i, "k"] = solow_update(simulation[i-1, "t"],
simulation[i-1, "k"])
}
ggplot(simulation, aes(x=t, y=k)) + geom_line() + theme_gdocs()
At \( t = t_0 \), \( k \) starts to increase and
it stops increasing when it attains the new steady state value.
Let \( A(0) = 1 \). Plot using a logarithmic scale for y-Axis.
simulation$KL = simulation$k * exp(g*simulation$t)
ggplot(simulation, aes(x=t, y=KL)) + geom_line() + scale_y_log10() + theme_gdocs()
Exercise: Reproduce the following graph.
Note that the growth rate, \( g_{K/L} \) of \( K/L \) satisfies
\[ g_{K/L}(t) = g + g_k(t), \]
where \( g \) is the exogenous growth rate of \( A \), \( g_k \) is the growth rate of \( k \).
After an increase in saving rate, we get \( \dot k > 0 \) and thereby \( g_k(t) > 0 \).
On the transition path, the growth of per capita capital is faster than on the BGP.
It seems to be consistent with observations about NICs.
Note that
\[ \begin{aligned} c^* &= (1 - s) f(k^*) \\ &= f(k^*) - (\delta + g + n) k^* \end{aligned} \]
When \( c^* \) is maximized, we should have (think of \( c^* \) as a function of \( k^* \))
\[ f'(k^*) = \delta + g + n \]
Let \( k^*_G \) be the unique stock level that satisfies the above equation. Golden rule capital stock.
Exercise: Reproduce the following graph.
Observe that the new steady state value, \( c^* \), for \( c = C/(AL) \) is smaller after the parameter change considered above.
\( Y = F(K, AL) \) implies that
\[ \begin{aligned} \dot Y &= \frac{\partial Y}{\partial K} \dot K + \frac{\partial Y}{\partial (AL)} \frac{d}{dt}(AL)\\ &= \frac{\partial Y}{\partial K} \dot K + \left[\frac{\partial Y}{\partial (AL)} A\right] \dot L + \left[\frac{\partial Y}{\partial (AL)} L\right] \dot A\\ &=: \frac{\partial Y}{\partial K} \dot K + \frac{\partial Y}{\partial L} \dot L + \left[\frac{\partial Y}{\partial (AL)} L\right] \dot A. \end{aligned} \]
We therefore have
\[ \frac{\dot Y}{Y} = \frac{\partial Y}{\partial K} \frac{K}{Y} \frac{\dot K}{K} + \frac{\partial Y}{\partial L} \frac{L}{Y} \frac{\dot L}{L} + \left[\frac{\partial Y}{\partial (AL)} \frac{AL}{Y} \right] g. \]
Define
\[ \begin{aligned} \alpha_K (t) &:= \frac{K}{Y} \frac{\partial Y}{\partial K} & (\text{capital elasticity of output}) \\ \alpha_L (t) &:= \frac{L}{Y} \frac{\partial Y}{\partial L} & (\text{labor elasticity of output}) \end{aligned} \]
Veryfy for the Cobb–Douglas family that they are constant: \( \alpha_K = \alpha \) and \( \alpha_L = 1-\alpha \).
By Euler's theorem on CRS functions,
\[ \alpha_K (t) + \alpha_L (t) = 1 \]
1% increase in capital input results in \( \alpha_K \)% increase in output.
\[ \alpha_K = \frac{K}{Y} \frac{\partial Y}{\partial K} \simeq \dfrac{ \dfrac{Y + \Delta Y}{Y} }{ \dfrac{K + \Delta K}{K} } \]
By employing this notation, the decomposition of \( \dot Y/Y \) becomes
\[ \frac{\dot Y}{Y} = \alpha_K \frac{\dot K}{K} + \alpha_L \frac{\dot L}{L} + R, \]
where
\[ R := \left[\frac{\partial Y}{\partial (AL)} \frac{AL}{Y} \right] g = \alpha_L g, \]
called the Solow residual.
All terms other than \( R \) can be obtained from data.
Equivalently,
\[ g_{Y/L} = \alpha_K g_{K/L} + R = \alpha_K g_{K/L} + a_L g. \]
In the steady state, \( \alpha_K \) fraction of growth in output per worker is attributable to capital accumulation. The rest is due to the technological progress.
After extended to incorporate human capital accumulation, the Solow model fits fairly well with data. See Mankiw, Romer and Weil (1992, QJE).