# Job Search II: Search and Separation¶

In addition to what’s in Anaconda, this lecture will need the following libraries:

```
!pip install --upgrade quantecon
```

## Overview¶

Previously we looked at the McCall job search model [McC70] as a way of understanding unemployment and worker decisions.

One unrealistic feature of the model is that every job is permanent.

In this lecture, we extend the McCall model by introducing job separation.

Once separation enters the picture, the agent comes to view

- the loss of a job as a capital loss, and
- a spell of unemployment as an
*investment*in searching for an acceptable job

The other minor addition is that a utility function will be included to make worker preferences slightly more sophisticated.

We’ll need the following imports

```
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
from numba import njit, jitclass, float64
from quantecon.distributions import BetaBinomial
```

## The Model¶

The model is similar to the baseline McCall job search model.

It concerns the life of an infinitely lived worker and

- the opportunities he or she (let’s say he to save one character) has to work at different wages
- exogenous events that destroy his current job
- his decision making process while unemployed

The worker can be in one of two states: employed or unemployed.

He wants to maximize

$$ {\mathbb E} \sum_{t=0}^\infty \beta^t u(y_t) \tag{1} $$

At this stage the only difference from the baseline model is that we’ve added some flexibility to preferences by introducing a utility function $ u $.

It satisfies $ u'> 0 $ and $ u'' < 0 $.

### The Wage Process¶

For now we will drop the separation of state process and wage process that we maintained for the baseline model.

In particular, we simply suppose that wage offers $ \{ w_t \} $ are IID with common distribution $ q $.

The set of possible wage values is denoted by $ \mathbb W $.

(Later we will go back to having a separate state process $ \{s_t\} $ driving random outcomes, since this formulation is usually convenient in more sophisticated models.)

### Timing and Decisions¶

At the start of each period, the agent can be either

- unemployed or
- employed at some existing wage level $ w_e $.

At the start of a given period, the current wage offer $ w_t $ is observed.

If currently *employed*, the worker

- receives utility $ u(w_e) $ and
- is fired with some (small) probability $ \alpha $.

If currently *unemployed*, the worker either accepts or rejects the current offer $ w_t $.

If he accepts, then he begins work immediately at wage $ w_t $.

If he rejects, then he receives unemployment compensation $ c $.

The process then repeats.

(Note: we do not allow for job search while employed—this topic is taken up in a later lecture)

## Solving the Model¶

We drop time subscripts in what follows and primes denote next period values.

Let

- $ v(w_e) $ be total lifetime value accruing to a worker who enters the current period
*employed*with existing wage $ w_e $ - $ h(w) $ be total lifetime value accruing to a worker who who enters the current period
*unemployed*and receives wage offer $ w $.

Here *value* means the value of the objective function (1) when the worker makes optimal decisions at all future points in time.

Our first aim is to obtain these functions.

### The Bellman Equations¶

Suppose for now that the worker can calculate the functions $ v $ and $ h $ and use them in his decision making.

Then $ v $ and $ h $ should satisfy

