@@ -465,8 +465,8 @@ def compute_reservation_wage(
465465 v = v_next
466466 i += 1
467467
468- res_wage = (1 - β) * (c + β * v @ q)
469- return v, res_wage
468+ w_bar = (1 - β) * (c + β * v @ q)
469+ return v, w_bar
470470```
471471
472472The cell computes the reservation wage at the default parameters
@@ -475,8 +475,8 @@ The cell computes the reservation wage at the default parameters
475475model = McCallModel()
476476c, β, w, q = model
477477v_init = w / (1 - β) # initial guess
478- v, res_wage = compute_reservation_wage(model, v_init)
479- print(res_wage )
478+ v, w_bar = compute_reservation_wage(model, v_init)
479+ print(w_bar )
480480```
481481
482482### Comparative Statics
@@ -569,8 +569,8 @@ def compute_res_wage_jitted(
569569 final_state = jax.lax.while_loop(cond, update, initial_state)
570570 v, i, error = final_state
571571
572- res_wage = (1 - β) * (c + β * v @ q)
573- return v, res_wage
572+ w_bar = (1 - β) * (c + β * v @ q)
573+ return v, w_bar
574574```
575575
576576Now we compute the reservation wage at each $c, \beta$ pair.
@@ -587,9 +587,9 @@ v_init = model.w / (1 - model.β)
587587for i, c in enumerate(c_vals):
588588 for j, β in enumerate(β_vals):
589589 model = McCallModel(c=c, β=β)
590- v, res_wage = compute_res_wage_jitted(model, v_init)
590+ v, w_bar = compute_res_wage_jitted(model, v_init)
591591 v_init = v
592- res_wage_matrix[i, j] = res_wage
592+ res_wage_matrix[i, j] = w_bar
593593
594594fig, ax = plt.subplots()
595595cs1 = ax.contourf(c_vals, β_vals, res_wage_matrix.T, alpha=0.75)
@@ -917,8 +917,8 @@ res_wages_volatility = []
917917for σ in σ_vals:
918918 μ = compute_μ_for_mean(σ, mean_wage)
919919 model = create_mccall_continuous(σ=float(σ), μ=float(μ))
920- res_wage = compute_reservation_wage_continuous(model)
921- res_wages_volatility.append(res_wage )
920+ w_bar = compute_reservation_wage_continuous(model)
921+ res_wages_volatility.append(w_bar )
922922
923923res_wages_volatility = jnp.array(res_wages_volatility)
924924```
0 commit comments