AJTMH Transactions of the Royal Society of Tropical Medicine and Hygiene
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am. J. Trop. Med. Hyg., 75(2 suppl), 2006, pp. 90-103
Copyright © 2006 by The American Society of Tropical Medicine and Hygiene

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by TEDIOSI, F.
Right arrow Articles by TANNER, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by TEDIOSI, F.
Right arrow Articles by TANNER, M.

AN APPROACH TO MODEL THE COSTS AND EFFECTS OF CASE MANAGEMENT OF PLASMODIUM FALCIPARUM MALARIA IN SUB-SAHARAN AFRICA

FABRIZIO TEDIOSI, NICOLAS MAIRE, THOMAS SMITH*, GUY HUTTON, JÜRG UTZINGER, AMANDA ROSS, AND MARCEL TANNER
Swiss Tropical Institute, Basel, Switzerland


ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
An important shortcoming of existing methods for estimating the cost-effectiveness of malaria control interventions is that the incidence of illness and transmission dynamics are assumed to be independent of the case management system. We have developed a model for case management and integrated it into a stochastic simulation of Plasmodium falciparum malaria dynamics. This allows us to predict the incidence of clinical episodes and of mortality while incorporating effects of case management on persistence of parasites and transmission. We make predictions for a range of different transmission intensities in sub-Saharan Africa and simulate a range of case management scenarios with different coverage rates. The model predicts that high treatment rates have a proportionately greater epidemiologic impact at low transmission levels. Further development is needed for models for health-seeking behavior and referral patterns. The current model is a first step towards useful predictions of the epidemiologic and economic consequences of introducing and/or scaling-up of malaria control interventions.


INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Models for the cost-effectiveness of different disease control strategies generally, and malaria control interventions in particular, usually assume the incidence of illness and transmission dynamics to be independent of the prevailing case management system.1,2 Decision trees typically start with the observation of an illness episode and consider how the episode is managed subsequently, but do not consider whether management itself affects the incidence of the disease. With infectious diseases such as malaria, this approach ignores the feedback that may arise from an effective case management system, which in turn reduces the frequency of infection, and thus impacts transmission dynamics.

Models for nosocomial infections and those recently developed to simulate the severe acute respiratory syndrome epidemic3,4 provide examples of approaches in which the operation of the health system interacts in a dynamic fashion with the biology of the infecting organism. In these models, prompt treatment of infections reduces epidemic spread, and thus results in both reductions in the subsequent burden of disease and in the requirements for treatment of secondary cases.

When the infection has a long time course, which may encompass several illness episodes, treatment may also reduce the subsequent burden of disease independently of its effect on transmission. In malaria, untreated Plasmodium falciparum infections can persist for many months, during which clinical attacks recur at irregular intervals.5 One of the current mainstays of malaria control is access to early diagnosis and effective treatment.6 Prompt and effective treatment not only reduces the reservoir hosts who are infective to mosquitoes, but also prevents recurrences. Longitudinal studies of malaria in endemic populations frequently record declines in incidence over time. A major reason for this is likely to be that study participants receive more frequent treatment and this reduces the incidence of subsequent malaria fever attacks, irrespective of any effect on transmission. Conventional cost-effectiveness analysis of treatment does not consider these effects.

Malaria control interventions, such as source reduction by means of environmental management, increasing the coverage of insecticide-treated nets (ITNs), indoor residual spraying, and (potentially) the introduction of a malaria vaccine, modify the demands on the health system, and thus affect both immediate direct impact and longer-term indirect effects of case management. This applies even when the intervention, such as vaccination, does not directly modify case management. It follows that prediction of the impact of preventative and curative interventions against malaria must take into account these dynamic effects.

This report presents a first attempt to develop a dynamic model including case management of P. falciparum malaria in a typical setting of sub-Saharan Africa. It has been integrated into a model for the clinical epidemiology and natural history of P. falciparum malaria.7 We compare the outcomes of different case management regimens in settings of different transmission intensities.


MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Epidemiologic model. The epidemiologic model is a stochastic individual-based simulation of P. falciparum malaria in endemic settings that uses a five-day time step. The primary input is the pattern of the entomologic inoculation rate (EIR) in the absence of malaria control interventions, with separate values of the EIR specified for each of the 73 five-day periods during the year.7,8 For the present analyses, we simulate populations of 100,000 individuals, with an approximately stationary age distribution matching that of the demographic surveillance site in Kilombero, Tanzania in 1997–1999.9

For every individual in the simulated population, each discrete P. falciparum infection is characterized by a simulated duration and parasite density at each five-day time point.9 The host acquires immunity as a function of exposure and this in turn modifies the parasite density, and infectivity to mosquitoes10,11 at subsequent time-points. At each time point, a clinical event, either uncomplicated clinical malaria, severe malaria, or death from either malaria or other causes, may occur. Probabilities for occurrence of these events depend on the parasite density, recent exposure, and age-dependent co-morbidity. They have been determined by functions that have been fitted to field data across a wide range of transmission settings.1214 In addition, the prevalence of anemia (hemoglobin levels less than 8 g/dL) is assigned at the population rather than the individual level, as a function of simulated age and parasite prevalence.15

Clinical events. There are five different entry points into the case management tree: no event, uncomplicated malaria, severe malaria, indirect malaria death, and non-malaria death or out-migration. They are defined as follows.

No event, includes asymptomatic malaria infections. In this case, the simulated individual continues to the next time point, with the natural history of P. falciparum infections unmodified by the case management model.

Uncomplicated clinical malaria comprises P. falciparum infections that may be treated either at home or in peripheral health facilities. The model assumes that the risk of uncomplicated clinical malaria depends on whether the parasite density exceeds a critical threshold, which in turn is a function of past exposure.13 The case management implications of uncomplicated clinical malaria further depend on whether the host has recently been treated for malaria. Two possibilities were considered, as follows. The first is uncomplicated clinical malaria in the absence of recent treatment. This is defined by no treatment over the previous 30 days (6 time points). The decision tree pathways for this scenario are shown in Figure 1aGo. They include entry into the formal health care system and receiving the first-line drug, self-treatment at home with the recommended first-line drug, or absence of seeking of malaria treatment. The second is an uncomplicated clinical malaria episode that occurs despite recent treatment history. Figure 1bGo shows the decision tree pathways for this scenario. An uncomplicated malaria case that was treated in the past 30 days is assumed to either seek care or not. If care is sought, it is assumed to take place in the formal health care system, with treatment being based on the second-line drug. We do not consider the possibility that patients self-treat after drug failure because this would very likely involve ineffective retreatment with the first-line drug, with no epidemiologic consequences.


