Bounded Rationality Alters the Dynamics of Paediatric Immunization Acceptance Tamer Oraby,, Chris T. Bauch Department of Mathematics, University of Texas Pan American, Edinburg, Texas, USA, Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada E-mail: orabytf@utpa.edu Appendix I: Table of parameters
Table. Model parameters and symbols with their definition. Family Parameters/ Definition Baseline Ref. Symbols values or ranges e Vaccine efficacy Assumed m Vaccine cost and accessibility. utils Assumed Vaccine V i Vaccine adverse event i (i =,..., (.99,.7,.,. utils Assumed q i Probability of V i s occurrence q = /, q = /4, [] q = / q = (q + q + q Disease β Transmission rate 85 year [] γ Recovery rate 65/ year [] J i Disease adverse event i (i =,..., (.9,,., utils Assumed p i Probability of J i s occurrence p = /8, p = /, [] p = /5, p = (p + p + p Demographic µ Birth/death rate /5 year [] Behavior κ Imitation rate 9 year Assumed c Slope of logit(p 6 Assumed Social δ Θ Social support of group Θ (Θ = V, N. utils Assumed α Θ Utility function exponent of group Θ (Θ = V, N [, 4] Cognitive λ Θ,l Loss aversion of group Θ [, 4] with l = c (commission, o (omission η Θ Parameter of weight function of group Θ [4, 5]
Appendix II: Model equilibria and stability analysis The model The model is given by ds = µ ( xe β S I µ S dt di = β S I (µ + γ I dt dx dt = κ x ( x ( P (π V N (I, x P (π N V (I, x ( where S and I are the proportions of susceptible and infected individuals in the population, µ is the birth/death rate, γ is the recovery rate, e is the vaccine efficacy, β is the disease transmission rate, x is the proportion of vaccinators (and rate of vaccination and κ is the imitation rate. The probability function P is defined by P (z = +exp( c z for some constant c >. The utilities are given by π N (N = π N (V = π V (N = π V (V = ω N ((I + ε p j u N,o (J j B ( j= ω N (q j u N,c (V j B m + j= ω N (( e (I + ε p j u N,o (J j B m ( j= ω V ((I + ε p j u V,o (J j B (4 j= ω V (q j u V,c (V j B m + j= j= ω V (( e (I + ε p j u V,o (J j B m (5 B = is the utility of the healthy status (the reference point; J, J, J, J are the utilities of mild, moderate, morbid, and death due to disease that occur with probabilities p, p, p, p to an infected child such that = B J > J > J > J. The quantities V, V, V, V are the utilities of no-side-effect, mild, morbid, and death due to adverse reaction to vaccination that occur with probabilities q, q, q, q to a vaccinated child such that = B V > V > V > V. The quantity m is the utility (normalized and given in utils of cost and effort to access vaccination.the parental perception of the probability of getting infected is proportional to the current prevalence I offset by a small probability < ε << making a total of I + ε to account for that at zero prevalence there might be a very faint fear
of infection due to immigration. Model equilibria Moving away from the rational decision model with social norms by adding a cognitive dimension to the model results in a plethora of dynamical behaviors of vaccine acceptance, by giving rise to new model equilibria relative to [6]. In addition to vaccine acceptance, incorporating vaccine efficacy and vaccine cost into the model has an effect on the disease dynamics. Social norms can rectify all deviations from the rational decision model of a partially-efficacious costly vaccine. Let M(y := ω V (q j u V,c (V j B m ω N (q j u N,c (V j B m j= + ω V ((y + ε p j u V,o (J j B + ω N ((y + ε p j u N,o (J j B j= j= j= ω V (( e (y + ε p j u V,o (J j B m j= j= ω N (( e (y + ε p j u N,o (J j B m be the simple average of both groups gain/loss of adopting the strategy of not-to-vaccinate when the disease prevalence is y. Here, when ε = (which is almost the case, M( is the average vaccine risk perceived by vaccinators and non-vaccinators. Thus, we will call M( is the average vaccine risk with immigration. The model equation ( has six fixed points. Three fixed points are disease-free equilibria: E = (S, I, x = ( e,,, E = (,,, and E ( e x,, x, where x solves the equation P (π V N (, x = P (π N V (, x or simply π V N (, x = π N V (, x according to the definition of P. Hence x = δ N +δ V (δ N + M( and so the equilibrium point E exists if and only if δ N > M( when M( < or δ V > M( when M( > (or simply, δ V > M( > δ N. Let the effective reproduction number R (x := β ( xe and let the basic reproduction number R := R (. The remaining three fixed points are disease-endemic equilibria: E 4 = ( R, E 5 = ( R, µ R ( R,, and E 6 = ( R, P ( π V N ( µ µ + γ µ R (x 6 R, x 6, where x 6 solves the equation ( ( R (x 6 µ, x 6 = P π N V R µ + γ R (x 6, x 6 R µ R ( R,, 4
which according to the definition of P is just the solution of π V N ( µ µ + γ ( R (x 6 µ, x 6 = π N V R µ + γ R (x 6, x 6 R or the intersection of the line (δ V + δ N x δ N and the function g defined by g(x := M( µ R (x R. Equilibrium point E 4 exists if and only if R ( > (or e < R, whereas E 5 exists if and only if R ( = R >. If < x 6 < then E 6 exists if and only if R (x 6 >. Dynamical regimes The Jacobian matrix of the system of differential equations (, in the main text, is given by J(S, I, x = β I µ β S µ e β I β S (µ + γ J, (I, x J, (I, x (6 where J, (I, x = κ x ( x ( P (I (π V N (I, x P (I (π N V (I, x and J, (I, x = κ x ( x ( P (x (π V N (I, x P (x (π N V (I, x +κ ( x ( P (π V N (I, x P (π N V (I, x The derivatives are given by P (I (π V N (I, x = cp (π V N (I, x( P (π V N (I, x p j ω N((I + ε p j u N,o (J j B+ j= ( e p j ω N (( e (I + ε p j u N,o (J j B m j= P (I (π N V (I, x = c P (π N V (I, x( P (π N V (I, x p j ω V ((I + ε p j u V,o (J j B+ j= ( e p j ω V (( e (I + ε p j u V,o (J j B m j= 5
and P (x (π V N (I, x = c (δ V + δ N P (π V N (I, x( P (π V N (I, x P (x (π N V (I, x = c (δ V + δ N P (π N V (I, x( P (π N V (I, x Equilibrium point E (pure vaccinator, disease-free is stable if the disease is not highly contagious and the vaccine is not very scary compared to the strength of social norms, exactly when R ( < (or e > R and δ V > M(. In contrast, E (non-vaccinator, disease-free is stable when R ( = R < and M( > δ N. The equilibrium point E (partial vaccinator, disease-free is not stable whenever it exists (since J, (, x > : for this equilibrium, social norms will always tend to move vaccine coverage up or down, hence the steady state is unstable. Equilibrium point E 4 (pure vaccinator, disease-endemic is stable if and only if vaccinator pressure is sufficiently large, δ V > M( µ within a particular range: R ( R, and if R ( falls < R γ ( < (7 γ + µ + γ µ + γ Notice that, the left hand side term of inequality (7 is strictly greater than (ensuring existence of the equilibrium point. The two inequalities in (7 can be rewritten in terms of vaccine efficacy (e as (µ + γ + γ (µ + γ µ R < e < (µ + γ γ (µ + γ (8 µ R The right hand side term is strictly less than (exactly by γ + ( γ µ. This region will R µ R appear in the simulation as a narrow gap between the two limits, see Figure (a in the main text. Equilibrium point E 5 (no vaccinator, disease-endemic is stable if and only if non-vaccinator pressure is sufficiently large, M( µ R ( R > δ N and if R ( falls within a certain range: < R γ ( < (9 γ + µ + γ µ + γ Again, the left hand side term is strictly greater than one and the right hand side term is very large that the upper inequality includes all of the the basic reproduction numbers of non-chronic childhood diseases. 6
γ (Note that is very close to one since lengths of incubation period of non-chronic childhood diseases µ + γ is much shorter than average human life length. Finally, equilibrium point E 6 (partial vaccinator, disease-endemic is stable if δ N + δ V < µ e µ µ + γ M ( µ + γ R (x 6 R where M is the derivative of M, and if R is sufficiently large, where P 6 = P ( π V N ( µ R (x 6 R, x 6. R (x 6 > κc(δ N + δ V x 6 ( x 6 P 6 ( P 6 µ Appendix III: Cumulative prospect theory Cumulative prospect theory (CPT amends the violation of stochastic dominance axiom by the prospect theory; c.f. [7]. It can also incorporate a multiple gains and losses [4]. In CPT, weights are given in terms of cumulative probabilities and not the probabilities. First, Θ = (z, p ; z, p ;... ; z k, p k is split into gains and losses say Θ + = (, (q + + q n ; x, q ; x, q ;... ; x n, q n and Θ = (y m, r m ;... ; y, r ; y, r ;, (r + + r m such that < x < x < < x n and y m < < y < y < where n + m = k. The rank-dependent utility is given by n n n m m m π(θ = [ω + ( q l ω + ( q l ] u(x i + [ω ( r l ω ( r l ] u(y j i= l=i l=i+ j= l=j l=j+ with the convention that k l=k+ a l =. The weighting functions ω + and ω have the same form as ω (see equation ( in the main text but need not to have the same parameters value. The new definition of utility works very well with the anomaly case in the PT of a prospect of only two possible outcomes with the same sign. In the model of vaccine acceptance, all the prospect s values are negative; and the prospect is of the kind Θ = (y m, r m ;... ; y, r ; y, r. Thus, π(θ = m j= [ω ( m l=j r l ω ( m l=j+ r l] u(y j. We use an approximation for ω ( m l=j r l ω ( m l=j+ r l that boils down to the form of the weight function of the PT, ω := ω, since ω (r j ω ( m l=j r l ω ( m l=j+ r l for j =,..., m as all r,..., r m 7
are small and ω ( m l=m r l ω ( m l=m+ r l = ω (r m. Appendix IV: The effect of the cognitive parameters To understand the effect of the cognitive parameters at different levels of vaccine cost on vaccine acceptance rates and disease eradication, we split them into three groups: the utility function exponents α V and α N ; the loss aversion indexes λ V,c, λ V,o, λ N,c and λ N,o ; and the weighting parameters η V and η N. We discuss each group separately and set the parameters of the other two groups equal to one. At vaccine efficacy e =.95 and for all the three groups, the more the vaccine is to cost, the larger the vaccinators group pressure is needed to have he same stability region that secures the possibility of full vaccine acceptance on condition of large initial acceptance. Moreover, the disease will be eradicated at that vaccine efficacy (Figures A, A and A. Compare those figures to Figures A4, A5 and A6 illustrated for vaccine efficacy e =.9. In the later cases, the full vaccine acceptance without diseases eradication E 4 appears in the place of E with slightly smaller group pressure needed to cover the same region in the parameter plane. Cognitive parameters, seemingly, at the fixed vaccintors group pressure can lead to full vaccine refusal (Figures A (a, A (a and A (a which warrant a careful consideration of subsidies and rewards when dealing with human perception to gain, loss and weighting of probabilities. As the cost of vaccination increases, a larger group pressure is required to maintain the same area of bistability in Figures A (a, A (a, A (a, A4 (a, A5 (a and A6 (a, compare panel (b in Figures A, A, and A to their respective panels (c and (d as well as in Figures A4, A5, and A6. In that case a large enough initial vaccine acceptance is required to reach a full vaccine acceptance and eradicate the disease. 8
(a.. (c.9.7..... (b.. (d...9.7....9.7... Figure A. (a α V α N plane at m = and δ V =.. The equilibrium points E (pure vaccinator, disease-free and E 5 (no vaccinator, disease-endemic are stable in the red region which includes the point (, corresponding to the rational decision model, given the values of the rest of the cognitive parameters below whereas in the blue region the only stable point is E 5. Contour plots for values of M( calculated at each pair of (α V, α N at (b m =, (c m =., and (d m =.. The full vaccination equilibrium point E is stable in the regions where δ V is larger than the value of M(. R ( R Thus, (a follows from (b at δ V =.. A similar contour plot of M( µ shows that it is always positive so E 5 is always stable. The rest of the parameters are κ = 9, c = 6, e =.95, η V =, η N =, λ V,c = λ V,o =, λ N,c = λ N,o =. (not shown here 9
(a (b Non vaccinators loss aversion (λ N,c λ N,c.. Vaccinators loss aversion (λ V,c λ V,c (c Non vaccinators loss aversion (λ N,c 5 5 5. (d Non vaccinators loss aversion (λ N,c..9.7 Vaccinators loss aversion (λ V,c Vaccinators loss aversion (λ V,c Figure A. (a λ V,c λ N,c plane with λ V,o = λ V,c and λ N,o = λ N,c at m = and δ V =.. The equilibrium points E (pure vaccinator, disease-free and E 5 (no vaccinator, disease-endemic are stable in the red region which includes the point (, corresponding to the rational decision model, given the values of the rest of the cognitive parameters below whereas in the blue region the only stable point is E 5. Contour plots for values of M( calculated at each pair of (λ V,c, λ N,c with λ V,o = λ V,c and λ N,o = λ N,c at (b m =, (c m =., and (d m =.. The full-vaccination equilibrium point E is stable in the regions where δ V is larger than the value of M(. Thus, (a follows from (b at δ V =.. A similar contour plot of M( µ R ( R (not shown here shows that it is always positive so E 5 is always stable. The rest of the parameters are κ = 9, c = 6, e =.95, α V =, α N =, η V =, η N =.
(a (b... (c (d...... Figure A. (a η V η N plane at m = and δ V =.. The equilibrium points E (pure vaccinator, disease-free and E 5 (no vaccinator, disease-endemic are stable in the red region which includes the point (, corresponding to the rational decision model, given the values of the rest of the cognitive parameters below whereas in the blue region the only stable point is E 5. Contour plot for values of M( calculated at each pair of (η V, η N at (b m =, (c m =., and (d m =.. The full vaccination equilibrium point E is stable in the regions where δ V is larger than the value of M(. R ( R Thus, (a follows from (b at δ V =.. A similar contour plot of M( µ shows that it is always positive so E 5 is always stable. The rest of the parameters are κ = 9, c = 6, e =.95, α V =, α N =, λ V,c = λ V,o =, λ N,c = λ N,o =. (not shown here
(a.. (b...9.7... (c (d...9.7......9.7... Figure A4. (a α V α N plane at m = and δ V =.. The equilibrium points E 4 (pure vaccinator, disease-endemic and E 5 (no vaccinator, disease-endemic are stable in the red region which includes the point (, corresponding to the rational decision model, given the values of the rest of the cognitive parameters below whereas in the blue region the only stable point is E 5. Contour plots for values of M( µ R ( calculated at each pair of (α V, α N at (b m =, (c m =., and (d m =.. The M( µ R full vaccination equilibrium point E 4 is stable in the regions where δ V is larger than the value of R ( R. Thus, (a follows from (b at δ V =.. A similar contour plot of M( µ R ( R (not shown here shows that it is always positive so E 5 is always stable. The rest of the parameters are κ = 9, c = 6, e =.9, η V =, η N =, λ V,c = λ V,o =, λ N,c = λ N,o =.
(a (b λ N,c Non vaccinators loss aversion (λ N,c.. λ V,c Vaccinators loss aversion (λ V,c (c Non vaccinators loss aversion (λ N,c 5 5 5. (d Non vaccinators loss aversion (λ N,c..9.7 Vaccinators loss aversion (λ V,c Vaccinators loss aversion (λ V,c Figure A5. (a λ V,c λ N,c plane with λ V,o = λ V,c and λ N,o = λ N,c at m = and δ V =.. The equilibrium points E 4 (pure vaccinator, disease-endemic and E 5 (no vaccinator, disease-endemic are stable in the red region which includes the point (, corresponding to the rational decision model, given the values of the rest of the cognitive parameters below whereas in the blue region the only stable point is E 5. Contour plots for values of M( µ R ( R calculated at each pair of (λ V,c, λ N,c with λ V,o = λ V,c and λ N,o = λ N,c at (b m =, (c m =., and (d m =.. The full-vaccination equilibrium point E 4 is stable in the regions where δ V is larger than the value of M( µ R ( R. R ( R Thus, (a follows from (b at δ V =.. A similar contour plot of M( µ shows that it is always positive so E 5 is always stable. The rest of the parameters are κ = 9, c = 6, e =.9, α V =, α N =, η V =, η N =. (not shown here
(a (b.. (c (d.... Figure A6. (a η V η N plane at m = and δ V =.. The equilibrium points E 4 (pure vaccinator, disease-endemic and E 5 (no vaccinator, disease-endemic are stable in the red region which includes the point (, corresponding to the rational decision model, given the values of the rest of the cognitive parameters below whereas in the blue region the only stable point is E 5. Contour plot for values of M( µ R ( calculated at each pair of (η V, η N at (b m =, (c m =., and (d m =.. The M( µ R full vaccination equilibrium point E 4 is stable in the regions where δ V is larger than the value of R ( R. Thus, (a follows from (b at δ V =.. A similar contour plot of M( µ R ( R (not shown here shows that it is always positive so E 5 is always stable. The rest of the parameters are κ = 9, c = 6, e =.9, α V =, α N =, λ V,c = λ V,o =, λ N,c = λ N,o =. References. Centers for Disease Control and Prevention, Vaccines: Vac-Gen/Some Misconceptions - Risk from disease versus risk from vaccines. ( Available at: http://www.cdc.gov/vaccines/vac- 4
gen/6mishome.htm. (Accessed: 5th March 4. Bauch, C. T. & Bhattacharyya, S. Evolutionary game theory and social learning can determine how vaccine scares unfold. PLoS Comput. Biol. 8, e45 (.. Trepel, C., Fox, C. R. & Poldrack, R. A. Prospect theory on the brain? Toward a cognitive neuroscience of decision under risk. Cognitive Brain Res., 4 5 (5. 4. Tversky, A. & Kahneman, D. Advances in prospect theory: Cumulative representation of uncertainty. J. Risk Uncertainty 5, 97 (99. 5. Wu, G. & Gonzalez, R. Curvature of the probability weighting function. Manage. Sci. 4, 676 69 (996. 6. Oraby, T., Thampi, V. & Bauch, C. T. The influence of social norms on the dynamics of vaccinating behaviour for paediatric infectious diseases. Proc. R. Soc. B: Behaviour 8, 7 (4. 7. Fennema, H. & Wakker, P. P. Original and cumulative prospect theory: A discussion of empirical differences. J. Behav. Decis. Making, 5 64 (997. 5