 Research
 Open Access
 Published:
Analysis of the spatial distribution of infant mortality by cause of death in Austria in 1984 to 2006
International Journal of Health Geographics volume 7, Article number: 21 (2008)
Abstract
Background
In Austria, over the last 20 years infant mortality declined from 11.2 per 1,000 life births (1985) to 4.7 per 1,000 in1997 but remained rather constant since then. In addition to this time trend we already reported a nonrandom spatial distribution of infant mortality rates in a recent study covering the time period 1984 to 2002.
This present study includes four additional years and now covers about 1.9 million individual birth certificates. It aimes to elucidate the observed nonrandom spatial distribution in more detail. We split up infant mortality into six groups according to the underlying cause of death. The underlying spatial distribution of standardized mortality ratios (SMR) is estimated by univariate models as well as by two models incorporating all six groups simultaneously.
Results
We observe strong correlations between the individual spatial patterns of SMR's except for "Sudden Infant Death Syndrome" and to some extent for "Peripartal Problems". The spatial distribution of SMR's is nonrandom with an area of decreased risk in the SouthEast of Austria. The group "Sudden Infant Death Syndrome" clearly and the group "Peripartal Problems" slightly show deviations from the common pattern. When comparing univariate and multivariate SMR estimates we observe that the resulting spatial distributions are very similar.
Conclusion
We observe different nonrandom spatial distributions of infant mortality rates when grouped by cause of death. The models applied were based on individual data thereby avoiding ecological regression bias. The estimated spatial distributions do not substantially depend on the employed estimation method. The observed nonrandom spatial patterns of Austrian infant mortality remain to appear ambiguous.
Background
Infant mortality rate in Austria was higher than average in the European Union until 1987 but now has reached the European Union means [1, 2] (Figure 1). Despite this very welcome temporal trend, in a recent study [3] including the years 1984–2002, we observed an explicit nonuniform spatial distribution indicating lower risks in the SouthEast of Austria. This study was based on about 1.6 million individual birth certificates allowing the adjustment of the infant mortality rates for a large number of covariates. The then observed nonuniform spatial distribution was surprising because the covariates comprised many important anthropometric as well as socioeconomic factors. To our knowledge, these covariates covered the most important prognostic factors, so only minor systematic differences in infant mortality rates due to variables not included in the model were expected. In order to gain more insight for the reasons of the observed spatial distribution for the present study we decided to split infant mortality into six groups defined by the underlying cause of death and to include the births of four more years so that now the study spans the time period of 1984 to 2006.
Methods
Analyses were performed stepwise: 1) Estimation of expected counts of infant deaths for each Austrian administrative district using logistic regression models and individual birth certificates and 2) univariate and multivariate modelling of indirectly standardized mortality rates (SMR's) using Markov Chain Monte Carlo (MCMC) methods in Winbugs [4] based on observed and expected counts.
Data
We analyzed 1,911,625 live births in 99 administrative districts, recorded between 1984 and 2006 in the Birth Certificate Registry managed by Statistics Austria [5]. Inclusion criteria for the data set were that the infants had been born as singletons between the 24^{th} and 44^{th} week of gestation, and to mothers between 13 to 50 years of age. Second or later born children were included if the period to the time of birth of a previous sibling was at least 40 weeks.
From birth certificates we extracted the following variables:

a)
data on the newborn:
survival status after one year, year of birth, gestational age, birth weight, length at birth, sex;