Figure 1
View larger version (27K):
[in this window]
[in a new window]
 
    FIGURE 1. Decision tree pathways. na = not applicable.

 
Severe malaria episodes are those clinical malaria episodes that are life-threatening if they are left untreated. We consider these as equivalent to those events that would have led to an admission diagnosis of severe malaria, had the patient presented to a health facility.12 The current model assumes that a severe malaria case can either be treated as an in-patient or not be treated at all. In the former case it is assumed that compliance is 100%. There are three possible clinical outcomes for treatment of a severe malaria episode, namely, death, recovery with neurologic sequelae, or full recovery (Figure 1cGo).

Indirect malaria death considers those deaths that would not have occurred in the absence of prior malaria exposure but which do not meet the criteria for severe malaria.12

Non-malaria death and out-migration correspond to events that occur independently of the parasitologic status of the host. These events are simulated to maintain the correct age-structure of the simulated population.

When more than one simulated clinical attack occurs within 30 days of another attack, these are counted as the same episode. Thus, there can be several treatments for one episode. The severity assigned to the episode is assigned to that of the most severe malaria attack within the 30-day period.

Each decision tree pathway predicts the outcome in terms of whether the parasites are cleared, and the clinical outcome (i.e., death, recovery with long-term sequelae or full recovery). The epidemiologic effects of the case management depend stochastically on the values of the joint probabilities of the clinical and parasitologic outcomes, conditional on the clinical event. These conditional probabilities are computed by calculating the probabilities for each branch of the decision tree pathways (Figure 1Go). For the model of uncomplicated malaria, the probabilities associated with each branch in the decision tree were obtained from the literature (Table 1Go).


View this table:
[in this window]
[in a new window]
 
TABLE 1
Model inputs used for efficacy and malaria treatment-seeking behavior*
 
For the severe malaria model, we used in-patient case fatality rates, and Qh(a), from a recent study in Tanzania,16 estimated corresponding community case fatality rates, Qc(a). Qh(a) varies with age a, taking values from 3% to 13%. Qc(a) takes values estimated previously using our model for severe malaria and mortality assuming the two risks to be related via a constant odds ratio, {varphi}1, taking a value of 2.09,12 i.e.,


Formula 1

We assume negligible drug-resistance to quinine, so that parasites are cleared in all hospitalized cases who survive. We assign a probability of sequelae, Rx, with a value independent of treatment (Table 1Go).

Disability adjusted life-years (DALYs). Years of life lived with disability are calculated using standard methods17 on the basis of the duration of disability, and respective disability weights (Table 2Go). These weights for different malaria-attributable disease conditions have been obtained from the Global Burden of Disease (GBD) study.18


View this table:
[in this window]
[in a new window]
 
TABLE 2
Disability weights and duration of disability used to calculate YLDs*
 
Years of life lost (YLLs) and DALYs (age-weighted) are calculated assuming age-specific life expectancies based on the life-table from Butajira, Ethiopia, with an average life expectancy of 46.6 years at birth.19 This life-table represents that of an east African setting, but is characterized by low malaria transmission. For example, it is very similar to that for Hai district, a high-altitude site in Tanzania.20

In a first step, YLLs and DALYs are presented with no discounting. Subsequently we compare the results with those obtained using a 3% discount rate, which is the one most commonly used in cost-effectiveness analyses.1,2

Malaria transmission intensity. The introduction of changes in case management (or of other interventions) leads to transient behavior, which may in principle modify the level of P. falciparum transmission. These effects on transmission are captured in the model by the effects on infectiousness of the human population resulting from clearing parasites. The simulation model predicts for time point t the proportion Formula 1m(t) of vectors that become infected at each feed on a human host.10,11 We adjust this to give {kappa}u(t) = 0.56 Formula 1m(t) to allow for the bias arising because Formula 1m(t) is estimated from artificial feed data.10 We record the value {kappa}u(0)(t) that {kappa}u(t) takes in the simulation of the reference scenario to which a change in the case management model has been applied, and compare this value to {kappa}u(1)(t), the prediction of {kappa}u(t) for the same time-point in the simulation with a change in case management. The effect of the change in case management on transmission is then modeled by a change in the EIR in adults at lv time units later (Emax(t + lv)), such that


Formula 2

where lv corresponds to the duration of the sporogonic cycle in the vector, and where Emax0(t + lv)/{kappa}u(0)(t) is the overall vectorial capacity. The consequences for the infection rates follow from details of the epidemiologic model.8

We considered four different intensities of transmission, each with the same seasonal pattern as that in Namawala, Tanzania21 (Table 3Go). For the reference scenario, we used an overall annual EIR of 21 infectious bites per year, which represents a typical level of transmission for a meso-endemic setting.22,23


View this table:
[in this window]
[in a new window]
 
TABLE 3
Scenarios modeled: health systems and transmission intensities*
 
Health systems. Uncomplicated malaria patients seeking formal health care in Tanzania are usually diagnosed during an out-patient visit either in a health center, a dispensary, or a hospital. A diagnostic test (normally light microscopy of finger prick blood smears) is performed on less than 10% of treated cases (National Malaria Control Program, 2004, unpublished data). Current Tanzanian national treatment guidelines recommend that an uncomplicated malaria episode be treated with sulfadoxine-pyrimethamine (SP) as the first-line drug. Amodiaquine serves as the second-line drug, and quinine is used for treatment of severe malaria.24 Our models assume that formal-sector treatment adheres to these practices (Tables 4Go and 5Go). We also assume that SP is used as the drug of choice for self-treatment, but the efficacy is assumed lower to account for lower quality of drugs purchased in the private sector. The levels of compliance and of drug resistance that we assume are given in Table 1Go.


View this table:
[in this window]
[in a new window]
 
TABLE 4
Sulfadoxine-pyrimethamine and amodiaquine doses and costs, by age and weight*
 

View this table:
[in this window]
[in a new window]
 
