Young and vulnerable: Spatial-temporal trends and risk factors for infant mortality in rural South Africa (Agincourt), 1992-2007
© Sartorius et al; licensee BioMed Central Ltd. 2010
Received: 15 December 2009
Accepted: 26 October 2010
Published: 26 October 2010
Infant mortality is an important indicator of population health in a country. It is associated with several health determinants, such as maternal health, access to high-quality health care, socioeconomic conditions, and public health policy and practices.
A spatial-temporal analysis was performed to assess changes in infant mortality patterns between 1992-2007 and to identify factors associated with infant mortality risk in the Agincourt sub-district, rural northeast South Africa. Period, sex, refugee status, maternal and fertility-related factors, household mortality experience, distance to nearest primary health care facility, and socio-economic status were examined as possible risk factors. All-cause and cause-specific mortality maps were developed to identify high risk areas within the study site. The analysis was carried out by fitting Bayesian hierarchical geostatistical negative binomial autoregressive models using Markov chain Monte Carlo simulation. Simulation-based Bayesian kriging was used to produce maps of all-cause and cause-specific mortality risk.
Infant mortality increased significantly over the study period, largely due to the impact of the HIV epidemic. There was a high burden of neonatal mortality (especially perinatal) with several hot spots observed in close proximity to health facilities. Significant risk factors for all-cause infant mortality were mother's death in first year (most commonly due to HIV), death of previous sibling and increasing number of household deaths. Being born to a Mozambican mother posed a significant risk for infectious and parasitic deaths, particularly acute diarrhoea and malnutrition.
This study demonstrates the use of Bayesian geostatistical models in assessing risk factors and producing smooth maps of infant mortality risk in a health and socio-demographic surveillance system. Results showed marked geographical differences in mortality risk across a relatively small area. Prevention of vertical transmission of HIV and survival of mothers during the infants' first year in high prevalence villages needs to be urgently addressed, including expanded antenatal testing, prevention of mother-to-child transmission, and improved access to antiretroviral therapy. There is also need to assess and improve the capacity of district hospitals for emergency obstetric and newborn care. Persisting risk factors, including inadequate provision of clean water and sanitation, are yet to be fully addressed.
Infant mortality is an important health indicator of a population given its strong link to socio-economic status (SES), health service access and quality, and maternal health.
In the absence of vital events registration, health and socio-demographic surveillance (HDSS) data provide a valuable source for estimating mortality rates, trends and risk factors. HDSS sites implementing the verbal autopsy (VA) to determine probable cause of death are often the only means in most developing and many middle-income countries to observe cause-specific mortality of a population on a longitudinal basis and are a valuable tool for assessing trends in burden of disease [1, 2].
Diarrhoea, pneumonia, malnutrition and malaria are the leading causes of death among infants in low income countries [3, 4]. Birth asphyxia and neonatal sepsis are responsible for most neonatal deaths . These diseases, that can be largely prevented or effectively treated at relatively low cost, cause almost 95% of preventable infant and child deaths . HIV/AIDS has emerged as a major cause of death among infants in recent years, though in few countries outside of Africa .
In 1990, there was a 20-fold difference in the rate of infant deaths between sub-Saharan African and industrialized countries (180 versus 9 deaths per 1000 live births). In 2000, this difference had increased to 29-fold with mortality rates of 175 and 6 per 1000 children respectively . This is because many sub-Saharan African countries have seen reversals in child mortality trends in recent years due to HIV/AIDS. In 2007, approximately 420 000 children became infected with HIV , mostly through mother-to-child transmission (MTCT) [8, 9] in resource poor settings particularly sub-Saharan Africa. Kahn et al showed a doubling of child mortality due to HIV in a rural South African population (Agincourt sub-district) between 1992 and 2003 from 39/1000 person-years to 77/1000 . Garrib et al in 2006 found very high levels of infant mortality in another rural area of South Africa, 67.5 per 1000 person-years, with HIV/AIDS estimated as the single largest cause of death in the under-5 age-group (41% of deaths) . Thus interventions to reduce infant and child mortality are urgently required. A study in Zambia estimated that the cost per averted infection was approximately US$890 . According to a study in Barbados the lifetime cost of treating an HIV infected child is US$ 8,665 . This is much lower than estimates from the US where the cost for perinatally infected infants was USD 113,476 for 9 years of survival, US$ 151,849 for 15 years, and US$ 228,155 for 25 years . According to a study in the Ivory Coast, the mean cost of treatment was € (euros) 254 per child-year for infected children, €108 more than the mean cost of treatment for HIV-negative children born to HIV-positive mothers (a 74% increase in treatment costs) . Thus despite the costs associated with HIV/AIDS prevention among young children [16, 17], lifetime treatments costs of HIV infected infants are far higher; hence preventive measures need to be prioritized and targeted to those at high risk in poor, resource limited settings.
Effective interventions such as prevention of mother to child transmission (PMTCT) are available. A comprehensive approach to PMTCT can reduce transmission rates to below 2% [18–20]. Yet health care access and inequity remain widespread problems in economically disadvantaged areas . Mpumalanga Province in northeast South Africa was an important destination for refugees fleeing the civil war in Mozambique from 1983 onwards. A formal peace agreement was signed in 1992, yet despite voluntary repatriation programmes, by 2000 it was estimated that more than 200,000 former Mozambican refugees were still inhabitants in the province . A study by Hargreaves et al  demonstrated higher mortality rates among children from former Mozambican refugee households when compared to those from South African-headed households in the Agincourt sub-district. They concluded that lack of legal status and poorer SES of Mozambican refugees partly explains this disparity.
Inequalities in health outcomes or access to services and benefits can occur across space and time. In some situations this can reflect a compositional effect with variations merely reflecting the different groups that inhabit different locations . However, certain inequalities in child health outcomes are avoidable and unjust. These may reflect underlying inequities in the distribution of wealth, resources and social privilege in a given society, rather than an individual's choice or behaviour. To be fair, society must strive to achieve equal opportunities for all children regardless of parental status (education, SES) and geographical location. High-quality services for children that bridge the social divide are an important means of achieving equity goals. If South Africa is to achieve the Millennium Development Goals by 2015, including MDG 4 to reduce child mortality, then there is need to scale-up coverage rapidly with access to high quality health care and social support, particularly in the most poor and marginalised communities . When population-wide intervention programmes are too costly to implement, it becomes necessary to target such efforts to high risk areas where adverse health events are the most likely to occur . To address inequity in child survival, service planners need to understand the underlying socio-demographic profile and other factors contributing to high risk. Spatial-temporal mapping of high risk communities identifies those with greatest need, rather than those that are easiest to reach . This provides evidence on where to target interventions for greatest impact  and generates hypotheses on the determinants of increased risk.
With the development of Markov Chain Monte Carlo (MCMC) methods and Bayesian software such as WinBUGS, geostatistical spatial-temporal modelling has increasingly been applied in epidemiological research , especially with regards to malaria risk and transmission. Gemperli et al (2004) carried out a Bayesian spatial analysis of infant mortality in Bali that confirmed well-known risk factors and found a spatial pattern of infant mortality that showed a clear relationship with established foci of malaria transmission [30–32]. Despite growing applications of spatial methodology in malaria research, fewer studies have analysed spatial variations in population dynamics including all-cause and cause-specific mortality, with little or no work on longitudinal data collected in relatively small geographic areas covered by health and socio-demographic surveillance.
Many individual and household level factors have been identified as key determinants of infant and child mortality. Since objects in close proximity are often more alike, common exposures (measured or unmeasured) may influence mortality similarly in households of the same geographical area, introducing spatial correlation in mortality outcomes. Repeated measurements on individuals and households are also expected to be correlated in time. Standard statistical methods assume independence of outcome measures, for example mortality data. Ignoring this correlation introduces bias in the risk analysis as the standard error of the risk factors is underestimated, thereby overestimating significance. Bayesian geostatistical models relax the assumption of independence by, for example, incorporating random effects to measure spatial correlation as a function of distance between locations.
The aim of this study was to assess changes in infant mortality patterns in rural northeast South Africa over time, determine mortality risk factors and produce cause-specific mortality maps to identify high risk areas. These insights can provide guidance on the best allocation of limited resources to reduce infant mortality in this and similar areas of the country.
Study area and population
A verbal autopsy (VA) was conducted on every death to determine its probable cause . Interviews administered by trained lay fieldworkers were assessed independently by two physicians to determine probably cause-of-death. Where consensus could not be reached, a third independent medical assessment was made. The VA was first validated in the mid-1990s  and again in 2006 with particular reference to HIV/AIDS related mortality. International Classification of Diseases (ICD-10) was used to classify main or underlying, immediate and contributory causes of death. For this study, cause-specific analysis was limited to main causes from 1992-2006 as VA's had not yet been assessed for 2007.
Covariates included: infant demographic variables (gender, nationality); 5-year time periods; maternal factors (former refugee status, age at pregnancy, death in first year of child's life, education); fertility factors (parity, birth intervals, sibling death); household mortality experience, socio-economic status (SES) and food security; distance to health facility; antenatal clinic attendance; and household elevation (climatic proxy). Every two years since 2001, an asset survey was conducted in all households within the HDSS . Information on living conditions and assets, building materials of main dwelling, water and energy supply, ownership of modern appliances and livestock, and means of transport available were recoded (one being higher SES and zero lower status), summed to give an overall score for a household, and then used to construct wealth quintiles for SES ranked by increasing score from most to least poor.
The negative binomial is an alternative for the commonly used Poisson distribution, often regarded as the default distribution for integer count data. The Poisson assumes that expected mean equals its variance. The negative binomial differs from the Poisson distribution in that it allows for the variance to exceed the mean. Since the negative binomial distribution has one more parameter than the Poisson distribution, the second parameter is used to adjust the variance independently of the mean. Our data displayed evidence of being highly overdispersed and thus the negative binomial model was chosen.
A preliminary negative binomial regression analysis was carried out to assess the relationship between infant mortality and each covariate. Covariates significant at the 10% level (without substantial missing values) were then incorporated into the multivariate model.
The multivariate Bayesian negative binomial model was fitted in WinBUGS to examine the association between the significant covariates and all-cause infant mortality. Observation dates were used to calculate the person-days contributed by each infant (offset). Spatial random effects were used at a village level to take into account spatial correlation. Temporal random effects were also used at yearly intervals to account for temporal correlation. Village specific random effects were modelled via a multivariate Gaussian process (multivariate Gaussian distribution with covariance matrix expressed as a parametric function of distance between pairs of village centroid points) . Standard Bayesian autoregressive (AR) approaches, with priors for the AR(1) and AR(2) processes defined by Schotman  and Zeller  respectively, as well as a Poisson generalized autoregressive moving average (GARMA) approach , were tested to model the temporal random effects. Various order models for the AR and MA terms were assessed and the one that best fitted the data was used. MCMC simulation was employed to estimate the model parameters . Further details of the statistical modelling approach are given in the appendix.
Model assessment and validation
Deviance Information Criterion (DIC)  was used as the first step in comparison of model fit and the one giving the lowest DIC was chosen. Models were then also validated by fitting the models for 1992-2006 and predicting outcomes for all infants in 2007. Credibility intervals were constructed and the model providing the best predictions (along with low DIC) were used as the final model. The negative binomial models, particularly the AR(1) and AR(2) to model the temporal random effect, provided the lowest DIC (8618.07 and 8617.34 respectively) by some margin when compared to other approaches such as GARMA. In Bayesian statistics, a credible interval is a posterior probability interval which is used for interval estimation, in contrast to point estimation (confidence intervals). In other words, the credibility interval refers to the distribution of parameter values while a confidence interval pertains to estimates of a single value. In this study the negative binomial AR(2) predicted the outcome much better than the AR(1) model based on these Bayesian credibility intervals. Thus the AR(2) process was used in the final multivariate model.
A baseline model was used that included no covariates but a constant and site-specific (village centroid) random effect. All identifying features (village centroids, geographic boundaries) were removed, and the prediction area expanded irregularly (~740 km2) to double the normal size, in order to ensure confidentiality and avoid stigmatizing of villages. The HIV/TB map is not shown for this reason. Simulation-based Bayesian kriging  at prediction points (regular grid) within the site was used to produce maps of mortality risk for the whole HDSS site. Model estimates were exponentiated to represent incidence rate ratios (IRR).
Data extraction and management was done using Microsoft SQL Server 2005. The analysis was carried out in STATA 10.0, WinBUGS and R. The predictions of the fitted spatial models were mapped in MapInfo Professional 9.5.
Demographic profile of study sample
Between 1992 and 2007 31,804 infants were either born or migrated into the Agincourt HDSS. Of these, 26,000 (81.8%) were born within the site and half (50.4%) were female. Just under two-thirds were South African citizens (20,375; 64.2%) and a little over one-third were born to Mozambicans (11,356; 35.8%). There were 737 infant deaths (2.3%) giving an overall mortality rate of 24.7 per 1,000 person years; of these, 175 deaths were within the perinatal period and 202 within the neonatal period.
Cause of death (1992-2006)
The top causes of death among infants, as assessed by verbal autopsy, were HIV/TB (n = 116), acute diarrhoea or malnutrition (n = 91), acute respiratory infection (ARI) or pneumonia (n = 82) and septicemia (n = 20). In total 300 infant deaths were attributed to infectious and/or parasitic causes. During 1992-2006, 230 infants (33.6%) had an unknown cause of death.
Temporal trends by cause
All-cause univariate risk factor analysis for infant mortality in the Agincourt sub-district, 1992-2007
1 year continuous
5 year period
Mother refugee status
South African citizen
Breast feeding and birth weight
Breast fed (diarrhoea & malnutrition)
Increasing birth weight (kilograms)
Mother Mozambican in-migrated post 1992 α
Mother in same household
Mother not in household
Mother residency status
Permanent (> = 6 months in site)
Increasing number of months resident during the previous 12 months
Mother died in child's first year
Mother age at pregnancy
None or primary
Father died in child's first year
Father died before birth
Household morbidity and mortality
Cumulative number of household deaths in year of birth (continuous)
Number of household admissions in year of birth (continuous)
Birth parity (continuous)
Preceding birth interval
Post birth interval
Previous birth stillborn
Previous sibling died
Preceding interval sibling death
Mother attended antenatal clinic
Number of antenatal clinic visits
Socio-economic status of household
SES absolute score (quintiles)
Food security status of household
Predicted food shortage in coming year
Same amount of food
Distance to nearest health facility
Minimum distance to health facility (straight-line) from household
< 5 km
> = 5 km
Minimum distance to health facility (network) from village centroid
< 5 km
> = 5 km
Elevation (meters) - rainfall proxy
A significant increase in the number of years of maternal education as well as antenatal clinic attendance was observed over the study period (both p < 0.001). However, significant increases in the number of mothers dying in their infants' first year, number of maternal deaths (< = 42 days after infants date of birth), as well as other household deaths over time was also observed (all p < 0.001).
Almost a third (30.2%, n = 91) of mothers who died in the infants' first year died of HIV/TB; this was a significant risk factor for infant mortality (IRR = 164.7, p < 0.001). Approximately 44% of the mothers that died in the infant's first year died of unknown causes, many probably unclassified HIV-related deaths.
Household water supply consisting of raw natural water (river, pond or dam) was a risk factor (IRR = 16.50, p = 0.010) for deaths due to acute diarrhoea and malnutrition, although numbers were small. Mother being of Mozambican origin also proved to be a significant risk factor for infant death due to diarrhoea or malnutrition (IRR = 1.66, p = 0.019).
All-cause multivariate risk factor analysis for infant mortality in Agincourt sub-district, 1992-2007, using Bayesian geostatistical and temporal models
Later year of birth
Cumulative household deaths
Mother died in infant's first year
Previous child died
Previous birth stillborn
The results indicate a worsening of infant mortality: year of birth was significantly associated with infant mortality and risk of death increased over the study period. The increase was particularly from 1998 onwards, and can be largely attributed to the HIV epidemic and its impact on mortality in the study area [9, 44], both direct (vertical transmission of HIV) and indirect (death of a caregiver). Mother's death in infant's first year was a major risk factor in this study, as was higher numbers of cumulative household deaths. Results confirmed the importance of other known risk factors [23, 30]. The protective association between increasing maternal education and infant mortality has been previously described [30, 45] and is possibly a result of better health awareness and utilization of health facilities , longer birth intervals , and higher income which improves infants' health through ability to purchase goods and services . A significant association of higher household SES was however not observed in this study. This has been shown elsewhere and may be explained by that fact that unlike endogenous maternal and demographic factors that substantially influence an infant's risk of death, the effects of SES factors on mortality increase as the child gets older due to exogenous factors which parents have more control over .
We examined health service access with respect to primary health care generally and antenatal care specifically. Distance to nearest primary health care facility was not a risk factor in this study. Antenatal clinic attendance and number of ANC visits was significantly protective, with no difference between South Africans and former Mozambican refugees. These finding suggest that factors other than geographic access may be key to understanding the risks associated with health care utilisation. These could include quality of care, level of available care (primary versus secondary), cost and social barriers. In South Africa, primary health care for children under the age of six is free, as is antenatal care. However, financial costs associated with transport and opportunity costs associated with lengthy waiting time  are some of the barriers described in this setting [51, 52]. Twine el al showed that the poorest households were less likely to apply for social support grants than those in higher socioeconomic strata due to barriers such as distance from government offices, lack of official documentation and education of caregiver and household head .
A recent study in Kenya found that, despite significant spatial variations in child mortality, these were not correlated with distance to health facilities . They concluded that geographic access to curative services did not influence population-level mortality given the density of health facilities in Kenya. They also suggest that when distance access targets are met, further improvements in child survival can only be achieved through renewed investigation of the social, behavioural and quality-of-care factors that may obstruct access to health care services. Similarly in rural South Africa, there is urgent need to evaluate and assure a high level of health service quality, assess and strengthen referral patterns for emergency obstetric, infant and child health care, and identify other barriers to accessing these and other government services.
Mothers' physical presence or absence had a significant impact on infant mortality: mother a temporary migrant (largely work-related) proved significantly protective, while conversely increasing number of months per year spent resident by the mother in the rural site was a risk. Brockerhoff  describes how maternal rural-urban migration may affect children through three types of living arrangement: children may remain in the village as foster-children in the care of their fathers or other relatives; children may accompany or follow their mothers to towns or cities; and children born after migrant mothers settle in an urban area may remain there through the first few years of life. (Note that in this study, infants born to mothers in urban areas would not be captured onto the HDSS database unless they later migrated into the rural household). Bledsoe et al,  reviewing evidence from West Africa, suggest that while fostered children may be disadvantaged compared to biological children (in terms of access to health care and nutrition), they may still be better off than if they had accompanied their migrant mothers. By staying home, these children avoid exposure to infectious diseases during a vulnerable period of their life, have continued access to economic resources of a non-migrant father, and benefit from remittances received from the migrant mother  - as well as better health care, nutrition and enhanced maternal health knowledge . In our study, migrant mothers had significantly higher education and came from households with significantly higher SES which may explain the protective effect of mothers' migration. According to Collinson et al,  since 1997 there has been an increasing trend in the number of temporary female labour migrants in the Agincourt sub-district, a poor area with limited employment opportunities with resulting pressures to migrate and remit wages back to the rural household.
The spatial distribution showed marked geographical differences in all-cause mortality risk, indicating variation even within a relatively small area. The highest infant mortality risk was in those villages on the eastern border of the site. Much of this spatial distribution can be explained by the migration patterns of former Mozambican refugees (who constitute about a third of the Agincourt HDSS population) who entered South Africa via the Kruger National Park, a wild game conservation area situated between the eastern border of the site and Southern Mozambique. They remained a vulnerable group, poorer in more isolated villages with less infrastructure and generally further away from health facilities, with poor access to water and sanitation as well as labour markets and legal rights . However, our study indicates that the all-cause infant mortality risk pattern is not being driven by former refugee status alone, a finding supported by Hargreaves et al  who found no difference in mortality rates between South African and former Mozambican refugee infants between 1992 and 2000, despite significant differences in the 1-4 year age group. Multiple factors are driving the observed all-cause spatial risk pattern, including Mozambican origin of mother for certain infectious causes, maternal death in first year of infant's life, lower maternal education, poor quality of and limited access to neonatal care, poor antenatal clinic attendance, and increased vulnerability of households with a high mortality burden. These factors should be better elucidated and quantified in order to contribute meaningfully to policy and programmes.
With regard to the geographical distribution of infectious infant deaths (particularly HIV/TB) there was a distinct spatial pattern of mortality with an increasing gradient towards the east of the site where communities appear to have increased risk and suitable interventions need to be directed accordingly. One village in particular had a significantly higher mortality rate (all-cause and HIV) when compared to all other villages. Diarrhoea and malnutrition-related mortality was clustered in the south east of the site suggesting greater problems with clean water and sanitation, services that need to be assessed and addressed by local government. Breastfeeding had a protective effect on all-cause as well as diarrhoea and malnutrition-related infant mortality (Table 1). Breastfeeding protects infants through decreased exposure to contaminated water and food, optimal nutrition, and improved resistance to infection however there is risk of HIV transmission through breast milk. In South Africa, Ministry of Health policy on breastfeeding by HIV positive mothers has evolved in response to emerging research ; current recommendations are to breastfeed exclusively during the first 6 months with administration of anti-retrovirals to HIV positive mothers , especially those with low CD4 counts. Mothers or infants receiving highly active anti-retroviral therapy (HAART) prophylaxis should continue prophylaxis for one week after breastfeeding has ended . Infant mortality due to diarrhoea, malnutrition and their interaction is a complex problem in poor, HIV prevalent African settings. Addressing this requires a multifaceted approach including provision of clean water and sanitation, promoting infant nutrition, and strengthened primary care services for mothers and infants to reduce the risk of HIV transmission through breast milk .
Addressing health inequities in populations is a major challenge , and research that documents and quantifies inequities is needed to inform policies to close health gaps in the developing world. Evidence on reducing inequities within countries is growing; successful approaches include those that improve geographic access to health interventions in poor communities, subsidize health care and health inputs for the poor, and empower poorer communities . The results of our study indicate the need for interventions in villages to the east of the site, many of which have a large proportion of former refugees, to reduce the higher burden of infant deaths due to infectious and parasitic causes. HAART for HIV began in 2007 in this district and its impact cannot thus be captured during the time frame of this study. This research does, however, provide useful insight into spatial-temporal mortality patterns before HAART rollout and will allow post-rollout assessment of its impact on infant mortality. Such evaluation has the potential to identify areas needing improved access to treatment, specifically prevention of mother-to-child transmission and anti-retroviral therapy.
Of concern is the high number of neonatal deaths (particularly in the perinatal period), their gradual increase over the study period, and the highest risk area being in close proximity to a health facility. This suggests problems of service quality rather than geographic access, and highlights the need to assess and improve the capacity of sub-district health facilities for antenatal, emergency obstetric and newborn care; improve coverage of deliveries by skilled birth attendants; and advise mothers on appropriate care-seeking for sick babies. Part of the perinatal mortality burden observed may relate to maternal HIV since the same village experienced highest risk for neonatal and infant mortality. A meta-analysis by Brocklehurst et al  in 1998 found an association between maternal HIV infection and adverse perinatal outcomes, including low birth weight and pre-term delivery.
A limitation of the study is the potential to miss infant deaths, particularly neonatal deaths, which would underestimate the overall infant mortality burden. Infants that are born and then die during the 12 months between HDSS census update rounds may not be reported, particularly if the mother migrated out of the household; similarly, death among in-migrant infants who die before they are enumerated in the annual household census may be missed. However, infant death ascertainment has improved in the study site , and the proportion of infants who were in-migrants decreased significantly over time, reducing the bias towards the end of the study period. Determination of cause of death through verbal autopsy is more problematic for diseases that have less specific symptoms such as HIV/AIDS . The prevalence of HIV infection in a population and the resulting rate of HIV-associated co-morbidity and death due to malnutrition in children, for example, may affect the performance (such as specificity) of the tool. Thus it is likely that the HIV burden is underestimated due to the misclassification of deaths as AIDS-related conditions such as malnutrition or diarrhoea, or there being placed in the "unknown cause" category. The significant increase in number of infant deaths attributed to unknown causes since the late 1990s (Figure 2) is concurrent with the rise in HIV-related mortality in the area. Levels of stigma associated with HIV are high in South Africa, particularly prior to the introduction of HAART. The ability to make a diagnosis on VA depends, in large part, on the quality of information provided by the respondent. This may have been compromised in some cases in an effort to disguise HIV as a likely cause of death, partly explaining the increase in unknown causes.
By estimating the true spatial distribution of the infant mortality burden in rural northeast South Africa, this study has shown variation across a relatively small geographical area. The approach used Bayesian geostatistical models in order to assess risk factors, correctly estimate the standard errors (significance) of these risk factors and produce smooth maps of infant mortality risk from spatially correlated longitudinal mortality data in a health and socio-demographic surveillance system. Findings indicate the need for interventions targeted at villages with excess infant mortality risk due to both a direct and indirect impact of HIV. Essential interventions include improved prevention of mother-to-child transmission programmes, and antiretroviral therapy for HIV positive mothers to ensure their survival during their infants' critical first year(s) of life. From our study, it is clearly inadequate to consider maternal health separately from infant and neonatal health. This is consistent with other studies which showed that maternal health directly affects infants' health . Policy should thus have greater emphasis on interventions targeting the mother-infant pair. We also conclude that the non-random clustering of infant mortality due to diarrhoea and malnutrition in the south-east part of the site represents a breakdown in basic services (or, indeed, their absence); there is hence need to assess and improve water and sanitation in these villages. The high levels of perinatal mortality, in some instances in close proximity to health facilities, is of concern, indicating need to strengthen the capacity of sub-district facilities for emergency obstetric and newborn care. Recommendations from this study will have applications to other similar rural settings within South Africa and potentially beyond.
Appendix: Statistical Model
logit (pit) = β0 + βX it + φ it (multivariate spatial model)
logit (pit) = β0 + βX it + α t (multivariate temporal model)
logit (pit) = β0 + βX it + φ it + α t (multivariate spatial-temporal model)
logit (pit) = β0 + φ it (spatial kriging model) i.e. constant and spatial random effect with no covariates
where β0 is the incidence rate where all covariates are zero (i.e. the constant), Xit denotes the covariates, β is the vector of regression coefficients, φ it the village-specific random effect, μ i the individual level random effect and α t the temporal random effect. Following a Bayesian model specification, noninformative normal prior distributions were adopted for the regression coefficients β and an informative (based on estimates from Stata) and non-informative gamma prior distribution for the over-dispersion parameter r were adopted and tested [lower DIC dictating which was used]. We assume that φ it has a multivariate normal distribution, φ it ~ MVN (0,Σ), with variance-covariance matrix Σ. We also assume an isotropic stationary spatial process, where Σkl = σw 2 exp(-φd kl ), d kl is the Euclidean distance between villages k and l, σw 2 is the geographical variability known as the sill, φ is a smoothing parameter that controls the rate of correlation decay with increasing distance and measures the range of geographical dependency. A noninformative gamma prior was adopted for phi [φ], which is the smoothing parameter that controls the rate of correlation decay, as well as uniform prior with a distribution between φ min and φ max . Both approaches were tested and the approach providing the best fit was then used. The range is defined as the minimum distance at which spatial correlation between locations is below 5%. This distance can be calculated as 3/u meters. The second order year level autoregressive temporal random effect (αt), for t = 1 to 16 years, was modelled as a normal distribution with mean αmean [t = 3,..,16] = ρ0 + ρ*α[t-1] + ρ*α[t-2] and a noninformative gamma distribution for the variance parameter. The first two autoregressive terms were specified as αmean  <- ρ0 + l and αmean  <- ρ0 + ρ*α + l. Noninformative normal prior distributions were adopted for the ρ and l coefficients .
MCMC simulation was applied to fit the models. We ran a single chain sampler with a burn-in of 5000 iterations. Convergence was assessed by running the simulation until the Monte Carlo error for each parameter of interest was less than 5% of the sample standard deviation. The chains thereafter were sampled every single iteration until a sample size of 10,000 had been attained.
This work was supported by a PhD fellowship from the South African Centre for Epidemiological Modeling and Analysis (SACEMA), a National Research Foundation (NRF) Centre of Excellence. Additional funding was provided by a travel grant from the NRF Knowledge Interchange and Collaboration (KIC), Wits Faculty of Health Sciences Medical Research Endowment Fund (MREF) (Grant number: SARB000) and a research fellowship from the MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt) through The Wellcome Trust, UK. Benjamin Clark Provided support in preparing the data. The Agincourt health and socio-demographic surveillance system and geographic information system was funded by the Wellcome Trust, UK (Grant numbers: 058893/Z/99/A, 069683/Z/02/Z, 069683/Z/08/Z), the University of the Witwatersrand, the South African Medical Research Council, and the Andrew W. Mellon Foundation, USA.
- Hammer G, Somé F, Müller O, Kynast-Wolf G, Kouyaté B, Becher H: Pattern of cause-specific childhood mortality in a malaria endemic area of Burkina Faso. Malaria Journal. 2006, 5: 47-10.1186/1475-2875-5-47.View ArticlePubMedPubMed CentralGoogle Scholar
- Adjuik M, Smith T, Clark S, Todd J, Garrib A, Kinfu Y, Kahn K, Mola M, Ashraf A, Masanja H, Kubaje A, Sacarlal J, Alam N, Marra A, Gbangou A, Mwageni E, Binka F: Cause specific mortality rates in developing countries. 2006, Bulletin of the World Health Organization, 84: 181-188.Google Scholar
- Black R, Morris S, Bryce J: Where and why are 10 million children dying every year?. Lancet. 2003, 361: 2226-2234. 10.1016/S0140-6736(03)13779-8.View ArticlePubMedGoogle Scholar
- Ezzati M, Lopez A, Rodgers A, Murray C: Comparative quantification of health risks: global and regional burden of disease attributable to selected major risk factors. 2004, Geneva: World Health OrganizationGoogle Scholar
- Walker N, Schwartländer B, Bryce J: Meeting international goals in child survival and HIV/AIDS. Lancet. 2002, 360: 284-89. 10.1016/S0140-6736(02)09550-8.View ArticlePubMedGoogle Scholar
- UNICEF: Progress since the world summit for children: a statistical review. 2001, New York: UNICEF, [http://www.unicef.org/pubsgen/wethechildren-stats/sgreport_adapted_stats_eng.pdf]Google Scholar
- UNICEF, UNAIDS, WHO: Children and AIDS; Second Stocktaking Report. 2008, NY: UNICEF, 11-Google Scholar
- Dabis F, Leroy V, Castetbon K, Spira R, Newell M, Salamon R: Preventing mother-to-child transmission of HIV-1 in Africa in the year 2000. AIDS. 2000, 14: 1017-1026. 10.1097/00002030-200005260-00014.View ArticlePubMedGoogle Scholar
- Newell M: Prevention of mother-to-child transmission of HIV: challenges for the current decade. Bulletin of World Health Organization. 2001, 79: 1138-1144.Google Scholar
- Kahn K, Garenne M, Collinson M, Tollman S: Mortality trends in a new South Africa: Hard to make a fresh start. Scand J Public Health Suppl. 2007, 69: 26-34. 10.1080/14034950701355668.View ArticlePubMedPubMed CentralGoogle Scholar
- Garrib A, Jaffar S, Knight S, Bradshaw D, Bennish M: Rates and causes of child mortality in an area of high HIV prevalence in rural South Africa. Trop Med Int Health. 2006, 11 (12): 1841-8. 10.1111/j.1365-3156.2006.01738.x.View ArticlePubMedGoogle Scholar
- Stringer E, Sinkala M, Ng'uni C, McFarlane YA, Chipimo M, Myzece E, et al: In an urban African setting, the cost-effectiveness of PMTCT services improves as services expand. International Conference on AIDS [15th: 2004: Bangkok, Thailand]. Int Conf AIDS. 2004, 11-16. 15: abstract no. TuPeC4933.Google Scholar
- Adomakoh S, Kumar A, St John A: Reducing mother to child transmission of HIV-1 in Barbados. Cost-effectiveness of the PACTG 076 protocol in a middle income, low prevalence setting. International Conference on AIDS. Int Conf AIDS. 2002, 14: abstract no. LbPp2210Google Scholar
- Sansom SL, Anderson JE, Farnham PG, Dominguez K, Soorapanth S, Clark J, Sukalac T, Earp MJ, Bohannon B, Fowler MG, PSD Consortium: Updated estimates of healthcare utilization and costs among perinatally HIV-infected children. J Acquir Immune Defic Syndr. 2006, 41 (4): 521-6. 10.1097/01.qai.0000191286.70331.7b.View ArticlePubMedGoogle Scholar
- Giraudon I, Leroy V, Msellati P, Elenga N, Ramon R, Welffens-Ekra C, Dabis F: The costs of treating HIV-infected children in Abidjan, Ivory Coast, 1996-1997. Article in French Sante. 1999, 9 (5): 277-81.PubMedGoogle Scholar
- Stover J, Walker N, Garnett G, Salomon JA, Stanecki K, Ghys P, et al: Can we reverse the HIV/AIDS pandemic with an expanded response?. Lancet. 2002, 360: 73-77. 10.1016/S0140-6736(02)09339-X.View ArticlePubMedGoogle Scholar
- Schwartländer B, Stover J, Walker N, Bollinger L, Gutierrez JP, McGreevey W, et al: Resource needs for HIV/AIDS. Science. 2001, 292: 2434-36. 10.1126/science.1062876.View ArticlePubMedGoogle Scholar
- European Collaborative Study: Mother-to-child transmission of HIV infection in the era of highly active antiretroviral therapy. CID. 2005, 40: 458-465. 10.1086/427287.View ArticleGoogle Scholar
- Navér L, Lindgren S, Belfrage E, Gyllensten K, Lidman K, Gisslén M, et al: Children born to HIV-1-infected women in Sweden in 1982-2003: trends in epidemiology and vertical transmission. Journal of Acquired Immune Deficiency Syndrome. 2006, 42: 484-489. 10.1097/01.qai.0000224571.30119.3d.View ArticleGoogle Scholar
- Newell M, Huang S, Fiore S, Thorne C, Mandelbrot L, Sullivan J, et al: Characteristics and management of HIV-1-infected pregnant women enrolled in a randomised trial: differences between Europe and the USA. BMC Infectious Diseases. 2007, 7: 60-10.1186/1471-2334-7-60.View ArticlePubMedPubMed CentralGoogle Scholar
- Victora C, Wagstaff A, Schellenberg J, Gwatkin D, Claeson M, Habicht JP: Applying an equity lens to child health and mortality: more of the same is not enough. Lancet. 2003, 362: 233-41. 10.1016/S0140-6736(03)13917-7.View ArticlePubMedGoogle Scholar
- Johnston N: The Regularisation of Former Mozambican Refugees in South Africa: Experiences and Lessons. 2000, Acornhoek, Refugee Research ProgrammeGoogle Scholar
- Hargreaves J, Collinson M, Kahn K, Clark S, Tollman SM: Childhood mortality among former Mozambican refugees and their hosts in rural South Africa. Int J Epidemiol. 2004, 33 (6): 1271-8. 10.1093/ije/dyh257.View ArticlePubMedGoogle Scholar
- WHO: Reaching the Poor: Challenges for Child Health in the Western Pacific Region. Chapter 7: Barriers to access to child health care. 2007Google Scholar
- Yates R: Women and children first: an appropriate first step towards universal coverage. Bull World Health Organ. 2010, 88 (6): 474-475. 10.2471/BLT.09.074401.View ArticlePubMedPubMed CentralGoogle Scholar
- Benzler J, Sauerborn R: Rapid risk household screening by neonatal arm circumference. Results from a cohort study in Burkina Faso. Tropical Medicine and International Health. 1998, 3: 62-74. 10.1046/j.1365-3156.1998.00340.x.View ArticleGoogle Scholar
- Mulholland E, Smith L, Carneiro I, Becher H, Lehmannd D: Equity and child-survival strategies. Bulletin of the World Health Organization. 2008, 86 (5): 10.2471/BLT.07.044545.
- Antunes J, Waldman E: Trends and spatial distribution of deaths of children aged 12-60 months in Sa~ o Paulo, Brazil, 1980-98. 2002, Bulletin of the World Health Organization, 80: 391-398.Google Scholar
- Elliot P, Wakefield J, Best N, Briggs D: Spatial Epidemiology-Methods and Applications. 2000, Oxford University Press LondonGoogle Scholar
- Gemperli A, Vounatsou P, Kleinschmidt I, Bagayoko M, Lengeler C, Smith T: Spatial patterns of infant mortality in Mali; the effect of malaria endemicity. Am J Epidemiol. 2004, 159: 64-72. 10.1093/aje/kwh001.View ArticlePubMedGoogle Scholar
- Kazembe L, Appleton C, Kleinschmidt I: Spatial analysis of the relationship between early childhood mortality and malaria endemicity in Malawi. Geospatial Health. 2007, 2 (1): 41-50.View ArticlePubMedGoogle Scholar
- Gosoniu L, Vounatsou P, Tami A, Nathan R, Grundmann H, Lengeler C: Spatial effects of mosquito bednets on child mortality. BMC Public Health. 2008, 8: 356-10.1186/1471-2458-8-356.View ArticlePubMedPubMed CentralGoogle Scholar
- Kahn K, Tollman SM, Collinson MA, Clark S, Twine R, Clarke B, et al: Research into health, population and social transitions in rural South Africa: Data and methods of the Agincourt Health and Demographic Surveillance System. Scandinavian Journal of Public Health. 2007, 35 (Suppl 69): 8-20. 10.1080/14034950701505031.View ArticleGoogle Scholar
- Kahn K, Tollman SM, Garenne M, Gear JS: Who dies from what? Determining cause of death in South Africa's rural north-east. Trop Med Int Health. 1999, 4 (6): 433-441. 10.1046/j.1365-3156.1999.00415.x.View ArticlePubMedGoogle Scholar
- Kahn K, Tollman SM, Garenne M, Gear JS: Validation and application of verbal autopsies in a rural area of South Africa. Trop Med Int Health. 2000, 5 (11): 824-831. 10.1046/j.1365-3156.2000.00638.x.View ArticlePubMedGoogle Scholar
- Kahn K: Dying to make a fresh start: mortality and health transition in a new South Africa [dissertation]. 2006, Umeå University, Faculty of Medicine, Public Health and Clinical MedicineGoogle Scholar
- Diggle P, Tawn J, Moyeed R: Model-based geostatistics. Appl Stat. 1998, 47: 299-350.Google Scholar
- Schotman P: Priors for the AR  model: parameterization issues and time series considerations. Econometric Theory. 1994, 10: 579-595. 10.1017/S0266466600008677.View ArticleGoogle Scholar
- Zeller A: An introduction to Bayesian inference and econometrics. 1996, John Wiley and Sons Ltd/Inc: New YorkGoogle Scholar
- Benjamin M, Rigby R, Stasinopoulos D: Generalized autoregressive moving average models. J Am Stat Association. 2003, 98 (461): 214-223. 10.1198/016214503388619238.View ArticleGoogle Scholar
- Gelfand A, Smith A: Sampling-based approaches to calculating marginal densities. J Am Stat Association. 1990, 85: 398-409. 10.2307/2289776.View ArticleGoogle Scholar
- Spiegelhalter D, Best N, Carlin B, van der Linde A: Bayesian measures of model complexity and fit. J R Stat Soc. 2002, 64 (B): 583-639. 10.1111/1467-9868.00353.View ArticleGoogle Scholar
- Gelfand A, Ravishanker N, Ecker M: Modeling and inference for point-referenced binary spatial data. Generalized Linear Models: A Bayesian Perspective. Edited by: Dey D, Ghosh S, Mallick B. 1999, Marcel Dekker Inc, 373-386.Google Scholar
- Tollman S, Kahn K, Sartorius B, Collinson M, Clark S, Garenne M: Implications of mortality transition for primary health care in rural South Africa: a population-based surveillance study. Lancet. 2008, 372: 893-901. 10.1016/S0140-6736(08)61399-9.View ArticlePubMedPubMed CentralGoogle Scholar
- Farah A, Preston S: Child mortality differentials in Sudan. Population and Development Review. 1982, 8: 365-383. 10.2307/1972992.View ArticleGoogle Scholar
- Jain A: Determinants of regional variation in infant mortality in rural India. Infant Mortality in India: Differentials and Determinants. Edited by: Jain A, Visaria L. 1988, Sage Publications, 127-167.Google Scholar
- Cleland J, van Ginneken J: Maternal education and child survival in developing countries: the search for pathways of influence. Soc Sci Med. 1989, 27: 1357-1368. 10.1016/0277-9536(88)90201-8.View ArticleGoogle Scholar
- Schultz T: Interpretation of relations among mortality, economics of the household and the health environment. Proceedings of the Meeting on Socio-Economic Determinants and Consequences of Mortality Mexico City. 1979, Geneva: World Health Organization, June 19-25Google Scholar
- Manda S: Birth intervals, breastfeeding and determinants of childhood mortality in Malawi. Soc Sci Med. 1999, 48 (3): 301-312. 10.1016/S0277-9536(98)00359-1.View ArticlePubMedGoogle Scholar
- Bigdeli M, Annear P: Barriers to access and the purchasing function of health equity funds: lessons from Cambodia. Bull World Health Organ. 2009, 87: 560-564. 10.2471/BLT.08.053058.View ArticlePubMedPubMed CentralGoogle Scholar
- Twine R, Collinson MA, Polzer TJ, Kahn K: Evaluating access to a child-oriented poverty alleviation intervention in rural South Africa. Scand J Public Health. 2007, 69: 118-27. 10.1080/14034950701430008.View ArticleGoogle Scholar
- Goudge J, Gilson L, Russell S, Gumede T, Mills A: The household costs of health care in rural South Africa with free public primary care and hospital exemptions for the poor. Trop Med Int Health. 2009, 14 (4): 458-67. 10.1111/j.1365-3156.2009.02256.x.View ArticlePubMedGoogle Scholar
- Moïsi J, Gatakaa H, Noor A, Williams T, Bauni E, Tsofa B, Levine O, Scott JA: Geographic access to care is not a determinant of child mortality in a rural Kenyan setting with high health facility density. BMC Public Health. 2010, 10: 142-View ArticlePubMedPubMed CentralGoogle Scholar
- Brockerhoff M: The impact of rural-urban migration on child survival. Health Transition Review. 1994, 4: 127-149.PubMedGoogle Scholar
- Bledsoe C, Anastasia B: Child fosterage and child mortality in sub-Saharan Africa: some preliminary questions and answers. Mortality and Society in Sub-Saharan Africa. Edited by: Etienne van de Walle, Gilles Pison and Mpembele Sala-Diakanda. 1992, Oxford: Clarendon Press, 279-302.Google Scholar
- Hildebrandt N, McKenzie D: The effects of migration on child health in Mexico. 2004, Mimeo. Stanford UniversityGoogle Scholar
- Collinson M, Tollman S, Garenne M, Kahn K: Temporary female migration and labour force participation in rural South Africa. 2001, A working paper of the Agincourt Health and Population Unit, University of the Witwatersrand. A working paper of Princeton University's, Center for Migration and Development, Working paper #01-10aGoogle Scholar
- Dolan C, Tollman S, Nkuna V, Gear J: The links between legal status and environmental health: A case study of Mozambican refugees and their hosts in the Mpumalanga (eastern Transvaal) Lowveld, South Africa. International Journal of Health and Human Rights. 1997, 2 (2): 62-84. 10.2307/4065272.View ArticleGoogle Scholar
- Coovadia H, Rollins N, Bland M, Little K, Coutsoudis A, Bennish M, Newell ML: Mother-to-child transmission of HIV-1 infection during exclusive breastfeeding in the first 6 months of life: an intervention cohort study. The Lancet. 2007, 369: 1107-1116. 10.1016/S0140-6736(07)60283-9.View ArticleGoogle Scholar
- WHO: Guidelines on HIV and infant feeding. Principles and recommendations for infant feeding in the context of HIV and a summary of evidence. 2010, WHO, accessed 4 Sep 2010, [http://www.who.int/child_adolescent_health/documents/9789241599535/en/index.html]Google Scholar
- Patwari AK: Diarrhoea and malnutrition interaction. Indian J Pediatr. 1999, 66 (1 Suppl): S124-34.PubMedGoogle Scholar
- Feachem R: Poverty and inequity: a proper focus for the new century. Bull World Health Organ. 2000, 78: 1-2.PubMedPubMed CentralGoogle Scholar
- Gwatkin D, Bhuiya A, Victora C: Making health systems more equitable. Lancet. 2004, 364: 1273-1280. 10.1016/S0140-6736(04)17145-6.View ArticlePubMedGoogle Scholar
- Brocklehurst P, French R: The association between maternal HIV infection and perinatal outcome: a systematic review of the literature and meta-analysis. Br J Obstet Gynaecol. 1998, 105 (8): 836-48.View ArticlePubMedGoogle Scholar
- Garenne M, Fauveau V: Potential and limits of verbal autopsies. 2006, Bulletin of the World Health Organization, 84: 161-256.Google Scholar
- Newell ML, Brahmbhatt H, Ghys PD: Child mortality and HIV infection in Africa: a review. AIDS. 2004, 18 (2): S27-S34. 10.1097/00002030-200406002-00004.View ArticlePubMedGoogle Scholar
- Gelfand A, Vounatsou P: Proper multivariate conditional autoregressive models for spatial data analysis. Biostat. 2003, 4: 11-25. 10.1093/biostatistics/4.1.11.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/10/645/prepub