b)
data on the mother:
age, level of education, marital status, time interval from previous birth, parity, place of living (administrative district), means of subsistence, citizenship.
Two adjacent districts were joined because of the small number of their deaths recorded, so that a total of 98 districts entered into analysis.
Infant mortality (within the first year of life) was split up into six groups defined by the underlying cause of death. The following groups were formed: a) "Infections, respiratory diseases"; b) "Peripartal problems"; c) "Immaturity"; d) "Malformations"; e) "Sudden Infant Death Syndrom (SIDS)"; f) "All other". The definition of the groups by ICD9 and ICD10 codes was done by clinical experience [see Additional file 1].
Calculation of expected counts
For each administrative district we calculated expected counts of infant deaths by means of logistic regression models using individual data in SAS [6]. All variables listed above were included except the administrative district. Based on the estimated regression model expected counts were calculated by summing up predicted probabilities for death by administrative district. This procedure was done for overall infant mortality including all six causesofdeath groups as well as for each of the six single groups. As the focus of this study is the description of the six groupspecific spatial distributions, the results of the logistic regression models are not shown. For a detailed description of the effects of the covariates on overall infant mortality we refer to Waldhoer, et al., 2006 [3].
Univariate and multivariate modelling of SMR's
For the univariate estimation of SMR by group of cause of death we used the conditional autoregressive model (CAR) also known as BYM or convolution model introduced by Besag et al. [7]. This model assumes, that the observed number of counts in spatial unit i is Poisson distributed with expectation μ_{i}. Log(μ_{i}) is assumed to be a sum of a spatially structured and unstructured random error S and U, respectively. S describes the information which is common to neighbouring units due to the spatial distribution of the underlying common risk factors. U stands for the spatially unstructured heterogeneity which may not be explained by factors in the model. This kind of model, also often called convolution model, is frequently used in spatial epidemiology. One reason for its popularity is its straightforward estimation by MCMC techniques e.g. in Winbugs [4].
The estimation of univariate CAR models for each causeofdeath group accounts for spatial autocorrelation of SMR's within the groups but ignores the fact that the 6 group specific spatial distributions may share information which could be used for a more precise estimation of the SMR's. In recent years some authors concentrated on the simultaneous analysis of two or more spatial distributions [8–10]. We used two different approaches, a multivariate CAR model and a shared component model, which also can be realized using Winbugs. In a multivariate CAR model (see e.g. Winbugs, Geobugs manual example 7) the univariate Gaussian conditional distribution for the random error terms S is replaced by a multivariate conditional distribution for the six random errors S_{1}S_{6} as in our case. This multivariate conditional distribution has a variancecovariance matrix with diagonal elements representing the conditional variances of the S_{i}'s and offdiagonal elements representing the covariance between the S_{i}'s.
KnorrHeld and Best [10] propose a shared component model where the spatial variation of two diseases is partitioned into a shared component and twodisease specific components (see e.g. Winbugs, Geobugs Manual example 8). Each of both components itself is modelled as a sum of a structured and unstructured random error. In our study we defined the expected number of deaths μ_{ij} for the ith district and kth group as log(μ_{ ik }) = log(E_{ ik }) + α_{ k }+ δ_{ k }* ϕ_{ i }+ ψ_{ ik }
E_{ik} is the expected number based on the logistic regression adjusted for the covariates, α_{k} is a group specific intercept, δ_{k} represents the strength with which the shared random error ϕ_{i} determines the rate of the kth group and Ψ_{ik} is the kth group specific random error of the ith district. Both ϕ and ψ are then partitioned into a structured and unstructured random error. In order to make the model identifiable, the sum of the log(δ_{k}) is set to zero.
The Deviance Information Criterion (DIC) [4] was used for model comparison. Moran's spatial autocorrelation coefficient I [11] was used for testing for a nonrandom uniform spatial distribution.
Results
In Table 1 the number of deaths pertaining to each of the six groups can be seen. The group "malformations" contributes most, while "infections and respiratory diseases" contributes least to the number of deaths.
Overall infant mortality
Figure 2 shows the spatial distribution of the SMR's of overall infant mortality (i.e. all six groups combined). In the SouthEast of Austria the risk of death within the first year of life clearly is lower than in the rest of Austria. The range of the SMR's is 0.83 to 1.21 showing differences in infant mortality of up to 45% within Austria.
Mortality by cause of death estimated by a multivariate CAR model
Figure 3 shows the spatial distributions for the six groups based on the multivariate CAR model. The groups "infections, respiratory diseases", "immaturity", "malformations" and "all other" exhibit a similar spatial distribution with clearly reduced SMR's in the SouthEast of Austria. In the group "peripartal problems" the area of reduced risk is shifted towards the West. "SIDS" has a different spatial distribution inasmuch as the area of lowest risk is clearly in the North of Austria.
In Table 2 the correlation coefficients between the structured components S of the six groups along with their 95%credibility intervals are shown. All correlations are positive apart from the associations with "SIDS". Eight out of 15 credibility intervals do not cover zero and two intervals overlap zero to just a very small amount so that only five of the 15 correlations are left which show an effect just by chance. The associations are strong with values going up to 0.83 and therefore pointing to similar spatial distribution. Correlations with group "SIDS" are all nonpositive pointing to a different spatial distribution of "SIDS" than the rest.
Mortality by cause of death estimated by a shared component model
Figure 4 shows the spatial distribution of the shared component ϕ. This shared component may be interpreted similar to a mean of all six spatial distributions while their importance for group k is described by the weight δ_{k}. Unsurprisingly, this spatial distribution resembles the map of overall infant mortality strongly. The values of δ_{k} in table 3 show that the groups "infections, respiratory diseases" and "all other" are strongly weighted, whereas "immaturity" and "SIDS" are least weighted with the shared component distribution. All credibility intervals are wide.
The ranking of the fractions of the total betweenarea variation in risk is in accordance to the ranking of the δ_{k} (Table 3). This fraction is the ratio of the empirical variance of the shared component, divided by the sum of the variance of shared and specific components. Clearly the most variability explained by the shared component occurs in the groups "infections, respiratory diseases" and "all other", the least fraction is explained in the group "SIDS". Credibility intervals for the fractions shared are very wide (Figure 5). Only the group "SIDS" exhibits a narrow interval which means that this group obviously does not have much in common with all other groups and therefore has a different spatial distribution.
Comparison of models
The DIC value for the univariate model is 3198.8, compared to a DIC of 3169.2 (multivariate CAR model) and 3160.6 (shared component model). Having the smallest DIC value, the shared component model seems to fit the data best.
All spatial autocorrelation coefficients (Moran's I) are significant except for group 6 and range between 0.025 and 0.421 pointing to spatial autocorrelation and nonhomogeneous risks.
Discussion
In a recent paper [3] we were able to show that in the last 2 decades infant mortality was not randomly spatially distributed in Austria. In this present study including four additional years we were able to show that the spatial distribution of mortality within first year of life depends on the cause of death. To our knowledge such remarkable nonrandom spatial distributions have not been described up to now in a small, well developed country with homogeneously spread highquality medical care.
Inequalities in quality of medical care within Austria should especially be perceivable in the groups "Infections, respiratory diseases", "Peripartal problems" and "Immaturity" [12–14]. Because of the observed nonrandom geographical distribution in these groups inequalities in quality of medical care could be reckoned as the underlying cause. However, this hypothesis contradicts the spatial distribution of the groups "Malformations" and "All other". These groups primarily comprise of diseases of which the outcome only weakly depends on quality of medical care as e.g. malformations, tumours, accidents and metabolic diseases and therefore should be distributed randomly in Austria.
The geographical distribution of SIDS nearly is inversely related to the other groups. Among other things, SIDS has been associated with prone and side positions for infant sleep, smoke exposure, soft bedding and bedding surfaces, overheating, use of pacifiers at sleeping time and room sharing without bed sharing [15]. Apart from these risk factors which are not available for our study, the regional distribution of SIDS may also depend on the quality of diagnosis which itself is linked to autopsy. In our study autopsy rates reached a high value of 80.4% but nevertheless autopsy rates also are nonrandomly distributed within Austria [16] and therefore may influence the spatial distribution of SIDS. We intend to examine the reasons for the observed deviant regional distribution of SIDS in the future.
Infant mortality depends on many risk factors as e.g. social status and age of mother, gestational age, etc[3]. Many of these factors have been included in our models so that the observed regional patterns of SMR's are not dependent on the spatial distribution of these risk factors within Austria.
This study differs from many other studies dealing with spatial epidemiology as we used individualbased data and therefore avoided the well known ecological regression bias [17]. The size of the data base (1.9 million observations) prohibited direct model estimation with spatially structured and unstructured random errors because to our knowledge all available statistical software packages are uncapable to handle such large data sets for these kinds of models. Therefore the effect of 13 variables which have been shown to have a significant effect on infant mortality [3] was modelled by logistic regression models without random effects. Subsequently, the spatial structure was analyzed via observed and estimated deaths per district in WinBugs by means of univariate and multivariate based models including spatially structured and unstructured random errors.
The results showed that the estimated SMR's to some extent differ in their absolute values but the resulting choropleth maps using quartiles of SMR's do not and thus maintain the underlying spatial structure. In our opinion these similarities of the six maps based on 3 different estimation methods comes primarily from two facts: Firstly, the expected numbers of counts are large and therewith random variability does not superimpose the underlying spatial effect to a large extent. Secondly, the spatial effect of infant mortality is strong, even when mortality is split up into groups of causes of death.
In this study we used multivariate models analysing six spatial distributions simultaneously since it was to be expected that the spatial distributions are similar due to underlying covariates. The purpose of these kinds of models for disease mapping is borrowing strength for SMR estimation from spatial distributions of other diseases in addition to borrowing strength from neighbouring spatial units of the same disease [10]. This leads to shrinking mortality rates to a global mean based on the different group specific spatial distributions as well as smoothing of rates within disease groups. Consequently, this multivariate approach may produce more precise estimates than their univariate counterparts do.
The use of the multivariate CAR model specifically gives insight into pairwise associations by providing correlation coefficients between groups by cause of death. The CAR model may be contrasted with the shared component model which estimates something similar to a grand mean which apart from δ_{k} is common to all disease groups. The relevance of this shared component may then be quantified by the fraction of variability explained by the shared component. This approach concentrates on the shared component and therefore may be useful when searching for the spatial distribution of underlying risk factors.
The results of our three different models are very similar in terms of the resulting disease maps so that a decision which model is to be preferred is not clear. However, the ranking of the DIC values shows that the model estimating six univariate BYM models has the largest DIC value and therefore seems to be less preferable. The multivariate CAR model exhibits a distinct lower DIC value than the univariate BYM model but a larger value than the shared component model. Apart from the DIC value, the choice between the multivariate CAR and the shared component model may be based on whether ones interest lies in pairwise associations or in detecting a spatial component common to all spatial disease distributions.
In our study we observed different nonrandom spatial distributions of infant mortality rates grouped by cause of death. The models were based on individual data and therefore excluded the possibility of ecological regression bias. In order to propagate the usefulness of two new methods, we used two estimation methods which have rarely been used in practice because of their recency. The results show that in our study the estimation method does not change the spatial distribution substantially.
The results of this as well as our previous study [3] show that the spatial distribution of the independent variables in the logistic regression model, like gestational age, age of mother etc. cannot explain the observed spatial disparities in infant mortality in Austria.
In order to gain more insight into the underlying causes for these disparities additional studies covering larger data sets are necessary. Adding ecological data based on the administrative district level requires rather limited efforts but is prone to the ecological regression bias and may thus only serve as a hypothesis generating approach. Certainly more efficient in terms of causal interpretations would be the addition of individual data which can be gathered in a retrospective or prospective way. Since retrospective or prospective acquisition of data for the total population of all newborns is not affordable, one of two feasible approaches would be a retrospective casecontrol study based on hospital based records comparing data from deceased and surviving newborns. Though gathering clinical data this way seems to be feasible, acquisition of further family related data (e.g. life style of parents etc.) could result in selection bias related problems since the response rate as well as the distribution of variables may differ between parents of deceased and surviving newborns. The second choice would be a prospective cohort study in which a group of newborns are followed to the end of the first year of life. This would allow the inclusion of new interesting variables which may be gathered soon after birth for all newborns. However, this would require a very large number of newborns (n ~200.000) because of the small infant mortality rate and therewith a 1year followup of all Austrian newborns for about 3 consecutive years being not affordable in practice.
Conclusion
In spite of a highly developed health care system we observed a clear spatial disparity in infant mortality in Austria. This is surprising since infant mortality is a sensitive parameter of health care systems and one should expect great efforts on the part of health policy to compensate for these inequalities. We were able to show that official birth certificates containing only limited clinical data give valuable but not sufficient information to explain the observed disparities. We conclude, that further studies are needed which involve additional individual as well as clinical data in order to explain the observed disparities in mortality rates in Austria.
References
 1.
European health for all database, WHO: 2007, Copenhagen, Denmark
 2.
Waldhoer T, Vutuc C, Haidinger G, Mittlböck M, Kirchner L, Wald M: Trends in infant mortality in Austria between 1984 and 2002. Wien Klin Wochenschr. 2005, 117 (15–16): 548553. 10.1007/s0050800504014.
 3.
Waldhoer T, Haidinger G, Wald M, Heinzl H: Nonrandom geographical distribution of infant mortality in Austria 1984–2002. Wien Klin Wochenschr. 2006, 118 (11–12): 3417. 10.1007/s0050800606105.
 4.
Winbugs Spiegelhalter D, Thomas A, Best N, Lunn D: 2004, MRC Biostatistics Unit Cambridge
 5.
Statistik Austria: Direktion Bevölkerung, Geburtsblattregister.
 6.
SAS. NC 27513: SAS Institute Inc., 2002–2003.
 7.
Besag J, York Y, Molli'e A: Bayesian Image Restoration with two applications in spatial statistics. Ann Inst Statist Math. 1991, 43: 159. 10.1007/BF00116466.
 8.
Held L, Natario I, Fenton SE, Rue H, Becker N: Towards joint disease mapping. Stat Methods Med Res. 2005, 14 (1): 6182. 10.1191/0962280205sm389oa.
 9.
Held L, Giusi G, Frank C, Rue H: Joint spatial analysis of gastrointestinal infectious diseases. Stat Methods Med Res. 2006, 15 (5): 46580.
 10.
KnorrHeld L, Best N: A shared component model for detecting joint and selective clustering of two diseases. J R Statist Soc. 2001, 164 (1): 7385. 10.1111/1467985X.00187.
 11.
Moran PAP: Notes on continuous stochastic phenomena. Biometrika. 1950, 37: 1723.
 12.
Armenian SH, Singh J, Arrieta AC: Risk factors for mortality resulting from bloodstream infections in a pediatric intensive care unit. Pediatr Infect Dis J. 2005, 24 (4): 30914. 10.1097/01.inf.0000157086.97503.bd.
 13.
Phibbs CS, Baker LC, Caughey AB, Danielsen B, Schmitt SK, Phibbs RH: Level and volume of neonatal intensive care and mortality in verylowbirthweight infants. N Engl J Med. 2007, 356 (21): 216575. 10.1056/NEJMsa065029.
 14.
Alderliesten ME, Stronks K, van Lith JM, Smit BJ, Wal van der MF, Bonsel GJ, Bleker OP: Ethnic differences in perinatal mortality A perinatal audit on the role of substandard care. Eur J Obstet Gynecol Reprod Biol. 2007
 15.
Moon RY, Horne RS, Hauck FR: Sudden infant death syndrome. Lancet. 2007, 370 (9598): 157887. 10.1016/S01406736(07)616626.
 16.
Waldhoer T, Berzlanovich A, Vutuc C, Haidinger G: Rates of postmortem examination in Austria: the effect of distance between location of death and site of examination. J Clin Epidemiol. 2003, 56 (9): 8915. 10.1016/S08954356(03)000763.
 17.
Richardson S, GuihenneucJouyaux C, Lasserre V: Re: Ecologic studies–biases, misconceptions, and counterexamples. Am J Epidemiol. 1996, 143 (5): 5223.
Acknowledgements
The authors thank Dr Elena Dreassi, University of Firenze, for providing Winbugs code and Dr. Gerald Haidinger, Medical University of Vienna, for editing the manuscript.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
TW had the main responsibility for carrying out the study and writing the manuscript.
MW participated in conceiving the study and writing the medical part of the manuscript.
HH participated in conceiving the study and writing the manuscript.
All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
Waldhoer, T., Wald, M. & Heinzl, H. Analysis of the spatial distribution of infant mortality by cause of death in Austria in 1984 to 2006. Int J Health Geogr 7, 21 (2008) doi:10.1186/1476072X721
Received
Accepted
Published
DOI
Keywords
 Infant Mortality
 Credibility Interval
 Standardize Mortality Ratio
 Birth Certificate
 Administrative District