TABLE 5
Intravenous quinine doses and costs, by age and weight
 
We evaluate the effects of different case management in our model by varying either malaria treatment seeking-behavior or availability of treatment. Table 3Go summarizes the set of scenarios considered that include seven different sets of assumptions for the level of treatment.

No treatment. In this model, we assume no access to anti-malarial treatments.

Reference case management. In this model, we took the same probability of seeking treatment of uncomplicated malaria as in our recent simulation of a malaria vaccine trial,25 but assume 20% of treatments to be self-treatment and 80% to use formal care. We use a value of 48% for the probability of seeking treatment of severe malaria (Table 1Go).26,27 The treatment rates for uncomplicated episodes are low because the model for clinical episodes was fitted to very intensive surveillance data from Senegal, which included minor fevers that would be very unlikely to lead to treatment seeking.13 The treatment rates were estimated by triangulating the predictions of this model for clinical episodes with health system attendance data from Manhiça, Mozambique.25

Moderate coverage. In this model, we use recent data from the Tanzanian National Malaria Control Program of the Ministry of Health.28 This report states that 27% of children less than five years of age were treated within 24 hours in health facilities, 13% at home, and 2% at traditional healers. The remaining 58% received no treatment within 24 hours from the onset of disease.29 We use these percentages of 27% of children receiving formal care, and 13% self-treatment of uncomplicated malaria episodes (Table 1Go). We use the same coverage of formal care for severe malaria as in the reference model. In the context of our five-day time-step, we do not distinguish in our simulations whether treatment is within 24 hours or not.

Abuja target coverage. Sixty percent of uncomplicated episodes are treated with the appropriate drug. The simulated coverage of hospital treatment of severe episodes remains at 48%.

Complete coverage. We assume that 100% of uncomplicated clinical episodes are treated with formal-sector care. We also assume that 100% of severe malaria episodes are treated.

Effective treatment. We assume that 100% of clinical episodes are treated with formal-sector care and that there are no treatment failures. Moreover, we assume that all severe episodes receive in-patient treatment.

Anti-pyretic treatment. We assume the same anti-malarial treatment coverage as in the reference scenario, but in addition we assume that 50% of episodes are treated at home by paracetamol. The paracetamol is assumed to provide only symptomatic relief and thus to have no effect on the epidemiologic outcomes.

This defined the baseline status of the simulated populations. We ran simulations over a 90-year period at different transmission intensities under the assumptions of the reference case management scenario. The rather low level of treatment in this scenario is intended to approximate conditions prevailing in many areas studied in Africa of limited drug availability, drug resistance, and non-treatment or under-treatment of minor febrile attacks.

To explore the dynamic impact of different case management options, we then simulated the transient behavior over the next 5-, 10-, and 20-year periods for different case management scenarios, assuming the vectorial capacity for P. falciparum transmission to follow the same seasonal pattern as during the baseline period. We compared outcomes with those of scenarios in which the reference case management regimen continued.

Costing. Both marginal and average costs of health care were computed. The marginal cost of treatment is the additional financial or opportunity costs that is incurred when treating each additional case, but does not include the fixed cost of the infrastructure. The average costs include all those costs involved in delivering the intervention, including the use of spare capacity, and those health care resources diverted from other uses.30

Uncomplicated malaria. Direct costs of an uncomplicated malaria case seeking care at formal-sector facilities, Cdo, comprise the cost of an out-patient visit, the cost of drug treatment, and other costs incurred by the patients, i.e.,


Formula 3

where Ho is the patient (household) cost when visiting formal-sector outpatient facilities (excluding fees), Do is the cost of out-patient drug treatment, and Vo is the non-drug costs of an out-patient visit. Do is computed as


Formula 4

where Dod is the cost of drug per day and Ld is the number of days of therapy. The drug regimens and thus price depend on patient age and weight (Table 4Go), with the prices, which include distribution costs to districts, corresponding to those in the medical store department catalog31 of the Tanzanian Ministry of Health. W, the % additional cost of drug wastage, takes a value of 25% throughout.32

The non-drug cost of an outpatient visit is computed from published data on proportions of out-patients reporting at different levels of the health system, on the proportion, pt, of cases undergoing diagnostic tests, and on unit costs after exclusion of drug costs (Table 6Go). In the average analysis, the non-drug cost, Vao, is thus given by


View this table:
[in this window]
[in a new window]
 
TABLE 6
Health-seeking behavior and unit cost assumptions
 

Formula 5

In the case of marginal costs there is an adjustment for the proportion of recurrent non-fixed costs, i.e.,


Formula 6

The patient (household) costs per outpatient visit, Ho, comprise travel expenses, expenses related to medical supplies Hm, and non-medical supplies, Hn, such as the purchases of food and drinks or costs of spending the night away from home while seeking care33 (Table 6Go) so that


Formula 7

In self-treatment it is assumed that patients do not incur in any additional costs to purchase the drug because drugs are likely to be purchased from a private shop close to the patient’s home.

Severe malaria. The direct health care costs of a severe malaria case Cdi are given by


Formula 8

where Vi is the non-drug cost of in-patient care, Di is the cost of drug treatment, and Hi is the patient (household) cost when visiting formal-sector in-patient facilities. Di is computed by multiplying the costs by the duration for which they are incurred. During the first day of treatment, the drug dosage and consequently the costs are different, so overall Di is given by


Formula 9

where Di1 is the cost for the first day, Di2 is the cost per day thereafter, and Lt is the length of treatment (in days) (Table 6Go). The non-drug cost of in-patient care in the average analysis is given by


Formula 10

where Li(o) is the average length of stay, which varies depending on the outcome o, and Ni is the in-patient cost (Table 6Go). Correspondingly, in the marginal analysis the non drug cost is


Formula 11

The costs incurred by patients are the same as for an outpa-tient visit for the first day. For subsequent days of stay, we include only the costs of medical and non-medical supplies (Hm and Hn respectively) so that


Formula 12


RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Reference scenario simulation. Simulated patterns of age-prevalence and age-incidence for the reference scenario (Figure 2Go) are similar to those for typical meso-endemic settings in Africa to which the models were fitted.9,12,13 The direct cost per capita is stable over time. The predicted infectiousness of the host population, Formula 12u(t), fluctuates seasonally around a value of approximately 3% (Figure 3aGo). Since the reference scenario uses the same transmission pattern and health system to construct the baseline population as are applied during the follow-up, there is no trend over the 20-year simulation period in these epidemiologic variables or in the treatment costs.


