GitHub

@@ -38,6 +38,7 @@ $$ (eq:hugo1)

3838

where $\beta \in (0,1)$ and $u(c)$ is strictly increasing, twice differentiable,

3939

and strictly concave.

404041+4142

We assume that $u(0)$ is well defined.

42434344

We require that $c_t \geq 0$ and $ a_t \geq 0$.

@@ -54,17 +55,27 @@ Furthermore, $a=0$ when the worker is

5455

employed.

55565657

The probability

57-

of finding a job is $p(a)$ where $p$ is an increasing, strictly concave,

58+

of finding a job is $p(a)$.

59+60+

$p$ is an increasing, strictly concave,

5861

and twice differentiable function of $a$ that satisfies

5962

$p(a) \in [0,1]$ for $a \geq 0$, $p(0)=0$.

606364+65+

```{note}

66+

When we compute examples below, we'll use assume the same $p(a)$ function as

67+

{cite}`Hopenhayn_Nicolini_97`, namely, $p(a) = 1 - \exp(- r a)$, where $r$ is a parameter

68+

that we'll calibrate to hit the same target that {cite}`Hopenhayn_Nicolini_97` did, namely,

69+

an empirical hazard rate of leaving unemployment.

70+

```

71+6172

The consumption good is nonstorable.

62736374

An unemployed worker has no savings and cannot borrow or lend.

6475657666-

An **insurance agency** or **planner** is the unemployed worker's only source of consumption

67-

smoothing over time and across states.

77+

The unemployed worker's only source of consumption

78+

smoothing over time and across states is an **insurance agency** or **planner**.

687969807081

Once a worker has found a job, he is beyond the planner's grasp.

@@ -104,7 +115,7 @@ Now let $V^u$ be the expected discounted present value of utility for an

104115

unemployed worker who chooses consumption, effort pair $(c,a)$

105116

optimally.

106117107-

It satisfies the Bellman equation

118+

Value $V^u$ satisfies the Bellman equation

108119109120

$$

110121

