Loading [MathJax]/jax/output/SVG/jax.js
Research article

Beyond Americas War on Drugs: Developing Public Policy to Navigate the Prevailing Pharmacological Revolution

  • This paper places America's “war on drugs” in perspective in order to develop a new metaphor for control of drug misuse. A brief and focused history of America's experience with substance use and substance use policy over the past several hundred years provides background and a framework to compare the current Pharmacological Revolution with America's Nineteenth Century Industrial Revolution. The paper concludes with cautions about growing challenges and provides suggestions for navigating this revolution and reducing its negative impact on individuals and society.

    Citation: Andrew Golub, Alex S. Bennett, Luther Elliott. Beyond Americas War on Drugs: Developing Public Policy to Navigate the Prevailing Pharmacological Revolution[J]. AIMS Public Health, 2015, 2(1): 142-160. doi: 10.3934/publichealth.2015.1.142

    Related Papers:

    [1] Liumei Wu, Baojun Song, Wen Du, Jie Lou . Mathematical modelling and control of echinococcus in Qinghai province, China. Mathematical Biosciences and Engineering, 2013, 10(2): 425-444. doi: 10.3934/mbe.2013.10.425
    [2] Xi-Chao Duan, Xue-Zhi Li, Maia Martcheva . Dynamics of an age-structured heroin transmission model with vaccination and treatment. Mathematical Biosciences and Engineering, 2019, 16(1): 397-420. doi: 10.3934/mbe.2019019
    [3] Hongbin Guo, Michael Yi Li . Global dynamics of a staged progression model for infectious diseases. Mathematical Biosciences and Engineering, 2006, 3(3): 513-525. doi: 10.3934/mbe.2006.3.513
    [4] Kai Wang, Zhidong Teng, Xueliang Zhang . Dynamical behaviors of an Echinococcosis epidemic model with distributed delays. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1425-1445. doi: 10.3934/mbe.2017074
    [5] Shuang-Hong Ma, Hai-Feng Huo . Global dynamics for a multi-group alcoholism model with public health education and alcoholism age. Mathematical Biosciences and Engineering, 2019, 16(3): 1683-1708. doi: 10.3934/mbe.2019080
    [6] Rundong Zhao, Qiming Liu, Huazong Zhang . Dynamical behaviors of a vector-borne diseases model with two time delays on bipartite networks. Mathematical Biosciences and Engineering, 2021, 18(4): 3073-3091. doi: 10.3934/mbe.2021154
    [7] Ruixia Zhang, Shuping Li . Analysis of a two-patch SIS model with saturating contact rate and one- directing population dispersal. Mathematical Biosciences and Engineering, 2022, 19(11): 11217-11231. doi: 10.3934/mbe.2022523
    [8] Tailei Zhang, Hui Li, Na Xie, Wenhui Fu, Kai Wang, Xiongjie Ding . Mathematical analysis and simulation of a Hepatitis B model with time delay: A case study for Xinjiang, China. Mathematical Biosciences and Engineering, 2020, 17(2): 1757-1775. doi: 10.3934/mbe.2020092
    [9] Mingtao Li, Xin Pei, Juan Zhang, Li Li . Asymptotic analysis of endemic equilibrium to a brucellosis model. Mathematical Biosciences and Engineering, 2019, 16(5): 5836-5850. doi: 10.3934/mbe.2019291
    [10] Shengqiang Liu, Lin Wang . Global stability of an HIV-1 model with distributed intracellular delays and a combination therapy. Mathematical Biosciences and Engineering, 2010, 7(3): 675-685. doi: 10.3934/mbe.2010.7.675
  • This paper places America's “war on drugs” in perspective in order to develop a new metaphor for control of drug misuse. A brief and focused history of America's experience with substance use and substance use policy over the past several hundred years provides background and a framework to compare the current Pharmacological Revolution with America's Nineteenth Century Industrial Revolution. The paper concludes with cautions about growing challenges and provides suggestions for navigating this revolution and reducing its negative impact on individuals and society.


    An estimated 15 million people are infected with clonorchiasis, predominantly in east Asia (ie, China, South Korea, northern Vietnam, and parts of Russia) [1,2]. Clonorchiasis infection is primarily related to liver and biliary disorders, especially cholangiocarcinoma.

    The first report on clonorchiasis in the medical literature was published in the Lancet on August 21, 1875 [3]. In 1910, Kobayashi first proved that freshwater fish act as the second intermediate host, and Muto discovered in 1918 that freshwater snails serve as the first intermediate host [4]. The life cycle of clonorchiasis is characterised by an alternation of sexual and asexual reproduction in different hosts [5,6,7]. Eggs laid by hermaphroditic adult worms, which can reach the intestine with bile fluids and are taken away with the faeces [8]. After freshwater snails ingest the eggs, miracidia hatch in the intestine and penetrate into the intestinal wall. By asexual reproduction, sporocysts, rediae, and then cercariae are produced. Cercariae escape from the snails, adhere to freshwater fish, and then develop into mature metacercariae in the subcutaneous tissues or muscles of the fish. When people eat raw or insufficiently cooked infected fish, metacercariae will enter the human body through gastric juice digestion as the final host (Figure 1).

    Figure 1.  Life cycle of Clonorchis sinensis.

    Mathematical models have played significant roles for understanding and controlling of infectious diseases since the publication [9]. This work was motivated by the recent work for mathematical modeling of clonorchiasis [10], which had considered the transmission dynamic of clonorchiasis among human beings, snails and fish. However, snails are infected only after ingested eggs in faeces produced by parasitic worms in the bile ducts of infectious human being, fish are infected only after cercariae escaped from infectious snails and adhere to freshwater fish. In order to investigate the efficient strategies to control clonorchiasis, we must set up a mathematical model which responds the basic life cycle of clonorchiasis, and consider eggs both in an infected individual and in freshwater, cercariae both in infectious snails and in freshwater. Praziquantel and albendazole are recommended drugs for treatment of clonorchiasis. The efficacy of those drugs depends on the treatment schedule and infection intensity. We will consider the treatment (or recovery) rate in our mathematical model of clonorchiasis such that we can compare the efficacy of different control strategy.

    This paper is organized as follows. In section 2, we will set up a mathematical model to describe the transmission dynamic of clonorchiasis, calculate the basic reproductive number, and analyze the relationship between the the basic reproductive number and the global property of the model. In section 3, we will fit the data of human clonorchiasis infection ratios of Guangzhou City of Guangdong Province in China from 2006-2012 by using the Genetic Algorithm and calculate the 95% confidence intervals of infection ratios of human. In section 4, we will give the prevention and control strategy of clonorchiasis.

    Both clonorchiasis and schistosomiasis have an important intermediate medium snail as vector. There are many papers to study the infectious mechanisms and control strategies of schistosomiasis [11,12,13,14,15,16,17]. They use the qualitative and the quantitative analysis to study the human$ - $snail transmission of schistosomiasis. Especially, Chen et. al [12] investigated the human$ - $cattle$ - $snail transmission of schistosomiasis.

    In order to investigate the control strategies of clonorchiasis, we set up the following $ S_hI_hR_h $ $ ES_sI_sCS_fI_f $ model by the the life cycle of clonorchiasis:

    $ {Sh(t)=λhβhSh(t)If(t)μhSh(t),Ih(t)=βhSh(t)If(t)γIh(t)μhIh(t),Rh(t)=γIh(t)μhRh(t),E(t)=θepeIh(t)δeβsSs(t)E(t)μeE(t),Ss(t)=λsβsSs(t)E(t)μsSs(t),Is(t)=βsSs(t)E(t)μsIs(t),C(t)=θcpcIs(t)δcβfSf(t)C(t)μcC(t),Sf(t)=λfβfSf(t)C(t)qSf(t)μfSf(t),If(t)=βfSf(t)C(t)qIf(t)μfIf(t), $ (2.1)

    where $ S_h(t) $, $ I_h(t) $ and $ R_h(t) $ denote the number of the susceptible, the infectious and the recovered for human beings, $ E(t) $ and $ C(t) $ denote the number of eggs and cercariae, $ S_s(t) $ and $ I_s(t) $ denote the number of the susceptible and the infectious for snails, $ S_f(t) $ and $ I_f(t) $ denote the number of the susceptible and the infectious for freshwater fish.

    Parameters in system (2.1) are summarized in Table 1 and all parameters are assumed to be nonnegative. $ \delta_e $ and $ \delta_c $ denote the consumption rate of eggs per snail and consumption rate of cercariae per fish, respectively. $ p_e $ denotes the egg proliferation rate by hermaphroditic adult worms in an infected individual, those eggs are attached with the faeces. Only a part of the faeces flows into fresh water. Parameter $ \theta_e $ denotes the ratio of eggs in fresh water. Similarly, $ p_c $ denotes the cercaria proliferation rate per infected snail in fresh water. $ \theta_c $ denotes the proportion of cercariae survival in fresh water.

    Table 1.  Description of model parameters and their values (unit:$ year^{-1} $).
    Parameter Value Interpretation Source
    $ \lambda_h $ $ 2.2404\times10^5 $ The recruitment rate of human population fitting
    $ \lambda_s $ $ 4.2741\times10^9 $ The recruitment rate of snail population fitting
    $ \lambda_f $ $ 3.5799\times10^9 $ The recruitment rate of fish population fitting
    $ \beta_h $ $ 3.6\times10^{-12} $ Transmission rate from infected fish to human fitting
    $ \beta_s $ $ 4.25 \times10^{-11} $ Transmission rate from egg to snail [13]
    $ \beta_f $ $ 9.784 \times10^{-9} $ Transmission rate from cercaria to fish fitting
    $ \mu_h $ 0.0128 Death rate of human being [14]
    $ \mu_s $ 1 Death rate of snail [12,13]
    $ \mu_f $ 0.5 Death rate of fish fitting
    $ \mu_e $ 0.3519 Death rate of egg fitting
    $ \mu_c $ 0.4616 Death rate of cercaria fitting
    $ \gamma $ 0.79 Recovery rate of human being [14]
    $ p_c $ $ 2.22\times10^3 $ The cercariae proliferation rate per infected snail [13]
    $ p_e $ $ 3.33\times10^6 $ The egg proliferation rate per infected individual [15]
    $ \theta_e $ 0.3499 The ratio of eggs in fresh water fitting
    $ \theta_c $ 0.9 The proportion of cercariae survive in fresh water fitting
    $ q $ 4.5 The harvest rate of fish fitting

     | Show Table
    DownLoad: CSV

    The number of eggs consumed by snails compared to the number of eggs in the environment is too small [17]. Thus, the term $ \delta_e\beta_sS_s(t)E(t) $ can be removed from the fourth equations of system (2.1). Similarly, the term $ \delta_c\beta_fS_f(t)C(t) $ can be removed from the seventh equations of system (2.1). System (2.1) can be written in the following form:

    $ {Sh(t)=λhβhSh(t)If(t)μhSh(t),Ih(t)=βhSh(t)If(t)γIh(t)μhIh(t),Rh(t)=γIh(t)μhRh(t),E(t)=θepeIh(t)μeE(t),Ss(t)=λsβsSs(t)E(t)μsSs(t),Is(t)=βsSs(t)E(t)μsIs(t),C(t)=θcpcIs(t)μcC(t),Sf(t)=λfβfSf(t)C(t)qSf(t)μfSf(t),If(t)=βfSf(t)C(t)qIf(t)μfIf(t). $ (2.2)

    In this paper, we assume that all solutions of system (2.2) satisfy the following initial conditions:

    $ Sh(0)0,Ih(0)0,Rh(0)0,E(0)0,Ss(0)0,Is(0)0,C(0)0,Sf(0)0,If(0)0. $

    Let

    $ Ω={(Sh(t),Ih(t),Rh(t),E(t),Ss(t),Is(t),C(t),Sf(t),If(t))R9+    Sh(t)+Rh(t)+Ih(t)λhμh, Ss(t)+Is(t)λsμs,Sf(t)+If(t)λfμf+q,     E(t)θepeλhμeμh, C(t)θcpcλsμcμs}. $

    Lemma 2.1. The solutions of system (2.2) are positive if all initial values are positive and $ \Omega $ is a positively invariant set of system (2.2).

    Proof. At first, we will prove the positivity of the solution of system (2.2). By system (2.2), we have

    $ {Sh(t)βhSh(t)If(t)μhSh(t),Ih(t)γIh(t)μhIh(t),Rh(t)μhRh(t),E(t)μeE(t),Ss(t)βsSs(t)E(t)μsSs(t),Is(t)μsIs(t),C(t)μcC(t),Sf(t)βfSf(t)C(t)qSf(t)μfSf(t),If(t)qIf(t)μfIf(t), $

    for $ t\geq 0. $ Considering the following auxiliary system:

    $ {˜Sh(t)=βh˜Sh(t)˜If(t)μh˜Sh(t),˜Ih(t)=γ˜Ih(t)μh˜Ih(t),˜Rh(t)=μh˜Ih(t),˜E(t)=μe˜E(t),˜Ss(t)=βs˜Ss(t)˜E(t)μs˜Ss(t),˜Is(t)=μs˜Is(t),˜C(t)=μc˜C(t),˜Sf(t)=βf˜Sf(t)˜C(t)q˜Sf(t)μf˜Sf(t),˜If(t)=q˜If(t)μf˜If(t). $

    We can obtain

    $ ˜Sh(t)=˜Sh(0)et0(βh˜If(s)μh)ds,˜Ih(t)=˜Ih(0)et0(γμh)ds,˜Rh(t)=˜Rh(0)et0μhds,˜E(t)=˜E(0)et0μeds,˜Ss(t)=˜Ss(0)et0(βs˜E(s)μs)ds,˜Is(t)=˜Is(0)et0μsds,˜C(t)=˜C(0)et0μcds,˜Sf(t)=˜Sf(0)et0(βf˜C(s)qμf)ds,˜If(t)=˜If(0)et0(qμf)ds. $

    By the comparison principle, we have obtained that all solutions of system (2.2) are positive when their initial values are positive. Therefore, $ R^9_+ = \{(S_h(t), I_h(t), R_h(t), E(t), S_s(t), I_s(t), $ $ C(t), S_f(t), I_f(t))\in R^9\mid S_h(t)\geq 0, I_h(t)\geq0, R_h(t)\geq 0, E(t)\geq0, S_s(t)\geq 0, I_s(t)\geq0, C(t)\geq0, S_f(t)\geq 0, I_f(t)\geq0\} $ is positive invariant with respect to system (2.2).

    And then, we will prove that $ \Omega $ is a positively invariant set of system (2.2). Let $ X(t) = (S_h(t), I_h(t), R_h(t), E(t), S_s(t), I_s(t), C(t), S_f(t), I_f(t)) $ be any positive solution of system (2.2) with initial condition $ X(0)\in \Omega $. Let $ N_h(t) = S_h(t)+R_h(t)+I_h(t) $, $ N_s(t) = S_s(t)+I_s(t) $, and $ N_f(t) = S_f(t)+I_f(t) $ denote the total population of human beings, the total population of snails, and the total population of fish, respectively. The total population of human beings satisfies the differential equation:

    $ Nh(t)=λhμhNh(t). $

    This implies that $ N_h(t)\rightarrow \frac{\lambda_h}{\mu_h} $ as $ t\rightarrow +\infty, $ or $ \lim_{t\rightarrow +\infty}(S_h(t)+I_h(t)+R_h(t)) = \frac{\lambda_h}{\mu_h} $. Similarly, we have

    $ Ns(t)=λsμsNs(t), $
    $ Nf(t)=λf(q+μf)Nf(t), $

    $ N_s(t)\rightarrow \frac{\lambda_s}{\mu_s} $ and $ N_f(t)\rightarrow \frac{\lambda_f}{\mu_f+q} $ as $ t\rightarrow +\infty, $ or $ \lim_{t\rightarrow +\infty}(S_s(t)+I_s(t)) = \frac{\lambda_s}{\mu_s} $ and $ \lim_{t\rightarrow +\infty} (S_f(t)+I_f(t)) = \frac{\lambda_f}{\mu_f+q}. $

    By the 4th and 7th equation of system (2.2), we have

    $ E(t)θepeλhμhμeE(t),C(t)θcpcλsμsμcC(t), $

    $ \lim_{t\rightarrow +\infty} E(t) \leq \frac{\theta_ep_e\lambda_h}{\mu_e\mu_h} $, $ \lim_{t\rightarrow +\infty} C(t) \leq \frac{\theta_cp_c\lambda_s}{\mu_c\mu_s}. $ Thus, $ \Omega $ being a positively invariant set of system (2.2) has been proven.

    The basic reproduction number of system (2.2) can be derived by the next generation matrix as in Diekmann et al. [18,19], Van den Driessche and Watmough [20]. Let $ x = (I_h, I_s, I_f, E, C, S_h, R_h, $ $ S_s, S_f)^T $. System (2.2) can be written as the following form:

    $ x=F(x)V(x), $

    where $ \mathscr{V}(x) = \mathscr{V}^-(x)-\mathscr{V}^+(x), $ which denotes the transfer rate of individuals into or out of each population set and $ \mathscr{F}(x) $ denotes the rate of occurrence of new infections and

    $ \mathscr{F}(x) = \left(βhSh(t)If(t)βsSs(t)E(t)βfSf(t)C(t)000000\right), \quad \mathscr{V}^+(x) = \left(000θepeIh(t)θcpcIs(t)λhγIh(t)λsλf\right), $
    $ \mathscr{V}^-(x) = \left((μh+γ)Ih(t)μsIs(t)(q+μf)If(t)μeE(t)μcC(t)βhSh(t)If(t)+μhSh(t)μhRh(t)βsSs(t)E(t)+μsSs(t)βfSf(t)C(t)+qSf(t)+μfSf(t)\right). $

    According to the method in Diekmann et al. [18], the derivatives $ D \mathscr{F}(x) $ and $ D \mathscr{V}(x) $ are partitioned as

    $ D \mathscr{F}(x) = \left(F000\right), \ \ D \mathscr{V}(x) = \left(V0J1J2\right), $

    where $ F, V, J_1 $ and $ J_2 $ are given by

    $ F=(00βhλhμh00000βsλsμs00000βfλfq+μf0000000000),  V=(μh+γ00000μs00000q+μf00θepe00μe00θcpc00μc), $ (2.3)
    $ J_1 = \left(00βhλhμh00γ0000000βsλsμs00000βfλfq+μf\right),\ \ J_2 = \left(μh0000μh0000μs0000q+μf\right). $

    We can obtain the inverse matrix of $ V $ as

    $ V1=(1μh+γ000001μs000001q+μf00θePeμe(μh+γ)001μe00θcPcμcμs001μc). $

    The basic reproduction number is defined as the spectral radius of the nonnegative matrix $ FV^{-1} $ [20]. Therefore, we have obtained the formula of the basic reproduction number:

    $ R_0 = \rho(FV^{-1}) = \sqrt[3]{\frac{\beta_h \beta_s\beta_f \lambda_h \lambda_s \lambda_f \theta_c p_c \theta_e p_e}{\mu^{2}_s (q+\mu_f)^{2} \mu_c \mu_e\mu_h(\mu_h+\gamma)}}. $

    At first, we will prove that the system (2.2) admits a unique endemic equilibrium $ x^* = (S_h^{*}, I_h^{*}, R_h^{*}, $ $ E^{*}, S_s^{*}, I_s^{*}, C^{*}, S_f^{*}, I_f^{*}) $. And then, we will prove that $ x^* $ is globally asymptotically stable in the region $ R^9_{+}. $ Here $ S_h^{*}, I_h^{*}, R_h^{*}, E^{*}, S_s^{*}, I_s^{*}, C^{*}, S_f^{*}, I_f^{*} $ satisfy the following equilibrium equations:

    $ {λh=βhShIf+μhSh,βhShIf=γIh+μhIh,μhRh=γIh,μeE=θepeIh,λs=βsSsE+μsSs,μsIs=βsSsE,μcC=θcpcIs,λf=βfSfC+(q+μf)Sf,(q+μf)If=βfSfC. $ (2.4)

    By solving equations (2.4), we have

    $ I_f^{*} = \frac{\mu_e\mu_c\mu^{2}_s(q+\mu_f)^2\mu_h(\gamma+\mu_h)(R_0^{3}-1)}{\theta_ep_e\lambda_h\beta_h\beta_s(\theta_cp_c\beta_f\lambda_s(\mu_f+q)+\mu_c\mu_s(\mu_f+q)^2)+\mu_e\mu_c\mu^{2}_s(\mu_f+q)^2\beta_h(\gamma+\mu_h)}, $
    $ S_h^{*} = \frac{\lambda_h}{\beta_h I_f^{*}+\mu_h},\ I_h^{*} = \frac{\lambda_h\beta_hI_f^{*}}{(\gamma+\mu_h)(\mu_h+\beta_hI_f^{*})},\ R_h^{*} = \frac{\gamma I_h^{*}}{\mu_h}, \ E^{*} = \frac{\theta_ep_eI_h^{*}}{\mu_e}, $
    $ S_s^{*} = \frac{\lambda_s}{\beta_s E^{*}+\mu_s},\ I_s^{*} = \frac{\beta_sS_s^{*}E^{*}}{\mu_s},\ C^{*} = \frac{\theta_cp_cI_s^{*}}{\mu_c},\ S_f^{*} = \frac{\lambda_f}{\beta_f C^{*}+q+\mu_f}. $

    It is obvious to follow that $ I_f^{*} > 0 $ if and only if $ R_0 > 1 $. By the equations (2.4) and $ R_0 > 1 $, there is only one $ x^* \in \Omega $ for system (2.2). We have the following result.

    Theorem 2.2. If $ R_0 > 1, $ then the endemic equilibrium $ x^* $ is globally asymptotically stable in the interior of $ R^9_{+}. $

    Proof. The existence of $ x^* $ of system (2.2) has been shown in above. We will only prove that $ x^* $ of system (2.2) is globally asymptotically stable in the interior of $ R^9_{+} $ when $ R_0 > 1 $. Let $ g(a) = 1-a+\ln{a} $. It is easy to verify that for $ \forall a > 0 $,

    $ g(a) = 1-a+\ln{a}\leq 0, $

    and the equality holds if and only if $ a = 1 $.

    Let $ L_1 = \int_{S_h^{*}}^{S_h}\frac{z-S_h^{*}}{z}dz+I_h-I_h^{*}-\ln{\frac{I_h}{I_h^{*}}}, L_2 = \int_{S_s^{*}}^{S_s}\frac{z-S_s^{*}}{z}dz+I_s-I_s^{*}-\ln{\frac{I_s}{I_s^{*}}}, L_3 = \int_{S_f^{*}}^{S_f}\frac{z-S_f^{*}}{z}dz+I_f-I_f^{*}-\ln{\frac{I_f}{I_f^{*}}}, L_4 = E-E^{*}-E^{*}\ln{\frac{E}{E^{*}}}, L_5 = C-C^{*}-C^{*}\ln{\frac{C}{C^{*}}}. $ Differentiating $ L_1, L_2, L_3, L_4 $ and $ L_5 $ with respect to t along solution curves of system (2.2), we have

    $ L1=(1ShSh)Sh+(1IhIh)Ih=(1ShSh)[μhSh(1ShSh)+βhShIf(1ShIfShIf)]+βhShIf(1IhIh)(ShIfShIfIhIh)βhShIf(1ShIfShIf)(1ShSh)+βhShIf(1IhIh)(ShIfShIfIhIh)=βhShIf(2IhIhShShShIfIhShIfIh+IfIf)βhShIf(IfIflnIfIf+lnIhIhIhIh), $

    where $ 1-\frac{S_h^{*}}{S_h}\leq -\ln{\frac{S_h^{*}}{S_h}} $, $ 1-\frac{S_hI_fI_h^{*}}{S_h^{*}I_f^{*}I_h} \leq -\ln{\frac{S_hI_fI_h^{*}}{S_h^{*}I_f^{*}I_h}} $.

    $ L2=(1SsSs)Ss+(1IsIs)Is=(1SsSs)(βsSsE(1SsESsE)+usSs(1SsSs))+βsSsE(1IsIs)(SsSsEIsIs)βsSsE(2SsSsIsIs+EESsEIsSsEIs)βsSsE(EElnEE+lnIsIsIsIs). $
    $ L3=(1SfSf)Sf+(1IfIf)If=(1SfSf)(βfSfC(1SfCSfC)+(q+μf)Sf(1SfSf))+βfSfC(1IfIf)(SfCSfCIfIf)βfSfC(CClnCC+lnIfIfIfIf). $
    $ L4=(1EE)E=θepeIh(1EE+IhIhIhIhEE)θepeIh(IhIhlnIhIh+lnEEEE). $
    $ L5=(1CC)C=θcpcIs(1CC+IsIsIsIsCC)θcpcIs(IsIslnIsIs+lnCCCC). $

    Hence, we can define the following Lyapunov function:

    $ L=L1βhShIf+L2βsSsE+L3βfSfC+L4θepeIh+L5θcpcIs. $

    It follows that

    $ L=1βhShIfL1+1βsSsEL2+1βfSfCL3+1θepeIhL4+1θcpcIsL5(IfIflnIfIf+lnIhIhIhIh)+(EElnEE+lnIsIsIsIs)+(CClnCC+lnIfIfIfIf)  +(IhIhlnIhIh+lnEEEE)+(IsIslnIsIs+lnCCCC)=0. $

    Moreover, the equality $ L^{'} = 0 $ holds if and only if $ I_h = I_h^{*}, I_s = I_s^{*}, I_f = I_f^{*}, E = E^{*}, C = C^{*}, S_h = S_h^{*}, R_h = R_h^{*}, S_s = S_s^{*} $ and $ S_f = S_f^{*}. $ Therefore, $ x^* $ of system (2.2) is globally asymptotically stable in the region $ R^9_{+} $ if $ R_0 > 1. $

    If $ R_0 < 1 $, the endemic equilibrium $ x^* $ of system (2.2) does not exist. There is only the disease-free equilibrium $ x^0 $, where $ x^{0} = (S_h^{0}, 0, 0, 0, S_s^{0}, 0, 0, S_f^{0}, 0) $. We have the following result.

    Theorem 2.3. If $ R_0 < 1, $ then the disease-free equilibrium $ x^0 $ of system (2.2) is globally asymptotically stable in the region $ R^9_{+} $.

    Proof. In order to prove the Theorem 2.3, at first we will prove that the origin of a subsystem of system (2.2) is globally asymptotically stable in $ R^5_+ $, where $ R^5_+ = \{(I_h(t), I_s(t), I_f(t), E(t), C(t))\in R^5 \mid I_h(t)\geq 0, I_s(t)\geq 0, I_f(t)\geq 0, E(t)\geq 0, C(t)\geq 0\} $. And then, we will prove that the $ x^0 $ of system (2.2) is globally asymptotically stable in the region $ R^9_{+} $.

    Step 1 Let $ x = (I_h, I_s, I_f, G, C)^T $. By system (2.2), we have the following subsystem.

    $ {Ih(t)=βhSh(t)If(t)γIh(t)μhIh(t),Is(t)=βsSs(t)E(t)μsIs(t),If(t)=βfSf(t)C(t)qIf(t)μfIf(t),E(t)=θepeIh(t)μeE(t),C(t)=θcpcIs(t)μcC(t). $ (2.5)

    The origin $ \overline{x_{0}} = (0, 0, 0, 0, 0) $ is an equilibrium of the subsystem (2.5). All solutions of system (2.5) remain nonnegative and

    $ Ω0={(Ih(t),Is(t),If(t),E(t),C(t))R5+0Ih(t)λhμh,0Is(t)λsμs,     0If(t)λfp+μf, 0E(t)θepeλhμeμh, 0C(t)θcpcλsμcμs} $

    is positively invariant for system (2.5). If the origin $ \overline{x_{0}} $ is globally stable for $ \forall x \in \Omega_0 $, then, $ \overline{x_{0}} $ is globally stable for $ \forall x \in R^5_+ $.

    By system (2.5), for $ \forall x \in \Omega_0 $ we can obtain the following inequalities

    $ x(FV)x, $ (2.6)

    where $ F $ and $ V $ are expressed as (2.3). We construct a new linear system

    $ x=(FV)x, $ (2.7)

    where $ x \in R^5_+ $. By the Perron-Frobenius Theorem [21], there is a nonnegative left eigenvector $ \nu $ of the nonnegative matrix $ V^{-1}F $ corresponding to the eigenvalue $ \rho(V^{-1}F) = R_0, $ that is $ \nu^TV^{-1}F = R_0\nu^T. $ Define a Lyapunov function

    $ V_0 = \nu^TV^{-1}x. $

    Then the derivative of $ V_0 $ along with system (2.7) is

    $ V0=νTV1x=νTV1(FV)x=νTV1FxνTx(R01)νTx. $

    If $ R_0 < 1, $ then $ V_0^{'}\leq0. $ Let

    $ Ψ={(Ih(t),Is(t),If(t),E(t),C(t))Ω0V0=0}. $

    If $ R_0 < 1, $ $ V_0^{'} = 0 $ implies that $ \nu^Tx = 0. $ Thus, $ I_h = 0, I_s = 0, I_f = 0, G = 0, C = 0. $ Therefore, the largest invariant set of $ \Psi $ is the point $ \{\overline{x_{0}}\}. $ By LaSalle Invariance Principle [22], $ \overline{x_{0}} $ of the new linear system (2.7) is globally asymptotically stable in $ R^5_+ $. By inequalities (2.6) and the comparison theorem of ordinary differential equation, each component of the system (2.5) must approach to zero when its initial value in $ \Omega_0 $. We have finished the proof that the origin of the system (2.5) is globally asymptotically stable in the region $ \Omega_0 $ when $ R_0 < 1. $

    Step 2 Now, we will prove that $ x^0 $ of system (2.2) is globally asymptotically stable in the region $ R^9_{+} $. By Step 1, for $ \forall \varepsilon > 0, $ there exists a $ T > 0 $ such that for $ t > T, $ we have

    $ 0Ih(t)<ε, 0Is(t)<ε, 0If(t)<ε, 0E(t)<ε,0C(t)<ε. $

    When $ t > T, $ by the first equation of system (2.2), we have

    $ λhβhSh(t)εμhSh(t)Sh(t)λhμhSh(t). $

    We construct the following equations:

    $ S1(t)=λhμhS1(t),  limt+S1(t)=λhμh.S2(t)=λhβhS2(t)εμhS2(t),  limt+S2(t)=λhμh+βhε. $

    By the comparison theorem of ordinary differential equation, we have

    $ λhμh+βhεlimt+Sh(t)λhμh. $

    Let $ \varepsilon > 0 $ be sufficiently small. We have

    $ limt+Sh(t)=λhμh. $

    Similarly, from the third, fifth and eighth equations of system (2.2), we can prove

    $ limt+Rh(t)=0,limt+Ss(t)=λsμs,limt+Sf(t)=λfμf+p. $

    It follows that the disease-free equilibrium $ x^{0} = (\frac{\lambda_h}{\mu_h}, 0, 0, 0, \frac{\lambda_s}{\mu_s}, 0, 0, \frac{\lambda_f}{p+\mu_f}, 0) $ of system (2.2) is globally asymptotically stable.

    The data are infection ratios of human in Guangzhou City of Guangdong Providence in China from 2006-2012 [23] (see Figure 2). The infection ratio of human in model (2.2) is determined by the formula $ \frac{I_h(t)}{N_h(t)} $. The initial values of system (2.2) are given by $ S_h(0) = 9.4652\times 10^6, I_h(0) = 1.968\times 10^4, R_h(0) = 2.33 \times 10^6, E(0) = 120, S_s(0) = 6.95\times 10^3, I_s(0) = 130, C(0) = 100, S_f(0) = 2.2628\times 10^5, I_f(0) = 1.87 \times 10^4 $. Some parameter values are cited from references, the other parameter values are estimated by calculating the minimum sum of square (MSS):

    $ MSS=2012i=2006(dataiIh(i)Nh(i))2. $
    Figure 2.  The solid red lines represent the data reported by the epidemiological research and investigation [23] while the solid blue lines are simulated by using the model (2.2), the vertical segments are shown the 95% confidence intervals of infection ratios of human.

    The minimum sum of square (MSS) can be regarded as an optimization problem which is solved using Genetic Algorithm [24,25]. The efficiency of Genetic Algorithm depends on the appropriate choice of the starting population along with other associated parameters. The initial population size plays an important role in the quality and efficiency of the algorithm, and small population size results in local convergence or requirement of large generations. To avoid this, the population size 1000 and the gene size 50 are chosen. A crossover probability of 0.55 and mutation probability of 0.15 are chosen to maitain diversity in the population. By the parameter values in Table 1, we can estimate that the basic reproduction number of human clonorchiasis is $ R_0 = 1.41 $. This result is shown that clonorchiasis will be epidemic in Guangzhou City. Figure 3 presents the prediction of human clonorchiasis infection ratios and of the 95% confidence intervals.

    Figure 3.  The prediction of human clonorchiasis infection ratios in Guangzhou City.

    By simulation of the data of human clonorchiasis infection ratios of Guangzhou City of Guangdong Province from 2006–2012 (see Figure 2), we have obtained that the basic reproduction number is $ R_0 = 1.41 $ which indicates that clonorchiasis will be epidemic in Guangzhou City. The wide distribution of intermediate hosts and reservoir hosts, human eating habits, the lack of the food safety supervision of freshwater fish, and relatively undeveloped techniques for detection and treatment are contributing to the prevalence of clonorchiasis [23]. Therefore, it is very difficult to eradicate clonorchiasis even if some chemotherapy and control programmes have been implemented over several years in a few endemic areas [26,27].

    In Table 2, we have proposed seven types of control strategies dealing with three parameters $ \beta_h $, $ \gamma $, $ \theta_e $ which mean decreasing the transmission rate from infected fish to human, increasing the recovery rate of human being, and decreasing the ratio of the eggs in fresh water. If we decrease by half of the transmission rate from infected fish to human, increase to one for the recovery rate of human being, and decrease by half of the ratio of the eggs in fresh water, respectively, it is still impossible to eradicate the clonorchiasis, which means that single measure is difficult to eradicate the clonorchiasis. However, increasing the recovery rate has small effect on the epidemic of the clonorchiasis (see Figure 4(a)), decreasing the transmission rate from infected fish to human or decreasing the ratio of the eggs in fresh water has obvious effect on the epidemic of the clonorchiasis (see Figure 4(b), (c)).

    Table 2.  Comparison of different control strategies (FTH:Transmission rate from infected fish to human, REFW: Ratio of the eggs in fresh water).
    Plan FTH Recovery rate REFW $ R_0 $ Extinction Permanence
    1 $ \beta_h\div2 $ $ \gamma \times 1 $ $ \theta_e \times 1 $ 1.1269 No Yes
    2 $ \beta_h\times 1 $ $ \gamma = 1 $ $ \theta_e \times 1 $ 1.3140 No Yes
    3 $ \beta_h\times 1 $ $ \gamma \times 1 $ $ \theta_e \div 2 $ 1.1162 No Yes
    4 $ \beta_h\div2 $ $ \gamma = 1 $ $ \theta_e \times 1 $ 1.0429 No Yes
    5 $ \beta_h\div2 $ $ \gamma \times 1 $ $ \theta_e \div 2 $ 0.8859 Yes No
    6 $ \beta_h\times 1 $ $ \gamma =1 $ $ \theta_e\div 2 $ 1.0330 No Yes
    7 $ \beta_h\div2 $ $ \gamma =1 $ $ \theta_e\div 2 $ 0.7050 Yes No

     | Show Table
    DownLoad: CSV
    Figure 4.  The relationship between $ R_0 $ and $ \gamma, \theta_e, \beta_h $.

    The prevention and control strategy of clonorchiasis must be a combination of two or more measures, including health education, health promotion, chemotherapy and environmental reconstruction [28,29]. Health education includes the broadcast of educational programmes on television, broadcast and VCDs, billboard/propaganda paintings, the distribution of health guide booklets, and the transmission of disease-related knowledge to residents and school children. The contents of health education are avoid to eat raw or undercooked freshwater fish. Health promotion includes the investigation of the infection rates and distribution of freshwater fish and snails in endemic areas, under surveillance of the infected ponds, rapid, convenient and accurate detection of metacercaria-tainted fish. Environment reconstruction includes removing toilets and pigsties from fishpond areas.

    This research was partially supported by the National Nature Science Foundation of China (No.11871238, No.11871235, No.11471133), by self-determined research funds of CCNU from the colleges'basic research and operation of MOE (Grant No.CCNU16JCZX10, No.CCNU16A02009). We thank the anonymous reviewers very much for their valuable opinions on the original manuscript.

    The authors declare no competing interests.

    [1] Johnston LD (1991) Toward a theory of drug epidemics. In: Donohew DH, Sypher H, Bukoski W, editors. Persuasive Communication and Drug Abuse Prevention. Hillsdale, NJ: Lawrence Erlbaum. pp. 93-132.
    [2] Musto DF (1987) The American disease: Origins of narcotic control. New York: Oxford.
    [3] Hunt LG, Chambers CD (1976) The heroin epidemics: A study of heroin use in the U.S., part 2, 1965-75.
    [4] Becker HS (1967) History, culture, and subjective experience: An exploration of the social bases of drug-induced experience. J Health Soc Behav 8: 163-176. doi: 10.2307/2948371
    [5] Hamid A (1992) The developmental cycle of a drug epidemic: The cocaine smoking epidemic of 1981-1991. J Psychoactive Drugs 24: 337-348. doi: 10.1080/02791072.1992.10471658
    [6] Golub A, Johnson BD (1999) Cohort changes in illegal drug use among arrestees in Manhattan: From the heroin injection generation to the blunts generation. Subst Use Misuse 34: 1733-1763. doi: 10.3109/10826089909039425
    [7] Johnson BD, Golub A (2002) Generational trends in heroin use and injection among arrestees in New York City. In: Musto D, editor. One Hundred Years of Heroin. Westport, CT: Auburn House. pp. 91-128.
    [8] Golub A, Johnson BD (1997) Crack's decline: Some surprises across U.S. cities. National Institute of Justice Research in Brief.
    [9] Golub A, Johnson BD (2001) The rise of marijuana as the drug of choice among youthful adult arrestees. National Institute of Justice Research in Brief, NCJ 187490. Washington, DC: National Institute of Justice
    [10] Golub A, Johnson BD, Dunlap E (2005) Subcultural evolution and substance use. Addict Res Theory 13: 217-229. doi: 10.1080/16066350500053497
    [11] Bennett AS, Golub A (2012) Sociological Factors and Addiction. In: Shaffer HJ, LaPlante DA, Nelson SE, editors. Addiction Syndrome Handbook. Washington, DC: American Psychological Association. pp. 195-210.
    [12] Jonnes J (1999) Hep cats, narcs, and pipe dreams: A history of America's romance with illegal drugs. Baltimore, MD: Johns Hopkins.
    [13] Malcolm X, Haley A (1966) The autobiography of Malcolm X. New York: Ballantine.
    [14] Sifaneck SJ, Kaplan CD, Dunlap E, et al. (2003) Blunts and blowtjes: Cannabis use practices in two cultural settings and their implications for secondary prevention. Free Inquiry in Creative Sociology 31: 1-11.
    [15] Rogers EM (2003) Diffusion of innovations. New York: Free Press.
    [16] Jenkins P (1999) Synthetic panics: The symbolic politics of designer drugs. New York: New York University.
    [17] Reinarman C, Levine HG (1997) Crack in context: America's latest demon drug. In: Reinarman C, Levine HG, editors. Crack in America: Demon Drugs and Social Justice. Berkeley: California. pp. 1-17.
    [18] Hartman D, Golub A (1999) The social construction of the crack epidemic in the print media. J Psychoactive Drugs 31: 423-433. doi: 10.1080/02791072.1999.10471772
    [19] Schivelbusch W (1993) Tastes of paradise: A social history of spices, stimulants, and intoxicants. New York: Vintage. xiv, 236 p. p.
    [20] Goodman J, Lovejoy PE, Sherratt A (1995) Consuming habits: Drugs in history and anthropology. Oxon, UK, and New York: Routledge. xii, 244 p. p.
    [21] Wiebe RH (1980) The search for order, 1877-1920. Westport, Conn.: Greenwood Press. xiv, 333 p. p.
    [22] Smith CS (2007) Urban disorder and the shape of belief: The great Chicago fire, the Haymarket bomb, and the model town of Pullman. Chicago: The University of Chicago Press. xix, 395 p., 324 p. of plates p.
    [23] Cowan BW (2005) The social life of coffee: The emergence of the British coffeehouse. New Haven Conn.: Yale University Press. xii, 364 p. p.
    [24] Musto DF (2002) Drugs in America: A documentary history. New York: New York University Press. xv, 574 p. p.
    [25] White WL (2004) The Lessons of Language: Historical Perspectives on the Rhetoric of Addiction. In: Acker SWTaCJ, editor. Altering American Consciousness: The Hisotry of Alcohol and Drug Use in the United States, 1800-2000.
    [26] Faupel CE, Horowitz AM, Weaver GS (2004) The Sociology of the American Drug Use. New York: Mcgraw-Hill.
    [27] Bergen-Cico DK (2011) War and Drugs: The Role of Military Conflict in the Development of Substance Abuse. Boulder, Co: Paradigm.
    [28] Courtwright DT (2001) Forces of habit : drugs and the making of the modern world. Cambridge, Mass.: Harvard University Press. viii, 277 p.
    [29] Spillane J (2001) Cocaine: From medical marvel to modern menace. Baltimore: Johns Hopkins.
    [30] Courtwright DT (2001) Dark paradise: A history of opiate addiction in America. Cambridge, Mass.: Harvard University Press. xiii, 326 p. p.
    [31] Hickman T (2004) The Double Meaning of Addiction: Habitual Narcotic Use and the Logic of Professionalizing Medical Authority. In: Acker SWTaCJ, editor. Altering American Consciousness: The History of Alcohol and Drug Use in the United States, 1800-2000.
    [32] Musto DF (1987) The American disease: origins of narcotic control. New York: Oxford University Press. xvi, 384 p. p.
    [33] Musto D (1989) America's First Cocaine Epidemic. Wilson Quarterly 13: 59-64.
    [34] Courtwright D, Joseph H, Des Jarlais D (1989) Addicts who survived: An oral history of narcotic use in America, 1923-1965. Knoxville: University of Tennessee Press. xvi, 399 p. p.
    [35] Acker CJ (2002) Creating the American junkie: Addiction research in the classic era of narcotic control. Baltimore: Johns Hopkins University Press. 276 p. p.
    [36] Schneider E (2008) Smack: Heroin and the American City. Philadelphia: University of Pennsylvania Press.
    [37] New York Times (1900) War on Opium. pp. 5.
    [38] Bennett AS (2009) The discovery and treatment of juvenile drug use in post-World-War II Los Angeles [Thesis (PhD)]. Pittsburgh, PA: Carnegie Mellon University. viii, 253 p. p.
    [39] Tracy SW, Acker CJ (2004) Altering American consciousness : the history of alcohol and drug use in the United States, 1800-2000. Amherst, Mass.: University of Massachusetts Press. vii, 414 p. p.
    [40] Miller M (1997) History and Epidemiology of Amphetamine Abuse in the United States. In: Klee H, editor. Amphetamine Misuse: International Perspectives on Current Trends. Amsterdam Harwood Academic Publishers. pp. 113-133.
    [41] Rasmussen N (2008) On Speed: The Many Lives of Amphetamine. New York: New York University Press.
    [42] New York Times (1955) Crackdown Aimed at Thrill Pills. pp. 1.
    [43] Campbell N (2000) Using Women: Gender, Drug Policy, and Social Justice. New Tork: Routledge. 270 p.
    [44] Kuzmarov J (2009) The myth of the addicted army: Vietnam and the modern war on drugs. Amherst, Mass.: University of Massachusetts Press. xii, 303 p. p.
    [45] Musto DF, Korsmeyer P (2002) The quest for drug control: Politics and federal policy in a period of increasing substance abuse, 1963–1981. New Haven: Yale University Press. xxiv, 312 p. p.
    [46] Robins LN (1974) The Vietnam Drug User Returns. Washington, DC: U.S. Government Printing Office.
    [47] Robins LN, Davis DH, Nurco DN (1974) How Permanent was Vietnam Drug Addiction? Am J Public Health 64: 38-43. doi: 10.2105/AJPH.64.12_Suppl.38
    [48] Zinberg NE (1984) Drug, set, and setting: The basis for controlled intoxicant use. New Haven: Yale University Press.
    [49] Schafer DP (1998) Culture: Beacon of the future. Westport, CT: Praeger.
    [50] Tylor EB (1871) Primitive Cultures. London: J. Murray.
    [51] Green MD (2000) The end of identify? The implications of postmodernity for political identification. Nationalism and Ethnic Politics 6: 68-90.
    [52] Allan K, Turner JH (2000) A formalization of postmodern theory. Sociological Perspectives 43:363-385. doi: 10.2307/1389533
    [53] Beck U (2000) The cosmopolitan perspective: Sociology of the second age of modernity. Brit J Sociol 51: 79-105. doi: 10.1080/000713100358444
    [54] Swidler A (1986) Culture in action: Symbols and strategies. Am Sociol Rev 51: 273-286. doi: 10.2307/2095521
    [55] Hebdige D (1979) Subculture: The meaning of style. London and New York: Routledge.
    [56] Preble EJ, Casey JJ (1969) Taking care of business: The heroin user's life on the street. Int J Addict 4: 1-24.
    [57] Williams T (1992) Crackhouse: Notes from the end of the line. Reading, MA: Addison-Wesley.
    [58] Dumit J (2002) Drugs for Life. Mol Interv 2: 124-127. doi: 10.1124/mi.2.3.124
    [59] DeGrandpre R (2000) Ritalin Nation: Rapid-Fire Culture and the Transformation of Human Consciousness. New York: Norton.
    [60] Kramer PD (1993) Listening to Prozac. New York: Viking.
    [61] Loe M (2004) The rise of viagra: How the little blue pill changed sex in America. New York: NYU.
    [62] Conrad P (2007) The medicalization of society. Baltimore: Johns Hopkins.
    [63] Crister G (2005) Generation Rx: How prescription drugs are altering American lives, minds and bodies. Boston: Houghton Mifflin.
    [64] Larance B, Degenhardt L, Lintzeris N, et al. (2011) Definitions related to the use of pharmaceutical opioids: Extramedical use, diversion, non-adherence and aberrant medication-related behaviours. Drug Alcohol Rev30: 236-245.
    [65] Fainaru-Wada M, Williams L (2006) Game of shadows: Barry Bonds, BALCO, and the steroids scandal that rocked professional sports. New York: Gotham.
    [66] Murray TH, Maschke KJ, Wasunna AA (2009) Performance–Enhancing Technologies in Sports: Ethical, Conceptual, and Scientific Issues. Baltimore: Johns Hopkins.
    [67] Finkel D (2009) The Good Soldiers. New York: Picador.
    [68] Senior J (2011) The Prozac, Paxil, Zoloft, Wellbutrin, Celexa, Effexor, Valium, Klonopin, Ativan, Restoril, Xanax, Adderall, Ritalin, Haldol, Resperdal, Seroquel, Ambien, Lunesta, Elavil, Trazodone War. New York Magazine.
    [69] Fields G (2009) White House Czar Calls for End to ‘War on Drugs'. Wall Street Journal.
    [70] ONDCP (2011) National Drug Control Strategy. Washington, DC: Office of National Drug Control Policy.
    [71] Hanson DJ (1995) Preventing alcohol abuse: Alcohol, culture and control. Westport, CT: Praeger.
    [72] Hanson DJ (1996) Alcohol education: What we must do. Westport, CT: Praeger.
    [73] Babor T (1986) Alcohol: Customs and rituals. New York: Chelsea House.
    [74] Mandelbaum DG (1965) Alcohol and culture. In: Marshal M, editor. Beliefs, behaviors and alcoholic beverages: A cross-cultural survey. Ann Arbor, MI: University of Michigan Press. pp.14-30.
    [75] World Health Organization (2010) European Status Report on Alcohol and Health 2010.
  • This article has been cited by:

    1. Martha Betson, Allen Jethro I. Alonte, Rico C. Ancog, Angelou Marie O. Aquino, Vicente Y. Belizario, Anna Monica D. Bordado, Jessica Clark, Ma. Christina G. Corales, Mary Grace Dacuma, Billy P. Divina, Matthew A. Dixon, Stephen A. Gourley, Jasmine Renette D. Jimenez, Ben P. Jones, Sheina Macy P. Manalo, Joaquin M. Prada, Arnoud H.M. van Vliet, Kezia C.L. Whatley, Vachel Gay V. Paller, 2020, 108, 9780128207505, 47, 10.1016/bs.apar.2020.01.036
    2. Xiao-Hong Huang, Men-Bao Qian, Guang-Hu Zhu, Yue-Yi Fang, Yuan-Tao Hao, Ying-Si Lai, jong-Yil Chai, Assessment of control strategies against Clonorchis sinensis infection based on a multi-group dynamic transmission model, 2020, 14, 1935-2735, e0008152, 10.1371/journal.pntd.0008152
    3. Alongkot Suvarnamani, Puntani Pongsumpun, 2021, Analyze of the Model for Cancer Transmission, 9781450390040, 77, 10.1145/3469951.3469965
    4. Yun-Ting He, Xiao-Hong Huang, Yue-Yi Fang, Qing-Sheng Zeng, Lai-De Li, Le Luo, Ying-Si Lai, jong-Yil Chai, Cost-effectiveness evaluation of different control strategies for Clonorchis sinensis infection in a high endemic area of China: A modelling study, 2022, 16, 1935-2735, e0010429, 10.1371/journal.pntd.0010429
    5. Lars Kamber, Christine Bürli, Helmut Harbrecht, Peter Odermatt, Somphou Sayasone, Nakul Chitnis, Krystyna Cwiklinski, Modeling the persistence of Opisthorchis viverrini worm burden after mass-drug administration and education campaigns with systematic adherence, 2024, 18, 1935-2735, e0011362, 10.1371/journal.pntd.0011362
    6. Yuhong Guo, Zhijun Liu, Lianwen Wang, Ronghua Tan, Modeling the role of information and optimal control on an SEIR epidemic model with infectivity in latent period, 2024, 47, 0170-4214, 1044, 10.1002/mma.9698
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(7671) PDF downloads(1361) Cited by(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog