Can we spend our way out of the AIDS epidemic? A world halting AIDS model
BMC Public Health volume 9, Article number: S15 (2009)
There has been a sudden increase in the amount of money donors are willing to spend on the worldwide HIV/AIDS epidemic. Present plans are to hold most of the money in reserve and spend it slowly. However, rapid spending may be the best strategy for halting this disease.
We develop a mathematical model that predicts eradication or persistence of HIV/AIDS on a world scale. Dividing the world into regions (continents, countries etc), we develop a linear differential equation model of infectives which has the same eradication properties as more complex models.
We show that, even if HIV/AIDS can be eradicated in each region independently, travel/immigration of infectives could still sustain the epidemic. We use a continent-level example to demonstrate that eradication is possible if preventive intervention methods (such as condoms or education) reduced the infection rate to two fifths of what it is currently. We show that, for HIV/AIDS to be eradicated within five years, the total cost would be ≈ $63 billion, which is within the existing $60 billion (plus interest) amount raised by the donor community. However, if this action is spread over a twenty year period, as currently planned, then eradication is no longer possible, due to population growth, and the costs would exceed $90 billion.
Eradication of AIDS is feasible, using the tools that we have currently to hand, but action needs to occur immediately. If not, then HIV/AIDS will race beyond our ability to afford it.
There has been a sudden increase in the amount of money donors are willing to spend on the worldwide HIV/AIDS epidemic. Present plans are to hold most of the money in reserve and spend it slowly. However, rapid spending may be the best strategy for halting this disease. The purpose of this paper is to estimate how and when to best spend the available money to achieve halting HIV/AIDS as soon as possible.
In the year 2000, leaders from around the globe derived a set of goals aimed at making the world safer, healthier and more equitable. One of these goals was to halt and begin to reverse the global HIV/AIDS epidemic by 2015. To realise this goal, in 2006, world leaders pledged to move towards universal access to HIV prevention, treatment, care and support by 2010. In the past decade, there has been a marked decrease in the number of AIDS-related deaths worldwide and the percentage of individuals living with HIV has begun to plateau. Indeed, a number of highly infected countries have seen declines in the number of new HIV infections . This success is owed largely to a six-fold increase in financing for HIV activities in low- and middle-income countries enabling better access to HIV treatments and care. However, the picture is not uniformly positive and some countries continue to see an increase in infections. Also, in spite of the improved access to HIV treatment and care, the majority of those who need antiretroviral therapies are not currently receiving them . As a result, the epidemic is outpacing the rate at which these drugs are being delivered. In fact, the number of new infections is significantly greater than the increase in the number of people on antiretrovirals. Eradication of HIV/AIDS is therefore not achievable using the current funding model. In the absence of additional funding, we contend that a significant improvement in results can be achieved by optimising implementation, i.e. how donated funds are spent on drug delivery, educational programs and the distribution of HIV/AIDS prevention measures.
While estimates on spending have been made for many scenarios [2–4], these are usually projected in decades or generations. World events could overtake the HIV/AIDS epidemic on this time scale, or even be precipitated by it if it continues to grow [5–9]. Furthermore, donor fatigue  could return, or other crises (tsunamis, wars, etc) could divert presently committed funds. Is there a quick fix that could avoid passing the problem on to the next generation?
Consider, for example, Bill and Melinda Gates, who have made stopping AIDS the top priority of their foundation , with $60 billion available . As an extreme case, for instance, ignoring the ability of economies to absorb so much cash , the Gates Foundation might spend all their money in two years, meeting the worldwide needs of $27 billion per year for all health needs of all developing countries  (including HIV/AIDS), with the hope that infrastructure would then be in place and the rest of the donors would pick up from there. The projected spending rate of the Gates Foundation is $3 billion per year , which is roughly comparable to the expected annual interest. It is being spread out over 20 years, more than doubling the duration of the epidemic, which has already lasted over 25 years. Something in between a jump start and a trickle would seem much more optimal than either of these extremes.
Furthermore, if we fail to spend at the optimal rate, the epidemic may blow up. Currently, the percentage of people living with HIV has stabilised; however, the overall number of people living with HIV has steadily increased . The result is that the epidemic could become totally out of control or cost immensely more to bring under control. Many modeling studies have shown a bigger bang for the buck if the money is spent earlier in an epidemic [16, 17], but none have investigated the optimal course of spending to bring an epidemic under control as rapidly as possible and avoid a growing gap .
A number of views have been expressed on how to tackle HIV/AIDS. One approach involves getting back to the fundamentals of epidemiology  and avoiding AIDS exceptionalism , i.e. treating HIV/AIDS just like other epidemics. Most people advocate a mixture of prevention and treatment options [17, 21]. One worldwide single (but not magic) bullet approach that has been proposed would use the fact that highly active antiretroviral therapy (HAART) reduces infectivity:
"This optimistic population-based model shows that, in 45 years, HIV prevalence could be reduced by more than 70 times from more than 7 cases per 1000 people to less than 0.1 case per 1000. The number of HIV- infected people could be reduced from 38 million to less than 1 million. The cost of therapy would be about US $7 billion per year, with costs declining from $15 to $1 billion. Such a programme would be expected to cost $338 billion over 45 years" .
The average cost over the 45 years would be $7.5 billion per year.
Condoms offer another single-mode approach. Above an estimable and possibly attainable threshold, condoms could halt an HIV/AIDS epidemic . If distributed to 1.5 billion men (half the male population of the world), at 100 condoms per person per year , and a cost of $0.02 each , total cost would be $3 billion per year. (Female condoms have their place, but are about ten times more expensive each.) For comparison, the USA spends $12 billion per year domestically on HIV/AIDS, plus $3 billion per year on HIV/AIDS research .
Condoms have the advantage of reducing other STIs (sexually transmitted infections), some of which enhance HIV infectivity. The threshold for decline of the epidemic is the product of compliance in use and condom effectiveness, and there is room for technical improvement of the latter . The problems to be tackled are not simply behavioural, but also involve design, manufacturing, marketing and subsidisation . Even in the behavioural realm, the emphasis on risk factors instead of collective behaviour may be counterproductive [28, 29].
The HIV/AIDS epidemic is often spoken of in terms of reducing the spread , or achieving sustainable financing . Decades-long investments have been made in vaccines  and more recently in microbicides . These approaches are slow in coming, and do not take into account the greatly increased amount of cash now available, at least for the short term. The donor community has created the potential of stepping ahead of the exploding HIV/AIDS epidemic by pooling money faster than the disease is spreading. Our aim is to use mathematical modelling to find ways to spend the available money as fast, efficiently and effectively as possible to bring this epidemic to a full halt, now, while this transient advantage is still available.
A large number of mathematical models have been used to describe the epidemiology of HIV [33–46]. These range from the simple to the complicated, and have used a variety of techniques to analyse and determine conclusions. One such technique is that of metapopulation models: dividing space into distinct regions, within which disease transmission occurs (often at different rates) and with "travel" between regions allowing different regions to be connected. Such models have been used to describe HIV in a spatial context without resorting to partial differential equations [47–53].
This paper presents a simple, linear, metapopulation model that describes the world regions at a number of different scales (continent, country, etc). We use this model to show that interventions can lead to eradication, but such interventions must be applied globally, not locally, and that time is of the essence if our interventions are to succeed.
An eradication model
In its simplest form, HIV infection can be described by the classical SI (susceptible-infected) model. The dynamics are not linear, depending, at a minimum, on the behaviour of susceptible (uninfected) individuals, as well as their interactions with those who are infected. While the SI model does not capture all the dynamics of HIV, it has been used to model the disease throughout its entire history [54–60].
Consider the two-dimensional SI model
Here, Λ is the rate of appearance of new susceptibles, β is the infection rate, μ is the background death rate and γ is the death rate due to disease.
This model has two equilibria,
The Jacobian matrix for this system is
with eigenvalues λ = -μ, - μ - γ. It follows that the threshold for eradication  in the SI model is
The SI model, while it captures the dynamics of HIV spread in a population, is difficult to manipulate in a metapopulation scenario. We propose instead a linear approximation of this model. The approximation of the epidemic as linear in the number of infectives is clearly a coarse one. However, although this model structure does not capture the transient dynamics of infection and interaction, it nevertheless serves as a predictor for eradication.
This is because the major stability analysis carried out on such models - finding the basic reproductive ratio - consists of linearisation of the nonlinear models. Thus, a linear version of the model has the advantage of having the same invasion threshold, although it has the disadvantage of only describing local stability of the disease-free equilibrium.
Since , it follows that the one-dimensional model
will always overestimate the epidemic. We shall refer to this as the I-only model.
This model has only the trivial equilibrium I = 0. If the trivial equilibrium is stable, then all trajectories will approach it. If the trivial equilibrium is unstable, then solutions will increase without bound.
Although the total population without infection is not constant, it should nevertheless be noted that the I-only model has the eradication threshold
which is the same eradication threshold as the SI model (as expected). It follows that there will be eradication in the I-only model if and only if there is eradication in the SI model.
To illustrate, we simulated two cases: R0,SI< 1 (Figure 1) and R0, SI > 1 (Figure 2). Parameters used in the simulations were Λ = 20 people·years-1, μ = years-1, γ = years-1, with β = 0.00007 people-1years-1 (Figure 1) and β = 0.0002 people-1years-1 (Figure 2). Despite the fact that the transient dynamics are vastly different in the two models, the linear approximation has the same eradication threshold as the more accurate SI model and, furthermore, always overestimates the epidemic. It follows that, for eradication purposes, the simple, linear model should determine whether our control methods will be sufficient.
We now extend the SI model to a metapopulation model with p regions. For 1 ≤ i ≤ p, we can write
S i = susceptible individuals in the ith region;
I i = infected individuals in the ith region;
Λ i = the rate of appearance of new susceptible individuals in the ith region;
μ i = background death rate in the ith region;
γ i = death rate due to disease in the ith region;
n ij = migration rate of susceptible individuals from jth region to ith region;
m ij = migration rate of infected individuals from jth region to ith region.
With these assumptions, we have n ii ≡ m ii ≡ 0. Then, for i = 1,..., p,
If the system (2) is well-posed, we can find a bound for the population of susceptibles in each patch. We can then extend the same idea to obtain the linear system of I-equations only. (See Appendix for details.)
From Theorem 1 in the Appendix, we have
Thus, the total population of susceptible individuals in the ith region is limited. So we can write
π i = β × (total population in the ith region without infection)
d i = μ i + γ i ,
where d i represents the total death rate in the ith region. Thus, for i = 1,..., p, the equations of the I-only model are
This can be rewritten as
where I(t) = (I1(t), I2(t),..., I p (t))Tand
A two-region example
The simplest nontrivial version of model (2) is the case when p = 2, which can be presented as follows:
The disease-free equilibrium for model (4) is with
By using the next-generation matrix method , we can calculate the basic reproductive ratio for model (4),
This is the eradication threshold condition for this simple case p = 2. The case for p = 2 can be further simplified under different assumption of the travel (or migration) between the two regions. Here we list three possibilities.
The two regions are isolated from each other. Then, for each of the two regions, we have, for i = 1, 2,
which has the eradication condition
We only allow the susceptible individuals to travel (or migrate) between the two regions. Then we have
The disease-free equilibrium is again . However, the basic reproductive ratio (or the eradication condition for the disease) is given by
The travel (or migration) between the two regions is unidirectional. Without loss of generality, we can assume that individuals only travel from Region 1 to Region 2. Thus, we have
The disease-free equilibrium for model (10) is , with
Using the next-generation matrix method again, we obtain the basic reproductive ratio
which is also the eradication condition for the disease in this case.
Comparing all of the above cases when allowing the travel (or migration) between the two regions, we know that all of the eradication conditions in (5), (9) and (11) are not as simple as (7), the case when the two regions are disconnected. All of them depend on the travel rate between the two regions to some extent. For example, for Case 2, denoting
calculation shows that
From this, one can derive conditions for the travel rate of susceptible individuals which can either help a disease which is otherwise dying out to persist locally or cause an otherwise partially persistent disease to persist globally in the two regions (see Appendix for the details). It follows that the disease must be reduced within each region and the effects of travel between regions must be accounted for. This explains why HIV eradication must be considered as a world problem, not just a problem for individual countries, or continents, to tackle independently.
Finding the eradication threshold condition T0
Since system (3) is linear, the stability of the disease-free equilibrium is determined by the sign of the real parts of the eigenvalues of the matrix K. The only equilibrium for the I-only model (3) is the trivial, disease-free equilibrium (0,...,0). The linear stability of the trivial equilibrium for (3) is determined by s(K), the stability modulus of K, which is defined as the maximal real part of all eigenvalues of the matrix K. Thus, we can define the overall eradication threshold as
The trivial equilibrium is globally asymptotically stable if T0 < 1 (for this linear system (3), the local stability is equivalent to global stability) and unstable if T0 > 1.
Noting that the right-hand side of I-equations in (2) is never bigger than the right-hand side of (3) and then by the standard comparison theorem (see [63, 64]), one can show that T0 is also the condition to guarantee the stability of the disease-free equilibrium of system (2). This implies that HIV can be eradicated under the condition that T0 < 1, but it can not be eradicated under the condition that T0 > 1.
A continent-level example
As a more realistic example, we divide the world into six continents: Africa, Asia, Europe, North America, Oceania and South America. Infection occurs within each continent and travel/immigration of infectives occurs between continents. See Figure 3.
In this case, model (2) is explicitly
Parameters for our model consist of death, birth, immigration and emigration rates at the continental level, as well as estimates for the number of sexual partners and infection rates. Immediate data describing these rates does not exist; however, we can determine many of these parameters at the continental level using country-specific data. The birth and death rates have been derived from crude birth and death rates per 1000 people. The continental birth and death rates are determined using weighted averages by country population (see Table 1). Data was taken from the CIA world factbook 2008 .
Net immigration data is also reported in the CIA world factbook; however, this does not enable us to determine separate immigration and emigration rates. We can, however, determine continent-specific immigration and emigration rates from the number of foreign-born individuals by country. These are individuals who now reside in a particular country who stated that they were born in a foreign country on their census form. The Global Migrant Origin Database  gives the number of foreign-born individuals by country of origin and country of destination in the 2000 round of country censuses. This data has been summarised to determine the number of foreign born by continent; that is, individuals who now reside in a continent who were not born in that continent (see Table 2). Dividing Table 2 by the total number of foreign born gives the fraction of foreign-born individuals who migrated from continent A to continent B. If the number of immigrants each year is known, these fractions can be used to approximate the number of immigrants and emigrants by continent.
We use the information in Tables 1 and 2 to determine the eradication threshold T0. The off-diagonal entries of our matrix are the off-diagonal entries from Table 2. The diagonal entries of our matrix are calculated according to the formula
where Λ i is the influx of new susceptibles, μ i is the death rate (Table 1), β is the infection rate, γ is the AIDS-specific death rate, η ij is the number of individuals born in continent j residing on continent i per year (Table 2) and P i is the population of continent i (Table 1).
As an example, we choose β = 0.0002 people-1years-1 and γ = 0.1 years-1, reflecting the fact that the typical time course of the disease is ten years. We calculate the influx of new susceptibles according to the formula
where p i is the continent's birth rate and c is an infected individual's annual number of (uninfected) sexual partners. We set c = 20; however, the factor modifies an individual's annual number of sexual partners according to the growth (or decline) of the continent they are in. Thus, for Europe, the annual number of sexual partners is 17, whereas for Oceania the number is 44. These numbers are likely overestimates, but our intention is to be conservative in our approach.
If we now intervene and reduce β to 0.0002 x 2/5, then the overall threshold is . It follows that reducing the infection rate to 0.00008 will result in eradication.
Reducing β by three fifths is, of course, a mammoth task. Furthermore, this model implicitly assumes that such a reduction is instantaneous. To compare this to slowly reducing β over a period of decades, we assumed that β would be reduced to two fifths of its current value in twenty years' time, but that reductions would occur evenly throughout that time. However, we also included a population increase of 3% per year. The results are shown in Figure 4. In this case, the population "overruns" our goal. The overall T0 is not reduced below one, because the population in twenty years' time will be too high for a 3/5 reduction in the infection rate to lead to eradication.
We now explore how much these intervention measures might cost. Interventions include changes in education, condoms, drugs and travel restrictions. The costs that can be applied to these different interventions vary greatly from region to region. For example, the cost of HAART in high-income regions is much greater than the cost of the same drug treatment regime in low-income regions. However, this lower cost is still expensive compared to the per-capita national health expenditures. In this section, we determine the costs that we are going to apply to our model so that we can determine an optimal way to spend donor money to eradicate HIV. We list some costs found in the literature associated with education, HIV testing, drug therapy and circumcision in Table 3. All costs are given in US dollars (USD).
Along with the cost of delivering drugs to a patient, it is also important that we determine what benchmark of viral load and cell count we will use for initiation of therapy. It is obvious that the costs of delivering therapy and monitoring will increase if these are started at an earlier time in HIV infection. For example, it is common practice to initiate therapy in the developed world when CD4 counts are less than 350 cells/mm3. However, in the developing world, it is advised that treatment start when CD4 counts are less than 200 cell/mm3. In a study by Vijayaraghavan et al. , it was found that treating patients with HIV according to developed- versus developing-world guidelines is highly cost-effective and may result in substantial long-term savings. We choose to initiate individuals earlier in life for our cost analysis. This means that we will include a cost of therapy in the upper range of that listed in Table 3. However, we will only choose this cost from those for developing countries since the donor funds are applied mainly to the epidemic in the developing world.
Note that we assume that the costs of delivering therapy and health care to HIV infected individuals will stay relatively constant over a short period of time. Since the donor money will be spent over a span of a few years rather than decades, this assumption can be expected to give realistic results. If we were to determine the cost of such interventons over a long period of time, we must also include variables describing a decrease in the costs of therapy and testing as drug manufacturing becomes more efficient, the discovery and development of new drugs which will affect the cost of therapy, the increase in education and other medical expenses (since administrative costs and the payment of healthcare workers will increase over time), etc.
To effectively determine the cost of treatment, we also must include the possibility of the evolution of drug resistance . It is difficult to put a dollar figure on the evolution of drug resistance, since it may occur during different periods of infection in different patients. In our cost analysis, we choose costs of treatment and testing in the upper range of that reported (see Table 3), since it is assumed that the evolution of resistance will increase the cost of treatment and require more tests to monitor a patient's viral load and cell count.
Referring back to model (3), it is evident that condoms and education will reduce π i , while treatment will reduce π i , but will also decrease d i . Education will include condom awareness, but may include other factors, such as abstinence and monogamy . Travel/immigration restrictions will reduce m ij and m ji , but we do not expect that these will be realistic; past attempts to restrict movement based on travel or immigration have simply driven the epidemic underground .
Note that we assume a population growth of 3% per year. This corresponds to the maximum population growth rate reported by the CIA World Factbook. We use the maximum growth rate so that cost estimates give an upper bound of what may be realistic for donor spending. We also assume that HIV-positive cases will continue to increase at a rate of 3% per year, if no interventions are undertaken.
Our cost formula is thus
where C is the cumulative cost, n is the proportion of men who must ultimately receive condoms, m is the fraction of infected individuals who receive treatment, r is the timescale, ϵ is the (fixed) cost of distribution and education, and θ is the cost of treatment (including the cost of paid health workers, testing etc). Thus, in the first term of (17), there are 3 billion men, whose numbers increase by 3% every year, given 100 condoms each, at a cost of $0.02, with a fixed condom education campaign of $60 million. In the second term of (17), we assume that one tenth of infected individuals require treatment; thus, currently there are 3.3 million individuals requiring treatment, whose numbers increase by 3% every year, at an average cost of $2500 per patient per year.
If we provided condoms to 3/5 of all men over a period of twenty years and treated nobody, then the cumulative cost is
Thus, by steadily reducing the infection rate each year, the costs would blow up, way beyond the available funds and we would not even reach the eradication threshold. See Figure 5A.
If we treated everyone who requires treatment for twenty years, then the cumulative cost is
assuming no fixed costs of condom education.
If instead we aimed to treat 50% of people requiring treatment and provided condoms to 1/5 of all men within twenty years, then the cumulative cost is
Thus, treatment is significantly more expensive than prevention (as expected).
Conversely, suppose we aimed to reduce infection to two fifths of its current value and treat everybody infected over the next five years. Then the cumulative cost is
Thus, we could reach our eradication threshold within five years, provide treatment to everyone who needs it and stay within our $60 billion budget (plus interest). See Figure 5B.
We have provided a method for predicting treatment and prevention levels necessary to eradicate HIV/AIDS, based on population and immigration data at the continental level. This method is easily applicable to a finer scale, such as country-level data. The mathematical model is linear, which allows it to be easily generalised. We have also developed a formula to estimate the cost of some of these interventions.
It should be stressed that our model does not attempt to quantify the prevalence of the disease or the time course of infection. The model is a predictor of eradication only and thus should not be used in other contexts. However, it provides us with insight into the degree of intervention measures that are needed to eradicate the disease: an across-the-board reduction in the infection rate by a factor of three fifths or more, and antiretroviral treatment for everyone who requires it.
These are steep requirements that will cost a lot of money. However, we have demonstrated that spending only some of the money now and saving the rest for the future will result in a loss of ability to eradicate the disease, using existing estimates. We have clear, specific goals and they need to be acted upon immediately.
Our model also demonstrates the folly of focusing locally. Even if the basic reproductive ratio in all continents were reduced below one, the disease could still be sustained by travel or immigration (see Appendix). We note that travellers/immigrants do not present a high risk factor at present, but they may be a significant obstacle when the disease is close to eradication. This is why we must consider HIV/AIDS from a global perspective.
It should be noted that, in addition to our model's obvious shortcomings, our parameter estimates also have limitations. In particular, the data we use currently assumes that the disease-specific death rate is negligible, compared to the background death rate. While this may be approximately true in some western countries, it is clearly not the case in Africa . However, while the numbers used here are illustrative, our calculations of T0 are based on analytical thresholds. Thus, the sensitivity of the results to our parameter choices is easily measured: variations in the data in Tables 1 and 2, as well as formulas (15) and (16) allows similar calculations to be performed at other scales than continental (e.g. country-level).
Another assumption that should be noted is that we assume that travellers/immigrants receive no intervention help, such as education or treatment. This may be likely if foreign nationals receive no health care and have limited access to education, but will not be universally true. However, in each case we have overestimated the severity of the disease, in order to provide a strict criterion: if our overall T0 - calculated from the linear model with the overestimated parameters - is less than one, then we will have eradication. Our immigration/emigration data did not include tourism, especially sex tourism  which potentially adds a huge effect to the impact of travellers.
Although the cost of a condom is relatively cheap compared to other intervention measures, it is necessary to acknowledge the fact that a condom taken does not mean a condom used. Condom promotion requires a well- developed program infrastructure with excellent logistic competence of its staff which may be difficult to achieve. Many areas, especially rural and remote villages, are not accessible; for example, the healthcare system covers less than 50% of the population in many sub-Saharan African countries. In addition, the effectiveness of condoms in HIV prevention is much less than 100%, perhaps mostly due to condom failure [24, 72]. In addition, there are civil-war zones, such as eastern Congo, where services cannot be provided.
This further limits the lack of accessibility to reach out universally to all parts of the populations. However, model (3) has the property that decreasing T0 will still be beneficial. Thus, even if intervention methods are imperfect, a large-scale attempt at eradication will still be enormously effective.
Similarly, although we have estimated treatment costs, we have used aggregate numbers that do not explicitly take into account many of the specific costs associated with treatment, such as healthcare worker training. It has been shown that 100% HAART coverage in sub-Saharan Africa could be achieved in 10 years, if the number of healthcare workers were doubled every year , but this becomes expensive quite quickly.
There are, of course, many other intervention methods that could be examined, such as circumcision, STD control etc. These methods could be applied together with the proposed interventions; e.g. condom promotion could be combined with syphilis treatment without adding lots of additional costs. As we have learned from family planning programs, the more methods are offered to a client and the greater the choice is, the more likely the client is using one. This is likely to apply for HIV/AIDS control as well; more methods are advantageous and more likely to lead to positive results. Recently, UNAIDS has also made this proposition to promote as many methods as are available and proven effective.
It is currently advised that countries with different epidemic patterns will require different national strategies for implementing effective HIV prevention programs. Countries classified as low-level, concentrated, generalised or hyperendemic should concentrate their efforts on different national strategies.
In our model, we have applied the same strategy to all regions and each region includes countries with one or more of these classifications. We note, however, that the education component could be specific to each continent, or to each country.
Furthermore, although model (3) predicts eradication if T0 < 1, the threshold also has the property that lower values are beneficial. Thus, even if actual eradication cannot be achieved (as may be the case, in reality), intervention methods such as those mentioned here that reduce T0 will nevertheless be beneficial.
Thus, for example, if some countries make significant inroads against the disease, it will still have a global benefit even if eradication does not occur.
We view this model as a starting point. It provides a criterion for eradication, but we believe that a more accurate model, one that reflected the time course of the disease in all regions, is achievable. However, such a model requires accurate input from all countries; ideally, it could be updated as more data becomes available, or as new intervention methods are developed. Disease eradication has occurred on a global scale before, as happened with smallpox and is currently underway for polio. While the mechanics of organising the international community in this way are beyond the scope of this paper, given the scale of the problem and the immediacy of the required solution, we believe that this model can provide us with an easy-to-grasp way of understanding what needs to be done and when to do it. For the former, our strategy involves harnessing all existing intervention techniques in significant strength. For the latter, the time is now.
Well-posedness of system (2)
Consider the system (2) with the following initial conditions:
The following theorem assures that the system is well-posed.
If (S10,⋯, Sp0, I10,⋯, Ip0) ∈ D L , then, for any L >L*, the set D L is positively invariant for solutions of (2) with (18).
We first show that the solutions of (2) with (18) are nonnegative. For this purpose, let us rewrite the system (2) as follows:
where S(t) = (S1(t),⋯, S p (t))T, I(t) = (I1(t),⋯, I p (t))T, Λ = (Λ1,⋯,Λ p )T,
with , for i = 1,..., p. Noting that the off-diagonal matrix elements of A(t) are nonnegative, we conclude that the entries of the matrix are all nonnegative. Indeed, let and rewrite A(t) as
where E denotes the p × p identity matrix. Then all entries of (t) are nonnegative, and hence so are entries of . We also have
Noting that the scalar matrix -G(t)E is commutative with any p × p matrix (hence with (t)), we have
implying that all entries of are nonnegative. Now, from (19), we have
Similarly, for any t ≥ 0, all entries of are nonnegative. Now, (20) leads to
implying I(t) ≥ 0 for t ≥ 0.
Finally, we show that S i (t) and I i (t) are bounded for t ≥ 0 and i = 1,..., p. Let N(t) = S1(t) + ⋯ + S p (t) + I1(t) + ⋯ + I p (t). By the nonnegativity of S i (t) and I i (t), for i = 1,..., p, we
This implies that N(t) is bounded with an upper bound ; hence, so are S i (t) and I i (t)
for t ≥ 0 and i = 1,..., p. This completes the proof. □
Comparison of Case 1 and 2 in Section 3.1
It is obvious from the formulas in (13) that and reflect the influence of travel of susceptible individuals between the two patches, and hence may be called the travel-modified basic reproductive ratios for Region 1 and Region 2, respectively. The following observations are direct consequences of (12)-(13).
Their verifications are straightforward and are thus omitted.
(A1) Assume < 1 and < 1. If n12 > 0 and n21 > 0 satisfy either
then > 1 and < 1. By symmetry, the conditions parallel to the above can lead to < 1 and > 1. Here, and in the following, we omit such parallel conditions.
(A2) Assume < 1 and > 1. If n12 > 0 and n21 > 0 satisfy either
then > 1 but > 1.
The biological meanings of (A1)-(A2) can be obtained in terms of the threshold condition meaning of and . For example, (A1) implies that travel of the susceptible individuals can help an otherwise dying out disease persist locally. Roughly speaking, a larger inflow of susceptible individuals favours the persistence of the disease in the region. (A2) shows that appropriate travel rates may cause an otherwise partially persistent disease to persist globally in both regions.
UNAIDS: Report on the global AIDS epidemic. 2008, [http://www.unaids.org/en/KnowledgeCentre/HIVData/GlobalReport/2008/2008_Global_report.asp]
Piot P: AIDS: from crisis management to sustained strategic response. Lancet. 2006, 368 (9534): 526-530. 10.1016/S0140-6736(06)69161-7.
Stover J, Bertozzi S, Gutierrez JP, Walker N, Stanecki KA, Greener R, Gouws E, Hankins C, Garnett GP, Salomon JA, Boerma JT, De Lay P, Ghys PD: The global impact of scaling up HIV/AIDS prevention programs in low- and middle- income countries. Science. 2006, 311 (5766): 1474-1476. 10.1126/science.1121176.
UNAIDS: AIDS in Africa: Three Scenarios to 2025. 2005, Geneva: Joint United Nations Programme on HIV/AIDS, [http://www.unaids.org/unaids_resources/images/AIDSScenarios/AIDS-scenarios-2025_report_en.pdf]
Barnett T, Whiteside A: AIDS in the Twenty-First Century: Disease and Globalization. 2006, Basingstoke, Hampshire: Palgrave Macmillan, [http://www.unaids.org/unaids_resources/images/AIDSScenarios/AIDS-scenarios-2025_report_en.pdf]2
Bloom DE, Rosenfield A: A moment in time: AIDS and business. AIDS Pat Care and STDs. 2000, 14 (9): 509-517. 10.1089/108729100438908.
Bloom DE, Mahal A: Does the AIDS epidemic threaten economic growth?. J Econometrics. 1997, 1-296.
Haaker M: The Macroeconomics of HIV/AIDS. 2004, Washington DC: International Monetary Fund
Tandon A: Macroeconomic Impact of HIV/AIDS in the Asian and Pacific Region. 2005, Manila: Asian Development Bank
Ewing T: AIDS programme faces donor fatigue. Nature. 1990, 346 (6285): 595-
Gates W, Gates M: Keynote Address, Opening Session, XVI International AIDS Conference, August 13 2006, Toronto. 2006, [http://www.kaisernetwork.org/health_cast/hcast_index.cfm?display=detail&hc=1791]
Okie S: Global health-the Gates-Buffett effect. N Engl J Med. 2006, 355 (11): 1084-1088. 10.1056/NEJMp068186.
Steinbrook R: Message from Toronto-deliver AIDS treatment and prevention. N Engl J Med. 2006, 355 (11): 1081-1084. 10.1056/NEJMp068167.
Shiffman J: Donor funding priorities for communicable disease control in the developing world. Health Policy Plan. 2006, 21 (6): 411-420. 10.1093/heapol/czl028.
McMahon R: Global Health and the Gates Foundation [Council on Foreign Relations.]. 2006, [http://www.cfr.org/publication/11004/global_health_and_the_gates_foundation.html]
Rowley JT, Anderson RM, Ng TW: Reducing the spread of HIV infection in sub-Saharan Africa: some demographic and economic implications. AIDS. 1990, 4: 47-56. 10.1097/00002030-199001000-00007.
Stover J: Influence of mathematical modeling of HIV and AIDS on policies and programs in the developing world. Sex Transm Dis. 2000, 27 (10): 572-578. 10.1097/00007435-200011000-00005.
Kates J: Financing the response to HIV/AIDS in low and middle income countries: funding for HIV/AIDS from the G7 and the European Commission. 2005, [http://www.sarpn.org.za/documents/d0001507/index.php]
Clark CF: AIDS and the Arrows of Pestilence. 1994, Golden, Colorado: Fulcrum Publishing
Bayer R, Fairchild AL: Changing the paradigm for HIV testing-the end of exceptionalism. N Engl J Med. 2006, 355 (7): 647-649. 10.1056/NEJMp068153.
Schwartländer B, Stover J, Walker N, Bollinger L, Gutierrez JP, McGreevey W, Opuni M, Forsythe S, Kumaranayake L, Watts C, Bertozzi S: Resource needs for HIV/AIDS. Science. 2001, 292 (5526): 2434-2436. 10.1126/science.1062876.
Montaner JS, Hogg R, Wood E, Kerr T, Tyndall M, Levy AR, Harrigan PR: The case for expanding access to highly active antiretroviral therapy to curb the growth of the HIV epidemic. Lancet. 2006, 368 (9534): 531-536. 10.1016/S0140-6736(06)69162-9.
Moghadas SM, Gumel AB, McLeod RG, Gordon R: Could condoms stop the AIDS epidemic?. J Theor Med. 2004, 5 (3-4): 171-181. 10.1080/10273660412331315147.
Gordon R: A critical review of the physics and statistics of condoms and their role in individual versus societal survival of the AIDS epidemic. J Sex Marital Ther. 1989, 15: 5-30.
Dugger CW: U.S. jobs shape condoms' role in foreign aid, New York Times (October 29). 2006, [http://www.nytimes.com/2006/10/29/world/29condoms.html]
UNAIDS: Joint United Nations Program on HIV/AIDS. 2004, [http://www.unaids.org/en]
Cohen DA, Farley TA: Social marketing of condoms is great, but we need more free condoms. Lancet. 2004, 364 (9428): 13-14. 10.1016/S0140-6736(04)16611-7.
Cohen DA, Farley TA, Mason K, Ridgeway G: The collectivity of sexual behaviour. Int J STD AIDS. 2006, 17 (3): 151-156. 10.1258/095646206775809213.
Trucco E: Studies in imitative behavior: A generalization of the Rashevsky model; Its mathematical properties. Bull Math Biol. 1954, 16 (4): 279-316.
United Nations General Assembly, Sixtieth session, Agenda item 45. Follow-up to the outcome of the twenty-sixth special session: implementation of the Declaration of Commitment on HIV/AIDS. Scaling up HIV prevention, treatment, care and support. 2006, [http://data.unaids.org/pub/InformationNote/2006/20060324_HLM_GA_A60737_en.pdf]
Kremer M, Snyder CM: Why is There No AIDS Vaccine? [Working Paper]. 2004, Cambridge, MA: National Bureau of Economic Research
Smith RJ, Bodine EN, Wilson DP, Blower SM: Evaluating the potential impact of vaginal microbicides to reduce the risk of acquiring HIV in female sex workers. AIDS. 2005, 19: 413-421. 10.1097/01.aids.0000161771.44276.92.
Hethcote HW, Van Ark JW: Modeling HIV Transmission and AIDS in the United States. Lect Notes in Biomath 95. 1992, Berlin-Heidelberg-New York: Springer-Verlag
Huang W, Cooke KL, Castillo-Chavez C: Stability and bifurcation for a multiple group model for the dynamics of HIV transmission. SIAM J Appl Math. 1992, 52: 835-854. 10.1137/0152047.
Mode CJ, Sleeman CK: Stochastic Processes in Epidemiology. HIV/AIDS, Other Infectious Diseases and Computers. 2000, Singapore, New Jersey: World Scientific
Korenromp EL, van Vliet C, Grosskurth H, Gavyole A, Ploeg van der CPB, Fransen L, Hayes RJ, Habbema JDF: Model-based evaluation of single-round mass treatment of sexually transmitted diseases for HIV control in a rural African population. AIDS. 2000, 14: 573-93. 10.1097/00002030-200003310-00013.
Blower SM, Aschenbach AN, Gershengorn HB, Kahn JO: Predicting the unpredictable: Transmission of drug-resistant HIV. Nat Med. 2001, 7: 1016-1020. 10.1038/nm0901-1016.
Boily MC, Bastos FI, Desai K, Mâsse B: Changes in the transmission dynamics of the HIV epidemic after the wide-scale use of antiretroviral therapy could explain increases in sexually transmitted infections-results from mathematical models. Sex Transm Dis. 2004, 31: 100-112. 10.1097/01.OLQ.0000112721.21285.A2.
Rapatski BL, Suppe F, Yorke JA: HIV epidemics driven by late disease stage transmission. J Acquir Immune Defic Syndr. 2005, 38: 241-253.
Vickerman P, Terris-Prestholt F, Delany S, Kumaranayake L, Rees H, Watts C: Are targeted HIV prevention activities still cost-effective in high prevalence settings? Results from an STI treatment intervention for sex workers in Hillbrow, South Africa. Sex Transm Dis. 2006, 33 (10 Suppl): S122-32.
Coffee M, Lurie MN, Garnett GP: Modelling the impact of migration on the HIV epidemic in South Africa. AIDS. 2007, 21: 343-350. 10.1097/QAD.0b013e328011dac9.
Nagelkerke NJD, Moses S, de Vlas SJ, Bailey RC: Modelling the public health impact of male circumcision for HIV prevention in high prevalence areas in Africa. BMC Infect Dis. 2007, 7: 16-10.1186/1471-2334-7-16.
Klein E, Laxminarayan R, Smith D, Gilligan C: Economic incentives and mathematical models of disease. Environment and Development Economics. 2007, 12: 707-732. 10.1017/S1355770X0700383X.
Garira Mukandavire: Effects of public health educational campaigns and the role of sex workers on the spread of HIV/AIDS among heterosexuals. Theor Pop Biol. 2007, 72 (3): 346-365. 10.1016/j.tpb.2007.07.002.
Xu Y, Xu GQ: Modeling of HIV/AIDS Transmission and Its Analysis. DCDIS Proceedings of the International Conference on Life System Modeling and Simulation. 2007, 488-492.
Gran JM, Wasmuth L, Amundsen EJ, Lindqvist BH, Aalen OO: Growth rates in epidemic models: Application to a model for HIV/AIDS progression. Statistics in Medicine. 2008, 27 (23): 4817-4834. 10.1002/sim.3219.
Lloyd AL, May RM: Spatial Heterogeneity in Epidemic Models. J Theor Biol. 1996, 179: 1-11. 10.1006/jtbi.1996.0042.
Arino J, Davis JR, Hartley D, Jordan R, Miller JM, van den Driessche P: A multi-species epidemic model with spatial dynamics. Math Med Biol. 2005, 22 (2): 129-42. 10.1093/imammb/dqi003.
Lloyd AL, Jansen VAA: Spatiotemporal dynamics of epidemics: synchrony in metapopulation models. Math Biosci. 2004, 188 (1-2): 1-16. 10.1016/j.mbs.2003.09.003.
Liebovitch LS, Schwartz IB: Migration induced epidemics: dynamics of flux-based multipatch models. Physics Letters A. 2004, 332: 256-267. 10.1016/j.physleta.2004.09.071.
Arino J, van den Driessche P: A multi-city epidemic model. Math Pop Stud. 2003, 10 (3): 175-193. 10.1080/08898480306720.
Fitzgibbon WE, Langlais M, Morgan JJ: A mathematical model for indirectly transmitted diseases. Math Biosci. 2007, 206 (2): 233-248. 10.1016/j.mbs.2005.07.005.
Hsieh YH, van den Driessche P, Wang L: Impact of Travel Between Patches for Spatial Spread of Disease. Bull Math Biol. 2007, 69: 1355-1375. 10.1007/s11538-006-9169-6.
Jacquez JA, Simon CP, Koopman J, Sattenspiel L, Perry T: Modeling and analyzing HIV transmission: the effect of contact patterns. Math Biosci. 1988, 92: 119-199. 10.1016/0025-5564(88)90031-4.
Castillo-Chavez CC, Cooke K, Huang W, Levin SA: On the role of long incubation periods in the dynamics of AIDS I: single population models. J Math Biol. 1989, 27: 373-398. 10.1007/BF00290636.
Busenberg S, Cooke KL, Thieme HR: Demographic change and persistence of HIV/AIDS in a heterosexual population. SIAM J Appl Math. 1991, 51: 1030-1051. 10.1137/0151052.
Thieme HR, Castillo-Chavez CC: How May Infection-Age-Dependent Infectivity Affect the Dynamics of HIV/AIDS?. SIAM J App Math. 1993, 53 (5): 1447-1479. 10.1137/0153068.
Brauer F, van den Driessche P: Models for transmission of disease with immigration of infectives. Math Biosci. 2001, 171 (2): 143-154. 10.1016/S0025-5564(01)00057-8.
McCluskey CC: A model of HIV/AIDS with staged progression and amelioration. Math Biosci. 2003, 181: 1-16. 10.1016/S0025-5564(02)00149-9.
Gumel AB, McCluskey CC, van den Driessche P: Mathematical Study of a Staged-Progression HIV Model with Imperfect Vaccine. Bull Math Biol. 2006, 68: 2105-2128. 10.1007/s11538-006-9095-7.
Heffernan JM, Smith RJ, Wahl LM: Perspectives on the Basic Reproductive Ratio. J R Soc Interface. 2005, 2 (4): 281-293. 10.1098/rsif.2005.0042.
van den Driessche P, Watmough J: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002, 180: 29-48. 10.1016/S0025-5564(02)00108-6.
Smith HL: Monotone Dynamical Systems. An Introduction to the Theory of Competitive and Cooperative Systems. Math Surveys Monogr. 1995, 41: 1-174. Amer. Math. Soc., Providence, R.I.
Smith HL, Waltman P: The Theory of the Chemostat. 1995, Cambrige, UK: Cambrige University Press
Central Intelligence Agency World Factbook. [https://www.cia.gov/library/publications/the-world-factbook]
The Global Migrant Origin Database. [http://www.migrationdrc.org/research/typesofmigration/global_migrant_origin_database.html]
Vijayaraghavan A, Efrusy MB, Mazonson PD, Ebrahim O, Sanne IM, Santas CC: Cost-effectiveness of alternative strategies for initiating and monitoring highly active antiretroviral therapy in the developing world. J Acquir Imme Defic Syndr. 2007, 46: 91-100.
Gordon R, Tyson JE: Fatal sexually transmitted diseases (FSTDs), such as AIDS, select for the evolution of monogamy and provide a model for background extinction. J Biol Systems. 1993, 1: 1-26. 10.1142/S0218339093000021.
Burris S: Driving the epidemic underground? A new look at law and the social risk of HIV testing. AIDS Public Policy Journal. 1997, 12 (2): 66-78.
Lopez AD, Mathers CD: Measuring the global burden of disease and epidemiological transitions: 2002-2030. Ann Trop Med Parasitol. 2006, 100 (5-6): 481-499. 10.1179/136485906X97417.
Abdullah AS, Ebrahim SH, Fielding R, Morisky DE: Sexually transmitted infections in travelers: implications for prevention and control. Clin Infect Dis. 2004, 39: 533-538. 10.1086/422721.
White ND, Hill DM, Bodemeier S: Male condoms that break in use do so mostly by a "blunt puncture" mechanism. Contraception. 2008, 77: 360-365. 10.1016/j.contraception.2008.01.014.
Barnighausen T, Bloom DE, Humair S: Human Resources for Treating HIV/AIDS: Needs, Capacities, and Gaps. AIDS Pat Care and STDs. 2007, 21 (11): 799-812. 10.1089/apc.2007.0193.
Balakrishnan P, Solomon S, Kumarasamy N, Mayer KH: Low-cost monitoring of HIV infected individuals on highly active antiretroviral therapy (HAART) in developing countries. Indian J Med Res. 2005, 121: 345-355.
Yazdanpanah Y: Costs associated with combination therapy in HIV-infected patients. J Antimicrob Chemother. 2004, 53: 558-561. 10.1093/jac/dkh142.
Kumarasamy N, Solomon S, Chaguturu SK, Cecelia AJ, Vallabhaneni S, Flanigan TP, Mayer KH: The changing natural history of HIV disease: Before and after the introduction of generic antiretroviral therapy in southern India. Clin Infect Dis. 2005, 41: 1525-1528. 10.1086/497267.
Meng X, Anderson AF, Hou X, Wang Y, Sun L, Zhang X, Li Z, Qui B, Lang Y, Zhang L: A pilot project for the effective delivery of HAART in rural China. AIDS Pat Care and STDs. 2006, 20 (3): 213-219. 10.1089/apc.2006.20.213.
Dunne M: Antiretroviral drug development: the challenge of cost and access. AIDS. 2007, 21 (suppl 4): S73-S79. 10.1097/01.aids.0000279709.47298.8a.
Hounton SH, Akonde A, Zannou DM, Bashi J, Meda N, Newlands D: Costing universal access of highly active antiretroviral therapy in Benin. AIDS Care. 2008, 20 (5): 582-587. 10.1080/09540120701868303.
Beaston-Blaakman A, Shepard DS, Stone N, Shevitz H: Cost-effectiveness of clinical interventions for AIDS wasting. AIDS Care. 2007, 19 (8): 996-1001. 10.1080/09540120701335238.
Shevitz AH, Wilson IB, McDermott AY, Spiegelman D, Skinner SC, Antonsson K, Layne JE, Beaston-Blaakman A, Shepard DS, Gorbach SL: A comparison of the clinical and cost-effectiveness of 3 intervention strategies for AIDS wasting. J Acquir Immune Defic Syndr. 2005, 38 (4): 399-406. 10.1097/01.qai.0000152647.89008.2b.
Fieno JV: Costing adult male circumcision in high HIV prevalence, low circumcision rate countries. AIDS Care. 2008, 20 (5): 515-520. 10.1080/09540120701867149.
We thank Shoshana Magnet, Andrea Brandon, Anthony Wilson, Avril Newman, Matt Davison, Frederick Chen and Beni Sahai for technical discussions. We are also grateful to two anonymous reviewers for comments that greatly improved the manuscript. RJS is supported by an NSERC Discovery Grant, an Early Researcher Award and funding from MITACS. JMH is supported by an NSERC Discovery Grant and funding from MITACS. Handling editor for this manuscript was Connell McCluskey.
This article has been published as part of BMC Public Health Volume 9 Supplement 1, 2009: The OptAIDS project: towards global halting of HIV/AIDS. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2458/9?issue=S1.
The authors declare that they have no competing interests
RJS wrote the manuscript and performed numerical simulations. JL wrote the appendix. RG wrote the introduction. JMH collected the data. All authors read and approved the final manuscript.
About this article
Cite this article
Smith, R.J., Li, J., Gordon, R. et al. Can we spend our way out of the AIDS epidemic? A world halting AIDS model. BMC Public Health 9 (Suppl 1), S15 (2009). https://doi.org/10.1186/1471-2458-9-S1-S15