Figure 2
View larger version (31K):
[in this window]
[in a new window]
 
    FIGURE 2. Predicted age-prevalence and age-incidence curves by transmission intensities The reference scenario is the transmission intensity observed at Namawala, Tanzania divided by a factor 16. Low transmission = Namawala transmission divided by 64; high transmission = one-fourth of Namawala transmission; very high transmission = transmission observed at Namawala. a, age-prevalence curve of patent parasit-emia; b, age-incidence curve of uncomplicated episodes; c, age-incidence curve of severe episodes; d, age-incidence curve of mortality.

 

Figure 3
View larger version (25K):
[in this window]
[in a new window]
 
    FIGURE 3. Infectivity of the human population. a, reference scenario; b, scenario assuming full treatment coverage.

 
Over the 20-year simulation period, the total number of undiscounted DALYs lost due to malaria in our population of 100,000 people is approximately 481,000, which corresponds to a rate of 0.24 DALYs per capita per year. Since most of these DALYs are due to mortality, the total number of YLLs is very close to that of DALYs (Table 7Go). If YLLs and DALYs are discounted at a standard rate of 3%, the total number of DALYs is considerably lower.


View this table:
[in this window]
[in a new window]
 
TABLE 7
Years of life lived with disabilities (YLLs),* disability adjusted life-years,* and direct costs{dagger}
 
The total undiscounted direct average costs to treat malaria episodes with the reference case management model amounts to US$ 485,793 over the 20-year simulation period. This corresponds to US$ 4.86 per capita and, on average, US$ 0.24 per capita per year (Table 7Go). Out-patient visits account for 32% of total direct costs, drug treatments of both uncomplicated and severe episodes for 7%, hospital admissions of severe episodes for 40%, and patient costs for 22%. The marginal cost, i.e., additional financial or opportunity costs that would be incurred when introducing a new control intervention, is approximately 58% of the average cost.

Effect of changing levels of access to case management. Comparison of the reference scenario with the extreme scenario with no treatment of malaria episodes (either uncomplicated or severe) showed noticeable epidemiologic effects of treatment despite the low attendance rates for uncomplicated episodes in the reference health system. In children less than 10 years of age, the no treatment scenario predicted higher prevalence of infection (Figure 4Go), a higher anemia prevalence (Figure 5Go), and a slight increase in the incidence of clinical episodes, with age patterns as shown in Figure 6a–cGo. There was only a small effect on the incidence of severe malaria (Figure 6d–fGo). Since the reference health system includes a relatively high treatment rate of severe episodes, the largest differences between the no treatment and reference scenarios were in the mortality rates, with a substantially higher mortality rate predicted with no treatment, especially in the second half of the first year of life (Figure 6g–iGo).


Figure 4
View larger version (20K):
[in this window]
[in a new window]
 
    FIGURE 4. Age-prevalence curves of parasitemia under different case management scenarios during a simulated 20-year follow-up period. {blacksquare} = no treatment; {triangleup} = reference; {blacktriangleup} = moderate coverage (first year of follow-up); * = complete coverage (first year of follow-up); {circ} = complete coverage (10th year of follow-up); • = complete coverage (20th year of follow-up).

 

Figure 5
View larger version (24K):
[in this window]
[in a new window]
 
    FIGURE 5. Age-prevalence curves of anemia (hemoglobin level < 8 g/dL) The reference scenario is the transmission intensity observed at Namawala, Tanzania divided by 16. Low transmission = Namawala transmission divided by 64; high transmission = one-fourth of Namawala transmission; very high transmission = transmission observed at Namawala.

 

Figure 6
View larger version (35K):
[in this window]
[in a new window]
 
    FIGURE 6. Age-incidence curves under different case management scenarios at different durations of simulation during a 20-year follow-up period. a–c, age-incidence of uncomplicated episodes; d–f, age-incidence of severe episodes; g–i, age-incidence of mortality.

 
The reference scenario is in a state of equilibrium throughout the simulation, although there is stochastic variation over time in the outputs. The no treatment scenario quickly reaches a new equilibrium. Therefore, only the average prevalence over 20 years is shown for this scenario (Figure 4Go), and the effects on incidence of morbidity and mortality do not change much during the follow-up period (Figure 6Go). The total number of DALYs lost over the simulated 20-year period was 13% higher with the no treatment regimen than with the reference. The effect on transmission, as measured by Formula 12u(t) is negligible.

The complete coverage scenario for treatment leads to a rapid decrease in transmission, as measured by Formula 12u(t) (Figure 3bGo). This stabilizes quickly at a value approximately 60% of that in the reference scenario, implying that treatment of all the clinical episodes (including minor episodes) can reduce the inoculation rate by about 40%. Complete coverage predicted very substantial decreases in prevalence of parasitemia (Figure 4Go), anemia (Figure 5Go), and incidence of uncomplicated episodes (Figure 6Go), but these outcomes, which reflect the dynamics of immunity, required an extended period to reach equilibrium. Although treatment of severe episodes in the complete coverage health system is no different from the reference, the effects on transmission and persistence of parasites result in substantial reductions in incidence of severe morbidity and mortality (Figure 6Go), leading to a total number of DALYs lost over 20 years of only approximately 45% of those in the reference scenario (Table 7Go).

The epidemiologic effects of high treatment coverage were concentrated in the youngest age groups, resulting in a substantial shift in the age of peak incidence of uncomplicated episodes (to older ages) and of mortality (to younger ages, in which a greater proportion of the mortality is contributed by indirect deaths). Since the changes in the age-prevalence and age-incidence curves caused by complete coverage were also time-dependent, transient effects can be seen throughout the 20-year follow-up period. As a result of the shifts in age-incidence, 10 years into the simulation an increase in incidence of clinical episodes above baseline levels is evident in individuals more than 10 years of age (Figure 6bGo). Since the shifts in the peak of the incidence curves to older age-groups accumulates over time, the benefit of the intensive treatment regimen decreases with time.

The distribution of direct costs is also changed as a function of treating all malaria episodes with a first-line drug. In the scenarios that we simulated, there would then be little need for second-line treatment; thus, many severe episodes could be prevented, which in turn reduces in-patient costs. Our model predicts that if all uncomplicated episodes were treated with the first-line drug, out-patient visit costs would account for 61% of total direct costs, drug treatments for 7%, inpatients admissions for 3%, and patient costs for 30%. The total direct costs to treat all malaria episodes would be approximately 7.4-times those of the reference scenario.

The moderate coverage scenario predicts effects on prevalence and on the clinical outcomes (Figures 4Go–6GoGo) more similar to those for complete coverage than to those for no treatment. This is despite the assumption in the moderate coverage health system of treatment rates for uncomplicated episodes much less than 50% and thus much closer to those for no treatment than those for complete coverage. This implies that within our models, there is a highly non-linear relationship between health outcomes and treatment coverage for uncomplicated malaria, with a very high marginal impact of increases in coverage when it starts from a low level. This conclusion is supported by the simulation of effective treatment, which gave very similar results to that of complete coverage for all the outcomes except mortality. Mortality rates and thus DALYs lost (Table 7Go) were reduced by the effective treatment health system to substantially below those with the complete coverage health system because all severe cases were assumed to be treated as in-patients.

The Abuja target coverage simulation is intended to simulate the impact of achievement of the coverage of 60% targeted in the Abuja Declaration on Roll Back Malaria in Africa, signed by African Heads of State and Government in April 2000. The YLLs, DALYs, and costs are intermediate between those predicted for the reference and those for full coverage.

The antipyretic treatment simulation is intended to indicate the sensitivity of the costs to allowance for the large number of people who self-treat with anti-pyretics, in addition to those who use anti-malarials. Unfortunately, it is difficult to obtain good estimates for actual levels of self-treatment in sub-Saharan Africa. The coverage of 50% for self-treating with anti-pyretics that we use is associated with an increase in drug costs of between 9% and 16% (Table 7Go), which indicates that although this is not likely to be a factor dominating costs, it could be of some importance because these are out-of-pocket expenses paid predominantly by poor patients. It would be important to include better estimates of drug choice and coverage for self-treatment.

Effect of transmission intensity. In higher transmission settings, simulated parasite prevalence was higher, and peaked at a younger age (Figure 2aGo). The incidence of uncomplicated malaria episodes also increased with transmission intensity in young children, but the reverse pattern was observed in older individuals (Figure 2bGo), matching the pattern to which the model was fitted.13 The incidence of severe episodes showed a similar pattern, but with a steeper decrease in incidence with age at high transmission, and consequently a crossing of the age-incidence curves at younger age (Figure 2cGo). The pattern for mortality was similar, with the mortality rate independent of transmission intensity at the age of approximately 3–4 years.

The average level of transmission to the vector, Formula 12u(t), was similar to that for the reference scenario for all values of EIR investigated, but the amplitude of the seasonal variation in Formula 12u(t) increased with the transmission intensity.

The total number of YLLs and DALYs lost over the simulated 20-year period are only slightly higher in the high transmission settings compared with areas of lower transmission intensity. The total direct costs are determined by the number of uncomplicated and severe episodes treated, which are higher than in the reference scenario in the low transmission setting and lower in the high transmission setting (Figure 7Go). However, these figures depend strongly on the model predictions for rates of severe malaria in adults, of which we are highly uncertain.12


Figure 7
View larger version (14K):
[in this window]
[in a new window]
 
    FIGURE 7. Direct costs in relation to transmission intensity. Upper line = undiscounted; lower line = discounted.

 
Changing access to case management in different malaria transmission intensity settings. The effects of an increased or decreased level of access to case management also vary according to the prevailing malaria transmission intensity (Figure 8Go). In a setting of moderate malaria transmission, treating all uncomplicated episodes over 20 years would lead to only a small difference in the total number of uncomplicated episodes, but reduces incidence of severe episodes by 49%, and the number of deaths and DALYs lost by 57% (Table 7Go).


Figure 8
View larger version (19K):
[in this window]
[in a new window]
 
    FIGURE 8. Effect of changing case management in different transmission settings. The status of the population at baseline was defined by simulating the reference health system for 90 years for each transmission intensity. The follow-up period then ran for 20 years with the chosen health system, and the results are the average for the 20 years. The incidence rates quickly settled to a new level. Because of the uncertainty in the rates of severe morbidity and mortality in adults, the predictions for only for children less than 10 years of age are shown. {blacksquare} = no treatment; {triangleup} = reference case management; x = full coverage.

 
In the low transmission setting, treating all uncomplicated cases over 20 years reduces incidence of uncomplicated episodes by 71%, severe episodes by 88%, and the number of deaths and DALYs lost by 83%. In the high transmission setting, complete treatment coverage would increase the total number of uncomplicated episodes over time, and would lead to a reduction of only 31% of severe episodes and DALYs lost and 34% of the number of deaths.

Treating everyone has a much greater effect on incidence at low transmission intensities. At high transmission, a high level of coverage always appears to be beneficial in terms of reducing incidence of severe episodes and of mortality, but may even lead to an increase in incidence of uncomplicated episodes. Within the model, this is because a very high treatment rate is associated with a reduction in exposure to asexual blood stage parasites and thus in acquired blood-stage immunity.

The economic implications of changing levels of access to case management also differ according to the malaria transmission intensity. Simulation over a 20-year period under the assumption of complete treatment coverage of uncomplicated malaria episodes would increase direct costs by a factor of almost 10 in a highly endemic setting, but the increase would be only 92% in low transmission settings (Table 7Go).


DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We present a first attempt to use a dynamic model for case-management of malaria in sub-Saharan Africa. For prediction of the effects of the case management, we considered a range of different transmission intensities characteristic for large parts of malaria-endemic Africa, and our simulations are thus likely to be of broad applicability. Lower EIRs than considered here (i.e., less than five infectious bites per person per year) are characteristic of highly urbanized settings, the African highlands, and areas located at the current distribution edges of P. falciparum transmission.34,35 Our model needs further development and validation to make meaningful predictions for such settings.

