The decline in paediatric malaria admissions on the coast of Kenya
© Okiro et al. 2007
Received: 27 July 2007
Accepted: 15 November 2007
Published: 15 November 2007
Skip to main content
© Okiro et al. 2007
Received: 27 July 2007
Accepted: 15 November 2007
Published: 15 November 2007
There is only limited information on the health impact of expanded coverage of malaria control and preventative strategies in Africa.
Paediatric admission data were assembled over 8.25 years from three District Hospitals; Kilifi, Msambweni and Malindi, situated along the Kenyan Coast. Trends in monthly malaria admissions between January 1999 and March 2007 were analysed using several time-series models that adjusted for monthly non-malaria admission rates and the seasonality and trends in rainfall.
Since January 1999 paediatric malaria admissions have significantly declined at all hospitals. This trend was observed against a background of rising or constant non-malaria admissions and unaffected by long-term rainfall throughout the surveillance period. By March 2007 the estimated proportional decline in malaria cases was 63% in Kilifi, 53% in Kwale and 28% in Malindi. Time-series models strongly suggest that the observed decline in malaria admissions was a result of malaria-specific control efforts in the hospital catchment areas.
This study provides evidence of a changing disease burden on the Kenyan coast and that the most parsimonious explanation is an expansion in the coverage of interventions such as the use of insecticide-treated nets and the availability of anti-malarial medicines. While specific attribution to intervention coverage cannot be computed what is clear is that this area of Kenya is experiencing a malaria epidemiological transition.
Since the inception of the Roll Back Malaria (RBM) movement in 1996 , billions of health dollars have been committed by the international donor community to reduce the burden of malaria in Africa , estimated to be over a million deaths directly due to Plasmodium falciparum annually . However there are remarkably few documentations of changes in disease burden associated with increases in access and use of interventions funded by new international donor agency money.
In Kenya the expansion of coverage of both ITN and effective ACT therapy (artemether-lumefanthrine) has occurred very recently. Between 2004 and 2005 ITN coverage among children aged less than five years rose from 7% to 24% and by the end of 2006 had risen to 67% coverage . Despite delays in implementing the revised drug policy supporting the use of AL , over 85% of rural clinics had AL in stock between August and December 2006 . Paediatric admission time-series from three Kenyan hospitals were assembled to explore the impact of parallel changes in intervention coverage and malaria disease burden.
Paediatric ward in-patient registers at Malindi and Msambweni were identified for all months from January 1999 to March 2007. These were arranged serially to check whether these represented a continuous, uninterrupted series. Each admission entry in the registers was recorded on a separate tally sheet indicating the month of admission, whether a primary working diagnosis of malaria had been defined for the child, whether the admission diagnosis was not malaria and whether the child had survived admission. Separately, hospital death certificates were reviewed for the same period to identify any paediatric deaths that may have not been recorded in the admission ward books. All paediatric admissions were assumed to be aged between birth and 15 years. Individual register entries were not reconciled with patient notes and thus we have assumed that the admission diagnosis remained the clinical management diagnosis and it is used here as the diagnosis of analysis. The reliability of slide confirmed malaria diagnosis at admission was not validated from hospital or laboratory records as these procedures are variously performed in most Kenyan district hospitals  and results are rarely used to refine a diagnosis [8–10].
At Kilifi district hospital a sophisticated paediatric ward surveillance system has been in situ since 1989 to form the basis of a series of clinical studies on the pathogenesis of malaria, pneumonia, malnutrition and neonatal illness [11–14]. Clinical and laboratory coverage is provided 24 hours each day throughout the year. Demographic details and clinical histories are recorded on every admission and a finger prick blood sample taken for malaria parasitology and basic haematology. Following a clinical examination at admission further clinical and laboratory investigations are undertaken as indicated. Clinical and laboratory findings, clinical progress and response to therapy were reviewed at discharge to derive a primary diagnosis and recorded on a standard admission proforma later entered onto a centralized database. Data for the present study were reassembled for the period January 1999 through to March 2007 and summarized by month as a discharge diagnosis of malaria or non-malaria for all admissions aged between birth and 15 years of age.
Seasonal patterns of malaria hospitalization are related to monthly rainfall precipitations along the Kenyan coast  and relate principally to the population dynamics of the dominant vector species . Monthly rainfall data in decimal mm were obtained from meteorological offices located 3 km (Malindi), 2 km (Kilifi) and 25 km (Msambweni) from the respective district hospitals and were complete across the same time-period as the paediatric ward surveillance (Figure 1). The Kenyan coast is an area where sporogony in the vector population and hence malaria transmission, is not limited by ambient temperature and therefore temperature data were not included in the analysis.
Natural population growth over the 8.25 observation years is likely to have affected the size of population accessing these three hospitals. Precise rates of admission adjusted for monthly denominator size were not possible for the hospitals as the catchment for these services cannot be rigorously defined over time . Non-malaria admissions data were also collected, however, to help calibrate for population and other issues affecting hospital usage.
Analysis was undertaken using STATA version 9.2 (Statacorp 2003, College Station, USA). Hospital admission data obtained over the 99 months of surveillance between January 1999 and March 2007. Malaria admission case totals were assembled chronologically by admission month and non-malaria admissions similarly assessed for comparison purposes. They were examined with time series analysis that used malaria admission cases as the main outcome.
Smoothing techniques (moving averages) were used to filter short-term annual fluctuations and thus highlight longer-term trends present in the each of the data series. Malaria, non-malaria and rainfall, time series data from Kilifi, Kwale and Malindi were subject to 13-point moving average to aid visual interpretation of trends. Deviations in monthly rainfall were also examined by comparing monthly values of these parameters obtained during the study period with synoptic mean values computed from the 8.25 years of surveillance. The disparity between long term values and current values was referred to as an anomaly and had either positive or negative values.
A regression model expressing malaria cases as a linear combination of non-malaria cases and rainfall was used to test for trend. The non-malaria cases and rainfall allow trend statistics to be presented that are "aware" of other potential longitudinal influences on malaria admissions. A continuous variable indicating the time in months from the start of the observation period was included in the regression model. The coefficient of time in the model estimates the trend in the series (i.e. the month-to-month change in the number of admission cases). A P-value of < 0.05 was considered significant. There were strong seasonal fluctuations observed for malaria admissions, with a higher proportion of admissions recorded in the rainy season than in the dry season. To control for the confounding effect of seasonality on the trend, seasonal effects were included in the regression model by the use of indicator terms (dummy variables). The 11 dummy variables are constructed as time series with the value one for observations falling in a given month and zero when not. The final month was used as the baseline for comparison. There are thus 11 potential intercepts generated depending on the choice of baseline month for comparison; we elected to plot and discuss the intercept and trend that had the maximum correlation with the 13-point moving average. These analyses were conducted separately for each district.
A key assumption when using ordinary least squares regression is that the model residuals are independent  and this is often violated by longitudinal data . In the presence of autocorrelation standard errors of parameter estimates are underestimated resulting in an overestimation of significance. The Durbin-Watson statistic was hence used to test for serial autocorrelation . Serial correlation was found in all malaria admission data series (results not shown) indicating that correction was necessary. An analysis of autocorrelation and the sample partial autocorrelation function of deseasonalised data [19, 21] showed that a correction of this serial correlation with a lag of two months was optimal. Newey-West standard errors [21, 22] with a lag of two months were therefore used to correct for autocorrelation and potential heteroscedasticity.
Longitudinal data was obtained from all the hospitals for the period starting January 1999 to March 2007. During the 8.25 years of surveillance there were a total of 76,101 paediatric admissions: 41,715 admissions in Kilifi, 13,492 in Kwale and 20,894 admissions in Malindi. A total of 34% of the total admissions had a diagnosis of malaria; 32% (13,919) in Kilifi, 40% (5,361) in Kwale, and 36% (7,498) in Malindi.
Trends in paediatric admissions during the period January1999 to March 2007 at three sites on the Kenyan coast.
Intercept (95% confidence interval)
Trend (95% confidence interval)
Total across sites
259.57 (155.64, 363.49)
-1.67† (-2.29, -1.05)
80.12 (58.88, 101.35)
-0.43† (-0.65, -0.21)
119.15 (84.00, 154.31)
-0.33† (-0.59, -0.08)
Total across sites
236.72 (194.95, 278.50)
73.47 (50.16, 96.78)
-0.01‡ (-0.31, 0.29)
79.45 (57.55, 101.35)
0.91† (0.60, 1.23)
81.98 (38.27, 125.58)
0.10‡ (-0.65, 0.85)
48.96 (4.55, 93.37)
-0.04‡ (-0.89, 0.80)
84.32 (52.82, 115.83)
0.02‡ (-0.40, 0.44)
Trend of monthly malaria admissions during period the period January 1999 to March 2007 at three sites on the Kenyan coast.
95% confidence interval
Significant (P < 0.001) upward trends in non-malaria admission cases were shown in both the Kilifi and Malindi hospitals over the study period (Table 1). There was no statistically significant difference in Kwale. By March 2007, non-malaria admission cases in Kilifi and Malindi were estimated to have increased by 56% and over 100% respectively (Figure 2 – right panel) while remaining relatively consistent in Kwale. Again including rainfall as an additional covariate had little effect on the intercept and slope of the trends observed (results not shown).
Between January 1999 and March 2007 (8.25 years) admissions due to malaria have systematically declined at three sites along the Kenyan coast. These declines were examined in relation to possible changes in overall hospital utilization by a naturally growing paediatric population size and within and between year variations in rainfall. The fact that the significant decline in malaria admissions occurred against a significant rise in non-malaria admissions in Kilifi and Malindi and remained unchanged at Kwale during the surveillance period (Figure 1, right panel) suggests that the declining admissions due to malaria were specific to malaria. It could be argued that an increase in non-malaria cases might be attributed to a changing diagnostic pattern over the surveillance period with fewer misclassifications or more doctors accepting as true malaria test results. This cannot be ruled out at Malindi or Kwale but at Kilifi standardized diagnostic practices have applied to all paediatric admissions throughout the entire surveillance period and thus changes in the pattern of diagnosis seem an unlikely explanation for the trends observed.
The monthly incidence of malaria is coupled to seasonal rainfall patterns on the Kenyan coast. The addition of rainfall as a covariate to the model had no effect on the slope or intercept of seasonally corrected trend lines because rainfall patterns were shown to have remained relatively constant during the observation period with no significant increasing or decreasing trends observed in any of the study sites. Attributing the long-term reduction in malaria cases to a constant change in rainfall is implausible as there was no evidence of a substantial decline in rainfall except for the anomaly recorded in 2004/5, where rainfall was lower on average than that expected (Figure 3, right panel). This was however followed by a considerable substantial increase in amount of rainfall recorded in 2006. Despite this there was no associated increase in malaria admissions observed in 2006 or early 2007.
The use of three different hospital admission series increases the external validity of a single observation from one hospital. We consider the data from Kilifi as the diagnostic gold standard based on the use of microscopy and a discharge diagnosis supported by review of clinical notes. Using this as our point of reference, external validity of this result is provided by data from Msambweni (Kwale) and Malindi. It is reassuring then that a consistent declining pattern in malaria admission is maintained across all study sites lending weight to the strength of the findings. By 2006, annual malaria admissions had decreased by an average of 57% compared to 1999 across all sites.
Data on the serial use of prompt, effective antimalarial treatment in each of the three districts was not possible to assemble. However a semi-qualitative, temporal description of drug efficacy is possible in accordance with national first line recommendations for treatment between 1999 and 2007. Between January 1999 and early 2000 chloroquine was still the only available antimalarial in government clinics despite wide-spread resistance documented on the coast . Between mid-2000 and April 2006 sulphadoxine-pyrimethamine (SP) was the nationally recommended first line drug available in most government clinics and drug efficacy declined rapidly over this period in Kilifi and Kwale until most studies showed at least 25% failure rates when undertaken after 2002 [26, 27]. The policy was changed in April 2004 to artemether-lumefathrine (AL) but was not effectively implemented until December 2006 . The efficacy of AL was defined in Kilifi as in excess of 94% in 2002–2003 . Patient access to antimalarials within 48 hours of onset of symptoms has been documented in Kwale district in 2001 (15%) and 2006 (17%) (; unpublished data), both representing very low estimates of prompt access and we assume not dissimilar to anticipated access figures for Kilifi and Malindi districts during the same period. However, the widespread availability of even sub-efficacious SP over-the-counter from 2000 onwards in Kwale [30, 31] and Kilifi  may have had a suppressive effect on clinical disease risks operating similar to strategies prompting intermittent presumptive treatment in young children .
The use of hospital data provides a useful indicator of the long-term and short-term impact of scaling-up malaria interventions. It is not possible to definitely attribute carefully controlled and adjusted changes in malaria admission rates to expanded coverage of preventative and curative interventions; however these seem to be the most parsimonious explanations for the observations reported here along the Kenyan coast.
Footnote : The intercept refers to the seasonally adjusted level of malaria and non-malaria cases and rainfall at the start of the observation period. Models for malaria and non-malaria included a covariate for rainfall. The regressions were performed with Newey-West standard errors with a lag of two months for models with malaria and non-malaria case outcomes and a lag of one month for models with rainfall.
Footnote : The intercept refers to the seasonally adjusted level of malaria cases at the start of the observation period. Models included a covariate for non-malaria cases and rainfall. The regression was performed with Newey-West standard errors with a lag of two months for models with malaria and non-malaria case outcomes.
Footnote : Information on net deliveries was assembled from a variety of sources for each district. We assumed that ITN coverage was less than 5% between January 1999 and December 2002. This position is supported by reviews of net use undertaken in Kenya during this period . In Kilifi district a large scale trial of ITN was completed in 1993  and by 1997 there were few net replacements or net re-treatments . From 2002 Population Services International (PSI) had launched a retail sector, social marketing campaign for ITN distribution . Annual sales figures for the coast were obtained from PSI for the period January 2002 through to December 2004. Toward the end of 2004 PSI delivered heavily subsidized ITN through clinics and monthly sales by geo-positioned clinics in each of the three districts were available for the period November 2004 to March 2007. In September 2006 the Ministry of Health launched a large-scale free distribution campaign of free ITN to children under the age of five years . Net distribution volumes were recorded per geo-located distribution point within each district. Finally, small-scale community distribution projects over the surveillance by the district heath management team, NGO's and philanthropic organizations were recorded through interviews with DHMT members and district stakeholders and recorded as volumes of distribution and month of delivery. Cumulative monthly ITN distribution volumes per district were computed per capita using population size estimates derived from national census data and annual growth rates. This was implemented with projected population growth rate curves from national inter-censal district-specific annual rates of net population increase derived in 1989 and 1999 .
This study received financial support from The Wellcome Trust, UK (#058922), the British Department for International Development Kenya Programme (#KEN/2004/085), Ministry of Health, Government of Kenya (DFID agreement number #AF/06/416254) and The Kenyan Medical Research Institute. RWS is supported by the Wellcome Trust as Principal Research Fellow (#079081). We are indebted to the medical records officers in Kwale and Malindi: Katana Ndegwa, Joseph Rimba, Priscilla Folonzi and Ambrose Munga for their help with assembling data from these hospitals. The invaluable support provided by Drs James Nokes, Jay Berkley, Antony Scott, Charles Newton and Evasius Bauni at the KEMRI programme in Kilifi for maintaining the long-term hospital surveillance. We also wish to thank Dr Willis Akhwale of the Ministry for Health, Veronica Musembi of PSI-Kenya, and Samuel Waweru of the Metrological department for their assistance in obtaining various data sets. We are grateful to Wendy Omeara for her comments on earlier versions of the manuscript. This paper is published with the permission of the Director KEMRI.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.