Le SIDA en Afrique subsaharienne (serveur d'exploration)

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

Childhood mortality in sub-Saharan Africa: cross-sectional insight into small-scale geographical inequalities from Census data

Identifieur interne : 002262 ( Pmc/Corpus ); précédent : 002261; suivant : 002263

Childhood mortality in sub-Saharan Africa: cross-sectional insight into small-scale geographical inequalities from Census data

Auteurs : Lawrence Kazembe ; Aileen Clarke ; Ngianga-Bakwin Kandala

Source :

RBID : PMC:3488715

Abstract

Objectives

To estimate and quantify childhood mortality, its spatial correlates and the impact of potential correlates using recent census data from three sub-Saharan African countries (Rwanda, Senegal and Uganda), where evidence is lacking.

Design

Cross-sectional.

Setting

Nation-wide census samples from three African countries participating in the 2010 African Census round. All three countries have conducted recent censuses and have information on mortality of children under 5 years.

Participants

111 288 children under the age of 5 years in three countries.

Primary and secondary outcome measures

Under-five mortality was assessed alongside potential correlates including geographical location (where children live), and environmental, bio-demographic and socioeconomic variables.

Results

Multivariate analysis indicates that in all three countries the overall risk of child death in the first 5 years of life has decreased in recent years (Rwanda: HR=0.04, 95% CI 0.02 to 0.09; Senegal: HR=0.02 (95% CI 0.02 to 0.05); Uganda: HR=0.011 (95% CI 0.006 to 0.018). In Rwanda, lower deaths were associated with living in urban areas (0.79, 0.73, 0.83), children with living mother (HR=0.16, 95% CI 0.15 to 0.17) or living father (HR=0.38, 95% CI 0.36 to 0.39). Higher death was associated with male children (HR=1.06, 95% CI 1.02 to 1.08) and Christian children (HR=1.14, 95% CI 1.05 to 1.27). Children less than 1 year were associated with higher risk of death compared to older children in the three countries. Also, there were significant spatial variations showing inequalities in children mortality by geographic location. In Uganda, for example, areas of high risk are in the south-west and north-west and Kampala district showed a significantly reduced risk.

Conclusions

We provide clear evidence of considerable geographical variation of under-five mortality which is unexplained by factors considered in the data. The resulting under-five mortality maps can be used as a practical tool for monitoring progress within countries for the Millennium Development Goal 4 to reduce under-five mortality in half by 2015.


Url:
DOI: 10.1136/bmjopen-2012-001421
PubMed: 23089207
PubMed Central: 3488715

Links to Exploration step

PMC:3488715

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Childhood mortality in sub-Saharan Africa: cross-sectional insight into small-scale geographical inequalities from Census data</title>
<author>
<name sortKey="Kazembe, Lawrence" sort="Kazembe, Lawrence" uniqKey="Kazembe L" first="Lawrence" last="Kazembe">Lawrence Kazembe</name>
<affiliation>
<nlm:aff id="af1">Department of Statistics, University of Namibia, Windhoek, Namibia</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Clarke, Aileen" sort="Clarke, Aileen" uniqKey="Clarke A" first="Aileen" last="Clarke">Aileen Clarke</name>
<affiliation>
<nlm:aff id="af2">Division of Health Sciences, University of Warwick, Warwick Medical School, Coventry, Warks, UK</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Kandala, Ngianga Bakwin" sort="Kandala, Ngianga Bakwin" uniqKey="Kandala N" first="Ngianga-Bakwin" last="Kandala">Ngianga-Bakwin Kandala</name>
<affiliation>
<nlm:aff id="af2">Division of Health Sciences, University of Warwick, Warwick Medical School, Coventry, Warks, UK</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="af3">Malaria Public Health and Epidemiology Group, Centre for Geographic Medicine, KEMRI-University of Oxford-Wellcome Trust Collaborative Programme, Nairobi, Kenya</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">23089207</idno>
<idno type="pmc">3488715</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3488715</idno>
<idno type="RBID">PMC:3488715</idno>
<idno type="doi">10.1136/bmjopen-2012-001421</idno>
<date when="2012">2012</date>
<idno type="wicri:Area/Pmc/Corpus">002262</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">002262</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Childhood mortality in sub-Saharan Africa: cross-sectional insight into small-scale geographical inequalities from Census data</title>
<author>
<name sortKey="Kazembe, Lawrence" sort="Kazembe, Lawrence" uniqKey="Kazembe L" first="Lawrence" last="Kazembe">Lawrence Kazembe</name>
<affiliation>
<nlm:aff id="af1">Department of Statistics, University of Namibia, Windhoek, Namibia</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Clarke, Aileen" sort="Clarke, Aileen" uniqKey="Clarke A" first="Aileen" last="Clarke">Aileen Clarke</name>
<affiliation>
<nlm:aff id="af2">Division of Health Sciences, University of Warwick, Warwick Medical School, Coventry, Warks, UK</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Kandala, Ngianga Bakwin" sort="Kandala, Ngianga Bakwin" uniqKey="Kandala N" first="Ngianga-Bakwin" last="Kandala">Ngianga-Bakwin Kandala</name>
<affiliation>
<nlm:aff id="af2">Division of Health Sciences, University of Warwick, Warwick Medical School, Coventry, Warks, UK</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="af3">Malaria Public Health and Epidemiology Group, Centre for Geographic Medicine, KEMRI-University of Oxford-Wellcome Trust Collaborative Programme, Nairobi, Kenya</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">BMJ Open</title>
<idno type="eISSN">2044-6055</idno>
<imprint>
<date when="2012">2012</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<sec>
<title>Objectives</title>
<p>To estimate and quantify childhood mortality, its spatial correlates and the impact of potential correlates using recent census data from three sub-Saharan African countries (Rwanda, Senegal and Uganda), where evidence is lacking.</p>
</sec>
<sec>
<title>Design</title>
<p>Cross-sectional.</p>
</sec>
<sec>
<title>Setting</title>
<p>Nation-wide census samples from three African countries participating in the 2010 African Census round. All three countries have conducted recent censuses and have information on mortality of children under 5 years.</p>
</sec>
<sec>
<title>Participants</title>
<p>111 288 children under the age of 5 years in three countries.</p>
</sec>
<sec>
<title>Primary and secondary outcome measures</title>
<p>Under-five mortality was assessed alongside potential correlates including geographical location (where children live), and environmental, bio-demographic and socioeconomic variables.</p>
</sec>
<sec>
<title>Results</title>
<p>Multivariate analysis indicates that in all three countries the overall risk of child death in the first 5 years of life has decreased in recent years (Rwanda: HR=0.04, 95% CI 0.02 to 0.09; Senegal: HR=0.02 (95% CI 0.02 to 0.05); Uganda: HR=0.011 (95% CI 0.006 to 0.018). In Rwanda, lower deaths were associated with living in urban areas (0.79, 0.73, 0.83), children with living mother (HR=0.16, 95% CI 0.15 to 0.17) or living father (HR=0.38, 95% CI 0.36 to 0.39). Higher death was associated with male children (HR=1.06, 95% CI 1.02 to 1.08) and Christian children (HR=1.14, 95% CI 1.05 to 1.27). Children less than 1 year were associated with higher risk of death compared to older children in the three countries. Also, there were significant spatial variations showing inequalities in children mortality by geographic location. In Uganda, for example, areas of high risk are in the south-west and north-west and Kampala district showed a significantly reduced risk.</p>
</sec>
<sec>
<title>Conclusions</title>
<p>We provide clear evidence of considerable geographical variation of under-five mortality which is unexplained by factors considered in the data. The resulting under-five mortality maps can be used as a practical tool for monitoring progress within countries for the Millennium Development Goal 4 to reduce under-five mortality in half by 2015.</p>
</sec>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Becher, H" uniqKey="Becher H">H Becher</name>
</author>
<author>
<name sortKey="Muller, O" uniqKey="Muller O">O Muller</name>
</author>
<author>
<name sortKey="Jahn, A" uniqKey="Jahn A">A Jahn</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Armstrong Schellenberg, Jrm" uniqKey="Armstrong Schellenberg J">JRM Armstrong-Schellenberg</name>
</author>
<author>
<name sortKey="Nathan, R" uniqKey="Nathan R">R Nathan</name>
</author>
<author>
<name sortKey="Abdulla, S" uniqKey="Abdulla S">S Abdulla</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Binka, Fn" uniqKey="Binka F">FN Binka</name>
</author>
<author>
<name sortKey="Indome, F" uniqKey="Indome F">F Indome</name>
</author>
<author>
<name sortKey="Smith, T" uniqKey="Smith T">T Smith</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Silverman, Jg" uniqKey="Silverman J">JG Silverman</name>
</author>
<author>
<name sortKey="Decker, Mr" uniqKey="Decker M">MR Decker</name>
</author>
<author>
<name sortKey="Wirth, K" uniqKey="Wirth K">K Wirth</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kazembe, L" uniqKey="Kazembe L">L Kazembe</name>
</author>
<author>
<name sortKey="Mpeketula, Pmg" uniqKey="Mpeketula P">PMG Mpeketula</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Johnson, Hl" uniqKey="Johnson H">HL Johnson</name>
</author>
<author>
<name sortKey="Liu, L" uniqKey="Liu L">L Liu</name>
</author>
<author>
<name sortKey="Fischer Walker, C" uniqKey="Fischer Walker C">C Fischer-Walker</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Magadi, Ma" uniqKey="Magadi M">MA Magadi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kandala, Nb" uniqKey="Kandala N">NB Kandala</name>
</author>
<author>
<name sortKey="Fahrmeir, L" uniqKey="Fahrmeir L">L Fahrmeir</name>
</author>
<author>
<name sortKey="Klasen, S" uniqKey="Klasen S">S Klasen</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="De Sherbinin, A" uniqKey="De Sherbinin A">A de Sherbinin</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Wang, L" uniqKey="Wang L">L Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Haines, Mr" uniqKey="Haines M">MR Haines</name>
</author>
<author>
<name sortKey="Preston, Sh" uniqKey="Preston S">SH Preston</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Box Steffensmeier, Jm" uniqKey="Box Steffensmeier J">JM Box-Steffensmeier</name>
</author>
<author>
<name sortKey="Jones, Bs" uniqKey="Jones B">BS Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Cox, Dr" uniqKey="Cox D">DR Cox</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Sastry, N" uniqKey="Sastry N">N Sastry</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bolstad, Wm" uniqKey="Bolstad W">WM Bolstad</name>
</author>
<author>
<name sortKey="Manda, So" uniqKey="Manda S">SO Manda</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kandala, Nb" uniqKey="Kandala N">NB Kandala</name>
</author>
<author>
<name sortKey="Ghilagaber, G" uniqKey="Ghilagaber G">G Ghilagaber</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Adebayo, Sb" uniqKey="Adebayo S">SB Adebayo</name>
</author>
<author>
<name sortKey="Fahrmeir, L" uniqKey="Fahrmeir L">L Fahrmeir</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Caselli, G" uniqKey="Caselli G">G Caselli</name>
</author>
<author>
<name sortKey="Pozzi, L" uniqKey="Pozzi L">L Pozzi</name>
</author>
<author>
<name sortKey="Vaupel, Jw" uniqKey="Vaupel J">JW Vaupel</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hennerfeind, A" uniqKey="Hennerfeind A">A Hennerfeind</name>
</author>
<author>
<name sortKey="Brezger, A" uniqKey="Brezger A">A Brezger</name>
</author>
<author>
<name sortKey="Fahrmeir, L" uniqKey="Fahrmeir L">L Fahrmeir</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hanagal, Dd" uniqKey="Hanagal D">DD Hanagal</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Vaupel, Jw" uniqKey="Vaupel J">JW Vaupel</name>
</author>
<author>
<name sortKey="Manton, Kg" uniqKey="Manton K">KG Manton</name>
</author>
<author>
<name sortKey="Stallard, E" uniqKey="Stallard E">E Stallard</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Banerjee, S" uniqKey="Banerjee S">S Banerjee</name>
</author>
<author>
<name sortKey="Wall, Mm" uniqKey="Wall M">MM Wall</name>
</author>
<author>
<name sortKey="Carlin, Bp" uniqKey="Carlin B">BP Carlin</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bastos, Ls" uniqKey="Bastos L">LS Bastos</name>
</author>
<author>
<name sortKey="Gamerman, D" uniqKey="Gamerman D">D Gamerman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Li, Y" uniqKey="Li Y">Y Li</name>
</author>
<author>
<name sortKey="Ryan, L" uniqKey="Ryan L">L Ryan</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gemperli, A" uniqKey="Gemperli A">A Gemperli</name>
</author>
<author>
<name sortKey="Vounatsou, P" uniqKey="Vounatsou P">P Vounatsou</name>
</author>
<author>
<name sortKey="Kleinschmidt, I" uniqKey="Kleinschmidt I">I Kleinschmidt</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Escarams, G" uniqKey="Escarams G">G Escarams</name>
</author>
<author>
<name sortKey="Carrasco, Jl" uniqKey="Carrasco J">JL Carrasco</name>
</author>
<author>
<name sortKey="Aponte, Jj" uniqKey="Aponte J">JJ Aponte</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kamman, Ee" uniqKey="Kamman E">EE Kamman</name>
</author>
<author>
<name sortKey="Wand, Mp" uniqKey="Wand M">MP Wand</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kneib, T" uniqKey="Kneib T">T Kneib</name>
</author>
<author>
<name sortKey="Fahrmeir, L" uniqKey="Fahrmeir L">L Fahrmeir</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Brezger, A" uniqKey="Brezger A">A Brezger</name>
</author>
<author>
<name sortKey="Kneib, T" uniqKey="Kneib T">T Kneib</name>
</author>
<author>
<name sortKey="Lang, S" uniqKey="Lang S">S Lang</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Brooks, Sp" uniqKey="Brooks S">SP Brooks</name>
</author>
<author>
<name sortKey="Gelman, A" uniqKey="Gelman A">A Gelman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gould, Jb" uniqKey="Gould J">JB Gould</name>
</author>
<author>
<name sortKey="Herrchen, B" uniqKey="Herrchen B">B Herrchen</name>
</author>
<author>
<name sortKey="Pham, T" uniqKey="Pham T">T Pham</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Manda, Som" uniqKey="Manda S">SOM Manda</name>
</author>
<author>
<name sortKey="Feltbower, Rg" uniqKey="Feltbower R">RG Feltbower</name>
</author>
<author>
<name sortKey="Gilhorpe, Ms" uniqKey="Gilhorpe M">MS Gilhorpe</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Browne, Wj" uniqKey="Browne W">WJ Browne</name>
</author>
<author>
<name sortKey="Goldstein, H" uniqKey="Goldstein H">H Goldstein</name>
</author>
<author>
<name sortKey="Rasbash, J" uniqKey="Rasbash J">J Rasbash</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zhou, L" uniqKey="Zhou L">L Zhou</name>
</author>
<author>
<name sortKey="Hanson, Te" uniqKey="Hanson T">TE Hanson</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Carlin, Bp" uniqKey="Carlin B">BP Carlin</name>
</author>
<author>
<name sortKey="Banerjee, S" uniqKey="Banerjee S">S Banerjee</name>
</author>
<author>
<name sortKey="Wall, M" uniqKey="Wall M">M Wall</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gastafson, P" uniqKey="Gastafson P">P Gastafson</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Spiegelhalter, Dj" uniqKey="Spiegelhalter D">DJ Spiegelhalter</name>
</author>
<author>
<name sortKey="Best, Ng" uniqKey="Best N">NG Best</name>
</author>
<author>
<name sortKey="Carlin, Bp" uniqKey="Carlin B">BP Carlin</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article">
<pmc-dir>properties open_access</pmc-dir>
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">BMJ Open</journal-id>
<journal-id journal-id-type="iso-abbrev">BMJ Open</journal-id>
<journal-id journal-id-type="hwp">bmjopen</journal-id>
<journal-id journal-id-type="publisher-id">bmjopen</journal-id>
<journal-title-group>
<journal-title>BMJ Open</journal-title>
</journal-title-group>
<issn pub-type="epub">2044-6055</issn>
<publisher>
<publisher-name>BMJ Publishing Group</publisher-name>
<publisher-loc>BMA House, Tavistock Square, London, WC1H 9JR</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">23089207</article-id>
<article-id pub-id-type="pmc">3488715</article-id>
<article-id pub-id-type="publisher-id">bmjopen-2012-001421</article-id>
<article-id pub-id-type="doi">10.1136/bmjopen-2012-001421</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Epidemiology</subject>
<subj-group>
<subject>Research</subject>
</subj-group>
</subj-group>
<subj-group subj-group-type="hwp-journal-coll">
<subject>1506</subject>
<subject>1692</subject>
<subject>1699</subject>
<subject>1719</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Childhood mortality in sub-Saharan Africa: cross-sectional insight into small-scale geographical inequalities from Census data</article-title>
<alt-title alt-title-type="short">Childhood mortality in sub-Saharan Africa</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Kazembe</surname>
<given-names>Lawrence</given-names>
</name>
<xref ref-type="aff" rid="af1">1</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Clarke</surname>
<given-names>Aileen</given-names>
</name>
<xref ref-type="aff" rid="af2">2</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kandala</surname>
<given-names>Ngianga-Bakwin</given-names>
</name>
<xref ref-type="aff" rid="af2">2</xref>
<xref ref-type="aff" rid="af3">3</xref>
</contrib>
</contrib-group>
<aff id="af1">
<label>1</label>
Department of Statistics, University of Namibia, Windhoek, Namibia</aff>
<aff id="af2">
<label>2</label>
Division of Health Sciences, University of Warwick, Warwick Medical School, Coventry, Warks, UK</aff>
<aff id="af3">
<label>3</label>
Malaria Public Health and Epidemiology Group, Centre for Geographic Medicine, KEMRI-University of Oxford-Wellcome Trust Collaborative Programme, Nairobi, Kenya</aff>
<author-notes>
<corresp>
<label>Correspondence to</label>
Dr Kandala Ngianga-Bakwin;
<email>N-B.Kandala@warwick.ac.uk</email>
</corresp>
</author-notes>
<pub-date pub-type="ppub">
<year>2012</year>
</pub-date>
<pub-date pub-type="collection">
<year>2012</year>
</pub-date>
<pub-date pub-type="epub">
<day>22</day>
<month>10</month>
<year>2012</year>
</pub-date>
<volume>2</volume>
<issue>5</issue>
<elocation-id>e001421</elocation-id>
<history>
<date date-type="received">
<day>13</day>
<month>6</month>
<year>2012</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>8</month>
<year>2012</year>
</date>
</history>
<permissions>
<copyright-statement>Published by the BMJ Publishing Group Limited. For permission to use (where not already granted under a licence) please go to http://group.bmj.com/group/rights-licensing/permissions</copyright-statement>
<copyright-year>2012</copyright-year>
<license license-type="open-access">
<license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See:
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/2.0/">http://creativecommons.org/licenses/by-nc/2.0/</ext-link>
and
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/2.0/legalcode">http://creativecommons.org/licenses/by-nc/2.0/legalcode</ext-link>
.</license-p>
</license>
</permissions>
<self-uri xlink:title="pdf" xlink:type="simple" xlink:href="bmjopen-2012-001421.pdf"></self-uri>
<self-uri content-type="draft-revisions-pdf" xlink:type="simple" xlink:href="bmjopen-2012-001421.draft_revisions.pdf"></self-uri>
<self-uri content-type="reviewers-comments-pdf" xlink:type="simple" xlink:href="bmjopen-2012-001421.reviewer_comments.pdf"></self-uri>
<abstract>
<sec>
<title>Objectives</title>
<p>To estimate and quantify childhood mortality, its spatial correlates and the impact of potential correlates using recent census data from three sub-Saharan African countries (Rwanda, Senegal and Uganda), where evidence is lacking.</p>
</sec>
<sec>
<title>Design</title>
<p>Cross-sectional.</p>
</sec>
<sec>
<title>Setting</title>
<p>Nation-wide census samples from three African countries participating in the 2010 African Census round. All three countries have conducted recent censuses and have information on mortality of children under 5 years.</p>
</sec>
<sec>
<title>Participants</title>
<p>111 288 children under the age of 5 years in three countries.</p>
</sec>
<sec>
<title>Primary and secondary outcome measures</title>
<p>Under-five mortality was assessed alongside potential correlates including geographical location (where children live), and environmental, bio-demographic and socioeconomic variables.</p>
</sec>
<sec>
<title>Results</title>
<p>Multivariate analysis indicates that in all three countries the overall risk of child death in the first 5 years of life has decreased in recent years (Rwanda: HR=0.04, 95% CI 0.02 to 0.09; Senegal: HR=0.02 (95% CI 0.02 to 0.05); Uganda: HR=0.011 (95% CI 0.006 to 0.018). In Rwanda, lower deaths were associated with living in urban areas (0.79, 0.73, 0.83), children with living mother (HR=0.16, 95% CI 0.15 to 0.17) or living father (HR=0.38, 95% CI 0.36 to 0.39). Higher death was associated with male children (HR=1.06, 95% CI 1.02 to 1.08) and Christian children (HR=1.14, 95% CI 1.05 to 1.27). Children less than 1 year were associated with higher risk of death compared to older children in the three countries. Also, there were significant spatial variations showing inequalities in children mortality by geographic location. In Uganda, for example, areas of high risk are in the south-west and north-west and Kampala district showed a significantly reduced risk.</p>
</sec>
<sec>
<title>Conclusions</title>
<p>We provide clear evidence of considerable geographical variation of under-five mortality which is unexplained by factors considered in the data. The resulting under-five mortality maps can be used as a practical tool for monitoring progress within countries for the Millennium Development Goal 4 to reduce under-five mortality in half by 2015.</p>
</sec>
</abstract>
</article-meta>
</front>
<body>
<boxed-text position="float">
<caption>
<title>Article summary</title>
</caption>
<sec>
<title>Article focus</title>
<list list-type="bullet">
<list-item>
<p>Census and household data contain small-area geographical information, such as the residence of a child at the time of death. The impact of such spatial effects on mortality is of substantive interest.</p>
</list-item>
<list-item>
<p>Mortality data in the census and household are also subject to recall bias and heaping effects. The event of deaths is also censored at the time of the survey. We use recent statistical techniques to account for the survival nature of the data in the analysis.</p>
</list-item>
<list-item>
<p>We use data from recent census/household data from three sub-Saharan Africa (Rwanda, Senegal and Uganda) to investigate the importance of country-specific geographical factors on under-five mortality.</p>
</list-item>
</list>
</sec>
<sec>
<title>Key messages</title>
<list list-type="bullet">
<list-item>
<p>Our results provide clear evidence of considerable geographical inequalities of under-five mortality that is unexplained by socio-economic factors considered in the data.</p>
</list-item>
<list-item>
<p>Our findings indicate that public health interventions and health promotion to reduce under-five mortality, should take into account both individual and area variation to account for the diversity of settings.</p>
</list-item>
<list-item>
<p>Planning and intervention measures could have different outcomes in terms of effectiveness in areas with a high degree of variability. Homogeneous policy intervention strategies may not give the required outcomes as suggested by large significant inequalities in of under-five mortality within and between countries in our study.</p>
</list-item>
</list>
</sec>
<sec>
<title>Strengths and limitations of this study</title>
<list list-type="bullet">
<list-item>
<p>Our study seems the only one that has attempted to investigate geographical inequalities of under-five mortality beyond individual and household factors using merged census and household data from sub-Saharan Africa countries.</p>
</list-item>
<list-item>
<p>The major strength is the use of census and nationally representative household survey to investigate and explain district-level inequalities in under-five mortality using a novel approach that accounts simultaneously for individual, household and area factors.</p>
</list-item>
<list-item>
<p>The major limitation of this study is the cross-sectional nature of the data, which does not permit one to draw causal association between under-five mortality and the associated spatial effects including individual and household factors.</p>
</list-item>
</list>
</sec>
</boxed-text>
<sec id="s1">
<title>Introduction</title>
<p>There have been considerable gains in child survival in the world over the past 10 years. Recent reports on the ‘State of the World's Children’ indicate an overall decline in child mortality from 100/1000 children in 1999 to 72/1000 children in 2010.
<xref ref-type="bibr" rid="R1">1</xref>
A number of countries in North Africa, Eastern Europe, south-east Asia and Latin America have reduced under-five child mortality by half in the period between 1990 and 2010.
<xref ref-type="bibr" rid="R1">1</xref>
In contrast, countries in sub-Saharan Africa (SSA) have remarkably high rates and were ranked the top worst performers in reduction of child mortality, with very few making progress and the majority experiencing no change or a reversal in gains made some 10 or so years previously. In Ethiopia, Malawi and Namibia the decline in under-five child mortality has been substantial despite meagre resources, while in other SSA countries it has remained the same. A case in point is DR Congo which posted an under-five mortality of 199/1000 children in 1999 and the same again in 2009. In countries such as Chad, there has been an increase from 201/1000 children born in 1999 to 209/1000 children in 2009. In short, even with known solutions and international assistance, the transition from high mortality to low mortality is highly uneven in the SSA region.</p>
<p>Several studies have shown that child survival in the first 5 years of life is influenced by a myriad of risk factors. For instance, Becher
<italic>et al</italic>
<xref ref-type="bibr" rid="R2">2</xref>
quantified the effect of risk factors for childhood mortality in a typical rural setting of Burkina Faso. They performed a survival analysis of births within a population from a demographic surveillance system in 39 villages. In another study in rural Tanzania, Armstrong-Schellenberg
<italic>et al</italic>
<xref ref-type="bibr" rid="R3">3</xref>
conducted a community-based nested case–control study of postneonatal deaths in children under 5 years, in which they investigated demographic and socio-economic factors, health-seeking behaviour, the household environment including accessibility to healthcare and individual child care factors. A similar population-based case–control study was carried out to investigate potential risk factors for postneonatal and child mortality in northern Ghana.
<xref ref-type="bibr" rid="R4">4</xref>
Child mortality demonstrated gender-based disparities,
<xref ref-type="bibr" rid="R5">5</xref>
varied with socio-economic inequalities
<xref ref-type="bibr" rid="R2">2</xref>
and was influenced by variation in coverage of interventions.
<xref ref-type="bibr" rid="R4">4</xref>
At times, living in either urban or rural areas can disadvantage under-five children's health.
<xref ref-type="bibr" rid="R6">6</xref>
The general picture is that major causes of childhood mortality, summarised as disease and malnutrition, are exacerbated by socio-economic differences and varied intervention coverage,
<xref ref-type="bibr" rid="R7">7</xref>
and these risk factors apply at both individual and community levels.
<xref ref-type="bibr" rid="R8">8</xref>
</p>
<p>Many of these studies have been conducted at subnational or national level. Our search showed that relatively few studies have considered between country or cross-national analyses of childhood health and associated determinants. Magadi
<xref ref-type="bibr" rid="R8">8</xref>
examined the risk factors of malnutrition among children whose mothers are infected with HIV in SSA. She applied a multilevel logistic regression to the Demographic and Health Survey (DHS) data from 18 countries for the period 2003–2008. Another across-countries study was conducted by Kandala
<italic>et al</italic>
<xref ref-type="bibr" rid="R9">9</xref>
in which they considered geographical and socio-economic determinants of child undernutrition in Malawi, Zambia and Tanzania. This identified regional patterns which transcend national boundaries. In a similar study, Sherbinin
<xref ref-type="bibr" rid="R10">10</xref>
reported on biophysical and geographical correlates of child malnutrition in Africa. Wang,
<xref ref-type="bibr" rid="R11">11</xref>
using data from 60 countries, explored the global pattern of child mortality and investigated the determinants both at national and subnational level.</p>
<p>While tackling the issue of determinants of child mortality, it is common that sample survey data have been used. Only a few studies have considered census data, mostly to analyse demographic indicators such as fertility, yet with very limited statistical modelling of childhood mortality.
<xref ref-type="bibr" rid="R12">12</xref>
Census data provide a large cross-sectional database that would allow investigation of the association between mortality or health outcomes and risk factors. Because census data are a complete enumeration of all individuals in a country, the statistical analysis is likely to have more power than data derived from a survey. Further the censuses provide a picture of the country at any given time, therefore allowing a better understanding of risk factors, critical in explaining variations and crucial for implementing interventions in child mortality, thus.</p>
</sec>
<sec id="s2">
<title>Overview of the analysis of mortality data</title>
<p>A number of statistical models have been proposed when analysing the risk of child mortality in the first 5 years of life and its determinants. Most popular have been logistic regression models which assume child survival as a binary response (either the child lived beyond 5 years or died before the fifth birthday). In such models one estimates the probability of a child surviving and can include risk factors. The coefficients of risk factors can be interpreted as ORs. These models, nevertheless, ignore the time to event (death), and therefore fail to capture exposure to the risk of dying or conceal the evolution of the subject's state over time.
<xref ref-type="bibr" rid="R13">13</xref>
More appropriately, survival models can be used to analyse the hazards of child survival. Both logistic and survival analysis can be implemented within the basic generalised linear models (GLM) framework.</p>
<p>Research on survival analysis in demography and related fields has increased since the seminal work by Cox in the 1970s,
<xref ref-type="bibr" rid="R14">14</xref>
with application in child mortality or survival appearing in the 1990s and 2000s,
<xref ref-type="bibr" rid="R15 R16 R17">15–17</xref>
involving both standard proportional hazard models
<xref ref-type="bibr" rid="R12">12</xref>
and complex models.
<xref ref-type="bibr" rid="R18">18</xref>
Precursor to the use of Cox regression, life tables have been used to estimate the probabilities of survival of a given cohort.
<xref ref-type="bibr" rid="R19">19</xref>
The Cox regression and recent modern survival techniques, in contrast, permit effects of risk factors (determinants) to be estimated, which is critical for designing interventions. Cox and other survival regression models have other advantages including analysis of censored and truncated response data, and time-varying effects.
<xref ref-type="bibr" rid="R20">20</xref>
</p>
<p>Two important extensions in survival regression models that have received considerable attention recently are the inclusion of random effects and flexible modelling through semiparametric and nonparametric approach.
<xref ref-type="bibr" rid="R21">21</xref>
Such analyses have an added advantage compared to ordinary GLM. Models that incorporate random effects are commonly called generalised linear mixed models (GLMM), and those that account for non-linearity are referred to as generalised additive models (GAM) and when extended to include random effects they are known as generalised additive mixed models (GAMM).</p>
<p>The inclusion of random effects permits modelling of unmeasured and unobserved contextual factors in the models. These may act at family, community, district, regional or national levels since the underlying causes of neonatal mortality are multisectoral and interwoven.
<xref ref-type="bibr" rid="R8">8</xref>
<xref ref-type="bibr" rid="R17">17</xref>
Those operating at individual, family, community and regional levels can have a direct or intermediary effect on the outcome. Regionally, expenditure on health services and cultural differences can also affect the survival status of children. In survival models, unobserved factors are considered as frailties, which adjust for hierarchical clustering of survival data. In essence, frailties are group specific factors acting on child survival, which together with individual factors may protect or accelerate death.
<xref ref-type="bibr" rid="R22">22</xref>
Studies on child mortality by Sastry
<xref ref-type="bibr" rid="R15">15</xref>
employed nested frailty models to analyse child survival data clustered at community and family levels. Another study by Bolstad and Manda
<xref ref-type="bibr" rid="R16">16</xref>
showed significant heterogeneity at community level, which can be attributed to differential availability of resources at community level.</p>
<p>Recent studies have assumed that such unobserved factors vary spatially to give spatial frailty survival models. Banerjee
<italic>et al</italic>
<xref ref-type="bibr" rid="R23">23</xref>
developed parametric frailty specifications based on both areal (lattice) and on point-referenced (geostatistical) spatial models, and compared them with traditional independent and identically distributed frailty and non-frailty approaches under a Weibull baseline hazard function in the context of county-level infant mortality data. Bastos and Gamerman
<xref ref-type="bibr" rid="R24">24</xref>
used a dynamic survival model with spatial frailty to handle time-varying covariates in the presence of spatial effects. In another study, Li and Ryan
<xref ref-type="bibr" rid="R25">25</xref>
used a semiparametric frailty model to analyse spatial survival models. In the above studies, a single modelling framework was used to model the spatial and time-varying effects simultaneously. Moreover, some of the risk factors may be geographically varying, leading child mortality to vary in space,
<xref ref-type="bibr" rid="R6">6</xref>
<xref ref-type="bibr" rid="R18">18</xref>
<xref ref-type="bibr" rid="R26">26</xref>
with well documented space–time interactions.
<xref ref-type="bibr" rid="R27">27</xref>
In another study, Kandala
<italic>et al</italic>
<xref ref-type="bibr" rid="R9">9</xref>
applied a geo-additive model where spatial, non-linear and fixed effects were simultaneously modelled in a single framework. Again, none of these considered the census data.</p>
<p>The main objective of this study is to analyse small-scale geographical variability in under-five mortality in the sub-Saharan region, by applying existing spatial statistical methodology. Our aim is to extend the standard Cox regression model to a random-effects model to permit spatial clustering and heterogeneity using census data from a number of countries. Specifically, we apply GLMM with spatially correlated random effects proposed by Hennerfeind
<italic>et al</italic>
,
<xref ref-type="bibr" rid="R20">20</xref>
and used it to analyse factors associated with child survival in the first 5 years of life. This modelling approach falls within a group termed structured additive regression (STAR) models, introduced by Kamman and Wand.
<xref ref-type="bibr" rid="R28">28</xref>
STAR models are a comprehensive class of models which permit simultaneous estimation of nonlinear effects of continuous covariates, with both spatially unstructured and structured components, together with the usual fixed effects in the predictor.
<xref ref-type="bibr" rid="R29">29</xref>
</p>
<p>The spatial analysis can be approached in various ways, one based on the stationary Gaussian random fields, which applies when the place of residence is known exactly, given by geographical x−y coordinates, and the principle originates from geostatistics.
<xref ref-type="bibr" rid="R23">23</xref>
These can also be interpreted as two-dimensional surface smoothers based on radial basis functions, and have been employed by Kammann and Wand
<xref ref-type="bibr" rid="R28">28</xref>
to model the spatial component in Gaussian regression models. Another option is to use two-dimensional P-splines described in more detail in Brezger
<italic>et al.</italic>
<xref ref-type="bibr" rid="R30">30</xref>
The advantage of these approaches is that they allow prediction of risk for locations where there are no data, thus allowing us to quantify small-scale variability. If observations are aggregated in geographical regions, spatial effects can be estimated using the Markov random field (MRF) approach, widely used in disease mapping.
<xref ref-type="bibr" rid="R23">23</xref>
Modelling and inference can use a fully Bayesian approach. However, the empirical Bayesian approach via penalised likelihood techniques is also possible.
<xref ref-type="bibr" rid="R29">29</xref>
</p>
<p>A detailed description of the statistical methodology used is in the appendix. In this study, we provide simulation studies and apply these techniques to the 2000–2010 Census data from selected SSA countries whose under-five death rates are at relatively similar ranking. Rwanda and Uganda are ranked 31st and 27th in the world respectively with regard to improvement in under-five mortality, while Senegal is ranked at 42nd. Rwanda and Uganda have an estimated under-five mortality of 91 and 99/1000 live births respectively, as of 2010, reduced from 163 and 175/1000, respectively, estimated in 1990. Senegal has a rate of 75, estimate of 2010 compared to 139, estimated in 1990. Infant mortality rates, estimated at 2010, are 59, 63 and 50 for Rwanda, Uganda and Senegal, respectively. More details on child survival can be found at UNICEF.
<xref ref-type="bibr" rid="R1">1</xref>
</p>
</sec>
<sec id="s3">
<title>Data and methods</title>
<sec id="s3a">
<title>Data</title>
<p>Data were analysed from three countries, Rwanda, Uganda and Senegal purposively selected because relevant census variables were available to carry out survival analyses for the first 5 years of life (under-five mortality: 5q0).</p>
<p>For Rwanda, we analysed census 2001 data, while for Uganda and Senegal we used data from the 2002 round of the census. For all countries, data analysed were limited to an approximate 1% random sample of the census data, mainly due the complexity of the models fitted. Total samples considered were 44 000 for Rwanda, 33 500 for Senegal and 33 400 for Uganda. While the censuses have limited numbers of variables, each child record, derived from self-reported information given by the household head, consisted of age at time of death, and whether a child was alive or dead at the time of census, as well as other covariates which may influence child mortality. Our analysis was restricted to children below the age of 5 years. Since no information was available as to whether the child was alive at the census prior to the current enumeration, the survival information was right-censored. Factors influencing child mortality varied from country to country, and questions were not uniform across the three countries. To enable comparability of the results we selected similar covariates. We considered the following individual socio-demographic variables in the analysis as determinants of child mortality: region and place of residence, education level of mother, occupation of father, number of under-five children in household, whether the previous child died, whether the father or mother was alive, and ownership of dwelling unit. We constructed two indices: (1) shelter, for the following characteristics: specifically the type of dwelling unit (permanent, semipermanent, traditional), type of roof, wall and floor and type of toilet; (2) electronics for the following electronic assets: radio, cell phone, television, phone, iron and fridge. For spatial analysis, we used provinces and districts as units of analysis.
<xref ref-type="table" rid="BMJOPEN2012001421TB1">Table 1</xref>
gives a summary of the variables used.</p>
<table-wrap id="BMJOPEN2012001421TB1" position="float">
<label>Table 1</label>
<caption>
<p>Summary of selected covariates used in the model</p>
</caption>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="char" char="(" span="1"></col>
<col align="char" char="(" span="1"></col>
<col align="char" char="(" span="1"></col>
</colgroup>
<thead valign="bottom">
<tr>
<th align="left" rowspan="1" colspan="1">Covariate</th>
<th align="left" rowspan="1" colspan="1">Rwanda</th>
<th align="left" rowspan="1" colspan="1">Senegal</th>
<th align="left" rowspan="1" colspan="1">Uganda</th>
</tr>
</thead>
<tbody>
<tr>
<td colspan="4" rowspan="1">Residence</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Rural</td>
<td rowspan="1" colspan="1">38 154 (8.1)</td>
<td rowspan="1" colspan="1">23 490 (10.6)</td>
<td rowspan="1" colspan="1">29 410 (4.4)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Urban</td>
<td rowspan="1" colspan="1">6215 (5.5)</td>
<td rowspan="1" colspan="1">10 036 (6.5)</td>
<td rowspan="1" colspan="1">4028 (3.3)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Religion</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Christian</td>
<td rowspan="1" colspan="1">40 650 (7.9)</td>
<td rowspan="1" colspan="1">977 (6.9)</td>
<td rowspan="1" colspan="1">28 567 (4.4)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Muslim</td>
<td rowspan="1" colspan="1">837 (8.2)</td>
<td rowspan="1" colspan="1">32 378 (9.4)</td>
<td rowspan="1" colspan="1">4167 (4.0)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Others</td>
<td rowspan="1" colspan="1">4369 (7.1)</td>
<td rowspan="1" colspan="1">171 (10.5)</td>
<td rowspan="1" colspan="1">689 (2.7)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Age of the child</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> <1 year</td>
<td rowspan="1" colspan="1">10 516 (16.1)</td>
<td rowspan="1" colspan="1">4806 (15.5)</td>
<td rowspan="1" colspan="1">9440 (10.1)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> 1 year</td>
<td rowspan="1" colspan="1">7710 (11.1)</td>
<td rowspan="1" colspan="1">5276 (16.1)</td>
<td rowspan="1" colspan="1">9307 (1.8)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> 2 years</td>
<td rowspan="1" colspan="1">7853 (5.0)</td>
<td rowspan="1" colspan="1">6278 (11.5)</td>
<td rowspan="1" colspan="1">6937 (0.8)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> 3 years</td>
<td rowspan="1" colspan="1">6128 (3.7)</td>
<td rowspan="1" colspan="1">5628 (6.4)</td>
<td rowspan="1" colspan="1">3613 (0.4)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> 4 years</td>
<td rowspan="1" colspan="1">5701 (2.6)</td>
<td rowspan="1" colspan="1">5780 (4.5)</td>
<td rowspan="1" colspan="1">2303 (0.5)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> 5 years</td>
<td rowspan="1" colspan="1">6461 (2.1)</td>
<td rowspan="1" colspan="1">5811 (3.4)</td>
<td rowspan="1" colspan="1">1553 (9.8)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Sex of child</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Male</td>
<td rowspan="1" colspan="1">22 189 (8.5)</td>
<td rowspan="1" colspan="1">16 722 (8.2)</td>
<td rowspan="1" colspan="1">17 056 (4.4)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Female</td>
<td rowspan="1" colspan="1">22 146 (6.9)</td>
<td rowspan="1" colspan="1">16 793 (10.2)</td>
<td rowspan="1" colspan="1">16 382 (4.2)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Electronic Index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Least</td>
<td rowspan="1" colspan="1">762 (3.0)</td>
<td rowspan="1" colspan="1">10 004 (9.8)</td>
<td rowspan="1" colspan="1">2373 (2.3)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Less</td>
<td rowspan="1" colspan="1">22 250 (8.7)</td>
<td rowspan="1" colspan="1">3261 (10.8)</td>
<td rowspan="1" colspan="1">15 881 (4.9)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">619 (2.9)</td>
<td rowspan="1" colspan="1">8789 (11.2)</td>
<td rowspan="1" colspan="1">173 (1.7)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> More</td>
<td rowspan="1" colspan="1">19 249 (7.3)</td>
<td rowspan="1" colspan="1">4753 (8.5)</td>
<td rowspan="1" colspan="1">14 104 (4.1)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Most</td>
<td rowspan="1" colspan="1">1489 (5.0)</td>
<td rowspan="1" colspan="1">6719 (6.2)</td>
<td rowspan="1" colspan="1">907 (3.0)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Shelter index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Lowest</td>
<td rowspan="1" colspan="1">8168 (7.0)</td>
<td rowspan="1" colspan="1">6688 (5.4)</td>
<td rowspan="1" colspan="1">6644 (3.2)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Low</td>
<td rowspan="1" colspan="1">10 904 (7.3)</td>
<td rowspan="1" colspan="1">6785 (9.2)</td>
<td rowspan="1" colspan="1">6647 (3.9)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">7706 (8.3)</td>
<td rowspan="1" colspan="1">6743 (10.5)</td>
<td rowspan="1" colspan="1">7461 (4.0)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> High</td>
<td rowspan="1" colspan="1">8434 (8.9)</td>
<td rowspan="1" colspan="1">6162 (10.9)</td>
<td rowspan="1" colspan="1">5571 (5.7)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Highest</td>
<td rowspan="1" colspan="1">9157 (7.5)</td>
<td rowspan="1" colspan="1">7248 (10.8)</td>
<td rowspan="1" colspan="1">7169 (4.9)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Mother alive</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">39 687 (1.0)</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">25 501 (4.1)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1247 (10.1)</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">7937 (5.1)</td>
</tr>
<tr>
<td colspan="4" rowspan="1">Father alive</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">37 133 (0.7)</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">19 637 (4.0)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">3200 (8.5)</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">13 801 (4.7)</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Given in the table are the counts (and proportion dead) across covariates, N(%).</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3b">
<title>Statistical analysis</title>
<p>We examined spatial variation in under-five mortality with a flexible geo-additive semiparametric mixed model while simultaneously controlling for spatial dependence and possibly non-linear effects of covariates within a simultaneous, coherent regression framework. Individual data records were constructed for children in each country.</p>
<p>A more general spatial Cox regression model, according to Hennerfeind
<italic>et al</italic>
<xref ref-type="bibr" rid="R20">20</xref>
which captures a wide range of issues including spatial frailties was adapted and applied to determine factors associated with the risk of early childhood mortality. We applied a fully Bayesian approach based on Markov priors and using Markov Chain Monte Carlo (MCMC) techniques for inference and model checking. For model choice, we used the Deviance Information Criterion (DIC) developed as a measure of fit and model complexity.</p>
<p>The analysis was carried out using V.1.4 of the BayesX software package,
<xref ref-type="bibr" rid="R30">30</xref>
which permits Bayesian inference based on MCMC simulation techniques. For all models, 25 000 iterations were run with the initial 5 000 discarded and every 20th sample stored to give a final sample of 1 000 for parameter estimation. Convergence was evaluated by inspecting trace and autocorrelation plots of samples for each chain, as well as through numerical summaries such as the
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq1.jpg"></inline-graphic>
</inline-formula>
diagnostic statistic of Brooks and Gelman.
<xref ref-type="bibr" rid="R31">31</xref>
After 5000 iterations, all parameters showed signs of convergence in the trace plots. The values of
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq2.jpg"></inline-graphic>
</inline-formula>
also quickly approached 1 and were all below the value of 1.12, which indicated convergence of both pooled and within-interval widths to stability. Statistical methods have also been discussed in more detail in the appendix.</p>
</sec>
</sec>
<sec sec-type="results" id="s4">
<title>Results</title>
<p>
<xref ref-type="table" rid="BMJOPEN2012001421TB1">Table 1</xref>
gives a summary of the selected covariates across the three countries. There are evident disparities by place of residence for all three countries, with rural children slightly disadvantaged in mortality. The same picture was observed by age, with children less than 1 year disadvantaged compared to older children, with the proportion dying diminishing with increasing age. Children without a living mother or father were likely to die in their first 5 years of life. However, there was no clear pattern in relation to the shelter or electronics indices, or with religion or sex of the child. Similar results were obtained in the bivariate analyses presented in
<xref ref-type="table" rid="BMJOPEN2012001421TB3">tables 3</xref>
,
<xref ref-type="table" rid="BMJOPEN2012001421TB4">4</xref>
,
<xref ref-type="table" rid="BMJOPEN2012001421TB5">5</xref>
.</p>
<table-wrap id="BMJOPEN2012001421TB3" position="float">
<label>Table 3</label>
<caption>
<p>Fixed effects for Rwanda child survival</p>
</caption>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" char="." span="1"></col>
</colgroup>
<thead valign="bottom">
<tr>
<th align="left" rowspan="1" colspan="1">Variable</th>
<th align="left" rowspan="1" colspan="1">Bivariate analysis</th>
<th align="left" rowspan="1" colspan="1">Multivariate analysis</th>
</tr>
<tr>
<th rowspan="1" colspan="1"></th>
<th align="left" rowspan="1" colspan="1">HR (95% CI)</th>
<th align="left" rowspan="1" colspan="1">HR (95% CI)</th>
</tr>
</thead>
<tbody>
<tr>
<td rowspan="1" colspan="1">Intercept</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">0.04 (0.02 to 0.09)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Place of residence</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Urban</td>
<td rowspan="1" colspan="1">0.67 (0.60 to 0.75)</td>
<td rowspan="1" colspan="1">0.79 (0.73 to 0.83)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Rural</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Dwelling ownership</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">0.99 (0.89 to 1.11)</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr>
<td colspan="3" rowspan="1">Religion</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Christian</td>
<td rowspan="1" colspan="1">1.19 (1.03 to 1.38)</td>
<td rowspan="1" colspan="1">1.14 (1.05 to 1.27)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Muslim</td>
<td rowspan="1" colspan="1">1.28 (0.97 to 1.68)</td>
<td rowspan="1" colspan="1">0.95 (0.81 to 1.12)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Others</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Sex of child</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Male</td>
<td rowspan="1" colspan="1">1.20 (1.12 to 1.28)</td>
<td rowspan="1" colspan="1">1.06 (1.02 to 1.08)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Female</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Electronic index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Least</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Less</td>
<td rowspan="1" colspan="1">2.96 (1.96 to 4.46)</td>
<td rowspan="1" colspan="1">1.29 (1.16 to 1.51)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">0.74 (0.38 to 1.44)</td>
<td rowspan="1" colspan="1">0.64 (0.46 to 0.89)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> More</td>
<td rowspan="1" colspan="1">2.45 (1.62 to 3.69)</td>
<td rowspan="1" colspan="1">1.31 (1.16 to 1.51)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Most</td>
<td rowspan="1" colspan="1">1.64 (1.03 to 2.62)</td>
<td rowspan="1" colspan="1">1.15 (0.99 to 1.35)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Shelter index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Lowest</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Low</td>
<td rowspan="1" colspan="1">1.06 (0.95 to 1.18)</td>
<td rowspan="1" colspan="1">0.96 (0.88 to 1.01)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">1.20 (1.07 to 1.34)</td>
<td rowspan="1" colspan="1">0.98 (0.92 to 1.05)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> High</td>
<td rowspan="1" colspan="1">1.28 (1.15 to 1.43)</td>
<td rowspan="1" colspan="1">1.05 (0.98 to 1.13)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Highest</td>
<td rowspan="1" colspan="1">1.09 (0.98 to 1.22)</td>
<td rowspan="1" colspan="1">1.09 (1.03 to 1.19)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Mother alive</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">0.012 (0.011 to 0.013)</td>
<td rowspan="1" colspan="1">0.16 (0.15 to 0.17)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Father alive</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">0.016 (0.014 to 0.018)</td>
<td rowspan="1" colspan="1">0.38 (0.36 to 0.39)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="BMJOPEN2012001421TB4" position="float">
<label>Table 4</label>
<caption>
<p>Fixed effects for Senegal child survival</p>
</caption>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="left" span="1"></col>
<col align="left" span="1"></col>
</colgroup>
<thead valign="bottom">
<tr>
<th rowspan="1" colspan="1"></th>
<th align="left" rowspan="1" colspan="1">Bivariate analysis</th>
<th align="left" rowspan="1" colspan="1">Multivariate analysis</th>
</tr>
<tr>
<th align="left" rowspan="1" colspan="1">Variable</th>
<th align="left" rowspan="1" colspan="1">HR (95% CI)</th>
<th align="left" rowspan="1" colspan="1">HR (95% CI)</th>
</tr>
</thead>
<tbody>
<tr>
<td rowspan="1" colspan="1">Intercept</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">0.024 (0.015 to 0.045)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Place of residence</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Urban</td>
<td rowspan="1" colspan="1">0.63 (0.58 to 0.69)</td>
<td rowspan="1" colspan="1">1.01 (0.97 to 1.05)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Rural</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">0</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Dwelling ownership</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">1.55 (1.40 to 1.73)</td>
<td rowspan="1" colspan="1">1.19 (1.09 to 1.28)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">0</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Religion</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Christian</td>
<td rowspan="1" colspan="1">0.61 (0.36 to 1.03)</td>
<td rowspan="1" colspan="1">0.84 (0.56 to 1.18)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Muslim</td>
<td rowspan="1" colspan="1">0.85 90.53 to 1.35)</td>
<td rowspan="1" colspan="1">0.99 (0.71 to 1.42)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Others</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Sex of child</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Male</td>
<td rowspan="1" colspan="1">0.78 (0.73 to 0.84)</td>
<td rowspan="1" colspan="1">0.88 (0.85 to 0.91)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Female</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Electronic index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Least</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Less</td>
<td rowspan="1" colspan="1">1.12 (0.99 to 1.27)</td>
<td rowspan="1" colspan="1">1.04 (0.99 to 1.12)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">1.16 (1.06 to 1.27)</td>
<td rowspan="1" colspan="1">1.02 (0.98 to 1.06)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> More</td>
<td rowspan="1" colspan="1">0.90 (0.80 to 1.01)</td>
<td rowspan="1" colspan="1">1.001 (0.94 to 1.05)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Most</td>
<td rowspan="1" colspan="1">0.66 (0.59 to 0.73)</td>
<td rowspan="1" colspan="1">0.93 (0.87 to 1.01)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Shelter index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Lowest</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Low</td>
<td rowspan="1" colspan="1">1.63 (1.43 to 1.86)</td>
<td rowspan="1" colspan="1">1.01 (0.95 to 1.07)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">1.96 (1.69 to 2.18)</td>
<td rowspan="1" colspan="1">1.01 (0.95 to 1.06)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> High</td>
<td rowspan="1" colspan="1">1.94 (1.71 to 2.21)</td>
<td rowspan="1" colspan="1">0.99 (0.93 to 1.04)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Highest</td>
<td rowspan="1" colspan="1">1.95 (1.72 to 2.21)</td>
<td rowspan="1" colspan="1">0.96 (0.96 to 1.09)</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="BMJOPEN2012001421TB5" position="float">
<label>Table 5</label>
<caption>
<p>Fixed effects for Uganda child survival</p>
</caption>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" char="." span="1"></col>
</colgroup>
<thead valign="bottom">
<tr>
<th rowspan="1" colspan="1"></th>
<th align="left" rowspan="1" colspan="1">Bivariate analysis</th>
<th align="left" rowspan="1" colspan="1">Multivariate analysis</th>
</tr>
<tr>
<th align="left" rowspan="1" colspan="1">Variable</th>
<th align="left" rowspan="1" colspan="1">HRs (95% CI)</th>
<th align="left" rowspan="1" colspan="1">HRs (95% CI)</th>
</tr>
</thead>
<tbody>
<tr>
<td rowspan="1" colspan="1">Intercept</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">0.011 (0.006 to 0.017)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Place of residence</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Urban</td>
<td rowspan="1" colspan="1">0.68 (0.57 to 0.81)</td>
<td rowspan="1" colspan="1">1.04 (0.92 to 1.17)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Rural</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Employed</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">0.94 (0.85 to 1.05)</td>
<td rowspan="1" colspan="1">1.01 (0.91 to 1.13)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr>
<td colspan="3" rowspan="1">Under 5</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> None</td>
<td rowspan="1" colspan="1">3.74 (1.41 to 10.40)</td>
<td rowspan="1" colspan="1">4.18 (2.92 to 5.81)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> 1–3 children</td>
<td rowspan="1" colspan="1">0.92 (0.33 to 2.41)</td>
<td rowspan="1" colspan="1">0.54 (0.39 to 0.73)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> ≥4</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Married</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">1.14 (1.04 to 1.23)</td>
<td rowspan="1" colspan="1">2.10 (1.78 to 2.42)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Polygamy</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Yes</td>
<td rowspan="1" colspan="1">1.003 (0.84 to 1.19)</td>
<td rowspan="1" colspan="1">1.02 (0.95 to 1.12)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> No</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Sex of last birth</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Male</td>
<td rowspan="1" colspan="1">0.98 (0.89 to 1.08)</td>
<td rowspan="1" colspan="1">1.02 (0.92 to 1.17)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Female</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Education</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> None</td>
<td rowspan="1" colspan="1">1.26 (1.10 to 1.45)</td>
<td rowspan="1" colspan="1">1.48 (1.28 to 1.70)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Lower primary</td>
<td rowspan="1" colspan="1">1.02 (0.89 to 1.16)</td>
<td rowspan="1" colspan="1">1.21 (1.05 to 1.45)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Upper primary</td>
<td rowspan="1" colspan="1">0.66 (0.54 to 0.81)</td>
<td rowspan="1" colspan="1">0.84 (0.70 to 0.95)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Secondary</td>
<td rowspan="1" colspan="1">0.38 (0.32 to 0.62)</td>
<td rowspan="1" colspan="1">0.57 (0.36 to 0.81)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Tertiary</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Electronic index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Least</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Less</td>
<td rowspan="1" colspan="1">2.36 (1.75 to 3.10)</td>
<td rowspan="1" colspan="1">1.58 (1.15 to 2.16)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">0.84 (0.36 to 2.65)</td>
<td rowspan="1" colspan="1">0.53 (0.16 to 1.27)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> More</td>
<td rowspan="1" colspan="1">1.99 (1.51 to 2.64)</td>
<td rowspan="1" colspan="1">1.36 (1.01 to 1.90)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Most</td>
<td rowspan="1" colspan="1">1.27 (0.80 to 2.16)</td>
<td rowspan="1" colspan="1">1.14 (0.84 to 1.68)</td>
</tr>
<tr>
<td colspan="3" rowspan="1">Shelter index</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Lowest</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1">1.00</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Low</td>
<td rowspan="1" colspan="1">1.27 (1.06 to 1.52)</td>
<td rowspan="1" colspan="1">0.92 (0.83 to 1.03)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Medium</td>
<td rowspan="1" colspan="1">1.31 (1.09 to 1.53)</td>
<td rowspan="1" colspan="1">0.88 (0.81 to 0.99)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> High</td>
<td rowspan="1" colspan="1">1.97 (1.67 to 2.36)</td>
<td rowspan="1" colspan="1">1.32 (1.15 to 1.51)</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> Highest</td>
<td rowspan="1" colspan="1">1.68 (1.42 to 1.99)</td>
<td rowspan="1" colspan="1">1.12 (0.97 to 1.23)</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In
<xref ref-type="table" rid="BMJOPEN2012001421TB2">table 2</xref>
, model selection values are given for the discrete-time survival models with different specifications of the covariates for the three countries. For all the three datasets the models which combined fixed and random effects were better than those that did not combine effects, indicating the importance of both sets of factors at explaining child survival. For Rwanda data, the best model was model M3b, which combines fixed effects at individual and household levels and random effects at district and provincial levels. The DIC for model M3b was 7756.1 compared to the nearest model, M1b, with DIC=10240.3. Moving to Senegal data, again the model that combined fixed and random effects produced the best fit (model M3b). Model M3b had a DIC=14711.1 which is smaller compared to model M3a (DIC=15361.5;
<xref ref-type="table" rid="BMJOPEN2012001421TB2">table 2</xref>
). Similar results are obtained for the Uganda data with model M3b emerging as best fit, although model M2d was indistinguishable (see
<xref ref-type="table" rid="BMJOPEN2012001421TB2">table 2</xref>
).</p>
<p>
<xref ref-type="table" rid="BMJOPEN2012001421TB3">Tables 3</xref>
,
<xref ref-type="table" rid="BMJOPEN2012001421TB4">4</xref>
,
<xref ref-type="table" rid="BMJOPEN2012001421TB5">5</xref>
present estimates of fixed risk factors resulting from the models with the best fit. For Rwanda (
<xref ref-type="table" rid="BMJOPEN2012001421TB3">table 3</xref>
), there was an overall decrease of risk of a child dying in the first 5 years of life (HR=0.04, 95% CI 0.02 to 0.09). Children in urban areas were less likely to die than those in rural areas (HR=0.79, 95% CI 0.73 to 0.83). The relationship of child dying and household electronic assets was nonlinear. At level 2 compared to level 1, the risk was higher with HR=1.29 (95% CI 1.16 to 1.51), while at level 3 we observed a lower risk with HR=0.64 (95% CI 0.46 to 0.89) and this is reversed at level 4 with HR=1.31, 95% CI: 1.16 to 1.51. For the shelter index, the risk was reduced at lower levels and increased at higher levels of the index, although this relationship was not significant at p<0.05. It is interesting to note that a child with a living mother and father had a reduced risk of dying (
<xref ref-type="table" rid="BMJOPEN2012001421TB3">table 3</xref>
). Children up to 1 year of age were at increased hazard relative to those aged 5 years or older. At less than 1 year of age the log hazard was 0.77 (95% CI: 0.71 to 0.83), while at 1 year the log hazard was 0.48 (95%CI 0.42 to 0.55). As age increased, the hazard reduced. For example, those aged 2–4 years the log hazard was −0.02, −0.22 and −0.46, respectively.</p>
<table-wrap id="BMJOPEN2012001421TB2" position="float">
<label>Table 2</label>
<caption>
<p>Model comparison values based on Deviance Information Criterion (DIC) for the models</p>
</caption>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="left" span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" char="." span="1"></col>
<col align="char" char="." span="1"></col>
</colgroup>
<thead valign="bottom">
<tr>
<th align="left" rowspan="1" colspan="1">Model</th>
<th align="left" rowspan="1" colspan="1">Description</th>
<th align="left" rowspan="1" colspan="1">
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq53.jpg"></inline-graphic>
</inline-formula>
</th>
<th rowspan="1" colspan="1">
<italic>P</italic>
<sub>D</sub>
</th>
<th align="left" rowspan="1" colspan="1">DIC</th>
</tr>
</thead>
<tbody>
<tr>
<td colspan="5" rowspan="1">Rwanda</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M0</td>
<td rowspan="1" colspan="1">Province (RR)</td>
<td rowspan="1" colspan="1">24 060.2</td>
<td rowspan="1" colspan="1">12.4</td>
<td rowspan="1" colspan="1">24 084.9</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M1a</td>
<td rowspan="1" colspan="1">Fixed effects only</td>
<td rowspan="1" colspan="1">10 256.7</td>
<td rowspan="1" colspan="1">21.9</td>
<td rowspan="1" colspan="1">10 300.7</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M1b</td>
<td rowspan="1" colspan="1">Fixed+ Province (RE)</td>
<td rowspan="1" colspan="1">10 184.8</td>
<td rowspan="1" colspan="1">29.3</td>
<td rowspan="1" colspan="1">10 240.3</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M2a</td>
<td rowspan="1" colspan="1">Unstructured random effects (District)</td>
<td rowspan="1" colspan="1">19 911.4</td>
<td rowspan="1" colspan="1">22.4</td>
<td rowspan="1" colspan="1">19 956.2</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M2b</td>
<td rowspan="1" colspan="1">Structured spatial effects (District) only</td>
<td rowspan="1" colspan="1">19 212.3</td>
<td rowspan="1" colspan="1">25.7</td>
<td rowspan="1" colspan="1">19 963.5</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M3a</td>
<td rowspan="1" colspan="1">Structured effects (District)+ Unstructured (Province)</td>
<td rowspan="1" colspan="1">19 870.3</td>
<td rowspan="1" colspan="1">24.9</td>
<td rowspan="1" colspan="1">19 920.3</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M3b</td>
<td rowspan="1" colspan="1">Fixed 1 structured effects (District) Unstructured (Province)</td>
<td rowspan="1" colspan="1">7667.2</td>
<td rowspan="1" colspan="1">44.9</td>
<td rowspan="1" colspan="1">7756.1</td>
</tr>
<tr>
<td colspan="5" rowspan="1">Senegal</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M0</td>
<td rowspan="1" colspan="1">Province (RE)</td>
<td rowspan="1" colspan="1">20 536.7</td>
<td rowspan="1" colspan="1">10.4</td>
<td rowspan="1" colspan="1">20 557.5</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M1a</td>
<td rowspan="1" colspan="1">Fixed effects only</td>
<td rowspan="1" colspan="1">19 361.8</td>
<td rowspan="1" colspan="1">28.4</td>
<td rowspan="1" colspan="1">19 418.7</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M1b</td>
<td rowspan="1" colspan="1">Fixed+Province (RE)</td>
<td rowspan="1" colspan="1">19 362.7</td>
<td rowspan="1" colspan="1">28.4</td>
<td rowspan="1" colspan="1">19 419.1</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M2a</td>
<td rowspan="1" colspan="1">Unstructured random effects (District)</td>
<td rowspan="1" colspan="1">15 284.5</td>
<td rowspan="1" colspan="1">54.2</td>
<td rowspan="1" colspan="1">15 356.8</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M3a</td>
<td rowspan="1" colspan="1">Structured effects (District)+ unstructured (Province)</td>
<td rowspan="1" colspan="1">15 284.6</td>
<td rowspan="1" colspan="1">56.4</td>
<td rowspan="1" colspan="1">15 361.5</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M3b</td>
<td rowspan="1" colspan="1">Fixed+structured effects (District)—unstructured (Province)</td>
<td rowspan="1" colspan="1">14 574.1.6</td>
<td rowspan="1" colspan="1">68.6</td>
<td rowspan="1" colspan="1">14 711.1</td>
</tr>
<tr>
<td colspan="5" rowspan="1">Uganda</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M0</td>
<td rowspan="1" colspan="1">Province (RE)</td>
<td rowspan="1" colspan="1">11 681.1</td>
<td rowspan="1" colspan="1">63.3</td>
<td rowspan="1" colspan="1">11 807.5</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M1a</td>
<td rowspan="1" colspan="1">Fixed effects only</td>
<td rowspan="1" colspan="1">10 613.0</td>
<td rowspan="1" colspan="1">20.1</td>
<td rowspan="1" colspan="1">10 653.1</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M1b</td>
<td rowspan="1" colspan="1">Fixed 1 Province (RE)</td>
<td rowspan="1" colspan="1">10 425.9</td>
<td rowspan="1" colspan="1">76.8</td>
<td rowspan="1" colspan="1">10 606.6</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M2a</td>
<td rowspan="1" colspan="1">Unstructured random effects (District)</td>
<td rowspan="1" colspan="1">11 735.7</td>
<td rowspan="1" colspan="1">30.0</td>
<td rowspan="1" colspan="1">11 795.7</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M2d</td>
<td rowspan="1" colspan="1">Fixed effects+struclured spatial effects (District)</td>
<td rowspan="1" colspan="1">10 484.6</td>
<td rowspan="1" colspan="1">49.5</td>
<td rowspan="1" colspan="1">10 588.5</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M3a</td>
<td rowspan="1" colspan="1">Structured effects (District)+ unstructured (Province)</td>
<td rowspan="1" colspan="1">11 691.7</td>
<td rowspan="1" colspan="1">48.5</td>
<td rowspan="1" colspan="1">11 788.7</td>
</tr>
<tr>
<td rowspan="1" colspan="1"> M3b</td>
<td rowspan="1" colspan="1">Fixed+structured effects (District)—unstructured (Province)</td>
<td rowspan="1" colspan="1">10 468.9</td>
<td rowspan="1" colspan="1">59.6</td>
<td rowspan="1" colspan="1">10 588.1</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The spatial variability of risk of dying is shown in
<xref ref-type="fig" rid="BMJOPEN2012001421F1">figure 1</xref>
, with log hazard ranging between −8.23 and 3.14. There were a number of areas that were associated with increased risk of death compared to the overall mean. These areas are identified by the right map, with a white colour and appear in the south, west and at the centre of the country. There are also areas of reduced risk shown by a black colour.</p>
<fig id="BMJOPEN2012001421F1" position="float">
<label>Figure 1</label>
<caption>
<p>Left: structured spatial effects, at district level in Rwanda, of child surviva (model M3b). Shown are the posterior means. Right: corresponding posterior probabilities at 80% nominal level, white denotes region regions with strictly positive credible intervals, black denotes regions with strictly positive credible intervals and grey depicts regions of non-significant effects.</p>
</caption>
<graphic xlink:href="bmjopen2012001421f01"></graphic>
</fig>
<p>In
<xref ref-type="table" rid="BMJOPEN2012001421TB4">table 4</xref>
we present results for Senegal. Overall the risk of death decreases with HR=0.024 (95% CI: 0.015 to 0.045). The risk significantly varied with ownership of dwelling unit, electronic assets, sex and age of the child. Ownership of a dwelling unit was associated with increased risk (HR=1.19, 95%CI 1.09 to 1.28) compared to those households without a dwelling unit. Male children were more likely to survive the first 5 years compared to female children (HR=0.88, 95% CI 0.85 to 0.91). The risk of dying was positively associated with all ages, however, this risk decreased with age, ranging from 2.48 at age less that 1 year to 1.13 at age of 4 compared to those aged 5 years or more. For ownership of electronic assets, the risk was higher for those at the lowest level (level 1) and decreased with increasing electronic assets, although the relationship was marginally significant at p<0.1 for levels 2, 3 and 4(results not shown). Nevertheless, the results were significant, at p<0.05 for the level 5 category when compared with those at level 1 (HR=0.93, 95% CI 0.87 to 0.99). Turning to the spatial distribution of risk in
<xref ref-type="fig" rid="BMJOPEN2012001421F2">figure 2</xref>
, there was substantial variation, with estimates of log hazard ranging from −3.44 to 5.87 (left map). The right map defined areas associated with significantly high risk (shaded white) as well as those of significantly low risk (black shading). We could not identify a clear pattern to the risk by region.</p>
<fig id="BMJOPEN2012001421F2" position="float">
<label>Figure 2</label>
<caption>
<p>Left: unstructured spatial effects, at district level in Senegal, of child survival (model M3b). Shown are the posterior means. Right: corresponding posterior probabilities at 80% nominal level, white denotes regions with strictly negative credible intervals, black denotes regions with strictly positive credible intervals and grey depicts regions of nonsignificant effects.</p>
</caption>
<graphic xlink:href="bmjopen2012001421f02"></graphic>
</fig>
<p>Results for Ugandan data are given in
<xref ref-type="table" rid="BMJOPEN2012001421TB5">table 5</xref>
. Again the overall risk of death decreases (HR=0.011, 95% CI 0.006 to 0.017). Risk factors associated with under-five mortality were identified to be number of under-five children in the household, marital status, education level of mother, ownership of electronic assets and shelter characteristics. Families with fewer under-five children predisposed children to a high mortality risk compared to those with 4 or more children (HR=4.18, 95% CI 2.92 to 5.81), while those with 1–3 children had a reduced risk (HR=0.54, 95% CI 0.39 to 0.73). Being married also appeared to confer increased risk of a child dying compared to those with children of single mothers (HR=2.10, 95% CI 1.78 to 2.42). Our results showed that education level of the mother matters when it comes to child survival. Children with mothers who had no formal education or only lower primary education were more likely to die than those with tertiary education (HR=1.48, 95% CI 1.28 to 1.70 and 1.21, 95% CI 1.05 to 1.45, respectively).For those with secondary education, the risk was lower relative to those with tertiary education (HR=0.57, 95% CI 0.36 to 0.81). In relation to electronic assets, the risk was non-linear with increasing risk at level 2, reduced risk at level 3 and increased risk again at levels 4 and 5, compared to level 1 (
<xref ref-type="table" rid="BMJOPEN2012001421TB5">table 5</xref>
). We observed that the risk was lower at levels 2 and 3 of the shelter index and increased at level 4 and 5 relative to level 1. Nevertheless, the only significant difference was observed at levels 3 and 4 (1.32 and 1.12, respectively).</p>
<p>The geographical variation in risk is shown in
<xref ref-type="fig" rid="BMJOPEN2012001421F3">figure 3</xref>
. Estimates ranged from −0.61 (low risk) to 0.73 (high risk). See left plot. However, the significance map (right map) indicates that areas of high risk are in the south-west and north-west while those of low risk are in the north-east and centre-east. Notably Kampala district showed a significantly reduced risk.</p>
<fig id="BMJOPEN2012001421F3" position="float">
<label>Figure 3</label>
<caption>
<p>Left: structured spatial effects, at district level in Uganda, of child survival (model M3b). Shown are the posterior means. Right: corresponding posterior probabilities at 80% nominal level, white denotes regions with strictly negative credible intervals, black denotes regions with strictly positive credible intervals and grey depicts regions of non-significant effects.</p>
</caption>
<graphic xlink:href="bmjopen2012001421f03"></graphic>
</fig>
<p>The unstructured spatial effects at provincial level were also fitted.
<xref ref-type="fig" rid="BMJOPEN2012001421F4">Figure 4</xref>
shows caterpillar plots for the three countries at province and county level. No single province or county residual was significantly above or below zero indicating no difference in risk of death between provinces or counties in the three countries. However, there was clear variation in the risk of death, for example, in Rwanda there are four provinces with an estimated lower risk of death while six provinces have an estimated risk in the higher direction. For Senegal, there were four provinces with a reduced risk, and eight with estimated high risk. In Uganda, about a 100 counties were estimated to have a lower risk of child mortality, while another 70 had a high risk (
<xref ref-type="fig" rid="BMJOPEN2012001421F4">figure 4</xref>
).</p>
<fig id="BMJOPEN2012001421F4" position="float">
<label>Figure 4</label>
<caption>
<p>Left: unstructured spatial effects, at province level in Rwanda and Senegal and at county level in Uganda, of child survival (model M3b). Shown are the posterior means and corresponding error bars at probabilities at 80% nominal level.</p>
</caption>
<graphic xlink:href="bmjopen2012001421f04"></graphic>
</fig>
</sec>
<sec id="s5">
<title>Discussion and conclusion</title>
<p>The central question of this study was to identify risk factors associated with child mortality, which go beyond individual factors, and extending to include other factors such as the geographic location. These factors were assumed to be best captured by assuming spatially varying processes. In doing so, we applied a novel Bayesian framework which permitted estimation of risk at individual, household and area level in a unified framework.</p>
<p>Our modelling approach can be considered as an extension of the generalised linear model and can be classified as a spatial generalised linear model. These types of models have a complex structure, which is easily exploited using the Bayesian approach.</p>
<p>Despite the complexity of our approach the results obtained from our approach are consistent with what has been reported previously. For instance, the decreased risk of mortality by age is well established in SSA.
<xref ref-type="bibr" rid="R6">6</xref>
<xref ref-type="bibr" rid="R16">16</xref>
However, we found that the degree of association of various factors varied by country. We also observed a rural-urban divide in under-five mortality, with rural children more likely to die in the first 5 years than their urban counterparts.</p>
<p>The significance of spatial effects is they have an important influence on child survival in the three countries. These spatial effects may represent many factors, and are likely surrogates of factors not captured by the census survey instruments. They may include distal factors such as access to health care, availability of health care centres, reproductive health behaviour,
<xref ref-type="bibr" rid="R32">32</xref>
cultural and religious practices (including nutrition habits) specific to certain areas which may either benefit children or put them at increased risk of mortality, or they may represent other factors such as disease prevalence or cost or quality. Understanding geographical variability of mortality is an increasingly important research approach.
<xref ref-type="bibr" rid="R17">17</xref>
However, this has often been done implicitly using categorical variables to measure geographic effects.
<xref ref-type="bibr" rid="R31">31</xref>
In our approach, we explicitly introduced spatial effects and modelled them at the local level.</p>
<p>In our analyses of spatial frailties, we have used the traditional approach that extends hierarchical exchangeable frailties to incorporate spatial autocorrelation and heterogeneity. However, recent developments in the analysis of spatially correlated survival data mean that a multivariate spatial dependence structure has been proposed.
<xref ref-type="bibr" rid="R33">33</xref>
This is ideal for multiple spatially dependent data, or data that arise in spatially arranged settings.
<xref ref-type="bibr" rid="R23">23</xref>
<xref ref-type="bibr" rid="R34">34</xref>
Therefore, a conventional conditional autoregressive (CAR) model may not be ideal to capture spatial frailties. Models designed through multiple memberships multiple classification,
<xref ref-type="bibr" rid="R34">34</xref>
or using a mixture of Polya trees,
<xref ref-type="bibr" rid="R35">35</xref>
or multivariate CAR
<xref ref-type="bibr" rid="R23">23</xref>
offer desirable properties. Indeed, further studies which consider the multilevel and multivariate structure of survival data are worth exploring.
<xref ref-type="bibr" rid="R33">33</xref>
<xref ref-type="bibr" rid="R36">36</xref>
<xref ref-type="bibr" rid="R37">37</xref>
</p>
</sec>
<sec id="s6">
<title>Strength and limitations</title>
<p>It should be noted that these samples from census are 10 times bigger than those used in national surveys like DHS. However, there are some limitations in the present study that deserve attention. First, census data are cross-sectional in nature and the present study does not allow establishing temporality and thus causality of the observed associations. Given the self-reporting of children death in census, we cannot disregard the likelihood that mortality outcomes may be influenced by the respondent's recall. In addition, there was limited or lack of information for variables such as income and diseases data in census compared to surveys data, which are relevant to mortality. Nevertheless, our findings corroborate the notion that childhood mortality is an increasing public health issue in these countries, with evidence of considerable spatial variation across different provinces in the three countries.</p>
<p>Another important issue in the use of this data is the issue of data quality because of the fact that national census in developing countries are prone to incomplete or partial reporting of responses. Moreover, the use of complex questionnaires inevitably allows scope for inconsistent responses to be recorded for different questions resulting in a further complication in the assessment of morality outcomes.</p>
<p>In summary, the primary objective of this article was to illustrate a novel application of a recently developed structured additive regression model to analyse census data in SSA. The approach is data driven and has rarely been applied in the SSA region. Census data is often not used, thus this paper exemplifies the use of such techniques in census data as opposed to survey samples, with results providing confirmatory data, and therefore, providing more confidence in censuses and national surveys. Indeed, the present paper offers incremental new information regarding child survival at the individual or geographical levels in Africa. Furthermore, the results emphasise the fact that there are complex social and demographic processes operating in under-five mortality which can be more clearly understood using adequate statistical modelling which analyses the outcome of mortality beyond the individual child's risk factors and which incorporates distal factors such as the area where the child lives at the time of the survey. This has implications for these countries in terms of policy and planning for the achievement of Millennium Development Goal (MDG 4) to reduce under-five mortality by half by 2015.</p>
</sec>
<sec sec-type="supplementary-material">
<title>Supplementary Material</title>
<supplementary-material id="d34e143" content-type="local-data">
<caption>
<title>Author's manuscript</title>
</caption>
<media mimetype="application" mime-subtype="pdf" xlink:href="bmjopen-2012-001421.draft_revisions.pdf"></media>
</supplementary-material>
<supplementary-material id="d34e144" content-type="local-data">
<caption>
<title>Reviewer comments</title>
</caption>
<media mimetype="application" mime-subtype="pdf" xlink:href="bmjopen-2012-001421.reviewer_comments.pdf"></media>
</supplementary-material>
</sec>
</body>
<back>
<fn-group>
<fn>
<p>
<bold>Contributors:</bold>
KL was involved in conception and design, literature review, data analysis and interpretation, drafting the article, critical revisions for important intellectual content and approval of final article for submission. CA was involved in interpretation of results, critical revisions for important intellectual content and approval of final article for submission. KN-B was involved in conception and design, interpretation of results, critical revisions for important intellectual content and approval of final article for submission.</p>
</fn>
<fn>
<p>
<bold>Competing interests:</bold>
None.</p>
</fn>
<fn>
<p>
<bold>Ethics approval:</bold>
National Statistical Office of each country.</p>
</fn>
<fn>
<p>
<bold>Provenance and peer review:</bold>
Not commissioned; externally peer reviewed.</p>
</fn>
<fn>
<p>
<bold>Data sharing statement:</bold>
This study uses census and household data available free of charge to all.</p>
</fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="R1">
<label>1</label>
<mixed-citation publication-type="other">
<collab>UNICEF.</collab>
<article-title>The State of the World's children, 2010</article-title>
.
<comment>Unicef Annual Report</comment>
</mixed-citation>
</ref>
<ref id="R2">
<label>2</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Becher</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Muller</surname>
<given-names>O</given-names>
</name>
<name>
<surname>Jahn</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Risk factors of infant and child mortality in rural Burkina Faso</article-title>
.
<source>Bull World Health Organ</source>
<year>2004</year>
;
<volume>82</volume>
:
<fpage>265</fpage>
<lpage>73</lpage>
<pub-id pub-id-type="pmid">15259255</pub-id>
</mixed-citation>
</ref>
<ref id="R3">
<label>3</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Armstrong-Schellenberg</surname>
<given-names>JRM</given-names>
</name>
<name>
<surname>Nathan</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Abdulla</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Risk factors of child mortality in rural Tanzania</article-title>
.
<source>Trop Med Int Health</source>
<year>2002</year>
;
<volume>7</volume>
:
<fpage>506</fpage>
<lpage>11</lpage>
<pub-id pub-id-type="pmid">12031072</pub-id>
</mixed-citation>
</ref>
<ref id="R4">
<label>4</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Binka</surname>
<given-names>FN</given-names>
</name>
<name>
<surname>Indome</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Smith</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Impact of spatial distribution of permethrin-impregnated bednets on childhood mortality in rural northern Ghana</article-title>
.
<source>Am J Tropical Med Hygiene</source>
<year>1998</year>
;
<volume>59</volume>
:
<fpage>80</fpage>
<lpage>5</lpage>
</mixed-citation>
</ref>
<ref id="R5">
<label>5</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Silverman</surname>
<given-names>JG</given-names>
</name>
<name>
<surname>Decker</surname>
<given-names>MR</given-names>
</name>
<name>
<surname>Wirth</surname>
<given-names>K</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Gender-based disparities in infant and child mortality based on maternal exposure to spousal violence: the heavy burden borne by Indian girls</article-title>
.
<source>Archiv Pediatric Adolesc Med</source>
<year>2011</year>
;
<volume>165</volume>
:
<fpage>22</fpage>
<lpage>7</lpage>
</mixed-citation>
</ref>
<ref id="R6">
<label>6</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kazembe</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Mpeketula</surname>
<given-names>PMG</given-names>
</name>
</person-group>
<article-title>Quantifying spatial disparities in neonatal mortality using a structured additive regression model</article-title>
.
<source>PLoS One</source>
<year>2010</year>
;
<volume>5</volume>
:
<fpage>e11180</fpage>
<pub-id pub-id-type="pmid">20567519</pub-id>
</mixed-citation>
</ref>
<ref id="R7">
<label>7</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Johnson</surname>
<given-names>HL</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Fischer-Walker</surname>
<given-names>C</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Estimating the distribution of cuases of death among children age 1–59 months in high mortality countries with incomplete death certification</article-title>
.
<source>Int J Epidemiol</source>
<year>2010</year>
;
<volume>39</volume>
:
<fpage>1103</fpage>
<lpage>14</lpage>
<pub-id pub-id-type="pmid">20519334</pub-id>
</mixed-citation>
</ref>
<ref id="R8">
<label>8</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Magadi</surname>
<given-names>MA</given-names>
</name>
</person-group>
<article-title>Household and community HIV/AIDS status and child malnutrition in sub-Saharan Africa: evidence from the demographic and health surveys</article-title>
.
<source>Soc Sc Med</source>
<year>2011</year>
;
<volume>73</volume>
:
<fpage>436</fpage>
<lpage>46</lpage>
<pub-id pub-id-type="pmid">21729821</pub-id>
</mixed-citation>
</ref>
<ref id="R9">
<label>9</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kandala</surname>
<given-names>NB</given-names>
</name>
<name>
<surname>Fahrmeir</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Klasen</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Geo-additive models of childhood undernutrition in three sub-Saharan African countries</article-title>
.
<source>Popul Space Place</source>
<year>2009</year>
;
<volume>15</volume>
:
<fpage>119</fpage>
<lpage>26</lpage>
</mixed-citation>
</ref>
<ref id="R10">
<label>10</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>de Sherbinin</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>The biophysical and geographical correlates of child malnutrition in Africa</article-title>
.
<source>Popul Space Place</source>
<year>2011</year>
;
<volume>17</volume>
:
<fpage>27</fpage>
<lpage>46</lpage>
</mixed-citation>
</ref>
<ref id="R11">
<label>11</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Determinants of child mortality in LDCs: empirical findings from demographic and health surveys</article-title>
.
<source>Health Policy</source>
<year>2003</year>
;
<volume>65</volume>
:
<fpage>277</fpage>
<lpage>99</lpage>
<pub-id pub-id-type="pmid">12941495</pub-id>
</mixed-citation>
</ref>
<ref id="R12">
<label>12</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Haines</surname>
<given-names>MR</given-names>
</name>
<name>
<surname>Preston</surname>
<given-names>SH</given-names>
</name>
</person-group>
<article-title>The use of the census to estimate childhood mortality: comparisons from the 1900 and 1910 United States census public use samples</article-title>
.
<source>Hist Methods</source>
<year>1997</year>
;
<volume>30</volume>
:
<fpage>2</fpage>
</mixed-citation>
</ref>
<ref id="R13">
<label>13</label>
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Box-Steffensmeier</surname>
<given-names>JM</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>BS</given-names>
</name>
</person-group>
<source>Event history modeling: a guide for social scientists</source>
.
<publisher-name>Cambridge, UK: Cambridge University Press</publisher-name>
,
<year>2004</year>
</mixed-citation>
</ref>
<ref id="R14">
<label>14</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cox</surname>
<given-names>DR</given-names>
</name>
</person-group>
<article-title>Regression models and life-tables</article-title>
.
<source>JR Stat Soc Series B (Methodol)</source>
<year>1972</year>
;
<volume>34</volume>
:
<fpage>187</fpage>
<lpage>220</lpage>
</mixed-citation>
</ref>
<ref id="R15">
<label>15</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sastry</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>A nested frailty model for survival data, with an application to the study of child survival in Northeast Brazil</article-title>
.
<source>J Am Stat Assoc</source>
<year>1997</year>
;
<volume>92</volume>
:
<fpage>426</fpage>
<lpage>35</lpage>
<pub-id pub-id-type="pmid">12155403</pub-id>
</mixed-citation>
</ref>
<ref id="R16">
<label>16</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolstad</surname>
<given-names>WM</given-names>
</name>
<name>
<surname>Manda</surname>
<given-names>SO</given-names>
</name>
</person-group>
<article-title>Investigating child mortality in Malawi using family and community random effects: a Bayesian analysis</article-title>
.
<source>J Am Stat Assoc</source>
<year>2001</year>
;
<volume>96</volume>
:
<fpage>12</fpage>
<lpage>19</lpage>
</mixed-citation>
</ref>
<ref id="R17">
<label>17</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kandala</surname>
<given-names>NB</given-names>
</name>
<name>
<surname>Ghilagaber</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>A geo-additive Bayesian discrete-time survival model and its application to spatial analysis of childhood mortality in Malawi. Quality & quantity–international</article-title>
<source>.
<italic>J Methodol</italic>
</source>
<year>2006</year>
;
<volume>40</volume>
:
<fpage>935</fpage>
<lpage>57</lpage>
</mixed-citation>
</ref>
<ref id="R18">
<label>18</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Adebayo</surname>
<given-names>SB</given-names>
</name>
<name>
<surname>Fahrmeir</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Analysing child mortality in Nigeria with geoadditive discrete-time survival models</article-title>
.
<source>Stat Med</source>
<year>2005</year>
;
<volume>24</volume>
:
<fpage>709</fpage>
<lpage>28</lpage>
<pub-id pub-id-type="pmid">15696506</pub-id>
</mixed-citation>
</ref>
<ref id="R19">
<label>19</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caselli</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Pozzi</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Vaupel</surname>
<given-names>JW</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Family clustering in Sardinian longevity: a genealogical approach</article-title>
.
<source>Exp Gerontol</source>
<year>2006</year>
;
<volume>41</volume>
:
<fpage>727</fpage>
<lpage>36</lpage>
<pub-id pub-id-type="pmid">16793232</pub-id>
</mixed-citation>
</ref>
<ref id="R20">
<label>20</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hennerfeind</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Brezger</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Fahrmeir</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Geoadditive survival models</article-title>
.
<source>J Am Stat Assoc</source>
<year>2006</year>
;
<volume>101</volume>
:
<fpage>1065</fpage>
<lpage>75</lpage>
</mixed-citation>
</ref>
<ref id="R21">
<label>21</label>
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Hanagal</surname>
<given-names>DD</given-names>
</name>
</person-group>
<source>Modeling survival data using frailty models</source>
.
<publisher-name>Florida, USA: Chapman and Hall/CRC</publisher-name>
,
<year>2011</year>
</mixed-citation>
</ref>
<ref id="R22">
<label>22</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vaupel</surname>
<given-names>JW</given-names>
</name>
<name>
<surname>Manton</surname>
<given-names>KG</given-names>
</name>
<name>
<surname>Stallard</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>The impact of heterogeneity in individual frailty on the dynamics of mortality</article-title>
.
<source>Demography</source>
<year>1979</year>
;
<volume>16</volume>
:
<fpage>439</fpage>
<lpage>54</lpage>
<pub-id pub-id-type="pmid">510638</pub-id>
</mixed-citation>
</ref>
<ref id="R23">
<label>23</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Banerjee</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Wall</surname>
<given-names>MM</given-names>
</name>
<name>
<surname>Carlin</surname>
<given-names>BP</given-names>
</name>
</person-group>
<article-title>Frailty modeling for spatially correlated survival data with application to infant mortality in Minnesota</article-title>
.
<source>Biostatistics</source>
<year>2003</year>
;
<volume>4</volume>
:
<fpage>123</fpage>
<lpage>42</lpage>
<pub-id pub-id-type="pmid">12925334</pub-id>
</mixed-citation>
</ref>
<ref id="R24">
<label>24</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bastos</surname>
<given-names>LS</given-names>
</name>
<name>
<surname>Gamerman</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Dynamic survival models with spatial frailty</article-title>
.
<source>Lifetime Data Anal</source>
<year>2006</year>
;
<volume>12</volume>
:
<fpage>441</fpage>
<lpage>60</lpage>
<pub-id pub-id-type="pmid">17031498</pub-id>
</mixed-citation>
</ref>
<ref id="R25">
<label>25</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Ryan</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Modeling spatial survival data using the semi-parametric frailty models</article-title>
.
<source>Biometrics</source>
<year>2002</year>
;
<volume>58</volume>
:
<fpage>287</fpage>
<lpage>97</lpage>
<pub-id pub-id-type="pmid">12071401</pub-id>
</mixed-citation>
</ref>
<ref id="R26">
<label>26</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gemperli</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Vounatsou</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Kleinschmidt</surname>
<given-names>I</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Spatial patterns of infant mortality in Mali: the effect of malaria endemicity.</article-title>
<source>Am J Epidemiol</source>
<year>2004</year>
;
<volume>159</volume>
:
<fpage>64</fpage>
<lpage>72</lpage>
<pub-id pub-id-type="pmid">14693661</pub-id>
</mixed-citation>
</ref>
<ref id="R27">
<label>27</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Escarams</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Carrasco</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Aponte</surname>
<given-names>JJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Spatio-temporal analysis of mortality among children under the age of five in Manhia (Mozambique) during the period 1997–2005</article-title>
.
<source>Int J Health Geogr</source>
<year>2011</year>
;
<volume>10</volume>
:
<fpage>14</fpage>
<pub-id pub-id-type="pmid">21332980</pub-id>
</mixed-citation>
</ref>
<ref id="R28">
<label>28</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kamman</surname>
<given-names>EE</given-names>
</name>
<name>
<surname>Wand</surname>
<given-names>MP</given-names>
</name>
</person-group>
<article-title>Geoadditive models</article-title>
.
<source>J R Statist Soc B</source>
<year>2003</year>
;
<volume>53</volume>
:
<fpage>118</fpage>
</mixed-citation>
</ref>
<ref id="R29">
<label>29</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kneib</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Fahrmeir</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>A mixed model approach for geoadditive hazard regression</article-title>
.
<source>Scandinavian J Statist</source>
<year>2007</year>
;
<volume>34</volume>
:
<fpage>207228</fpage>
</mixed-citation>
</ref>
<ref id="R30">
<label>30</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brezger</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Kneib</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Lang</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Software for Bayesian inference based on Markov Chain Monte Carlo simulation techniques</article-title>
.
<source>J Statist Software</source>
<year>2005</year>
;
<volume>14</volume>
:
<fpage>11</fpage>
</mixed-citation>
</ref>
<ref id="R31">
<label>31</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brooks</surname>
<given-names>SP</given-names>
</name>
<name>
<surname>Gelman</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>General methods for monitoring convergence of iterative simulations</article-title>
.
<source>J Comp Graph Statist</source>
<year>1998</year>
;
<volume>7</volume>
:
<fpage>434</fpage>
<lpage>55</lpage>
</mixed-citation>
</ref>
<ref id="R32">
<label>32</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gould</surname>
<given-names>JB</given-names>
</name>
<name>
<surname>Herrchen</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Pham</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Small-area analysis: targeting high-risk areas for adolescent pregnancy prevention programs</article-title>
.
<source>Fam Plan Perspect</source>
<year>1998</year>
;
<volume>30</volume>
:
<fpage>173176</fpage>
</mixed-citation>
</ref>
<ref id="R33">
<label>33</label>
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Manda</surname>
<given-names>SOM</given-names>
</name>
<name>
<surname>Feltbower</surname>
<given-names>RG</given-names>
</name>
<name>
<surname>Gilhorpe</surname>
<given-names>MS</given-names>
</name>
</person-group>
<article-title>A multivariate random effects model for multiple spatially dependent survival data</article-title>
. In:
<person-group person-group-type="editor">
<name>
<surname>Tu</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Greenwood</surname>
<given-names>D</given-names>
</name>
</person-group>
, eds.
<source>Modern methods for epidemiology</source>
.
<publisher-loc>Dordrecht</publisher-loc>
:
<publisher-name>Springer</publisher-name>
,
<year>2012</year>
:
<fpage>157</fpage>
<lpage>72</lpage>
</mixed-citation>
</ref>
<ref id="R34">
<label>34</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Browne</surname>
<given-names>WJ</given-names>
</name>
<name>
<surname>Goldstein</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Rasbash</surname>
<given-names>J</given-names>
</name>
</person-group>
.
<article-title>Multiple membership multiple classification models</article-title>
.
<source>Stat Model</source>
<year>2001</year>
;
<volume>1</volume>
:
<fpage>103</fpage>
<lpage>24</lpage>
</mixed-citation>
</ref>
<ref id="R35">
<label>35</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Hanson</surname>
<given-names>TE</given-names>
</name>
</person-group>
<article-title>Spatially dependent Polya tree modeling of survival data</article-title>
.
<source>Biometrics</source>
<year>2011</year>
;
<volume>67</volume>
:
<fpage>391</fpage>
<lpage>403</lpage>
<pub-id pub-id-type="pmid">20731644</pub-id>
</mixed-citation>
</ref>
<ref id="R36">
<label>36</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Carlin</surname>
<given-names>BP</given-names>
</name>
<name>
<surname>Banerjee</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Wall</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Frailty modelling for spatially correlated survival data with application to infant mortality in Minnesota</article-title>
.
<source>Biostatistics</source>
<year>2003</year>
;
<volume>4</volume>
:
<fpage>123</fpage>
<lpage>42</lpage>
<pub-id pub-id-type="pmid">12925334</pub-id>
</mixed-citation>
</ref>
<ref id="R37">
<label>37</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gastafson</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Large hierarchical Bayesian analysis of multivariate survival data</article-title>
.
<source>Biometrics</source>
<year>1997</year>
;
<volume>53</volume>
:
<fpage>230</fpage>
<lpage>42</lpage>
<pub-id pub-id-type="pmid">9147593</pub-id>
</mixed-citation>
</ref>
<ref id="R38">
<label>38</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Spiegelhalter</surname>
<given-names>DJ</given-names>
</name>
<name>
<surname>Best</surname>
<given-names>NG</given-names>
</name>
<name>
<surname>Carlin</surname>
<given-names>BP</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Bayesian measures of model complexity and fit (with discussion)</article-title>
.
<source>J R Stat Soc B</source>
<year>2002</year>
;
<volume>64</volume>
:
<fpage>1</fpage>
<lpage>34</lpage>
</mixed-citation>
</ref>
</ref-list>
<app-group>
<app id="app1">
<title>Appendix</title>
<sec id="s7">
<title>Statistical methodology</title>
<p>In studying survival of the child, we assume T as time to event or survival time with t as the actual realisation. The probability that a survival time T is less than or equal to some value t is given by
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq3.jpg"></inline-graphic>
</inline-formula>
. In our context, F(t) is the cumulative probability that a child dies at or before some given time t, thus F(t)=P(Child dies at time ≤t). The instantaneous probability that an event will occur at any given small interval is defined as F’(t)=f(t). The proportion of children surviving to time t or beyond is derived as S(t)=1−F(t), which is also known as the survivor function. An important approach is to consider the duration analysis through the hazard rate. The hazard rate, which links the survival and failure functions, is of the form
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq4.jpg"></inline-graphic>
</inline-formula>
or equivalent to
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq5.jpg"></inline-graphic>
</inline-formula>
. The hazard rate, unlike the survivor function, describes the risk or event of ‘failure’ given that the individual has survived all along up to point t.</p>
<p>In the analysis of child mortality, our interest is to answer this question: Given that the child has survived up to month t, what is the likelihood it will survive in the subsequent months? Further to this, one is interested in how the hazard rate varies with respect to some covariates, for instance, will the hazard be the same for children living in urban and rural areas? One way to analyse such data is to use Kaplan–Meier survival curves and log-rank test. This is an exploratory analysis that permits assessment of any differences in child survival by various covariates.</p>
<p>An alternative model which captures the effect of covariates is to use Cox regression models or commonly referred to as the proportional hazard model. It should be pointed out, however, that various statistical models may be constructed, see Box-Steffensmeier and Jones
<xref ref-type="bibr" rid="R13">13</xref>
for an overview on the topic. In contrast to other data, several issues are considered when analysing survival data. Central are censoring and truncation of survival data, existence of time-varying covariates, occurrence of multiple causes of death, whether events occurrences were recorded in discrete-time, and the possibility of group-risk factors and confounders acting on the hazard. Thus, a more general model that incorporates all these issues if they are present in the data is needed.</p>
<p>We propose using a more general Cox model that captures a wide range of issues including spatial frailties. Thus, a spatial Cox regression model
<xref ref-type="bibr" rid="R20">20</xref>
was applied to determine factors associated with the risk of early childhood mortality. Assume that
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq6.jpg"></inline-graphic>
</inline-formula>
is the observed number of months lived or the censoring time for jth child in area i. Under Cox's model, the hazard function at time T=t is given by
<disp-formula id="BMJOPEN-2012-001421M1">
<label>1</label>
<graphic xlink:href="bmjopen2012001421eq1"></graphic>
</disp-formula>
where
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq7.jpg"></inline-graphic>
</inline-formula>
is the baseline hazard at time t, and the βs are a vector of regression coefficients for the fixed and time-invariant variables
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq8.jpg"></inline-graphic>
</inline-formula>
. The exponent of a coefficient, that is, exp(β), is interpreted as HR, that is, the ratio of instantaneous risks which is assumed to be constant over time. The HR compares rates of deaths in one group to some reference group, for a categorical variable, and to the mean for a continuous variable.</p>
<p>Since individuals are clustered in geographical regions, group-specific random frailty term,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq9.jpg"></inline-graphic>
</inline-formula>
, was introduced to augment the Cox model, that is,
<disp-formula id="BMJOPEN-2012-001421M2">
<label>2</label>
<graphic xlink:href="bmjopen2012001421eq2"></graphic>
</disp-formula>
</p>
<p>The above model indicated that childhood survival was influenced by both individual-specific factors
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq10.jpg"></inline-graphic>
</inline-formula>
and group-specific environmental factors
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq11.jpg"></inline-graphic>
</inline-formula>
. Here it was assumed that the environmental factors were approximated by geographical locations. In the case of geographical regions, spatially distributed random effects
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq12.jpg"></inline-graphic>
</inline-formula>
were assumed, while for the other unstructured heterogeneity a random effect,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq13.jpg"></inline-graphic>
</inline-formula>
, was specified such that
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq14.jpg"></inline-graphic>
</inline-formula>
. Fitting model (2) assumed a semiparametric additive predictor, which is known as geoadditive survival model,
<xref ref-type="bibr" rid="R20">20</xref>
<disp-formula id="BMJOPEN-2012-001421M3">
<label>3</label>
<graphic xlink:href="bmjopen2012001421eq3"></graphic>
</disp-formula>
where
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq15.jpg"></inline-graphic>
</inline-formula>
is the log-additive predictor at time t for child j in area i. The term
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq16.jpg"></inline-graphic>
</inline-formula>
is the log baseline hazard effect at time t. The other terms are as defined above.</p>
</sec>
<sec id="s8">
<title>Estimation: fully Bayesian approach</title>
<sec id="s8a">
<title>Prior distributions for covariate effects</title>
<p>Modelling and inference uses the fully Bayesian approach. In the Bayesian formulation, the specification of the proposed model (
<xref ref-type="disp-formula" rid="BMJOPEN-2012-001421M3">equation 3</xref>
) is complete by assigning priors to all unknown parameters. For the fixed regression parameters, a suitable choice is the diffuse prior, that is, p(γ)∝const, but a weakly informative Gaussian prior is also possible. The baseline hazard effect,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq17.jpg"></inline-graphic>
</inline-formula>
, was assigned a penalised spline with a second order random walk prior. Similarly, the time and continuous covariates were estimated non-parametrically through smoothness priors. We use the second-order Gaussian random walk prior to allow enough flexibility, while penalising abrupt changes in the function, as suggested by Brezger
<italic>et al.</italic>
<xref ref-type="bibr" rid="R30">30</xref>
The prior can be expressed in the pairwise difference form as
<disp-formula id="BMJOPEN-2012-001421M4">
<label>4</label>
<graphic xlink:href="bmjopen2012001421eq4"></graphic>
</disp-formula>
where
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq18.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq19.jpg"></inline-graphic>
</inline-formula>
is the variance, with diffuse priors
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq20.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq21.jpg"></inline-graphic>
</inline-formula>
for initial values.</p>
<p>For the unstructured spatial heterogeneity term,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq22.jpg"></inline-graphic>
</inline-formula>
is assumed to follow an exchangeable Gaussian prior with zero mean and variance,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq23.jpg"></inline-graphic>
</inline-formula>
, that is,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq24.jpg"></inline-graphic>
</inline-formula>
. Finally, for the spatial components
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq25.jpg"></inline-graphic>
</inline-formula>
, we assign a MRF prior.
<xref ref-type="bibr" rid="R30">30</xref>
This is analogous to random walk models. The conditional distribution of
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq26.jpg"></inline-graphic>
</inline-formula>
, given adjacent areas
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq27.jpg"></inline-graphic>
</inline-formula>
, is a univariate normal distribution with mean equal the average
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq28.jpg"></inline-graphic>
</inline-formula>
values of
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq29.jpg"></inline-graphic>
</inline-formula>
's neighbouring areas and variance equal to
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq30.jpg"></inline-graphic>
</inline-formula>
divided by the number of adjacent areas. This leads to a joint density of the form
<disp-formula id="BMJOPEN-2012-001421M5">
<label>5</label>
<graphic xlink:href="bmjopen2012001421eq5"></graphic>
</disp-formula>
where i∼j denotes that area i is adjacent to j, and assumes that parameter values
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq31.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq32.jpg"></inline-graphic>
</inline-formula>
in adjacent areas are similar. The degree of similarity is determined by the unknown precision parameter
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq33.jpg"></inline-graphic>
</inline-formula>
.</p>
<p>By writing
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq34.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq135.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq35.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq36.jpg"></inline-graphic>
</inline-formula>
, for a well-defined design matrix Z and a (possibly high-dimensional) vector of regression parameters β, all different priors (
<xref ref-type="disp-formula" rid="BMJOPEN-2012-001421M4">equations 4</xref>
and
<xref ref-type="disp-formula" rid="BMJOPEN-2012-001421M5">5</xref>
) can be expressed in a general Gaussian form
<disp-formula id="BMJOPEN-2012-001421M6">
<label>6</label>
<graphic xlink:href="bmjopen2012001421eq6"></graphic>
</disp-formula>
with an appropriate penalty matrix
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq37.jpg"></inline-graphic>
</inline-formula>
. Its structure depends on the covariate and smoothness of the function. In most cases,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq38.jpg"></inline-graphic>
</inline-formula>
is rank deficient and hence the prior for
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq39.jpg"></inline-graphic>
</inline-formula>
is improper. For the variances
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq40.jpg"></inline-graphic>
</inline-formula>
we a
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq41.jpg"></inline-graphic>
</inline-formula>
assume inverse Gamma priors
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq42.jpg"></inline-graphic>
</inline-formula>
, with hyperparameters
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq43.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq44.jpg"></inline-graphic>
</inline-formula>
chosen such that this prior is weakly informative.</p>
</sec>
<sec id="s8b">
<title>Posterior distribution</title>
<p>Fully Bayesian inference is based on the analysis of posterior distribution of the model parameters. In general, the posterior is highly dimensional and analytically intractable, which makes direct inference almost impossible. This problem is circumvented by using MCMC simulation techniques, whereby samples are drawn from the full conditional of parameters given the rest of the data. Under conditional independence assumptions the posterior distribution for the Bernoulli model is given by Bayes Theorem</p>
<p>
<disp-formula id="BMJOPEN-2012-001421M7">
<label>7</label>
<graphic xlink:href="bmjopen2012001421eq7"></graphic>
</disp-formula>
where the quantity
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq46.jpg"></inline-graphic>
</inline-formula>
is the prior density function, and
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq47.jpg"></inline-graphic>
</inline-formula>
denotes the likelihood of the data. More specifically, the posterior is given by
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq1"></graphic>
</disp-formula>
</p>
<p>where
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq48.jpg"></inline-graphic>
</inline-formula>
is a binary indicator coded 1 if an event occurs and 0 if an event does not occur at time t. For updating the full conditionals of parameters, we use a hybrid MCMC sampling scheme of the iteratively weighted least squares proposals, developed for GLMM by Brezger,
<xref ref-type="bibr" rid="R30">30</xref>
and Metropolis-Hastings algorithm. Full details are presented elsewhere.
<xref ref-type="bibr" rid="R20">20</xref>
<xref ref-type="bibr" rid="R36">36</xref>
</p>
</sec>
</sec>
<sec id="s9">
<title>Data analysis</title>
<p>A number of models were explored. The first model (M0) explored unstructured variation in child i at provincial level k</p>
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq2"></graphic>
</disp-formula>
<p>The second set of models estimated fixed effects only (M1a) and then we adjusted for unstructured random effects at province level (M1b).
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq3"></graphic>
</disp-formula>
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq4"></graphic>
</disp-formula>
</p>
<p>We also investigated geographical variation at district level. We fitted both unstructured (M2a)and structured random effects (M2b) using districts as variables.
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq5"></graphic>
</disp-formula>
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq6"></graphic>
</disp-formula>
</p>
<p>The last set of models combined fixed and random effects at district and province levels. In model M3a we estimated structured spatial effects at district level and unstructured effects at province level, and model (M3b) improved model M3a by combining with fixed effects.
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq7"></graphic>
</disp-formula>
<disp-formula>
<graphic xlink:href="bmjopen2012001421ueq8"></graphic>
</disp-formula>
</p>
<p>Model comparison was based on the DIC.
<xref ref-type="bibr" rid="R38">38</xref>
This is given by, where
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq49.jpg"></inline-graphic>
</inline-formula>
is the deviance of the model evaluated at the posterior mean of the parameters, and represents the fit of the model to the data. The component
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq50.jpg"></inline-graphic>
</inline-formula>
is the effective number of parameters, which assessed the complexity of the model. Since small values of
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq51.jpg"></inline-graphic>
</inline-formula>
indicate good fit while small values of
<inline-formula>
<inline-graphic xlink:href="bmjopen-2012-001421ileq52.jpg"></inline-graphic>
</inline-formula>
indicate a parsimonious model, small values of DIC indicate a better model. Models with differences in DIC of <3 compared with the best model cannot be distinguished, while those between 3–7 can be weakly differentiated.
<xref ref-type="bibr" rid="R38">38</xref>
</p>
</sec>
</app>
</app-group>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/Sante/explor/SidaSubSaharaV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 002262 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 002262 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/Sante
   |area=    SidaSubSaharaV1
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     PMC:3488715
   |texte=   Childhood mortality in sub-Saharan Africa: cross-sectional insight into small-scale geographical inequalities from Census data
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:23089207" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd   \
       | NlmPubMed2Wicri -a SidaSubSaharaV1 

Wicri

This area was generated with Dilib version V0.6.32.
Data generation: Mon Nov 13 19:31:10 2017. Site generation: Wed Mar 6 19:14:32 2024