- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Computational and Mathematical Methods in Medicine
Volume 2013 (2013), Article ID 753424, 12 pages
Mathematical Modeling of the HIV/Kaposi’s Sarcoma Coinfection Dynamics in Areas of High HIV Prevalence
1Department of Mathematics, University of Botswana, Gaborone, Botswana
2Department of Mathematics, University of Tennessee, TN, USA
3Department of Mathematics, University of Dar es Salaam, Dar es Salaam, Tanzania
4Department of Applied Mathematics, National University of Science and Technology, Bulawayo, Zimbabwe
Received 12 May 2013; Revised 29 August 2013; Accepted 12 September 2013
Academic Editor: Gabriel Turinici
Copyright © 2013 E. Lungu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
We formulate a deterministic system of ordinary differential equations to quantify HAART treatment levels for patients co-infected with HIV and Kaposi's Sarcoma in a high HIV prevalence setting. A qualitative stability analysis of the equilibrium states is carried out and we find that the disease-free equilibrium is globally attracting whenever the reproductive number . A unique endemic equilibrium exists and is locally stable whenever . Therefore, reducing to below unity should be the goal for disease eradication. Provision of HAART is shown to provide dual benefit of reducing HIV spread and the risk of acquiring another fatal disease for HIV/AIDS patients. By providing treatment to 10% of the HIV population, about 87% of the AIDS population acquire protection against coinfection with HIV and Kaposi's Sarcoma (KS). Most sub-Sahara African countries already have programmes in place to screen HIV. Our recommendation is that these programmes should be expanded to include testing for HHV-8 and KS counseling.
Kaposi's Sarcoma is a cancer that occurs mostly in humans with suppressed immune systems . The development of this cancer depends upon prior infection to the human herpesvirus-8 (HHV-8) , a virus which is usually transmitted either sexually or via saliva . For HIV-related Kaposi's Sarcoma development, immunosuppression is a necessary causal factor . Because HIV is an immunosuppressive virus, it promotes the development of Kaposi's Sarcoma in individuals dually infected with both viruses (HIV and HHV-8) and this combination has proved to be fatal and has made Kaposi's Sarcoma the fourth largest killer of people living with HIV/AIDS in sub-Sahara Africa .
In competent immune systems, acquisition of HHV-8 does not guarantee the development of KS; in fact, most individuals with a strong immune response could remain latently infected with HHV-8 throughout their lifetime . The HIV-1 growth factors stimulate the immune cells including the healthy and infected B-cells to proliferate. The activation of latently infected B-cells to proliferate only leads to production of more HHV-8 that, according to the theory proposed by Foreman et al. , may be responsible for infection of progenitor cells of endothelial origin, which once infected with HHV-8 develop into Kaposi's Sarcoma cells .
Recent research findings  show that highly active antiretroviral therapy (HAART) for HIV significantly decreases (and, in some instances, completely forces the cancer into remission) KS activity in a patient, but such treatment is only effective against persons who have seroconverted. Lungu et al. studied the within-host dynamics of HIV and KS and came to the conclusion that HIV infection accentuates the potential for infected individuals to develop KS conditions and that administration of HAART on KS individuals who have seroconverted results in the reversal of KS conditions. However, the same result could not be demonstrated if HAART was administered to HIV negative individuals.
From what has been described above, we formulate a mathematical model which includes the following classes: a susceptible class, , two infected classes of individuals infected with HIV-1 only, , individuals coinfected with HHV-8 and HIV-1, , two classes of asymptomatic infectives with HIV only, , asymptomatic coinfected, , a class of treated individuals, , and two AIDS classes, namely, individuals with full-blown AIDS only, , and individuals coinfected with full-blown AIDS and Kaposi's Sarcoma.
The structure of the paper is as follows. In Section 2 we formulate the model that describes the coinfection transmission dynamics. Existence of solutions is proved in Section 3, but our main interest is to simulate and study the dynamic behavior of the steady state solution. In particular, we want to quantify the proportion of HAART-treated individuals who never develop KS at various treatment rates and .
AIDS is the progressive stage at which an HIV-infected individual loses competency of his immune system and becomes prone to opportunistic infections. For this reason, HIV coinfection models have found a lot of space in HIV epidemiology studies. HIV and tuberculosis coinfection have been studied by Cohen et al. , Roeger et al. , and Ramkissoon et al.  while Barley et al. , Chiyaka et al. , and Mukandavire et al.  studied HIV-malaria coinfection dynamics with each study giving results pertinent to the coinfection under investigation. The study of HIV-KS coinfection dynamics is still in infancy and, to the best of our knowledge, no mathematical modeling study has been carried out to assess the coinfection dynamics of HIV and KS at the population level in sub-Sahara Africa. Because our model incorporates HAART administration to all infective classes, it can provide insights into treatment strategies and, in particular, decide whether the current policy based on a CD4 count threshold to access treatment should be continued as a strategy.
2. Model Formulation
We begin with a human population of susceptibles which is free of both KS and HIV denoted by . This population is replenished at constant rate through sexual maturity or immigration. Upon effective contact with individuals infected with the HIV virus, the new infectives progress into the infected class at rate , where denotes the force of infection.
For simplicity of the model, we assume that all individuals in the class are latently infected with the HHV-8 virus which causes the KS infection. Therefore, in our model every class, except for the susceptible class, is assumed to be at risk of developing KS. HIV infection is known to promote or enhance the development of KS condition [18, 19], and so we assume that individuals in the class can develop mild KS at rate and move into the class of coinfected individuals KS denoted by .
HIV-only and coinfected individuals are assumed to progress to the asymptomatic pre-AIDS classes ( and ) at the same constant rate . Furthermore, persons in the pre-AIDS class are deemed to be sexually interacting and that individuals with mild KS in the class can develop acute KS, which manifests in the form of visible lesions  and severe debilitation. Due to the nature of these symptoms, we assume that individuals who develop acute KS will no longer be sexually interacting. Those individuals in the class, who develop acute KS, die at a disease induced rate .
We assume that all infected individuals in the classes and are tested and if they are tested positive for HIV-1 they receive treatment at rate and move into the class of treated individuals, . Individuals in the pre-AIDS classes and have higher viral loads and therefore possess less competent immune systems. These individuals will then present for medical attention at a higher rate . HAART is known to reverse KS conditions for people with HIV , and so in the model we assume that the therapy is perfect in reversing mild KS conditions. Upon accessing HAART, individuals in the infective class and pre-AIDS class achieve full KS recovery and progress to the treated class . Additionally, individuals in the pre-AIDS class progress to acute KS at rate . Individuals in the pre-AIDS class are at risk of developing KS at a constant rate . Individuals in the pre-AIDS classes and develop clinical symptoms and progress to full-blown AIDS and , respectively, at rate . Progression to full-blown AIDS by persons in the class represents treatment failure, and they progress to the class at the same rate . Full-blown AIDS individuals can also develop KS at rate due to weakened immune systems. Persons in the full-blown AIDS classes have additional AIDS-induced mortality . Progression rate to acute KS for persons with full-blown AIDS is , which we assume to be greater than and . Individuals in both AIDS classes are also subject to a natural mortality rate of .
The force of infection depends on the probability of transmission per contact , the proportion of infected individuals in each category ( and ), the proportion of infected individuals in receipt of HAART , and the pre-AIDS classes ( and ). Individuals in receipt of HAART have reduced viral load  and are therefore assumed to be less infectious relative to infectives not in receipt of HAART. This reduced infectiousness is modeled by . Coinfected infectives and pre-AIDS individuals ( and ) have weaker immune systems and are therefore likely to carry higher viral loads than their counterparts in the ( and ) classes. This added infectiousness is modeled by and . Pre-AIDS infectives are more infectious because of increased viral load, and this increased infectiousness is modeled by where . Persons with full-blown AIDS exhibit symptoms related to HIV and therefore are assumed to be noninteracting. The total sexually active variable population at time is given by . Assuming homogeneous mixing, the time dependent force of infection for HIV is given by
The model flow diagram depicting this biological system is illustrated in Figure 1.
The above assumptions and formulations lead us to this nonlinear system of differential equations:
All parameters for the model system (2) are assumed to be nonnegative for all time .
3. Basic Properties
3.1. Positivity and Boundedness of Solutions
We denote by the set of points in with positive coordinates and consider the system (2) with initial values
In this section, we prove the following lemma.
Lemma 1. The system (2) can be written as a system of differential inequalities where If then for all and .
Without loss of generality we may assume that , since the case follows by approximating the system with a sequence .
Proof. Suppose that the assertion , for , is not true. Then there exists the smallest number , such that
Then is a decreasing function at , so that
From the differential inequality (4) for we get which is a contradiction.
For the state variables in our model, we always take
From Lemma 1 we conclude that
Thus, in the region the model is epidemiologically and mathematically well posed and we can use it to study the dynamics of HIV-KS coinfection.
3.1.1. Equilibrium States, Reproductive Number, and Stability Analysis
The model (2) possesses a disease-free equilibrium, , given by and at least one endemic equilibrium state whose existence is discussed in Section 3.3. Following Van Den Driessche and Watmough , the KS-induced reproductive number of system (2), , is given by the spectral radius of the matrix where the matrices and are given by where
From (12) we have calculated the reproduction number where
This number is a threshold such that if the disease clears from the population. If the steady state becomes unstable and the disease establishes itself into the population. This number is comprehensively analyzed further in Section 3.4 to reveal the impact of treatment.
3.2. Global Stability of the Disease-Free Equilibrium
If then fixed point is locally asymptotically stable. We now determine conditions which guarantee global asymptotic stability of the disease-free state . Rewrite model system (2) as where denotes the number of uninfected individuals and denotes the number of infected individuals including those latently infected and those who are infectious. The disease-free equilibrium state can now be written as where . To guarantee global stability, the following conditions must be satisfied: (i)for , is globally asymptotically stable; (ii) for ; (iii) is an -matrix (the off diagonal elements of are nonnegative) and is the region where the model makes biological sense.
Lemma 2. The fixed point is a globally asymptotically stable point of (2) provided that and conditions stated above are satisfied.
The biological interpretation of Lemma 2 is that this population is observing the one partner policy or abstinence from unprotected sex resulting in HIV being eliminated from the population regardless of the size of the initial subpopulations .
3.3. Existence of Endemic Equilibria and Bifurcation Analysis
In this section following the approach in , we establish conditions for the existence of endemic equilibria and investigate their stability. Denote the arbitrary endemic equilibrium of model (2) by where
Let be the force of infection. By setting the right hand side of model (2) to zero, we obtain, in terms of the force of infection,
From the solution , the condition is necessary for existence of the nontrivial endemic equilibrium and we summarize this as follows.
Lemma 3. The model (2) has a unique endemic equilibrium whenever .
Having proved the existence of the endemic equilibrium point we now investigate its stability using the Centre Manifold Theory  as described by Castillo-Chavez and Song . Consider the following general system of ordinary differential equations with a parameter
Without loss of generality, it is assumed that 0 is an equilibrium for system (27) for all values of the parameter , (i.e., for all .
Assume the following.(A1): is the linearise matrix of system (27) around the equilibrium 0 with evaluated at . Zero is a simple eigenvalue of and all other eigenvalues of have negative real parts. (A2): Matrix has a nonnegative right eigenvector and a left eigenvector corresponding to the zero eigenvalue.
Let be the th component of and
The local dynamics of system (27) around are totally determined by and . (i), . When with , is locally asymptotically stable and there exists a positive unstable equilibrium; when , is unstable and there exists a negative and locally asymptotically stable equilibrium. (ii), . When with , is unstable; when , is locally asymptotically stable, and there exists a positive unstable equilibrium. (iii), . When with , 0 is unstable, and there exists a locally asymptotically stable negative equilibrium; when , is stable, and a positive unstable equilibrium appears. (iv), . When changes from negative to positive, changes its stability from stable to unstable. Correspondingly, a negative unstable equilibrium becomes positive and locally asymptotically stable. Particularly, if and , then, a backward bifurcation occurs at .
We make the following change of variables: , , , , , , , and so that . We now use the vector notation . Then, model system (2) can be written in the form , such that where
The Jacobian matrix of system (29) at the disease-free equilibrium is given by from which it can be shown that is the same as in (14). If is taken as a bifurcation point and if we consider the case and solve for gives
The linearised system of the transformed equations (29) at has a simple zero eigenvalue. Hence, the Centre Manifold Theory  can be used to analyse the dynamics of system (29) near . It can be shown that the Jacobian of (29) at has a right eigenvector associated with the zero eigenvalue given by , where
The left eigenvector of associated with the zero eigenvalue at is given by where
For the model system (2), the associated nonzero partial derivatives of at the disease-free equilibrium are given by
From (28), it follows that
For the sign of , it is associated with the nonvanishing partial derivatives of ,
From (28), it follows that
3.4. Analysis of the Reproductive Number
To analyze the coinfection dynamics of KS and HIV, we investigate the KS-induced reproductive number in (14). In (14), we note that the partial reproductive numbers , , , , and represent the contribution to the reproduction number from the following classes of infectives HIV only (), coinfected (), Pre-AIDS (), coinfected pre-AIDS (), and those receiving HAART treatment (). The partial reproductive number representing the contribution of groups in receipt of ART is split into , , , , and and represents the contribution of HIV infectives, infectives dually infected with KS, Pre-AIDS infectives, and Pre-AIDS infectives dually infected with KS, respectively, who receive anti-retroviral therapy.
If we define to be the basic HIV reproductive number in the absence of KS and HIV treatment and set parameters related to KS and treatment to zero in (14) we obtain
If we define the partial reproductive number without treatment to be and set parameters related to treatment to zero then we obtain
If , then not providing treatment to all infective classes worsens the HIV in this population. Note that even if , the option not to provide treatment to all infective classes could increase the reproduction number above 1. Treatment to all infective classes is necessary to effectively control the HIV disease spread. Using sensitivity analysis on , we have reinforced some of these findings as shown in Figure 2.
From Figures 2(a) and 2(b) wem note that an increase in the rate of infected people who develop KS will lead to an increase in and consequently an increase in the HIV epidemic. However, an increase in the rate of progress to dual infection for people in the pre-AIDS class will result in a marginal increase of and hence the HIV epidemic.
Figures 2(c) and 2(d) show that treatment of infected individuals and dually infected individuals will have a much higher positive impact on and the epidemic than treatment at the pre-AIDS stage and this leads us to conclude that early administration of HAART on HIV infectives to curtail the growth of opportunistic infections, such as KS, will have a more positive impact than delayed therapy at the pre-AIDS stage.
4. Numerical Simulations
Using the R programming environment, we ran numerical simulations of the model. The following data were input as initial conditions:
Parameter values used in the numerical simulations of model system (2) are shown in Table 1. HIV/KS coinfection study is still in its infancy. A number of numeric values for the parameters shown in Table 1 are reasonable estimates. Some of the parameters were estimated by modifying baseline values from published literature and these are denoted with an asterisk in Table 1.
Our value for the rate of HIV transmission is the average of the minimum (0.011) and maximum values (0.95) for the same parameter in Baggaley et al.  and Boily et al. . We assumed that a death rate due to AIDS and KS coinfection would be ; in our model, and are summed. Our recruitment rate is calculated by incorporating an annual recruitment rate (0.029) from Malunguza et al. . In this study, we take the population of sub-Sahara Africa to be 767 million as in Mukandawire et al.  and Barley et al. .
Figure 3(a) shows that in the absence of treatment, the number of individuals with the AIDS/KS coinfection exceeds the number of those with just AIDS; that is, 54.0% of the total AIDS population has KS by the end of the model. Merely affecting a treatment rate of , is enough to reverse this relation as depicted in Figure 3(b) where of all AIDS patients are devoid of KS. If the treatment rate is increased to , the AIDS population without the coinfection increases to as depicted in Figure 3(c). Using parameter values in Table 1 causes the model to tend towards some endemic equilibrium over time. In other words, the DFE is unstable, and so we expect . In fact, carrying out these computations with , we obtain . This value is slightly higher than estimated values of for HIV/AIDS in European countries, which is consistent with our expectations .
Figure 3(d) provides an interesting insight into the predictions of our model. By providing treatment to of all infectives and pre-AIDS individuals the model fails to eradicate the coinfection, and approximately an prevalence of the coinfection among the full-blown AIDS cohort exists. This result means that we need to be cognisant of two objectives, namely, that we should have benchmark target goals in the provision of treatment and that work should be done to decrease the number of patients entering the pre-AIDS coinfection class. It is entirely unrealistic to presume that the coinfection can be eliminated entirely based on providing treatment to infectives and pre-AIDS individuals alone. In fact, Figure 3(d) shows that providing a treatment rate of results in about of the AIDS population losing the coinfection. Providing any further increases to and results in only a increase.
Our model makes no attempt to consider the cost associated with providing the kind of treatment we have discussed. However, it is noteworthy to consider that the gains associated with providing treatment to any more than of all infectives and pre-AIDS individuals are minimal. It is clear that there is an gap between where we would like to see the percentage of the full-blown AIDS cohort without coinfection and where our model currently estimates that percentage to be. Assuming that it is unrealistic to completely eradicate the coinfection, steps should be taken to shift the curve in Figure 3(d) closer to . We believe that a good way to do this is by screening patients for the coinfection once the HIV infection is first diagnosed. In general, this will serve to decrease the total number of patients who become pre-AIDS with the coinfection and, in turn, this will decrease the total number of individuals presenting with the coinfection.
In numerous HIV-positive cohorts, susceptibility to KS development is incredibly high. Recent findings that antiretroviral treatment for HIV clinical symptoms can reverse KS in most patients suggest that the only obstacle in preventing KS-related complications (or even deaths) is one's ability to access the treatment itself. This is especially the case in sub-Saharan Africa, where KS prevalence is higher than anywhere else in the world, yet access to proper treatment remains low. We have shown that providing treatment to just of HIV-infected individuals, regardless of their KS status, can offer a significant reduction in the overall number of individuals who end up with full-blown AIDS and KS. Certainly, we expect that, as we increase the treatment rate, the disparity will become more and more favorable. However, our results only apply to situations in which individuals develop KS on their own—that is, in the absence of interaction with HHV-8-infected individuals. While such a generalization may be perfectly valid in sub-Saharan Africa, it will not apply to other regions where HHV-8 needs to be acquired sexually or through saliva. In other words, this model is very specific to the sub-Saharan region.
In the future, we hope to extend this model to consider all modes of HHV-8 transmission. One way to do this is by introducing separate classes of HHV-8 infectious individuals. This will increase the overall reach of the model and hopefully provide a positive influence on policy making in areas affected by KS. For now, supporting KS education and awareness needs to become as important as providing universal access to antiretroviral treatment for HIV patients. Most southern African countries already have programs in place to encourage HIV screening; this is the perfect time to also test for HHV-8 and provide KS counseling. Regardless of how long the HIV/AIDS epidemic lasts in Africa, there is no reason people should have to fight KS as well.
This work was done under SAMSA Masamu Program with support from the National Science Foundation under Grant no. DMS-1050259. Lungu acknowledges the support of the University of Botswana, Massaro would like to acknowledge partial support from the NSF SCALE-IT graduate fellowship program at the University of Tennessee, and Malunguza acknowledges the support of the National University of Science and Technology.
- Biotrin International, “HHV-8: questions and answers,” 2005, http://www.diasorin.com/HHV8/FAQ.html.
- K. E. Foreman, J. Friborg Jr., W.-P. Kong et al., “Propagation of a human herpesvirus from AIDS-associated Kaposi's sarcoma,” New England Journal of Medicine, vol. 336, no. 3, pp. 163–171, 1997.
- D. C. Dugdale III, J. M. Vyas, D. Zieve, et al., “Kaposi’s sarcoma,” PubMed, 2012, http://www.ncbi.nlm.nih.gov/pubmedhealth/PMH0001682.
- F. Sitas, M. Parkin, Z. Chirenje, L. Stein, N. Mqoqi, and H. Wabinga, “Cancers,” in Disease and Mortality in Sub-Saharan Africa, chapter 20, 2nd edition, 2013.
- R. F. Baggaley, R. G. White, and M. C. Boily, “Systematic review of orogenital HIV-1 transmission probabilities,” International Journal of Epidemiology, vol. 37, no. 6, pp. 1255–1265, 2008.
- M.-C. Boily, R. F. Baggaley, L. Wang et al., “Heterosexual risk of HIV-1 infection per sexual act: systematic review and meta-analysis of observational studies,” The Lancet Infectious Diseases, vol. 9, no. 2, pp. 118–129, 2009.
- N. J. Malunguza, S. D. Hove-Musekwa, G. Musuka, and Z. Mukandavire, “Investigating al-cohol consumption as a risk factor for HIV transmission in heterosexual settings in sub-saharan african communities,” Bulletin of Mathematical Biology, vol. 74, no. 9, pp. 2094–2124, 2012.
- Z. Mukandavire, C. Chiyaka, G. Magombedze, G. Musuka, and N. J. Malunguza, “Assessing the effects of homosexuals and bisexuals on the intrinsic dynamics of HIV/AIDS in heterosexual settings,” Mathematical and Computer Modelling, vol. 49, no. 9-10, pp. 1869–1882, 2009.
- K. Barley, D. Murillo, S. Roudenko, A. M. Tameru, and S. Tatum, “A mathematical model of HIV and malaria co-infection in sub-saharan Africa,” Journal of AIDS and Clinical Research, vol. 3, no. 7, pp. 1–7, 2012.
- Z. Mukandavire and W. Garira, “Sex-structured HIV/AIDS model to analyse the effects of condom use with application to Zimbabwe,” Journal of Mathematical Biology, vol. 54, no. 5, pp. 669–699, 2007.
- Z. Mukandavire, N. J. Malunguza, C. Chiyaka, G. Musuka, and J. M. Tchuenche, “HIV/AIDS model assessing the effects of gender-inequality affecting women in african heterosexual settings,” International Journal of Biomathematics, vol. 3, no. 1, pp. 43–67, 2009.
- E. M. Lungu, “Stability of non-autonomous co-infection models,” in Proceedings of the 7th ICIAM Congress, 2011, https://www.masamu.auburn.edu/ConMan/ConMan.
- T. Cohen, M. Lipsitch, R. P. Walensky, and M. Murray, “Beneficial and perverse effects of isoniazid preventive therapy for latent tuberculosis infection in HIV-tuberculosis coinfected populations,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 18, pp. 7042–7047, 2006.
- L.-I. W. Roeger, Z. Feng, and C. Castillo-Chavez, “Modeling TB and HIV co-infections,” Mathematical Biosciences and Engineering, vol. 6, no. 4, pp. 815–837, 2009.
- S. Ramkissoon, H. G. Mwambi, and A. P. Matthews, “Modelling HIV and MTB co-infection including combined treatment strategies,” PLoS One, vol. 7, no. 11, Article ID e49492, 2012.
- C. Chiyaka, J. M. Tchuenche, W. Garira, and S. Dube, “A mathematical analysis of the effects of control strategies on the transmission dynamics of malaria,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 641–662, 2008.
- Z. Mukandavire, A. B. Gumel, W. Garira, and J. M. Tchuenche, “Mathematical analysis of a model for HIV-malaria co-infection,” Mathematical Biosciences and Engineering, vol. 6, no. 2, pp. 333–362, 2009.
- J. H. Von Roenn, “Treatment of HIV-associated kaposi sarcoma,” 2003, http://hivinsite.ucsf.edu/InSite?page=kb-06-02-04.
- D. P. Wilson, G. L. Matthew, A. E. Grulich, D. A. Cooper, and J. A. Kaldor, “Relation between HIV viral load and infectiousness: a model-based anal-ysis,” The Lancet, vol. 372, no. 9635, pp. 314–320, 2013.
- P. Van Den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.
- C. Castillo-Chavez, Z. Feng, and W. Huang, “On the computation of ℛ0 and its role on global stability,” http://math.la.asu.edu/chavez/2002/JB276.pdf.
- O. Sharomi and A. B. Gumel, “Re-infection-induced backward bifurcation in the transmission dynamics of chlamydia trachomatis,” Journal of Mathematical Analysis and Applications, vol. 356, no. 1, pp. 96–118, 2009.
- J. Carr, Applications Centre Manifold Theory, Springer, New York, NY, USA, 1981.
- C. Castillo-Chavez and B. Song, “Dynamical models of tuberculosis and their applications,” Mathematical Biosciences and Engineering, vol. 1, no. 2, pp. 361–404, 2004.
- H. Nishiura, “Correcting the actual reproduction number: a simple method to estimate R0 from early epidemic growth data,” International Journal of Environmental Research and Public Health, vol. 7, no. 1, pp. 291–302, 2010.