Our modeling approach expands the scope for predictions of the epidemiologic and economic consequences of malaria interventions as a direct function of the case management. In a first step, we have simulated different rates of treatment coverage, including the most extreme scenarios of either complete lack of treatment or full coverage. These two scenarios, together with a reference scenario largely constructed from real data obtained from Tanzania, were used for simulations up to 20 years. Costs were also built into our dynamic models, which will ultimately make it possible to predict the cost-effectiveness of the case management.

Our immediate purpose is to integrate effects of case management into our dynamic models of the clinical epidemiology and natural history of P. falciparum malaria in a typical setting of sub-Saharan Africa. This approach could readily be adapted to assess the costs of scaling up malaria treatment but this would entail more detailed analysis of the activities involved in changing treatment practices.

Our model can be used to make predictions of the effect of introducing a new malaria control intervention (e.g., malaria vaccine) or scaling-up of existing control measures (e.g., ITNs). The former motivated the development of our modeling approach. The current model is probably better at fulfilling this objective than in capturing the impact of changes in case management, including different levels of treatment and changes in national antimalarial drug policies. This is justified because the epidemiologic model was fitted mainly to cross-sectional data from various settings across sub-Saharan Africa. Since it is less able to capture longitudinal patterns within hosts, it does not claim to incorporate realistic patterns of treatment-seeking behavior or of referral patterns. In African settings where patients have limited resources, care-seeking patterns in general and malaria treatment-seeking behavior in particular are complex.3638

There is a need to build on simple models of the referral system32 so that we become more confident of the likely impacts of changes in national antimalarial drug policies. For example, parasite resistance to SP has reached critical levels in many parts of central and east Africa, including Tanzania,39,40 and this places a question mark against our longer-term predictions of cost-effectiveness, which assume that treatment with SP remains efficacious. In view of the public health, social, and economic significance of SP resistance, efforts are underway in Tanzania and elsewhere to change national policies towards artemisinin-based combination therapy (ACT). The Global Fund to Fight AIDS, Tuberculosis, and Malaria has recently approved a project to switch from SP (and amodiaquine) to ACT.41 Such a shift in drug policy is of considerable public health significance and will directly affect the case management of P. falciparum malaria. There is a need to adapt our model to field data from a range of different settings, and that will make it possible to explore the epidemiologic and economic consequences of the case management system under different scenarios. These should include simulation of the shift from SP to ACT conditional on various levels of SP resistance. All of these simulations need to take into account patterns of self-treatment, for which we currently have only weak data.

Regarding the epidemiologic model, we chose to define the seasonal pattern of P. falciparum transmission by using data from the village of Namawala in Tanzania. The high level of transmission measured there, even in comparison to other sites in Tanzania,22,42 made it possible to measure the intensity of transmission during the dry season. Multiplication of the Namawala rate by a constant therefore provides us with a reasonable estimate of the seasonal pattern, even for lower transmission areas where dry season transmission usually cannot not be measured.

The national malaria control program of Tanzania reports that 75% of people live in areas of stable malaria transmission, 17% in areas of unstable transmission (duration of transmission less than one month per year), and the remaining 8% in areas of unstable transmission (highly seasonal).29 However, historic and contemporary maps of malaria endemicity for Tanzania43,44 do not provide estimates of the inoculation rate with which to determine the distribution of EIR levels among the people within areas of stable transmission.

Our approach makes it possible to look at how such variations in transmission intensity might affect the impact of changes in the health system. However, our confidence in the present results is limited by uncertainties in our epidemiologic models (especially that for severe malaria and mortality, for which we had no data for older age groups12).

Increasing levels of treatment generally shift the age-prevalence and age-incidence curves so that the peaks are in older age-groups. We developed our parasitologic model mainly using archive data that predated the widespread use of anti-malarial chemotherapy. A delayed peak in the age-prevalence curve that we attribute to this effect was already apparent in the dataset from Navrongo, Ghana45 which was from a more recent period than the other studies.9 The effect of treatment on reducing acquired blood-stage immunity is very uncertain because asexual blood stage immunity is modeled as a function of both the number of distinct infections and of the cumulative parasite load, and we do not know what should be the relative contributions of these two different components of acquired immunity. The effects of cumulative parasite load are intended to simulate acquisition of immunity to antigenic variants that arise during the course of the infection.

We agree with other models for cost-effectiveness of malaria interventions in attributing most of the burden of disease to mortality rather than to disability associated with acute illness, sequelae, or anemia. Important methodologic issues requiring further investigation arise in the computation of the DALYs. In African populations with high infant mortality, there is generally an increase in life expectancy during the first few years of life, and this leads to perverse outcomes in the computation of DALYs if the effect of an intervention is partly to shift mortality to older ages in childhood. For example, if life expectancy at one year of age is 64 years, and life expectancy at five years of age is 72 years, then the number of life years gained by shifting age at death would be negative: 64–72 = –8. This effect is accentuated by discounting or age-weighting of the DALYs, but does not arise with the Japanese life tables used to compute DALYs in the GBD calculations used by the World Health Organization.17 We did not use these life tables because it is generally considered that local life tables should be used to compute DALYs in cost-effectiveness analyses,46 but there is a strong case for using death rates that exclude the health effect under investigation. For the present analyses, we used a life table from an east African site with low malaria incidence. In these life tables, where malaria plays only a small role, there is an increase in life expectancy over the first few years of life.

In conclusion, we have made a first attempt to develop a modeling framework to simulate the dynamic effects of the case management of P. falciparum malaria across a wide range of transmission intensities in sub-Saharan Africa. We discovered several deficiencies in our understanding of the relevant health systems. Our simulations of a range of scenarios indicate which of these uncertainties are most likely to be important for the prediction of cost-effectiveness of malaria interventions. Further development of our modeling approach offers more realistic evaluation of the epidemiologic and economic consequences of malaria interventions. This in turn will create a sound foundation for measuring the effects of introducing new antimalarial interventions (e.g., malaria vaccines), or scaling up those that are already known to be efficacious and cost-effective.


Received September 18, 2005. Accepted for publication March 27, 2006.

Acknowledgments: We thank Dan Anderegg for editorial assistance, and members of the Technical Advisory Group (Michael Alpers, Paul Coleman, David Evans, Brian Greenwood, Carol Levin, Kevin Marsh, F. Ellis McKenzie, Mark Miller, and Brian Sharp), the Project Management Team at the Program for Appropriate Technology in Health (PATH) Malaria Vaccine Initiative, and GlaxoSmithKline Biologicals S.A for assistance in this study.

