11. Present Values#
11.1. Overview#
This lecture describes the present value model that is a starting point of much asset pricing theory.
Asset pricing theory is a component of theories about many economic decisions including
consumption
labor supply
education choice
demand for money
In asset pricing theory, and in economic dynamics more generally, a basic topic is the relationship among different time series.
A time series is a sequence indexed by time.
In this lecture, we’ll represent a sequence as a vector.
So our analysis will typically boil down to studying relationships among vectors.
Our main tools in this lecture will be
matrix multiplication, and
matrix inversion.
We’ll use the calculations described here in subsequent lectures, including consumption smoothing, equalizing difference model, and monetarist theory of price levels.
Let’s dive in.
11.2. Analysis#
Let
\(\{d_t\}_{t=0}^T \) be a sequence of dividends or “payouts”
\(\{p_t\}_{t=0}^T \) be a sequence of prices of a claim on the continuation of the asset’s payout stream from date \(t\) on, namely, \(\{d_s\}_{s=t}^T \)
\( \delta \in (0,1) \) be a one-period “discount factor”
\(p_{T+1}^*\) be a terminal price of the asset at time \(T+1\)
We assume that the dividend stream \(\{d_t\}_{t=0}^T \) and the terminal price \(p_{T+1}^*\) are both exogenous.
This means that they are determined outside the model.
Assume the sequence of asset pricing equations
We say equations, plural, because there are \(T+1\) equations, one for each \(t =0, 1, \ldots, T\).
Equations (11.1) assert that price paid to purchase the asset at time \(t\) equals the payout \(d_t\) plus the price at time \(t+1\) multiplied by a time discount factor \(\delta\).
Discounting tomorrow’s price by multiplying it by \(\delta\) accounts for the “value of waiting one period”.
We want to solve the system of \(T+1\) equations (11.1) for the asset price sequence \(\{p_t\}_{t=0}^T \) as a function of the dividend sequence \(\{d_t\}_{t=0}^T \) and the exogenous terminal price \(p_{T+1}^*\).
A system of equations like (11.1) is an example of a linear difference equation.
There are powerful mathematical methods available for solving such systems and they are well worth studying in their own right, being the foundation for the analysis of many interesting economic models.
For an example, see Samuelson multiplier-accelerator
In this lecture, we’ll solve system (11.1) using matrix multiplication and matrix inversion, basic tools from linear algebra introduced in linear equations and matrix algebra.
We will use the following imports
import numpy as np
import matplotlib.pyplot as plt
11.3. Representing sequences as vectors#
The equations in system (11.1) can be arranged as follows:
Write the system (11.2) of \(T+1\) asset pricing equations as the single matrix equation
Carry out the matrix multiplication in (11.3) by hand and confirm that you recover the equations in (11.2).
In vector-matrix notation, we can write system (11.3) as
Here \(A\) is the matrix on the left side of equation (11.3), while
The solution for the vector of prices is
For example, suppose that the dividend stream is
Let’s write Python code to compute and plot the dividend stream.
T = 6
current_d = 1.0
d = []
for t in range(T+1):
d.append(current_d)
current_d = current_d * 1.05
fig, ax = plt.subplots()
ax.plot(d, 'o', label='dividends')
ax.legend()
ax.set_xlabel('time')
plt.show()
Now let’s compute and plot the asset price.
We set \(\delta\) and \(p_{T+1}^*\) to
δ = 0.99
p_star = 10.0
Let’s build the matrix \(A\)
A = np.zeros((T+1, T+1))
for i in range(T+1):
for j in range(T+1):
if i == j:
A[i, j] = 1
if j < T:
A[i, j+1] = -δ
Let’s inspect \(A\)
A
array([[ 1. , -0.99, 0. , 0. , 0. , 0. , 0. ],
[ 0. , 1. , -0.99, 0. , 0. , 0. , 0. ],
[ 0. , 0. , 1. , -0.99, 0. , 0. , 0. ],
[ 0. , 0. , 0. , 1. , -0.99, 0. , 0. ],
[ 0. , 0. , 0. , 0. , 1. , -0.99, 0. ],
[ 0. , 0. , 0. , 0. , 0. , 1. , -0.99],
[ 0. , 0. , 0. , 0. , 0. , 0. , 1. ]])
Now let’s solve for prices using (11.5).
b = np.zeros(T+1)
b[-1] = δ * p_star
p = np.linalg.solve(A, d + b)
fig, ax = plt.subplots()
ax.plot(p, 'o', label='asset price')
ax.legend()
ax.set_xlabel('time')
plt.show()
Now let’s consider a cyclically growing dividend sequence:
T = 100
current_d = 1.0
d = []
for t in range(T+1):
d.append(current_d)
current_d = current_d * 1.01 + 0.1 * np.sin(t)
fig, ax = plt.subplots()
ax.plot(d, 'o-', ms=4, alpha=0.8, label='dividends')
ax.legend()
ax.set_xlabel('time')
plt.show()
Compute the corresponding asset price sequence when \(p^*_{T+1} = 0\) and \(\delta = 0.98\).
Solution to Exercise 11.2
We proceed as above after modifying parameters and consequently the matrix \(A\).
δ = 0.98
p_star = 0.0
A = np.zeros((T+1, T+1))
for i in range(T+1):
for j in range(T+1):
if i == j:
A[i, j] = 1
if j < T:
A[i, j+1] = -δ
b = np.zeros(T+1)
b[-1] = δ * p_star
p = np.linalg.solve(A, d + b)
fig, ax = plt.subplots()
ax.plot(p, 'o-', ms=4, alpha=0.8, label='asset price')
ax.legend()
ax.set_xlabel('time')
plt.show()
The weighted averaging associated with the present value calculation largely eliminates the cycles.
11.4. Analytical expressions#
By the inverse matrix theorem, a matrix \(B\) is the inverse of \(A\) whenever \(A B\) is the identity.
It can be verified that the inverse of the matrix \(A\) in (11.3) is
Check this by showing that \(A A^{-1}\) is equal to the identity matrix.
If we use the expression (11.6) in (11.5) and perform the indicated matrix multiplication, we shall find that
Pricing formula (11.7) asserts that two components sum to the asset price \(p_t\):
a fundamental component \(\sum_{s=t}^T \delta^{s-t} d_s\) that equals the discounted present value of prospective dividends
a bubble component \(\delta^{T+1-t} p_{T+1}^*\)
The fundamental component is pinned down by the discount factor \(\delta\) and the payout of the asset (in this case, dividends).
The bubble component is the part of the price that is not pinned down by fundamentals.
It is sometimes convenient to rewrite the bubble component as
where
11.5. More about bubbles#
For a few moments, let’s focus on the special case of an asset that never pays dividends, in which case
In this case system (11.1) of our \(T+1\) asset pricing equations takes the form of the single matrix equation
Evidently, if \(p_{T+1}^* = 0\), a price vector \(p\) of all entries zero solves this equation and the only the fundamental component of our pricing formula (11.7) is present.
But let’s activate the bubble component by setting
for some positive constant \(c\).
In this case, when we multiply both sides of (11.8) by the matrix \(A^{-1}\) presented in equation (11.6), we find that
11.6. Gross rate of return#
Define the gross rate of return on holding the asset from period \(t\) to period \(t+1\) as
Substituting equation (11.10) into equation (11.11) confirms that an asset whose sole source of value is a bubble earns a gross rate of return
11.7. Exercises#
Give analytical expressions for an asset price \(p_t\) under the following settings for \(d\) and \(p_{T+1}^*\):
\(p_{T+1}^* = 0, d_t = g^t d_0\) (a modified version of the Gordon growth formula)
\(p_{T+1}^* = g^{T+1} d_0, d_t = g^t d_0\) (the plain vanilla Gordon growth formula)
\(p_{T+1}^* = 0, d_t = 0\) (price of a worthless stock)
\(p_{T+1}^* = c \delta^{-(T+1)}, d_t = 0\) (price of a pure bubble stock)
Solution to Exercise 11.4
Plugging each pair of the above \(p_{T+1}^*, d_t\) into Equation (11.7) yields:
\(p_t = \sum^T_{s=t} \delta^{s-t} g^s d_0\)
\(p_t = \sum^T_{s=t} \delta^{s-t} g^s d_0 + \delta^{T+1-t} g^{T+1} d_0\)
\(p_t = 0\)
\(p_t = c \delta^{-t}\)