@@ -394,9 +394,9 @@ def compute_V_fre(L0_arr, L1_arr, π_star, wf):
394394π_star = 0.5
395395V_fre_arr, PFA_arr, PD_arr = compute_V_fre(L0_arr, L1_arr, π_star, wf)
396396
397- plt.plot(range(T), V_fre_arr, label='$\min_{d} \overline{V}_{fre}(t,d)$')
397+ plt.plot(range(T), V_fre_arr, label=r '$\min_{d} \overline{V}_{fre}(t,d)$')
398398plt.xlabel('t')
399- plt.title('$\pi^*=0.5$')
399+ plt.title(r '$\pi^*=0.5$')
400400plt.legend()
401401plt.show()
402402```
@@ -434,8 +434,8 @@ for i, π_star in enumerate(π_star_arr):
434434
435435``` {code-cell} python3
436436plt.plot(π_star_arr, V_fre_bar_arr)
437- plt.xlabel('$\pi^*$')
438- plt.title('$\overline{V}_{fre}$')
437+ plt.xlabel(r '$\pi^*$')
438+ plt.title(r '$\overline{V}_{fre}$')
439439
440440plt.show()
441441```
@@ -447,12 +447,12 @@ $\left(PFA,PD\right)$ change as $\pi^{*}$ varies.
447447fig, axs = plt.subplots(1, 2, figsize=(14, 5))
448448
449449axs[0].plot(π_star_arr, t_optimal_arr)
450- axs[0].set_xlabel('$\pi^*$')
450+ axs[0].set_xlabel(r '$\pi^*$')
451451axs[0].set_title('optimal sample size given $\pi^*$')
452452
453453axs[1].plot(π_star_arr, PFA_optimal_arr, label='$PFA^*(\pi^*)$')
454454axs[1].plot(π_star_arr, PD_optimal_arr, label='$PD^*(\pi^*)$')
455- axs[1].set_xlabel('$\pi^*$')
455+ axs[1].set_xlabel(r '$\pi^*$')
456456axs[1].legend()
457457axs[1].set_title('optimal PFA and PD given $\pi^*$')
458458
@@ -587,7 +587,7 @@ plt.vlines(β, 0, β * wf.L0, linestyle="--")
587587plt.vlines(α, 0, (1 - α) * wf.L1, linestyle="--")
588588
589589ax.set(xlim=(0, 1), ylim=(0, 0.5 * max(wf.L0, wf.L1)), ylabel="cost",
590- xlabel="$\pi$", title="Value function")
590+ xlabel=r "$\pi$", title="Value function")
591591
592592plt.legend(borderpad=1.1)
593593plt.show()
@@ -689,7 +689,7 @@ plt.vlines(β, 0, wf.L0, linestyle='--')
689689plt.text(β+0.01, wf.L0/2, 'β')
690690plt.vlines(α, 0, wf.L0, linestyle='--')
691691plt.text(α+0.01, wf.L0/2, 'α')
692- plt.xlabel('$\pi$')
692+ plt.xlabel(r '$\pi$')
693693plt.title('Objective value function $V(\pi)$')
694694plt.legend()
695695plt.show()
@@ -734,11 +734,11 @@ for i, π_star in enumerate(π_star_arr):
734734 axs[row_i, col_i].vlines(π_optimal, V_baye_bar, V_baye.max(), linestyle='--')
735735 axs[row_i, col_i].text(π_optimal+0.05, (V_baye_bar + V_baye.max()) / 2,
736736 '${\pi_0^*}=$'+f'{π_optimal:0.2f}')
737- axs[row_i, col_i].set_xlabel('$\pi$')
738- axs[row_i, col_i].set_ylabel('$\overline{V}_{baye}(\pi)$')
739- axs[row_i, col_i].set_title('$\pi^*=$' + f'{π_star}')
737+ axs[row_i, col_i].set_xlabel(r '$\pi$')
738+ axs[row_i, col_i].set_ylabel(r '$\overline{V}_{baye}(\pi)$')
739+ axs[row_i, col_i].set_title(r '$\pi^*=$' + f'{π_star}')
740740
741- fig.suptitle('$\overline{V}_{baye}(\pi)=\pi^*V^0(\pi) + (1-\pi^*)V^1(\pi)$', fontsize=16)
741+ fig.suptitle(r '$\overline{V}_{baye}(\pi)=\pi^*V^0(\pi) + (1-\pi^*)V^1(\pi)$', fontsize=16)
742742plt.show()
743743```
744744
@@ -763,14 +763,14 @@ for i, π_star in enumerate(π_star_arr):
763763fig, axs = plt.subplots(1, 2, figsize=(14, 5))
764764
765765axs[0].plot(π_star_arr, V_baye_bar_arr)
766- axs[0].set_xlabel('$\pi^*$')
767- axs[0].set_title('$\overline{V}_{baye}$')
766+ axs[0].set_xlabel(r '$\pi^*$')
767+ axs[0].set_title(r '$\overline{V}_{baye}$')
768768
769769axs[1].plot(π_star_arr, π_optimal_arr, label='optimal prior')
770770axs[1].plot([π_star_arr.min(), π_star_arr.max()],
771771 [π_star_arr.min(), π_star_arr.max()],
772772 c='k', linestyle='--', label='45 degree line')
773- axs[1].set_xlabel('$\pi^*$')
773+ axs[1].set_xlabel(r '$\pi^*$')
774774axs[1].set_title('optimal prior given $\pi^*$')
775775axs[1].legend()
776776
@@ -799,10 +799,10 @@ V_baye_bar = V_baye.min()
799799```
800800
801801``` {code-cell} python3
802- plt.plot(range(T), V_fre_arr, label='$\min_{d} \overline{V}_{fre}(t,d)$')
803- plt.plot([0, T], [V_baye_bar, V_baye_bar], label='$\overline{V}_{baye}$')
802+ plt.plot(range(T), V_fre_arr, label=r '$\min_{d} \overline{V}_{fre}(t,d)$')
803+ plt.plot([0, T], [V_baye_bar, V_baye_bar], label=r '$\overline{V}_{baye}$')
804804plt.xlabel('t')
805- plt.title('$\pi^*=0.5$')
805+ plt.title(r '$\pi^*=0.5$')
806806plt.legend()
807807plt.show()
808808```
@@ -816,14 +816,14 @@ rule does better on average for all values of $\pi^{*}$.
816816``` {code-cell} python3
817817fig, axs = plt.subplots(1, 2, figsize=(14, 5))
818818
819- axs[0].plot(π_star_arr, V_fre_bar_arr, label='$\overline{V}_{fre}$')
820- axs[0].plot(π_star_arr, V_baye_bar_arr, label='$\overline{V}_{baye}$')
819+ axs[0].plot(π_star_arr, V_fre_bar_arr, label=r '$\overline{V}_{fre}$')
820+ axs[0].plot(π_star_arr, V_baye_bar_arr, label=r '$\overline{V}_{baye}$')
821821axs[0].legend()
822- axs[0].set_xlabel('$\pi^*$')
822+ axs[0].set_xlabel(r '$\pi^*$')
823823
824824axs[1].plot(π_star_arr, V_fre_bar_arr - V_baye_bar_arr, label='$diff$')
825825axs[1].legend()
826- axs[1].set_xlabel('$\pi^*$')
826+ axs[1].set_xlabel(r '$\pi^*$')
827827
828828plt.show()
829829```
0 commit comments