Financial support: The mathematical modeling study was supported by the PATH Malaria Vaccine Initiative and GlaxoSmithKline Biologicals S.A.

Disclaimer: This publication and the contents hereof do not necessarily reflect the endorsement, opinion, or view points of the PATH Malaria Vaccine Initiative or of GlaxoSmithKline Biologicals S.A.

* Address correspondence to Thomas Smith, Swiss Tropical Institute, PO Box, CH-4002 Basel, Switzerland. E-mail: Thomas-A.Smith{at}unibas.ch Back

Authors’ address: Fabrizio Tediosi, Nicolas Maire, Thomas Smith, Guy Hutton, Jürg Utzinger, Amanda Ross, and Marcel Tanner, Swiss Tropical Institute, PO Box, CH-4002 Basel, Switzerland, Telephone: 41-61-284-8273, Fax: 41-61-284-8105, E-mails: fabrizio.tediosi{at}unibas.ch, nicolas.maire{at}unibas.ch, Thomas-A.Smith{at}unibas.ch, guy.hutton{at}unibas.ch, juerg.utzinger{at}unibas.ch, amanda.ross{at}unibas.ch, and marcel.tanner{at}unibas.ch.

Reprint requests: Thomas Smith, Swiss Tropical Institute, Postfach, CH-4002 Basel, Switzerland.


REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Gold MR, Gold SR, Weinstein MC, 1996. Cost-Effectiveness in Health and Medicine. Oxford, United Kingdom: Oxford University Press.
  2. Drummond M, O’Brien BJ, Stoddart GL, Torrance G, 1997. Methods for the Economic Evaluation of Health Care Programmes. Second edition. New York: Oxford University Press.
  3. Cooper B, Lipsitch M, 2004. The analysis of hospital infection data using hidden Markov models. Biostatistics 5: 223–237.[Abstract]
  4. Lipsitch M, Cohen T, Cooper B, Robins JM, Ma S, James L, Gopalakrishna G, Chew SK, Tan CC, Samore MH, Fisman D, Murray M, 2003. Transmission dynamics and control of severe acute respiratory syndrome. Science 300: 1966–1970.[Abstract/Free Full Text]
  5. Collins WE, Jeffery GM, 1999. A retrospective examination of the patterns of recrudescence in patients infected with Plasmodium falciparum. Am J Trop Med Hyg 61: 44–48.[Abstract]
  6. World Health Organization, 1993. A Global Strategy for Malaria Control. Geneva: World Health Organization.
  7. Smith T, Killeen G, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M, 2006. Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: overview. Am J Trop Med Hyg 75 (Suppl 2): 1–10.[Free Full Text]
  8. Smith T, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, Tanner M, 2006. Relationships between the entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg 75 (Suppl 2): 11–18.[Abstract/Free Full Text]
  9. Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L, 2006. A model for natural immunity to asexual blood stages of Plasmodium falciparum in malaria-endemic areas. Am J Trop Med Hyg 75 (Suppl 2): 19–31.[Abstract/Free Full Text]
  10. Killeen GF, Ross A, Smith T, 2006. Infectiousness of malaria-endemic human populations to vectors. Am J Trop Med Hyg 75 (Suppl 2): 38–45.[Abstract/Free Full Text]
  11. Ross A, Killeen GF, Smith T, 2006. Relationships of host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg 75 (Suppl 2): 32–37.[Abstract/Free Full Text]
  12. Ross A, Maire N, Molineaux L, Smith T, 2006. An epidemiologic model of severe morbidity and mortality caused by Plasmodium falciparum. Am J Trop Med Hyg 75 (Suppl 2): 63–73.[Abstract/Free Full Text]
  13. Smith T, Ross A, Maire N, Rogier C, Trape JF, Molineaux L, 2006. An epidemiologic model of the incidence of acute illness in Plasmodium falciparum malaria. Am J Trop Med Hyg 75 (Suppl 2): 56–62.[Abstract/Free Full Text]
  14. Ross A, Smith T, 2006. The effect of malaria transmission intensity on neonatal mortality in malaria-endemic areas. Am J Trop Med Hyg 75 (Suppl 2): 74–81.[Abstract/Free Full Text]
  15. Carneiro I, Smith T, Lusingu J, Malima R, Utzinger J, Drakeley C, 2006. Modeling the relationship between the population prevalence of Plasmodium falciparum malaria and anemia. Am J Trop Med Hyg 75 (Suppl 2): 82–89.[Abstract/Free Full Text]
  16. Reyburn H, Drakeley C, Carneiro I, Jones C, Cox J, Bruce J, Riley E, Greenwood B, Whitty C, 2004. The epidemiology of severe malaria due to Plasmodium falciparum at different transmission intensities in NE Tanzania. LSHTM Malaria Centre Report 2002–2003: 6–7.
  17. Murray CJL, Lopez AD, 1996. Estimating causes of death: new methods and global and regional applications for 1990. Murray CJL, Lopez AD, eds. The Global Burden of Disease. Geneva: World Health Organization.
  18. Murray CJL, Lopez AD, 1996. The Global Burden of Disease: A Comprehensive Assessment of Mortality and Disability from Diseases, Injuries, and Risk Factors in 1990 and Projected to 2020. Cambridge, MA: Harvard University Press.
  19. Network INDEPTH, 2004. Model Life Tables for Sub-Saharan Africa. Aldershot, United Kingdom: Ashgate.
  20. Network INDEPTH, 2002. Population, Health and Survival at INDEPTH Sites. Ottawa, Ontario, Canada: International Development Research Centre.
  21. Smith T, Charlwood JD, Kihonda J, Mwankusye S, Billingsley P, Meuwissen J, Lyimo E, Takken W, Teuscher T, Tanner M, 1993. Absence of seasonal variation in malaria parasitaemia in an area of intense seasonal transmission. Acta Trop 54: 55–72.[ISI][Medline]
  22. Hay SI, Rogers DJ, Toomer JF, Snow R, 2000. Annual Plasmodium falciparum entomological inoculation rates (EIR) across Africa: literature survey, Internet access and review. Trans R Soc Trop Med Hyg 94: 113–127.[ISI][Medline]
  23. Robert V, Macintyre K, Keating J, Trape JF, Duchemin JB, Warren M, Beier JC, 2003. Malaria transmission in urban sub-Saharan Africa. Am J Trop Med Hyg 68: 169–176.[Abstract/Free Full Text]
  24. Ministry of Health, 2000. National Guidelines for Malaria Diagnosis and Treatment. Malaria Control. Series 2000. Dar es Salaam, Tanzania: Ministry of Health.
  25. Maire N, Aponte J, Ross A, Thompson R, Alonso P, Utzinger J, Smith T, 2006. Modeling a field trial of the RTS,S/ASO2A malaria vaccine. Am J Trop Med Hyg 75 (Suppl 2): 104–110.[Abstract/Free Full Text]
  26. McCombie SC, 1996. Treatment seeking for malaria: a review of recent research. Soc Sci Med 43: 933–945.[ISI][Medline]
  27. McCombie SC, 2002. Self-treatment for malaria: the evidence and methodological issues. Health Policy Plann 17: 333–344.[Abstract/Free Full Text]
  28. National Malaria Control Program, 2003. National Malaria Medium Term Strategic Plan 2002–2007. Dar es Salaam, Tanzania: National Malaria Control Program.
  29. National Malaria Control Program, 2004. Monitoring Malaria Situation and Control Activities in Tanzania 2001–2003. Health Facility and Community Survey. Dar es Salaam, Tanzania: National Malaria Control Programme Monitoring and Evaluation Unit. National Malaria Control Program.
  30. Hutton G, Tediosi F, 2006. The costs of introducing a malaria vaccine through the expanded program on immunization in Tanzania. Am J Trop Med Hyg 75 (Suppl 2): 119–130.[Abstract/Free Full Text]
  31. MSD, 2004. Medical Store Department price catalogue 2004.
  32. Goodman C, Coleman P, Mills A, 2000. Economic Analysis of Malaria Control in Sub-Saharan Africa. Geneva: Global Forum for Health Research.
  33. Adam T, Kakundwa C, Manzi F, Schellenberg JA, Mgalula L, de Savigny D, Mbuya C, Wilzynska K and the MCE team in Tanzania, 2004. Analysis Report on the Costs of IMCI in Tanzania. Multi-Country Evaluation of the Integrated Management of Childhood Illness (IMCI). World Health Organization. Geneva: Department of Child and Adolescent Health and Development, World Health Organization.
  34. Keiser J, Utzinger J, De Castro MC, Smith TA, Tanner M, Singer B, 2004. Urbanization in sub-Saharan Africa and implication for malaria control. Am J Trop Med Hyg 71: 118–127.[Abstract/Free Full Text]
  35. Zhou G, Minakawa N, Githeko AK, Yan G, 2004. Association between climate variability and malaria epidemics in the East African highlands. Proc Natl Acad Sci U S A 101: 2375–2380.[Abstract/Free Full Text]
  36. Tanner M, Vlassoff C, 1998. Treatment-seeking behaviour for malaria: a typology based on endemicity and gender. Soc Sci Med 46: 523–532.[ISI][Medline]
  37. Muela SH, Ribera JM, Mushi AK, Tanner M, 2002. Medical syncretism with reference to malaria in a Tanzanian community. Soc Sci Med 55: 403–413.[ISI][Medline]
  38. de Savigny D, Mayombana C, Mwageni E, Masanja H, Minhaj A, Mkilindi Y, Mbuya C, Kasale H, Reid G, 2004. Care-seeking patterns for fatal malaria in Tanzania. Malar J 3: 27.[Medline]
  39. Mugittu K, Ndejembi M, Malisa A, Lemnge M, Premji Z, Mwita A, Nkya W, Kataraihya J, Abdulla S, Beck HP, Mshinda H, 2004. Therapeutic efficacy of sulfadoxine-pyrimethamine and prevalence of resistance markers in Tanzania prior to revision of malaria treatment policy: Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase mutations in monitoring in vivo resistance. Am J Trop Med Hyg 71: 696–702.[Abstract/Free Full Text]
  40. Talisuna AO, Bloland P, D’Alessandro U, 2004. History, dynamics, and public health importance of malaria parasite resistance. Clin Microbiol Rev 17: 235–254.[Abstract/Free Full Text]
  41. Global Fund to Fight AIDS, Tuberculosis, and Malaria, 2004. Available from http://www.theglobalfund.org/search/portfolio.aspx?countryID=TNZ.
  42. Drakeley C, Schellenberg D, Kihonda J, Sousa CA, Arez AP, Lopes D, Lines J, Mshinda H, Lengeler C, Armstrong SJ, Tanner M, Alonso P, 2003. An estimation of the entomological inoculation rate for Ifakara: a semi-urban area in a region of intense malaria transmission in Tanzania. Trop Med Int Health 8: 767–774.[ISI][Medline]
  43. Craig M, Snow R, Le Sueur D, 1999. A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today 15: 105–111.[ISI][Medline]
  44. Clyde DF, 1967. Malaria in Tanzania. Oxford, United Kingdom: Oxford University Press.
  45. Owusu-Agyei S, Smith T, Beck H-P, Amenga-Etego L, Felger I, 2002. Molecular epidemiology of Plasmodium falciparum infections among asymptomatic inhabitants of a holoendemic malarious area in northern Ghana. Trop Med Int Health 7: 421–428.[ISI][Medline]
  46. Fox-Rushby J, Hanson K, 2001. Calculating and presenting disability adjusted life years (DALYs) in cost-effectiveness analysis. Health Policy Plann 16: 326–331.[Abstract/Free Full Text]
  47. Brewster DR, Kwiatkowski D, White NJ, 1990. Neurological sequelae of celebral malaria in children. Lancet 336: 1039–1043.[ISI][Medline]
  48. Health Research for Action, 1999. Health Care Financing in Tanzania: Costing Study of Health Services. Final Report. Volume 1. Laarstraat, Belgium.
  49. Alonso-Gonzalez M, Menendez C, Font F, Kahigwa E, Kimario J, Mshinda H, Tanner M, Bosch-Capblanch X, Alonso PL, 2000. Cost-effectiveness of iron supplementation and malaria chemoprophylaxis in the prevention of anaemia and malaria among Tanzanian infants. Bull World Health Organ 78: 97–107.[ISI][Medline]