$$ v(w_e) = u(w_e) + \beta \left[ (1-\alpha)v(w_e) + \alpha \sum_{w' \in \mathbb W} h(w') q(w') \right] \tag{2} $$

and

$$ h(w) = \max \left\{ v(w), \, u(c) + \beta \sum_{w' \in \mathbb W} h(w') q(w') \right\} \tag{3} $$

Equation (2) expresses the value of being employed at wage $ w_e $ in terms of

- current reward $ u(w_e) $ plus
- discounted expected reward tomorrow, given the $ \alpha $ probability of being fired

Equation (3) expresses the value of being unemployed with offer $ w $ in hand as a maximum over the value of two options: accept or reject the current offer.

Accepting transitions the worker to employment and hence yields reward $ v(w) $.

Rejecting leads to unemployment compensation and unemployment tomorrow.

Equations (2) and (3) are the Bellman equations for this model.

They provide enough information to solve for both $ v $ and $ h $.

### A Simplifying Transformation¶

Rather than jumping straight into solving these equations, let’s see if we can simplify them somewhat.

(This process will be analogous to our second pass at the plain vanilla McCall model, where we simplified the Bellman equation.)

First, let

$$ d := \sum_{w' \in \mathbb W} h(w') q(w') \tag{4} $$

be the expected value of unemployment tomorrow.

We can now write (3) as

$$ h(w) = \max \left\{ v(w), \, u(c) + \beta d \right\} $$or, shifting time forward one period

$$ \sum_{w' \in \mathbb W} h(w') q(w') = \sum_{w' \in \mathbb W} \max \left\{ v(w'), \, u(c) + \beta d \right\} q(w') $$Using (4) again now gives

$$ d = \sum_{w' \in \mathbb W} \max \left\{ v(w'), \, u(c) + \beta d \right\} q(w') \tag{5} $$

Finally, (2) can now be rewritten as

$$ v(w) = u(w) + \beta \left[ (1-\alpha)v(w) + \alpha d \right] \tag{6} $$

In the last expression, we wrote $ w_e $ as $ w $ to make the notation simpler.

### The Reservation Wage¶

Suppose we can use (5) and (6) to solve for $ d $ and $ v $.

(We will do this soon.)

We can then determine optimal behavior for the worker.

From (3), we see that an unemployed agent accepts current offer $ w $ if $ v(w) \geq u(c) + \beta d $.

This means precisely that the value of accepting is higher than the expected value of rejecting.

It is clear that $ v $ is (at least weakly) increasing in $ w $, since the agent is never made worse off by a higher wage offer.

Hence, we can express the optimal choice as accepting wage offer $ w $ if and only if

$$ w \geq \bar w \quad \text{where} \quad \bar w \text{ solves } v(\bar w) = u(c) + \beta d $$### Solving the Bellman Equations¶

We’ll use the same iterative approach to solving the Bellman equations that we adopted in the first job search lecture.

Here this amounts to

- make guesses for $ d $ and $ v $
- plug these guesses into the right-hand sides of (5) and (6)
- update the left-hand sides from this rule and then repeat

In other words, we are iterating using the rules

$$ d_{n+1} = \sum_{w' \in \mathbb W} \max \left\{ v_n(w'), \, u(c) + \beta d_n \right\} q(w') \tag{7} $$

$$ v_{n+1}(w) = u(w) + \beta \left[ (1-\alpha)v_n(w) + \alpha d_n \right] \tag{8} $$

starting from some initial conditions $ d_0, v_0 $.

As before, the system always converges to the true solutions—in this case, the $ v $ and $ d $ that solve (5) and (6).

(A proof can be obtained via the Banach contraction mapping theorem.)

## Implementation¶

Let’s implement this iterative process.

In the code, you’ll see that we use a class to store the various parameters and other objects associated with a given model.

This helps to tidy up the code and provides an object that’s easy to pass to functions.

The default utility function is a CRRA utility function

```
@njit
def u(c, σ=2.0):
return (c**(1 - σ) - 1) / (1 - σ)
```

Also, here’s a default wage distribution, based around the BetaBinomial distribution:

```
n = 60 # n possible outcomes for w
w_default = np.linspace(10, 20, n) # wages between 10 and 20
a, b = 600, 400 # shape parameters
dist = BetaBinomial(n-1, a, b)
q_default = dist.pdf()
```

Here’s our jitted class for the McCall model with separation.

```
mccall_data = [
('α', float64), # job separation rate
('β', float64), # discount factor
('c', float64), # unemployment compensation
('w', float64[:]), # list of wage values
('q', float64[:]) # pmf of random variable w
]
@jitclass(mccall_data)
class McCallModel:
"""
Stores the parameters and functions associated with a given model.
"""
def __init__(self, α=0.2, β=0.98, c=6.0, w=w_default, q=q_default):
self.α, self.β, self.c, self.w, self.q = α, β, c, w, q
def update(self, v, d):
α, β, c, w, q = self.α, self.β, self.c, self.w, self.q
v_new = np.empty_like(v)
for i in range(len(w)):
v_new[i] = u(w[i]) + β * ((1 - α) * v[i] + α * d)
d_new = np.sum(np.maximum(v, u(c) + β * d) * q)
return v_new, d_new
```

Now we iterate until successive realizations are closer together than some small tolerance level.

We then return the current iterate as an approximate solution.

```
@njit
def solve_model(mcm, tol=1e-5, max_iter=2000):
"""
Iterates to convergence on the Bellman equations
* mcm is an instance of McCallModel
"""
v = np.ones_like(mcm.w) # Initial guess of v
d = 1 # Initial guess of d
i = 0
error = tol + 1
while error > tol and i < max_iter:
v_new, d_new = mcm.update(v, d)
error_1 = np.max(np.abs(v_new - v))
error_2 = np.abs(d_new - d)
error = max(error_1, error_2)
v = v_new
d = d_new
i += 1
return v, d
```

### The Reservation Wage: First Pass¶

The optimal choice of the agent is summarized by the reservation wage.

As discussed above, the reservation wage is the $ \bar w $ that solves $ v(\bar w) = h $ where $ h := u(c) + \beta d $ is the continuation value.

Let’s compare $ v $ and $ h $ to see what they look like.

We’ll use the default parameterizations found in the code above.

```
mcm = McCallModel()
v, d = solve_model(mcm)
h = u(mcm.c) + mcm.β * d
fig, ax = plt.subplots()
ax.plot(mcm.w, v, 'b-', lw=2, alpha=0.7, label='$v$')
ax.plot(mcm.w, [h] * len(mcm.w),
'g-', lw=2, alpha=0.7, label='$h$')
ax.set_xlim(min(mcm.w), max(mcm.w))
ax.legend()
plt.show()
```

The value $ v $ is increasing because higher $ w $ generates a higher wage flow conditional on staying employed.

### The Reservation Wage: Computation¶

Here’s a function `compute_reservation_wage`

that takes an instance of `McCallModel`

and returns the associated reservation wage.

```
@njit
def compute_reservation_wage(mcm):
"""
Computes the reservation wage of an instance of the McCall model
by finding the smallest w such that v(w) >= h.
If no such w exists, then w_bar is set to np.inf.
"""
v, d = solve_model(mcm)
h = u(mcm.c) + mcm.β * d
w_bar = np.inf
for i, wage in enumerate(mcm.w):
if v[i] > h:
w_bar = wage
break
return w_bar
```

Next we will investigate how the reservation wage varies with parameters.

## Impact of Parameters¶

In each instance below, we’ll show you a figure and then ask you to reproduce it in the exercises.

### The Reservation Wage and Unemployment Compensation¶

First, let’s look at how $ \bar w $ varies with unemployment compensation.

In the figure below, we use the default parameters in the `McCallModel`

class, apart from
c (which takes the values given on the horizontal axis)

As expected, higher unemployment compensation causes the worker to hold out for higher wages.

In effect, the cost of continuing job search is reduced.

### The Reservation Wage and Discounting¶

Next, let’s investigate how $ \bar w $ varies with the discount factor.

The next figure plots the reservation wage associated with different values of $ \beta $

Again, the results are intuitive: More patient workers will hold out for higher wages.

### The Reservation Wage and Job Destruction¶

Finally, let’s look at how $ \bar w $ varies with the job separation rate $ \alpha $.

Higher $ \alpha $ translates to a greater chance that a worker will face termination in each period once employed.

Once more, the results are in line with our intuition.

If the separation rate is high, then the benefit of holding out for a higher wage falls.

Hence the reservation wage is lower.

## Exercises¶

### Exercise 1¶

Reproduce all the reservation wage figures shown above.

Regarding the values on the horizontal axis, use

```
grid_size = 25
c_vals = np.linspace(2, 12, grid_size) # unemployment compensation
beta_vals = np.linspace(0.8, 0.99, grid_size) # discount factors
alpha_vals = np.linspace(0.05, 0.5, grid_size) # separation rate
```

## Solutions¶

### Exercise 1¶

Here’s the first figure.

```
mcm = McCallModel()
w_bar_vals = np.empty_like(c_vals)
fig, ax = plt.subplots()
for i, c in enumerate(c_vals):
mcm.c = c
w_bar = compute_reservation_wage(mcm)
w_bar_vals[i] = w_bar
ax.set(xlabel='unemployment compensation',
ylabel='reservation wage')
ax.plot(c_vals, w_bar_vals, label=r'$\bar w$ as a function of $c$')
ax.legend()
plt.show()
```

Here’s the second one.

```
fig, ax = plt.subplots()
for i, β in enumerate(beta_vals):
mcm.β = β
w_bar = compute_reservation_wage(mcm)
w_bar_vals[i] = w_bar
ax.set(xlabel='discount factor', ylabel='reservation wage')
ax.plot(beta_vals, w_bar_vals, label=r'$\bar w$ as a function of $\beta$')
ax.legend()
plt.show()
```

Here’s the third.

```
fig, ax = plt.subplots()
for i, α in enumerate(alpha_vals):
mcm.α = α
w_bar = compute_reservation_wage(mcm)
w_bar_vals[i] = w_bar
ax.set(xlabel='separation rate', ylabel='reservation wage')
ax.plot(alpha_vals, w_bar_vals, label=r'$\bar w$ as a function of $\alpha$')
ax.legend()
plt.show()
```