Optimal Implementation of Intervention Strategies for Elderly People with Ludomania
Article information
Osong Public Health Res Perspect. 2014;5(5):266-273
Department of Mathematics, Kyungpook National University, Daegu, Korea
∗Corresponding author.
yongkuk@knu.ac.kr
Received 2014 July 21; Revised 2014 August 18; Accepted 2014 August 25.
Abstract
Objectives
Now-a-days gambling is growing especially fast among older adults. To control the gratuitous growth of gambling, well-analyzed scientific strategies are necessary. We tried to analyze the adequacy of the health of society mathematically through immediate treatment of patients with early prevention.
Methods
The model from Lee and Do was modified and control parameters were introduced. Pontryagin's Maximum Principle was used to obtain an optimal control strategy.
Results
Optimal control can be achieved through simultaneous use of the control parameters, though it varies from society to society. The control corresponding to prevention needed to be implemented in full almost all the time for all types of societies. In the case of the other two controls, the scenario was greatly affected depending on the types of societies.
Conclusion
Prevention and treatment for elderly people with ludomania are the main intervention strategies. We found that optimal timely implementation of the intervention strategies was more effective. The optimal control strategy varied with the initial number of gamblers. However, three intervention strategies were considered, among which, preventing people from engaging in all types of gambling proved to be the most crucial.
Problem gambling or ludomania is a type of disorder that consists of an urge to continuously gamble despite harmful negative consequences or a desire to stop and that is associated with both social and family costs. Problem gambling and wider gambling-related harm constitute a significant health and social issue [1]. To study problems associated with gambling, Shaffer and Korn [2] used the classic public health model for communicable disease, which examines the interaction among host, agent, environment, and vector. Moreover, some sociologists [3–6] have shown that a significant predictor of the occurrence of ludomania is peer pressure; in the sense that the occurrence depends on the number of individuals involved, the number of individuals who might be involved, as well as the frequency, duration, priority, and intensity of association with peers. Therefore, ludomania might be considered as a contagious disease. Recently, from the point of view of a communicable disease, Lee and Do [7,8] used a mathematical modeling approach to study the dynamics of problem gambling.
In this study, we adopted the optimal control theory to their model and tried to find optimal strategies for intervention. A variety of policies and services have been developed with the intent of preventing and reducing problem gambling and related harm. The prevalence and consequences of problem gambling as well as approaches to treatment can be found in the book by Petry [9]. We considered a basic model [7] to incorporate some important epidemiological features, such as time-dependent control functions. The extended model can then be used to determine cost-effective strategies for combating the spread of problem gambling in a given population; a mathematical modeling approach to study the dynamics of problem gambling.
2.1 Basic model
We considered the model of Lee and Do [7] without demographic effect as follows:
(1)dSdt=−αSL+PNdLdt=αSL+PN−βLPN−ϕL+ψPdPdt=βLPN+ϕL−ψP−γPHN−θP+τHdHdt=γPHN+θP−τH
The whole population N(t)=S(t)+L(t)+P(t)+H(t) consisting of older adults aged 65–80 years was divided into four classes: susceptible S(t), latent gamblers L(t), pathological gamblers P(t), and treated gamblers H(t). The susceptible population S(t) was a class of individuals who had never gambled more than five times in a single year in their life time. Using the per capita transition rate α, susceptible people entered the second compartment L(t), which was composed of individuals who gamble frequently but had two or less symptoms of problem gambling in the previous year. The transition rate α might be understood as the peer pressure from people in L(t) and P(t). Latent people were pathological gamblers, with the peer pressure transition rate β from people in P(t), or with the natural progression rate ϕ. The class of excessive gamblers P(t) consisted of problem and pathological gamblers. When problem or pathological gamblers sought help, they transited to class H(t) of individuals who were in treatment, with the peer pressure rate γ from people inH(t), or with the voluntary transition rate θ. By attending several types of psychotherapy, including Gamblers Anonymous, cognitive behavioral therapy, behavioral therapy, psychodynamic therapy, and family therapy [10], people in H(t) may have returned to P(t) with the transition rate τ. The rate τ was closely related to the efficacy of a cognitive–behavioral treatment package for pathological gambling [11].
2.2 Optimal control
Using sensitivity analysis, Lee and Do [7] showed that the best way to reduce gambling problems among elderly people is to minimize the value of α, which is similar to the claim of Shaffer and Korn [2] that primary prevention is most important. We considered three interventions to reduce gambling problems among elderly people: reducing α and β, and urging the pathological gamblers to take medical services, which resulted in increasing θ. Although we may have gained some insights from such constant controlling of the parameters, it is unrealistic to have constant controls to α, β, and θ over time. The goal was to show that it was possible to implement time-dependent control techniques while minimizing the cost of implementation of such control measures.
We formulated an optimal control problem for the transmission dynamics of gambling by adding control terms to the basic model (1) as follows:
(2)dSdt=−α(1−u1(t))SL+PNdLdt=α(1−u1(t))SL+PN−β(1−u2(t))LPN−ϕL+ψPdPdt=β(1−u2(t))LPN+ϕL−ψP−γPHN−(θ+ρu3(t))P+τHdHdt=γPHN+(θ+ρu3(t))P−τH
Here, we noted that N(t)=S(t)+L(t)+P(t)+H(t) was constant.
The control variables u1(t), u2(t), and u3(t) represent the amount of intervention related to the parameters α, β, and θ at time t, respectively. The factor of 1−u1(t) and 1−u2(t) reduced the per capita transition rate α from S to L and β from L toP, respectively. It was also assumed that the per capita transition rate θ from P to H increased at a rate proportional to u3(t); where ρ>0 was a rate constant.
We defined our control set to be:
U={(u1(t),u2(t),u3(t)):ui(t)isLebesguemeasurableon[0,T],0≤ui(t)≤1,i=1,2,3}.
An optimal control problem with the objective cost functional can be given by
(3)J(u1,u2,u3)=∫0T(ALL(t)+APP(t)+B12u12(t)+B22u22(t)+B32u32(t))dt
subject to the state system given by (2). In the objective cost functional, the quantities
AL,AP,B1,B2 and
B3 represented the weight constants. The costs associated with controls of transition rates were described in the terms
B1u12(t),B2u22(t) and
B3u32(t). The goal was to minimize the populations
L(t) and
P(t) of problem gamblers and the cost of implementing the controls. Therefore, optimal control functions
(u1∗,u2∗,u3∗) needed to be found such that:
(4)J(u1∗,u2∗,u3∗)=min{J(u1,u2,u3):(u1,u2,u3)∈U}
subject to the system of equations given by (2). In order to find an optimal solution, first we should define the Hamiltonian function
H for the problems (2) and (3), and then use Pontryagin's Maximum Principle
[12] to derive the characterization for an optimal control. The principle converts (2) and (3) into a problem of minimizing point wise a Hamiltonian,
H, with respect to
u1,
u2 and
u3. The Hamiltonian for our problem was the integrand of the objective functional coupled with the four right-hand sides of the state equations, where
X(t)=(S(t),L(t),P(t),H(t)),u(t)=(u1(t),u2(t),u3(t))
and
Λ(t)=(λ1(t),λ2(t),λ3(t),λ4(t)). Then
(5)H(X(t),u(t),Λ(t))=ALL(t)+APP(t)+B12u12(t)+B22u22(t)+B32u32(t)+Λ(t)(dX(t)dt)T
λ1′(t)=(λ1(t)−λ2(t))α(1−u1∗(t))L∗(t)+P∗(t)N(t)λ2′(t)=−AL+(λ1(t)−λ2(t))α(1−u1∗(t))S∗(t)N∗(t)+(λ2(t)−λ3(t))(β(1−u2∗(t))P∗(t)N∗(t)+ϕ)λ3′(t)=−AP+(λ1(t)−λ2(t))α(1−u1∗(t))S∗(t)N∗(t)+(λ2(t)−λ3(t))(β(1−u2∗(t))L∗(t)N∗(t)−ψ)+(λ3(t)−λ4(t))(γH∗(t)N∗(t)+θ+ρu3∗(t))λ4′(t)=(λ3(t)−λ4(t))(γP∗(t)N∗(t)−τ)
Let S∗(t),L∗(t),P∗(t),H∗(t) be optimal state solutions with associated optimal control variables u1∗(t),u2∗(t) and u3∗(t) for the optimal control problem (2) and (3). Then, there were adjoint variables λ1(t),λ2(t),λ3(t),λ4(t) that satisfied
(6)=ALL(t)+APP(t)+B12u12(t)+B22u22(t)+B32u32(t)+λ1(t)(−α(1−u1(t))SL(t)+P(t)N(t))+λ2(t)(α(1−u1(t))SL(t)+P(t)N(t)−β(1−u2(t))L(t)P(t)N(t)−ϕL(t)+ψP(t))+λ3(t)(β(1−u2(t))L(t)P(t)N(t)+ϕL(t)−ψP(t)−γP(t)H(t)N(t)−(θ+ρu3(t))P(t)+τH(t))+λ4(t)(γP(t)H(t)N(t)+(θ+ρu3(t))P(t)−τH(t))
with the transversality condition (or boundary condition)
(7)λj(T)=0,j=1,2,3,4.
Furthermore, the optimal controls u1∗(t),u2∗(t) and u3∗(t) were given by
(8)u1∗(t)=min{1,max{0,1B1(αS∗(L∗+P∗)(λ2−λ1)N)}}u2∗(t)=min{1,max{0,1B2(βL∗P∗(λ3−λ2)N)}}u3∗(t)=min{1,max{0,ρP∗(λ3−λ4)B3}}
(See Appendix 1 for the detailed derivation)
For numerical simulation, the forward–backward sweep method [13] based on 4th order Runge–Kutta algorithm was used to treat the problem. The problem consisted of eight ordinary differential equations describing states and adjoint variables along with three controls. The parameters in were adopted from a previous study [7] and used for our simulation. A natural shortcoming was that the controls were not 100% effective, so the upper boundary of the controls u1, u2, and u3 was chosen to be 0.6. The rate constant ρ for u3 was chosen to be 1.
Parameter values for the model.
depicts the variation of the maximum for the three controls, which illustrates that the control u2 became useless for AL≥AP, and the control u2 came into action only for AP>AL.
Variation of maximum controls subject to social structure. AL = 1, AP = 20, B1 = 5000, B2 = 500, and B3 = 50,000.
depicts the numerical simulation that was carried out in the time interval [0,20](years) with initial conditions S(0)=68500,L(0)=21000,P(0)=9000,H(0)=1500, so that N(0)=100000 with the weight values AL=1,AP=20,B1=5000,B2=500,B3=50000. The solid lines in the four graphs on the left show populations in different compartments in the absence of control efforts and the dotted line shows the states with implementation of the optimal controls. These graphs reveal the impact of control by the reduced number of gamblers and pathological gamblers, and increased number of susceptible gamblers. The rightmost graph shows the control profile, which says that we need a full three controls almost all the time.
Optimal control scenario with AL = 1, AP = 20, B1 = 5000, B2 = 500, and B3 = 50,000.
On the other hand, the control scenario would not be similar in all societies. The control scenario might be greatly affected by the number of gamblers and pathological gamblers, that is to say, the control scenario may vary depending on the initial conditions. To analyze the effect of the number of gamblers in society, keeping the total population unchanged, we varied the total number of gamblers and pathological gamblers from 5% to 35%, among which gamblers and pathological gamblers were in the ratio 7:3, and the proportion of treated gamblers was 5% of the total gambling population. Simulation results have been plotted in , which illustrates that the control u1 is implemented in full for almost all the time in all types of societies. In the case of u2 and u3, the scenario was more dramatic. Both of the controls had maximum implementation for a long time in a highly-gambling society only. As the percentage of gamblers fell, maximum implementation of u2 shrank gradually. However, in the case of u3, it reduced slowly up to ∼20%, after which it fell abruptly. For u2 and u3, if the gambling populations were <11% and <10%, respectively, maximum implementation was not necessary at all. However, in societies with a low percentage of gamblers u3 is used more than u2.
Variation of optimal controls subject to social structure for B1 = 5000, B2 = 500 and B3 = 50,000.
An optimal control strategy was analyzed with the help of Pontryagin's Maximum Principle for three control factors. The control scenario would not be similar in all societies. The control scenario might have been affected by its impact on society, and the impact of gamblers on society was introduced into the model through the coefficients AL and AP in the cost functional. The inequality AP>AL means that pathological gamblers are more detrimental than gamblers. The controls u1 and u3 were not affected by this, but the control u2 showed an important response to AL,AP, and the control u2 became useless for AL≥AP, and the control u2 came into action only for AP>AL. According to , we need a full three controls during almost all the time. However, the control scenario might also be greatly affected by the total number of gamblers and pathological gamblers in society. If the control u1 needs to be implemented in full for almost all the time for all types of societies, in the case of u2 and u3, the scenario is greatly dependent on types of societies.
In conclusion, it was conspicuous that simultaneous implementation of all the controls gave the most effective result. However, the control u1 corresponding to peer pressure on the susceptible gamblers was more crucial than the control u2 corresponding to peer pressure on the gamblers and u3 to pressure towards an urge for medical services. In addition, for AL≥AP, the control u2 became totally ineffective. Therefore, strategies should be taken to keep people away not only from problem gambling, but rather from gambling altogether.
All contributing authors declare no conflicts of interest.
References
1. Eades J.M.. Gambling addiction: the problem, the pain and the path to recovery. 2003. Vine Books.
2. Shaffer H.J., Korn D.A.. Gambling and related mental disorders: a public health analysis. Ann Rev Pub Health 232002;:171–212.
11910060.
3. Crane J.. The epidemic theory of ghettos and neighborhood effects on dropping out and teenage childbearing. Am J Sociol 96(5)1991 Mar;:1226–1259.
4. McGuire W.J.. Inducing resistance to persuasion. Adv Soc Psychol 11964;:191–229.
5. Sullivan C.J.. Early adolescent delinquency: assessing the role of childhood problems, family environment, and peer pressure. Natl Criminal Justice Ref Service 4(4)2006 Oct;:291–313.
6. Voss H.L.. Differential association and containment theory: a theoretical convergence. Social Forces 47(4)1969 Jun;:381–391.
7. Lee Y.S., Do T.S.. A mathematical modeling approach to gambling among older adults. Appl Math Comput 221(2)2013 Sep;:403–410.
8. Do T.S., Lee Y.S.. A differential equation model for the dynamics of youth gambling. Osong Public Health Res Perspect 5(4)2014;:1–9.
24955305.
9. Petry N.M.. Pathological gambling: etiology, comorbidity, and treatment. 2005. American Psychological Association. Washington, DC:
10. Fong T.W.. Types of psychotherapy for pathological gamblers. Psychiatry (Edgmont) 2(5)2005 May;:32–39.
21152147.
11. Sylvain C., Ladouceur R., Boisvert J.M.. Cognitive and behavioral treatment of pathological gambling: a controlled study. J Consult Clin Psychol 65(5)1997 Oct;:727–732.
9337491.
12. Agrachev A., Sachkov Y.. Control Theory from the geometric viewpoint. 2004. Springer. Berlin: p. 165.
13. Lenhart S., Workman J.T.. Optimal control applied to biological models. 2007. Chapman and Hall/CRC. p. 49.
Theorem 1
Let S∗(t),L∗(t),P∗(t),H∗(t) be optimal state solutions with associated optimal control variables u1∗(t),u2∗(t), and u3∗(t) for the optimal control problem (2) and (3). Then, there were adjoint variables λ1(t),λ2(t),λ3(t),λ4(t) that satisfied
λ1′(t)=(λ1(t)−λ2(t))α(1−u1∗(t))L∗(t)+P∗(t)N(t)λ2′(t)=−AL+(λ1(t)−λ2(t))α(1−u1∗(t))S∗(t)N∗(t)+(λ2(t)−λ3(t))(β(1−u2∗(t))P∗(t)N∗(t)+ϕ)λ3′(t)=−AP+(λ1(t)−λ2(t))α(1−u1∗(t))S∗(t)N∗(t)+(λ2(t)−λ3(t))(β(1−u2∗(t))L∗(t)N∗(t)−ψ)+(λ3(t)−λ4(t))(γH∗(t)N∗(t)+θ+ρu3∗(t))λ4′(t)=(λ3(t)−λ4(t))(γP∗(t)N∗(t)−τ)
with the transversality condition (or boundary condition)
(7)λj(T)=0,j=1,2,3,4.
Furthermore, the optimal controls u1∗(t),u2∗(t) and u3∗(t) were given by
(8)u1∗(t)=min{1,max{0,1B1(αS∗(L∗+P∗)(λ2−λ1)N)}}u2∗(t)=min{1,max{0,1B2(βL∗P∗(λ3−λ2)N)}}u3∗(t)=min{1,max{0,ρP∗(λ3−λ4)B3}}
Proof
To determine the adjoint equations and the transversality conditions, we used the Hamiltonian (6). By Pontryagin's Maximum Principle [12], setting S(t)=S∗(t),L(t)=L∗(t),P(t)=P∗(t),H(t)=H∗(t) and also differentiating the Hamiltonian (6) with respect to S(t),L(t),P(t),H(t), we obtained:
λ1′(t)=−∂H∂S=λ1(t)α(1−u1∗(t))L∗(t)+P∗(t)N(t)−λ2(t)α(1−u1∗(t))L∗(t)+P∗(t)N(t)λ2′(t)=−∂H∂L=−AL+λ1(t)α(1−u1∗(t))S∗(t)N∗(t)−λ2(t)α(1−u1∗(t))S∗(t)N∗(t)+λ2(t)β(1−u2∗(t))P∗(t)N∗(t)+λ2(t)ϕ−λ3(t)β(1−u2∗(t))P∗(t)N∗(t)−λ3(t)ϕλ3′(t)=−∂H∂P=−AP+λ1(t)α(1−u1∗(t))S∗(t)N∗(t)−λ2(t)α(1−u1∗(t))S∗(t)N∗(t)+λ2(t)β(1−u2∗(t))L∗(t)N∗(t)−λ2(t)ψ−λ3(t)β(1−u2∗(t))L∗(t)N∗(t)+λ3(t)ψ+λ3(t)(γH∗(t)N∗(t)+θ+ρu3∗(t))−λ4(t)(γH∗(t)N∗(t)+θ+ρu3∗(t))λ4′(t)=−∂H∂H=λ3(t)(γP∗(t)N∗(t)−τ)−λ4(t)(γP∗(t)N∗(t)−τ)
To obtain the optimality conditions (8), we also differentiated the Hamiltonian H, with respect to u1(t),u2(t),u3(t) and set it equal to zero.
0=∂H∂u1=B1u1∗(t)+λ1(t)αS∗(t)L∗(t)+P∗(t)N(t)−λ2(t)αS∗(t)L∗(t)+P∗(t)N(t)0=∂H∂u2=B2u2∗(t)+λ2(t)βL∗(t)P∗(t)N∗(t)−λ3(t)βL∗(t)P∗(t)N∗(t)0=∂H∂u3=B3u3∗(t)−λ3(t)ρP∗(t)+λ4(t)ρP∗(t)
Solving for optimal controls, we obtained:
u1∗(t)=αS∗(t)(L∗(t)+P∗(t))(λ2(t)−λ1(t))B1N(t)u2∗(t)=βL∗(t)P∗(t)(λ3(t)−λ2(t))B2N(t)u3∗(t)=ρP∗(t)(λ3(t)−λ4(t))B3
To determine an explicit expression for the optimal controls for 0≤ui∗(t)≤1,(i=1,2,3), a standard optimality technique was utilized. We considered the following three cases.
On the set {t:0<u1∗(t)<1}, we had ∂H/∂u1=0. Hence, the optimal control was:
u1∗(t)=αS∗(t)(L∗(t)+P∗(t))(λ2(t)−λ1(t))B1N(t)
On the set {t:u1∗(t)=0}, we had ∂H/∂u1≥0. This implies that:
λ1(t)αS∗(t)L∗(t)+P∗(t)N(t)−λ2(t)αS∗(t)L∗(t)+P∗(t)N(t)≥0
when we had that
αS∗(t)(L∗(t)+P∗(t))(λ2(t)−λ1(t))/B1N(t)≤0=u1∗(t)
On the set {t:u1∗(t)=1}, we had ∂H/∂u1≤0. This implied that
λ1(t)αS∗(t)L∗(t)+P∗(t)N(t)−λ2(t)αS∗(t)L∗(t)+P∗(t)N(t)≤−B1
when we had that
αS∗(t)(L∗(t)+P∗(t))(λ2(t)−λ1(t))/B1N(t)≥1=u1∗(t)
Combining these three cases above, we found a characterization of u1∗
u1∗(t)=min{1,max{0,1B1(αS∗(L∗+P∗)(λ2−λ1)N)}}
Using similar arguments, we also obtained the second and third optimal control function
u2∗(t)=min{1,max{0,1B2(βL∗P∗(λ3−λ2)N)}}u3∗(t)=min{1,max{0,ρP∗(λ3−λ4)B3}}
This work was supported by Kyungpook National University Research Fund, 2012.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Article information Continued
© 2014 Published by Elsevier B.V. on behalf of Korea Centers for Disease Control and Prevention.