V^u = \max_{a \geq 0} \biggl\{ u(0) - a + \beta \left[

@@ -147,21 +158,15 @@ to compute $V^u_{j+1}$.

147158148159

Another benchmark model helps set the stage for the model with private information that we ultimately want to study.

149160150-

In this model, the unemployment agency has full information about the unemployed work.

151-152-

We study optimal provision of insurance with

153-

full information.

161+

We temporarily assume that an unemployment insurance agency has full information about the unemployed worker.

154162155-

An insurance agency can set both

156-

the consumption and search effort of an unemployed person.

163+

We assume that the insurance agency can control both the consumption and the search effort of an unemployed worker.

157164158-

The

159-

agency wants to design an unemployment insurance contract to give

165+

The agency wants to design an unemployment insurance contract to give

160166

the unemployed worker expected discounted utility $V > V_{\rm aut}$.

161167162-

The planner wants to deliver value $V$ efficiently,

163-

meaning in a way that minimizes expected

164-

discounted cost, using $\beta$ as the discount factor.

168+

The agency, i.e., the planner, wants to deliver value $V$ efficiently,

169+

meaning in a way that minimizes an expected present value discounted costs, using $\beta$ as the discount factor.

165170166171

We formulate the optimal insurance problem

167172

recursively.

@@ -177,11 +182,11 @@ only at an increasing marginal cost in terms of the consumption good.

177182178183179184

Given $V$, the planner assigns first-period pair $(c,a)$ and promised

180-

continuation value $V^u$, should the worker be unlucky

181-

and not find a job.

185+

continuation value $V^u$ next period if the worker is unlucky

186+

and does not find a job this period.

182187183-

$(c, a, V^u)$ are chosen to be functions of $V$ and to

184-

satisfy the Bellman equation

188+

The planner sets $(c, a, V^u)$ as functions of $V$ and to

189+

satisfy the following Bellman equation for associated cost function $C(V)$:

185190186191

$$

187192

C(V) = \min_{c, a, V^u} \biggl\{ c + \beta [1 - p(a)] C(V^u) \biggr\} ,

@@ -211,7 +216,7 @@ at least $V$.

211216

Let $\theta$ be a Lagrange multiplier

212217

on constraint {eq}`eq:hugo6`.

213218214-

At an interior solution, the first-order

219+

At an interior solution, first-order

215220

conditions with

216221

respect to $c, a$, and $V^u$, respectively, are

217222

@@ -227,7 +232,7 @@ The envelope condition $C'(V) = \theta$ and the third equation

227232

of {eq}`eq:hugo7` imply that $C'(V^u) =C'(V)$.

228233229234

Strict convexity of $C$ then

230-

implies that $V^u =V$

235+

implies that $V^u =V$.

231236232237

Applied repeatedly over time,

233238

$V^u=V$ makes

@@ -250,11 +255,10 @@ employment and unemployment unless $V=V^e$.

250255251256

### Incentive Problem

252257253-

The preceding efficient insurance scheme requires that the insurance agency

254-

control both $c$ and $a$.

258+

The preceding efficient insurance scheme assumes that the insurance agency

259+

controls both $c$ and $a$.

255260256-

It will not do for the insurance agency

257-

simply to announce $c$ and then allow the worker to choose $a$.

261+

The insurance agency cannot simply provide $c$ and then allow the worker to choose $a$.

258262259263

Here is why.

260264

@@ -264,7 +268,7 @@ the autarky value $V_{\rm aut}$ by doing two things.

264268

It **increases** the unemployed worker's consumption $c$ and **decreases** his search

265269

effort $a$.

266270267-

But the prescribed

271+

The prescribed

268272

search effort is **higher** than what the worker would choose

269273

if he were to be guaranteed consumption level $c$ while he

270274

remains unemployed.

@@ -273,23 +277,22 @@ This follows from the first two equations of {eq}`eq:hugo7` and the

273277

fact that the insurance scheme is costly, $C(V^u)>0$, which imply

274278

$[ \beta p'(a) ]^{-1} > (V^e - V^u)$.

275279276-

But look at the worker's

280+

Now look at the worker's

277281

first-order condition {eq}`eq:hugo4` under autarky.

278282279283

It implies that if search effort $a>0$, then

280284

$[\beta p'(a)]^{-1} = [V^e - V^u]$, which is inconsistent

281-

with the preceding inequality

282-

$[ \beta p'(a) ]^{-1} > (V^e - V^u)$ that prevails when $a >0$ under

283-

the social

284-

insurance arrangement.

285+

with the inequality

286+

$[ \beta p'(a) ]^{-1} > (V^e - V^u)$ that prevails when $a >0$ when the agency controls

287+

both $a$ and $c$.

285288286289

If he were free to choose $a$, the worker would therefore want to

287290

fulfill {eq}`eq:hugo4`, either at equality so long as $a >0$, or by setting

288291

$a=0$ otherwise.

289292290293

Starting from the $a$ associated with

291-

the social insurance scheme,

292-

he would establish the desired equality

294+

the full-information social insurance scheme in which the agency controls both $c$ and $a$,

295+

the worker would establish the desired equality

293296

in {eq}`eq:hugo4` by *lowering* $a$, thereby decreasing

294297

the term $[ \beta p'(a) ]^{-1}$ (which also lowers $(V^e - V^u)$

295298

when the value of being

@@ -303,17 +306,17 @@ payment, set $a=0$, and never work again.

303306304307

Thus, since the worker does not take the

305308

cost of the insurance scheme into account, he would choose a search

306-

effort below the socially optimal one.

309+

effort below the socially optimal, full-information level.

307310308-

The efficient contract

311+

The full-information contract thus

309312

relies on the agency's ability to control *both* the unemployed

310313

worker's consumption *and* his search effort.

311314312315

## Private Information

313316314-

Following Shavell and Weiss (1979) {cite}`Shavell_Weiss_79` and

315-

Hopenhayn and Nicolini (1997) {cite}`Hopenhayn_Nicolini_97`, now assume that the unemployment insurance agency cannot

316-

observe or enforce $a$, though it can observe and control $c$.

317+

Following {cite}`Shavell_Weiss_79` and

318+

{cite}`Hopenhayn_Nicolini_97`, now assume that the unemployment insurance agency cannot

319+

observe or control $a$, though it can observe and control $c$.

317320318321

The worker is free to choose $a$, which puts expression {eq}`eq:hugo4`, the worker's first-order condition under autarky,

319322

back in the picture.

@@ -326,30 +329,30 @@ an instance of the so-called first-order approach to incentive problems.

326329

Given a contract, the individual will choose search effort according to

327330

first-order condition {eq}`eq:hugo4`.

328331329-

This fact leads the insurance agency

330-

to design the unemployment insurance contract to respect this restriction.

332+

This fact motivates the insurance agency

333+

to design an unemployment insurance contract that respects this restriction.

331334332-

Thus, the recursive contract design problem is now to minimize the right side of equation

335+

Thus, the contract design problem is now to minimize the right side of equation

333336

{eq}`eq:hugo5` subject to expression {eq}`eq:hugo6` and the incentive constraint {eq}`eq:hugo4`.

334337335338

Since the restrictions {eq}`eq:hugo4` and {eq}`eq:hugo6` are not linear

336-

and generally do not define a convex set, it becomes difficult

339+

and generally do not define a convex set, it becomes challenging

337340

to provide conditions under which the solution to the dynamic

338341

programming problem results in a convex function $C(V)$.

339342340343

* Sometimes this complication can be handled by convexifying

341-

the constraint set through the introduction of lotteries.

344+

the constraint set by introducing lotteries.

342345

* A common finding is that optimal plans do not involve

343346

lotteries, because convexity of the constraint set is a sufficient

344347

but not necessary condition for convexity of the cost function.

345-

* Following Hopenhayn and Nicolini (1997) {cite}`Hopenhayn_Nicolini_97`, we therefore proceed under the assumption that $C(V)$ is strictly convex in order to characterize the optimal solution.

348+

* In order to characterize the optimal solution, we follow Hopenhayn and Nicolini (1997) {cite}`Hopenhayn_Nicolini_97` by hopefully proceeding under the assumption that $C(V)$ is strictly convex.

346349347350

Let $\eta$ be the multiplier on constraint {eq}`eq:hugo4`, while

348351

$\theta$ continues to denote the multiplier on constraint {eq}`eq:hugo6`.

349352350353

But now we replace the weak inequality in {eq}`eq:hugo6` by an equality.

351354352-

The unemployment insurance agency cannot award a higher utility than

355+

* We do this because the unemployment insurance agency cannot award a higher utility than

353356

$V$ because that might violate an incentive-compatibility constraint

354357

for exerting the proper search effort in earlier periods.

355358

@@ -369,7 +372,7 @@ where the second equality in the second equation in {eq}`eq:hugo8` follows from

369372

of the incentive constraint {eq}`eq:hugo4` when $a>0$.

370373371374

As long as the

372-

insurance scheme is associated with costs, so that $C(V^u)>0$, first-order

375+

insurance scheme is associated with costs, so that $C(V^u)>0$, the first-order

373376

condition in the second equation of {eq}`eq:hugo8` implies that the multiplier $\eta$ is strictly

374377

positive.

375378

@@ -389,10 +392,10 @@ It also follows from {eq}`eq:hugo4` at equality that

389392

search effort $a$ rises as $V^u$ falls, i.e., it rises with the duration

390393

of unemployment.

391394392-

The duration dependence of benefits is designed to provide

393-

incentives to search.

395+

The of benefits on the duration of unemployment is designed to provide the worker an

396+

incentive to search.

394397395-

To see this, from the third equation of {eq}`eq:hugo8`, notice how

398+

To understand this, from the third equation of {eq}`eq:hugo8`, notice how

396399

the conclusion that consumption falls with the duration of

397400

unemployment depends on the assumption that more search effort

398401

raises the prospect of finding a job, i.e., that $p'(a) > 0$.

@@ -422,7 +425,7 @@ The lower bound is

422425

the expected lifetime utility in autarky,

423426

$V_{\rm aut}$.

424427425-

To compute the upper bound,

428+

To compute an upper bound,

426429

represent condition {eq}`eq:hugo4` as

427430428431

$$

@@ -524,7 +527,8 @@ class params_instance:

524527

### Parameter Values

525528526529527-

For the other parameters we have just loaded in the above Python code, we'll set brate the net interest rate $r$ to match the hazard rate -- the probability of finding a job in one period -- in US data.

530+

For the other parameters appearing in the above Python code, we'll calibrate parameter $r$

531+

that pins down the function $p(a) = 1 - \exp(- r a)$ to match an observerd hazard rate -- the probability that an unemployed worker finds a job each -- in US data.

528532529533

In particular, we seek an $r$ so that in autarky `p(a(r)) = 0.1`, where `a` is the optimal search effort.

530534

@@ -533,7 +537,7 @@ In particular, we seek an $r$ so that in autarky `p(a(r)) = 0.1`, where `a` is

533537

First, we create some helper functions.

534538535539

```{code-cell} ipython3

536-

# The probability of finding a job given search effort, a and interest rate r.

540+

# The probability of finding a job given search effort, a and parameter r.

537541

def p(a,r):

538542

return 1-np.exp(-r*a)

539543

@@ -558,15 +562,15 @@ $$

558562

V^u = \max_{a} \{u(0) - a + \beta\left[p_{r}(a)V^e + (1-p_{r}(a))V^u\right]\}

559563

$$ (eq:yad1)

560564561-

At the optimal choice of $a$, we have the first order condition for this problem as:

565+

At the optimal choice of $a$, we have first-order necessary condition:

562566563567

$$

564568

\beta p_{r}'(a)[V^e - V^u] \leq 1

565569

$$ (eq:yad2)

566570567571

with equality when a >0.

568572569-

Given an interest rate $\bar{r}$, we can solve the autarky problem as follows:

573+

Given a value of parameter $\bar{r}$, we can solve the autarky problem as follows:

570574571575

1. Guess $V^u \in \mathbb{R}^{+}$

572576

2. Given $V^u$, use the FOC {eq}`eq:yad2` to calculate the implied optimal search effort $a$

@@ -590,11 +594,11 @@ def Vu_error(self,Vu,r):

590594

return error

591595

```

592596593-

Since the calibration exercise is to match the hazard rate under autarky to the data, we must find an interest rate $r$ to match `p(a,r) = 0.1`.

597+

Since the calibration exercise is to match the hazard rate under autarky to the data, we must find a parameter $r$ to match `p(a,r) = 0.1`.

594598595599

The function below `r_error` calculates, for a given guess of $r$ the difference between the model implied equilibrium hazard rate and 0.1.

596600597-

This will be used to solve for the a calibrated $r^*$.

601+

We'll use this to compute a calibrated $r^*$.

598602599603

```{code-cell} ipython3

600604

# The error of our p(a^*) relative to our calibration target

@@ -622,25 +626,29 @@ To do so, we will use a bisection strategy.

622626623627

```{code-cell} ipython3

624628

r_calibrated = sp.optimize.brentq(r_error_Λ,1e-10,1-1e-10)

625-

print(f"Interest rate to match 0.1 hazard rate: r = {r_calibrated}")

629+

print(f"Parameter to match 0.1 hazard rate: r = {r_calibrated}")

626630627631

Vu_aut = sp.optimize.fsolve(Vu_error_Λ,15000,args = (r_calibrated))[0]

628632

a_aut = invp_prime(1/(params.β*(params.Ve-Vu_aut)),r_calibrated)

629633630634

print(f"Check p at r: {p(a_aut,r_calibrated)}")

631635

```

632636633-

Now that we have calibrated our interest rate $r$, we can continue with solving the model with private information.

637+

Now that we have calibrated our the parameter $r$, we can continue with solving the model with private information.

634638635639

+++

636640637641

### Computation under Private Information

638642639643

+++

640644641-

Our approach to solving the full model is a variant on Judd (1998) {cite}`Judd1998`, who uses a polynomial to approximate the value function and a numerical optimizer to perform the optimization at each iteration.

645+

Our approach to solving the full model follows ideas of Judd (1998) {cite}`Judd1998`, who uses a polynomial to approximate the value function and a numerical optimizer to perform the optimization at each iteration.

646+647+

```{note}

648+

For further details of the Judd (1998) {cite}`Judd1998` method, see {cite}`Ljungqvist2012`, Section 5.7.

649+

```

642650643-

In contrast, we will use cubic splines to interpolate across a pre-set grid of points to approximate the value function. For further details of the Judd (1998) {cite}`Judd1998` method, see {cite}`Ljungqvist2012`, Section 5.7.

651+

We will use cubic splines to interpolate across a pre-set grid of points to approximate the value function.

644652645653

+++

646654

@@ -650,9 +658,9 @@ $$

650658

C(V) = \min_{c,a,V^u} \{c + \beta\left[1-p(a)\right]C(V^u)\}

651659

$$ (eq:yad3)

652660653-

To solve this model, notice that in equations {eq}`eq:hugo21` and {eq}`eq:hugo22`, we have analytical solutions of $c$ and $a$ in terms of (at most) promised value $V$ and $V^u$ (and other parameters).

661+

Notice that in equations {eq}`eq:hugo21` and {eq}`eq:hugo22`, we have analytical solutions of $c$ and $a$ in terms of promised value $V$ and $V^u$ (and other parameters).

654662655-

We can substitute these equations for $c$ and $a$ and obtain the functional equation {eq}`eq:hugo23` that we want to solve.

663+

We can substitute these equations for $c$ and $a$ and obtain the functional equation {eq}`eq:hugo23`.

656664657665658666

```{code-cell} ipython3

@@ -740,7 +748,7 @@ def iterate_C(self,C_old,Vu_grid):

740748

return C_new,V_star,cons_star,a_star

741749

```

742750743-

The below code executes steps 4 and 5 in the Algorithm until convergence to a function $C^*(V)$.

751+

The following code executes steps 4 and 5 in the Algorithm until convergence to a function $C^*(V)$.

744752745753

```{code-cell} ipython3

746754

def solve_incomplete_info_model(self,Vu_grid,Vu_aut,tol = 1e-6,max_iter = 10000):

@@ -767,7 +775,7 @@ def solve_incomplete_info_model(self,Vu_grid,Vu_aut,tol = 1e-6,max_iter = 10000)

767775768776

+++

769777770-

Using the above functions, we create another instance of the parameters with the correctly calibrated interest rate, $r$.

778+

Using the above functions, we create another instance of the parameters with our calibrated parameter $r$.

771779772780

```{code-cell} ipython3

773781

##? Create another instance with the correct r now

@@ -792,7 +800,7 @@ V_star_interp = sp.interpolate.interp1d(Vu_grid,V_star)

792800793801

+++

794802795-

We want to graph the replacement ratio ($c/w$) and search effort $a$ as functions of the duration of unemployment.

803+

Let's graph the replacement ratio ($c/w$) and search effort $a$ as functions of the duration of unemployment.

796804797805

We'll do this for three levels of $V_0$, the lowest being the autarky value $V_{\rm aut}$.

798806

@@ -857,7 +865,7 @@ But for $V^u > V_{\rm aut}$, the planner makes the replacement ratio decline an

857865

### Interpretations

858866859867

The downward slope of the replacement ratio when $V^u > V_{\rm aut}$ is a consequence of the

860-

the planner's limited information about the worker's search effort.

868+

planner's limited information about the worker's search effort.

861869862870

By providing the worker with a duration-dependent schedule of replacement ratios, the planner induces the worker in effect to reveal

863871

his/her search effort to the planner.

@@ -891,7 +899,7 @@ unemployed worker with proper incentives, not to punish an unlucky worker

891899

who has been unemployed for a long time.

892900893901

The planner believes that a worker who has been unemployed a long time is unlucky, not that he has

894-

done anything wrong (i.e., has not lived up to the contract).

902+

done anything wrong (e.g.,that he has not lived up to the contract).

895903896904

Indeed, the

897905

contract is designed to induce the unemployed workers to search in

Read the original on github.com ↗