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.

A Bayesian Two Part Model Applied to Analyze Risk Factors of Adult Mortality with Application to Data from Namibia

Identifieur interne : 002A53 ( Pmc/Corpus ); précédent : 002A52; suivant : 002A54

A Bayesian Two Part Model Applied to Analyze Risk Factors of Adult Mortality with Application to Data from Namibia

Auteurs : Lawrence N. Kazembe

Source :

RBID : PMC:3774685

Abstract

Despite remarkable gains in life expectancy and declining mortality in the 21st century, in many places mostly in developing countries, adult mortality has increased in part due to HIV/AIDS or continued abject poverty levels. Moreover many factors including behavioural, socio-economic and demographic variables work simultaneously to impact on risk of mortality. Understanding risk factors of adult mortality is crucial towards designing appropriate public health interventions. In this paper we proposed a structured additive two-part random effects regression model for adult mortality data. Our proposal assumed two processes: (i) whether death occurred in the household (prevalence part), and (ii) number of reported deaths, if death did occur (severity part). The proposed model specification therefore consisted of two generalized linear mixed models (GLMM) with correlated random effects that permitted structured and unstructured spatial components at regional level. Specifically, the first part assumed a GLMM with a logistic link and the second part explored a count model following either a Poisson or negative binomial distribution. The model was used to analyse adult mortality data of 25,793 individuals from the 2006/2007 Namibian DHS data. Inference is based on the Bayesian framework with appropriate priors discussed.


Url:
DOI: 10.1371/journal.pone.0073500
PubMed: 24066052
PubMed Central: 3774685

Links to Exploration step

PMC:3774685

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">A Bayesian Two Part Model Applied to Analyze Risk Factors of Adult Mortality with Application to Data from Namibia</title>
<author>
<name sortKey="Kazembe, Lawrence N" sort="Kazembe, Lawrence N" uniqKey="Kazembe L" first="Lawrence N." last="Kazembe">Lawrence N. Kazembe</name>
<affiliation>
<nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">24066052</idno>
<idno type="pmc">3774685</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3774685</idno>
<idno type="RBID">PMC:3774685</idno>
<idno type="doi">10.1371/journal.pone.0073500</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">002A53</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">002A53</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">A Bayesian Two Part Model Applied to Analyze Risk Factors of Adult Mortality with Application to Data from Namibia</title>
<author>
<name sortKey="Kazembe, Lawrence N" sort="Kazembe, Lawrence N" uniqKey="Kazembe L" first="Lawrence N." last="Kazembe">Lawrence N. Kazembe</name>
<affiliation>
<nlm:aff id="aff1"></nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint>
<date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<p>Despite remarkable gains in life expectancy and declining mortality in the 21st century, in many places mostly in developing countries, adult mortality has increased in part due to HIV/AIDS or continued abject poverty levels. Moreover many factors including behavioural, socio-economic and demographic variables work simultaneously to impact on risk of mortality. Understanding risk factors of adult mortality is crucial towards designing appropriate public health interventions. In this paper we proposed a structured additive two-part random effects regression model for adult mortality data. Our proposal assumed two processes: (i) whether death occurred in the household (prevalence part), and (ii) number of reported deaths, if death did occur (severity part). The proposed model specification therefore consisted of two generalized linear mixed models (GLMM) with correlated random effects that permitted structured and unstructured spatial components at regional level. Specifically, the first part assumed a GLMM with a logistic link and the second part explored a count model following either a Poisson or negative binomial distribution. The model was used to analyse adult mortality data of 25,793 individuals from the 2006/2007 Namibian DHS data. Inference is based on the Bayesian framework with appropriate priors discussed.</p>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bendavid, E" uniqKey="Bendavid E">E Bendavid</name>
</author>
<author>
<name sortKey="Holmes, Cb" uniqKey="Holmes C">CB Holmes</name>
</author>
<author>
<name sortKey="Bhattacharya, I" uniqKey="Bhattacharya I">I Bhattacharya</name>
</author>
<author>
<name sortKey="Muller, G" uniqKey="Muller G">G Muller</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bendavid, E" uniqKey="Bendavid E">E Bendavid</name>
</author>
<author>
<name sortKey="Seligman, B" uniqKey="Seligman B">B Seligman</name>
</author>
<author>
<name sortKey="Kubo, J" uniqKey="Kubo J">J Kubo</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Obermeyer, Z" uniqKey="Obermeyer Z">Z Obermeyer</name>
</author>
<author>
<name sortKey="Rajaratnam, Jk" uniqKey="Rajaratnam J">JK Rajaratnam</name>
</author>
<author>
<name sortKey="Park, Ch" uniqKey="Park C">CH Park</name>
</author>
<author>
<name sortKey="Gakidou, E" uniqKey="Gakidou E">E Gakidou</name>
</author>
<author>
<name sortKey="Hogan, Mc" uniqKey="Hogan M">MC Hogan</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Ali, M" uniqKey="Ali M">M Ali</name>
</author>
<author>
<name sortKey="Jin, Y" uniqKey="Jin Y">Y Jin</name>
</author>
<author>
<name sortKey="Kim, Dr" uniqKey="Kim D">DR Kim</name>
</author>
<author>
<name sortKey="De, Zb" uniqKey="De Z">ZB De</name>
</author>
<author>
<name sortKey="Park, Jk" uniqKey="Park J">JK Park</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Turra, Cm" uniqKey="Turra C">CM Turra</name>
</author>
<author>
<name sortKey="Goldman, N" uniqKey="Goldman N">N Goldman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Fantahun, M" uniqKey="Fantahun M">M Fantahun</name>
</author>
<author>
<name sortKey="Berhane, Y" uniqKey="Berhane Y">Y Berhane</name>
</author>
<author>
<name sortKey="Hogberg, U" uniqKey="Hogberg U">U Högberg</name>
</author>
<author>
<name sortKey="Wall, S" uniqKey="Wall S">S Wall</name>
</author>
<author>
<name sortKey="Byass, P" uniqKey="Byass P">P Byass</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tucker Seeley, Rd" uniqKey="Tucker Seeley R">RD Tucker-Seeley</name>
</author>
<author>
<name sortKey="Li, Y" uniqKey="Li Y">Y Li</name>
</author>
<author>
<name sortKey="Sorensen, G" uniqKey="Sorensen G">G Sorensen</name>
</author>
<author>
<name sortKey="Subramanian, Sv" uniqKey="Subramanian S">SV Subramanian</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Neelon, B" uniqKey="Neelon B">B Neelon</name>
</author>
<author>
<name sortKey="Ghosh, P" uniqKey="Ghosh P">P Ghosh</name>
</author>
<author>
<name sortKey="Loebs, Pf" uniqKey="Loebs P">PF Loebs</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Winkelmann, R" uniqKey="Winkelmann R">R Winkelmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Olsen, Mk" uniqKey="Olsen M">MK Olsen</name>
</author>
<author>
<name sortKey="Schafer, J" uniqKey="Schafer J">J Schafer</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Su, L" uniqKey="Su L">L Su</name>
</author>
<author>
<name sortKey="Tom, Bdm" uniqKey="Tom B">BDM Tom</name>
</author>
<author>
<name sortKey="Farewell, Vt" uniqKey="Farewell V">VT Farewell</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Jin, X" uniqKey="Jin X">X Jin</name>
</author>
<author>
<name sortKey="Carlin, Bp" uniqKey="Carlin B">BP Carlin</name>
</author>
<author>
<name sortKey="Banerjee, S" uniqKey="Banerjee S">S Banerjee</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gelman, A" uniqKey="Gelman A">A Gelman</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>
<author>
<name sortKey="Van Der Linde, A" uniqKey="Van Der Linde A">A van der Linde</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Andersen, Rm" uniqKey="Andersen R">RM Andersen</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tanser, F" uniqKey="Tanser F">F Tanser</name>
</author>
<author>
<name sortKey="Hosegood, V" uniqKey="Hosegood V">V Hosegood</name>
</author>
<author>
<name sortKey="Benzler, J" uniqKey="Benzler J">J Benzler</name>
</author>
<author>
<name sortKey="Solarsh, G" uniqKey="Solarsh G">G Solarsh</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Stock, R" uniqKey="Stock R">R Stock</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Fahrmeir, L" uniqKey="Fahrmeir L">L Fahrmeir</name>
</author>
<author>
<name sortKey="Lang, S" uniqKey="Lang S">S Lang</name>
</author>
</analytic>
</biblStruct>
<biblStruct></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">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group>
<journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">24066052</article-id>
<article-id pub-id-type="pmc">3774685</article-id>
<article-id pub-id-type="publisher-id">PONE-D-13-14683</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0073500</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>A Bayesian Two Part Model Applied to Analyze Risk Factors of Adult Mortality with Application to Data from Namibia</article-title>
<alt-title alt-title-type="running-head">Bayesian Two Part Model for Adult Mortality</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Kazembe</surname>
<given-names>Lawrence N.</given-names>
</name>
<xref ref-type="aff" rid="aff1"></xref>
<xref ref-type="corresp" rid="cor1">
<sup>*</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<addr-line>Department of Statistics and Population Studies, University of Namibia, Windhoek, Namibia</addr-line>
</aff>
<contrib-group>
<contrib contrib-type="editor">
<name>
<surname>Hsiao</surname>
<given-names>Chuhsing Kate</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1">
<addr-line>National Taiwan University, Taiwan</addr-line>
</aff>
<author-notes>
<corresp id="cor1">* E-mail:
<email>lkazembe@unam.na</email>
</corresp>
<fn fn-type="conflict">
<p>
<bold>Competing Interests: </bold>
The author has declared that no competing interests exist.</p>
</fn>
<fn fn-type="con">
<p>Conceived and designed the experiments: LNK. Analyzed the data: LNK. Wrote the paper: LNK.</p>
</fn>
</author-notes>
<pub-date pub-type="collection">
<year>2013</year>
</pub-date>
<pub-date pub-type="epub">
<day>16</day>
<month>9</month>
<year>2013</year>
</pub-date>
<volume>8</volume>
<issue>9</issue>
<elocation-id>e73500</elocation-id>
<history>
<date date-type="received">
<day>10</day>
<month>4</month>
<year>2013</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>7</month>
<year>2013</year>
</date>
</history>
<permissions>
<copyright-year>2013</copyright-year>
<copyright-holder>Lawrence N</copyright-holder>
<license>
<license-p>Kazembe. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<abstract>
<p>Despite remarkable gains in life expectancy and declining mortality in the 21st century, in many places mostly in developing countries, adult mortality has increased in part due to HIV/AIDS or continued abject poverty levels. Moreover many factors including behavioural, socio-economic and demographic variables work simultaneously to impact on risk of mortality. Understanding risk factors of adult mortality is crucial towards designing appropriate public health interventions. In this paper we proposed a structured additive two-part random effects regression model for adult mortality data. Our proposal assumed two processes: (i) whether death occurred in the household (prevalence part), and (ii) number of reported deaths, if death did occur (severity part). The proposed model specification therefore consisted of two generalized linear mixed models (GLMM) with correlated random effects that permitted structured and unstructured spatial components at regional level. Specifically, the first part assumed a GLMM with a logistic link and the second part explored a count model following either a Poisson or negative binomial distribution. The model was used to analyse adult mortality data of 25,793 individuals from the 2006/2007 Namibian DHS data. Inference is based on the Bayesian framework with appropriate priors discussed.</p>
</abstract>
<funding-group>
<funding-statement>The author acknowledges financial support from the University of Namibia. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts>
<page-count count="10"></page-count>
</counts>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The improvement of health of the population has been a long-term focus in most developing countries, mostly with the bid to meet the Millennium Development Goals
<xref ref-type="bibr" rid="pone.0073500-Jamison1">[1]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Bendavid1">[2]</xref>
. African countries have channeled considerable resources aimed at boosting child and maternal health. National governments and development agencies have proposed strategies of advancing service delivery to effectively combat poor health in children and mothers through disease prevention and control
<xref ref-type="bibr" rid="pone.0073500-Jamison1">[1]</xref>
. Little focus, though, has been drawn on adult health, specifically on adult survival and mortality. In some countries adult mortality has slumped considerably, despite the remarkable gains in life expectancy and declining mortality in the 21st century
<xref ref-type="bibr" rid="pone.0073500-Bendavid2">[3]</xref>
. The strong impact of HIV on mortality, and the survival of adults has to some extent substantially affected the population structure of African communities
<xref ref-type="bibr" rid="pone.0073500-Bendavid1">[2]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Ngom1">[4]</xref>
, Ngom. The life-course effects of abject poverty have created structural spheres which have negatively impacted on the quality of life manifesting at old age, further exacerbating adult mortality. Moreover, many factors including behavioural, socio-economic and demographic variables work simultaneously to impact on risk of mortality. Understanding risk factors of adult mortality is crucial towards designing appropriate public health interventions.</p>
<p>Notwithstanding, the paucity of vital statistics in most sub-Saharan Africa makes it difficult to study patterns and trends, as well as risk factors of adult mortality. In most sub-Saharan Africa, mortality data is based on census data, however, census are many years apart to inform meaningful trends and patterns, and assess the effect of any mitigatory interventions. To fill the gap, estimates or projections from the WHO/UNDP have been used
<xref ref-type="bibr" rid="pone.0073500-Jamison1">[1]</xref>
. In some countries, adult mortality estimates are derived from demographic surveillance system (DSS), where these exist INDEPTH or through the vital registration system (VRS). Ideally, the VRS is not functional in many African countries. An alternative to this is to use an increasing cross-sectional body of data publicly available from Demographic and Health Surveys (DHS), which collects data on survivorship/widowhood/siblings. Recently, the DHS has collected data in eight African countries, including Namibia, under the module “Support For Those Who Have Died” which captured the household census mortality data
<xref ref-type="bibr" rid="pone.0073500-Ministry1">[6]</xref>
.</p>
<p>This study was aimed at understanding adult mortality by developing models that explain risk factors of adult mortality in Namibia. Following the conceptual framework as outlined in Rogers et al.
<xref ref-type="bibr" rid="pone.0073500-Rogers1">[7]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Rogers2">[8]</xref>
, we incorporated a number of relevant explanatory variables, in particular, behavioural and health factors, socio-economic and demographic variables that shape mortality levels in a community. Within this framework, behavioural and health factors include alcohol drinking and smoking, availability and access to health care. The economic and social variables encompass level of living, marriage and family characteristics, and social ties to communities e.g. religion. Demographic variables consist of age, sex and ethnicity. These variables are further grouped into individual, household and community factors. Contextual factors, which measure unobserved or unmeasured determinants, are introduced as random effects. These are assumed to have a structured and unstructured pattern, where the structured effects permit similarities across areas, while the unstructured effects allow within area heterogeneity.</p>
<p>The need for geographical analysis is critical as it may assist in social planning by highlighting areas of excess mortality. To assume homogeneity in terms of regional mortality would be an understatement. The 2010 poverty estimates in Namibia, the country of focus in this study, indicated huge disparities across regions
<xref ref-type="bibr" rid="pone.0073500-Namibia1">[9]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Kazembe1">[10]</xref>
. Moreover, the thirteen regions of Namibia remain disparate geographically, economically, culturally and socially and undoubtedly may explain the expected differentials in adult and old-age mortality (AMOA). Even so in these regions, there are considerable heterogeneity in health outcomes and socio-economic factors
<xref ref-type="bibr" rid="pone.0073500-Kazembe1">[10]</xref>
.</p>
<p>Past work on the epidemiology of adult mortality have considered either prevalence alone (i.e. whether death occurred in a household) and used logit models to analyze such
<xref ref-type="bibr" rid="pone.0073500-Obermeyer1">[11]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Houle1">[12]</xref>
, or severity only (i.e. compared the number of reported deaths) and applied count regression models
<xref ref-type="bibr" rid="pone.0073500-Ali1">[13]</xref>
<xref ref-type="bibr" rid="pone.0073500-Fantahun1">[15]</xref>
. In some cases where count data have excess zeros, zero-inflation (ZI) regression model have been used, particularly in general adult health
<xref ref-type="bibr" rid="pone.0073500-TuckerSeeley1">[16]</xref>
. Use of such ZI model in adult mortality is rare despite many examples in biostatistics, epidemiology and econometrics that deal with zero-inflated data. See for example Winkelmann
<xref ref-type="bibr" rid="pone.0073500-Winkelmann1">[17]</xref>
and references therein.</p>
<p>In this article, however, our proposal was that there are two processes: (i) whether death occurred in the household (prevalence part), and (ii) number of reported deaths, if death did occur (severity part). Wickelman
<xref ref-type="bibr" rid="pone.0073500-Winkelmann1">[17]</xref>
called such two-part processes as extensive (the zeros) and intensive (the positives) margins in a multi-index count model, representing the case that death did not occur and that death did occur, respectively. It is also referred to as a zero-hurdle model because it allows for a systematic difference in the statistical process governing individuals (observations) below the hurdle and individuals above the hurdle set at zero
<xref ref-type="bibr" rid="pone.0073500-Winkelmann1">[17]</xref>
<xref ref-type="bibr" rid="pone.0073500-Winkelmann2">[19]</xref>
. An alternative approach to two-part process is to use finite mixtures, which is a combination of zeros point mass distribution and the nonzero distribution
<xref ref-type="bibr" rid="pone.0073500-Winkelmann1">[17]</xref>
. Most work in this area, however, do not address concerns that occurrence of death (prevalence) and number of deaths (severity) reported in a household are joint processes, and that failing to account for the joint nature of these processes can bias estimates of risk factors on adult mortality
<xref ref-type="bibr" rid="pone.0073500-Olsen1">[20]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Su1">[21]</xref>
. Furthermore, several challenges arise in analyzing these data including how to: (i) handle correlation in the data between the prevalence and severity processes; (ii) model excess zero counts in the data; (iii) permit possible spatial correlation; and (iv) fit nonlinearity in metrical variables.</p>
<p>In this paper we proposed a structured additive two-part random effects model for adult mortality data to address the above issues. Therefore following Neelon et al.
<xref ref-type="bibr" rid="pone.0073500-Neelon1">[18]</xref>
, Winkelmann
<xref ref-type="bibr" rid="pone.0073500-Winkelmann2">[19]</xref>
, Olsen and Schafer
<xref ref-type="bibr" rid="pone.0073500-Olsen1">[20]</xref>
and Su et al.
<xref ref-type="bibr" rid="pone.0073500-Su1">[21]</xref>
, the proposed model specification consisted of two generalized linear mixed models (GLMM) with correlated random effects, where the first part assumed a GLMM with a logistic link and the second part explored a count model following either a Poisson or negative binomial distribution. Furthermore, the random effects in the model permitted structured and unstructured spatial components to account for spatial correlation at area level. The proposed model was then used to analyse adult mortality data of 25,793 individuals from the 2006/2007 Namibian DHS data. Inference was based on the Bayesian framework with appropriate priors discussed.</p>
<p>The remainder of this paper is organized as follows. Section 2 presents the data, as well as outlines the two-part model and its application to the Namibian DHS data. Results of the analysis are given in Section 3. This is followed by the discussion in Section 4 and we conclude.</p>
</sec>
<sec sec-type="methods" id="s2">
<title>Methods</title>
<sec id="s2a">
<title>Data</title>
<p>This study used data from the 2006/2007 Namibia DHS
<xref ref-type="bibr" rid="pone.0073500-Ministry1">[6]</xref>
, which included a household mortality module that has been implemented in 8 African countries since 2005 under the sub-theme of “Support for Persons Who Have Died”. The DHS are periodic cross-sectional health surveys funded by the U.S. Agency for International Development (USAID) Bureau for Global Health. The DHS have been conducted in over 70 countries since 1980s, and the data are publicly available at MEASURE DHS website. The data are completely anonymized before they are released for general public use. The DHS follows a cross-sectional multi-stage stratified survey design. In the Namibian context, a two-stage stratified sampling design was implemented to collect the data and provide direct estimates of demographic and health indicators at national and regional level. At first stage, a total of 500 enumeration areas (EA) from a sampling frame of 3750 EAs as defined in the 2001 Population and Housing Census, were selected stratified by urban-rural status with sampling probability proportional to the population of the region. At the second stage, a fixed number of 20 households were randomly drawn from each selected EA. From the selected households, all women of age 15–49 years were eligible for interview. A final representative probability sample of 10,000 households was selected. The sample allocation by region is given in the Appendix A of the Namibian DHS report
<xref ref-type="bibr" rid="pone.0073500-Ministry1">[6]</xref>
. The overall women response rates for urban and rural areas were 90.0% and 94.8% respectively, while for the 13 regions ranged from 88.4% (Khomas region) to 96.6% (Omusati region). The household response rate was 96.8% and 98.5% for urban and rural areas respectively, with a regional variation of 95.1% (Khomas region) to 99.1% (Oshikoto region).</p>
<p>The household module enumerated all usual household members and collected information on their age, gender and economic activities (education level, or employment status if not in school). All deaths in the household that occurred in the past 12 months preceding the survey were collected with full information on age at death and gender. The questions used to collect mortality data were: “Has any usual member of your household died in the last 12 months”, and if yes this was followed by “How many members died in the last 12 months”. Information were also collected on whether any household member was sick, whether any medical services, emotional or psychological, material or social support were available to the sick or deceased or any household member. In addition, approximate time to nearest health facility and type of health facility where health care was obtained were equally captured.</p>
</sec>
<sec id="s2b">
<title>Statistical Analysis</title>
<p>The number of household members that died was assumed to follow two process, which was better envisaged as a two-part mixture consisting of a point mass at zeros followed by truncated count data distribution for the non-zero observations. The model for the count can be either a Poisson or a negative binomial. These two-part models are also referred to as zero-hurdle models
<xref ref-type="bibr" rid="pone.0073500-Neelon1">[18]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Winkelmann2">[19]</xref>
. In particular, according to Winkelman
<xref ref-type="bibr" rid="pone.0073500-Winkelmann1">[17]</xref>
, a hurdle model combines a dichotomous model for the binary outcome of the count being below or above the hurdle (the selection variable), with a truncated model for outcomes above the hurdle. For the Poisson hurdle model we have
<disp-formula id="pone.0073500.e001">
<graphic xlink:href="pone.0073500.e001"></graphic>
<label>(1)</label>
</disp-formula>
<disp-formula id="pone.0073500.e002">
<graphic xlink:href="pone.0073500.e002"></graphic>
<label>(2)</label>
</disp-formula>
where
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e003.jpg"></inline-graphic>
</inline-formula>
indicates the response for household member
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e004.jpg"></inline-graphic>
</inline-formula>
in area
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e005.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e006.jpg"></inline-graphic>
</inline-formula>
is the mean for the truncated Poisson distribution. An alternative to the Poisson hurdle is a negative binomial which replacing
<xref ref-type="disp-formula" rid="pone.0073500.e002">equation (2</xref>
) is given by</p>
<p>
<disp-formula id="pone.0073500.e007">
<graphic xlink:href="pone.0073500.e007"></graphic>
<label>(3)</label>
</disp-formula>
<disp-formula id="pone.0073500.e008">
<graphic xlink:href="pone.0073500.e008"></graphic>
<label>(4)</label>
</disp-formula>
with parameters
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e009.jpg"></inline-graphic>
</inline-formula>
for the mean and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e010.jpg"></inline-graphic>
</inline-formula>
for over-dispersion. The component
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e011.jpg"></inline-graphic>
</inline-formula>
is a mortality propensity, while the count part models severity of mortality in the reported household. If one takes a less strict view of the two-part process then a finite mixture of two distributions can be assumed. A less strict view is assumed where zeros are generated from the same underlying process as positives, thus single-index models e.g. Poisson or negative binomial may apply. Otherwise, if the zero-generating process is not subject to such a constraint, then multi-index models are appropriate
<xref ref-type="bibr" rid="pone.0073500-Winkelmann1">[17]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Winkelmann2">[19]</xref>
. In this case the factors that explain whether death occured in a household within past 12 months may differ from the factors that determine the number of deaths following the first death in the same household.</p>
<p>A special case of the models above is the zero-inflated Poisson or zero-inflated Negative binomial. These models have a degenerate distribution at zero with untruncated Poisson or Negative binomial distribution. A zero-inflated Poisson is denoted by
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e012.jpg"></inline-graphic>
</inline-formula>
. Combining zero inflation and over-dispersion gives a zero inflated negative binomial defined as
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e013.jpg"></inline-graphic>
</inline-formula>
, where
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e014.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e015.jpg"></inline-graphic>
</inline-formula>
are the predictor and over-dispersion parameters respectively.</p>
<p>The zero-hurdle model can be extended to accommodate covariates and random effects. Since we have two parts, we introduced two GLMMs. For the prevalence component, we assumed a logit link while for the severity component we proposed a log link:
<disp-formula id="pone.0073500.e016">
<graphic xlink:href="pone.0073500.e016"></graphic>
<label>(5)</label>
</disp-formula>
where
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e017.jpg"></inline-graphic>
</inline-formula>
is the intercept for process
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e018.jpg"></inline-graphic>
</inline-formula>
, the terms
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e019.jpg"></inline-graphic>
</inline-formula>
 = (
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e020.jpg"></inline-graphic>
</inline-formula>
are vectors of regression parameters corresponding to the set of covariates,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e021.jpg"></inline-graphic>
</inline-formula>
(
<xref ref-type="table" rid="pone-0073500-t001">Table 1</xref>
). The non-linear components of continuous covariates (
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e022.jpg"></inline-graphic>
</inline-formula>
) were captured through terms
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e023.jpg"></inline-graphic>
</inline-formula>
 = (
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e024.jpg"></inline-graphic>
</inline-formula>
. The components
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e025.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e026.jpg"></inline-graphic>
</inline-formula>
are the unstructured heterogeneity and spatially structured variation terms, respectively, at regional level. Because of dependence in the binary and count outcomes, the random effects were correlated and were modelled using multivariate distributions as explained below.</p>
<table-wrap id="pone-0073500-t001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.t001</object-id>
<label>Table 1</label>
<caption>
<title>Summary of household members who died by selected characteristics.</title>
</caption>
<alternatives>
<graphic id="pone-0073500-t001-1" xlink:href="pone.0073500.t001"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Variable</td>
<td align="left" rowspan="1" colspan="1">Percentage died</td>
<td align="left" rowspan="1" colspan="1">Total</td>
<td align="left" rowspan="1" colspan="1">χ
<sup>2</sup>
test (p-value)</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>Residence</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Urban</td>
<td align="left" rowspan="1" colspan="1">5.4</td>
<td align="left" rowspan="1" colspan="1">10829</td>
<td align="left" rowspan="1" colspan="1">157.7 (<0.001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Rural</td>
<td align="left" rowspan="1" colspan="1">9.7</td>
<td align="left" rowspan="1" colspan="1">14964</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>
<sup>1</sup>
Nearest health facility</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Hospital</td>
<td align="left" rowspan="1" colspan="1">6.3</td>
<td align="left" rowspan="1" colspan="1">5531</td>
<td align="left" rowspan="1" colspan="1">54.5 (<0.001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Health centre</td>
<td align="left" rowspan="1" colspan="1">8.6</td>
<td align="left" rowspan="1" colspan="1">1880</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Clinic</td>
<td align="left" rowspan="1" colspan="1">8.5</td>
<td align="left" rowspan="1" colspan="1">17903</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>
<sup>2</sup>
Means to nearest health facilit</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Car/motorcycle</td>
<td align="left" rowspan="1" colspan="1">4.8</td>
<td align="left" rowspan="1" colspan="1">4013</td>
<td align="left" rowspan="1" colspan="1">75.6 (<0.001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Public transport/Animal cart</td>
<td align="left" rowspan="1" colspan="1">8.0</td>
<td align="left" rowspan="1" colspan="1">4877</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Walking</td>
<td align="left" rowspan="1" colspan="1">8.7</td>
<td align="left" rowspan="1" colspan="1">16021</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>
<sup>3</sup>
ime to nearest health facility</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Minutes</td>
<td align="left" rowspan="1" colspan="1">6.6</td>
<td align="left" rowspan="1" colspan="1">16730</td>
<td align="left" rowspan="1" colspan="1">135.9 (<0.001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Hours</td>
<td align="left" rowspan="1" colspan="1">10.6</td>
<td align="left" rowspan="1" colspan="1">8174</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Day</td>
<td align="left" rowspan="1" colspan="1">12.1</td>
<td align="left" rowspan="1" colspan="1">612</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>Sex of household head</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Male</td>
<td align="left" rowspan="1" colspan="1">6.0</td>
<td align="left" rowspan="1" colspan="1">15019</td>
<td align="left" rowspan="1" colspan="1">16.6 (<0.002)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Female</td>
<td align="left" rowspan="1" colspan="1">10.6</td>
<td align="left" rowspan="1" colspan="1">10774</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>Sex of household member</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Male</td>
<td align="left" rowspan="1" colspan="1">7.3</td>
<td align="left" rowspan="1" colspan="1">12020</td>
<td align="left" rowspan="1" colspan="1">12.8 (<0.001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Female</td>
<td align="left" rowspan="1" colspan="1">8.5</td>
<td align="left" rowspan="1" colspan="1">13773</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>
<sup>4</sup>
ducation of household member</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">None</td>
<td align="left" rowspan="1" colspan="1">8.5</td>
<td align="left" rowspan="1" colspan="1">3876</td>
<td align="left" rowspan="1" colspan="1">115.4 (<0.001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Primary</td>
<td align="left" rowspan="1" colspan="1">9.9</td>
<td align="left" rowspan="1" colspan="1">8230</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Secondary and higher</td>
<td align="left" rowspan="1" colspan="1">6.5</td>
<td align="left" rowspan="1" colspan="1">13283</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>
<sup>5</sup>
Marital status of household member</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Married</td>
<td align="left" rowspan="1" colspan="1">5.4</td>
<td align="left" rowspan="1" colspan="1">9464</td>
<td align="left" rowspan="1" colspan="1">203.7 (<0.0001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Never married</td>
<td align="left" rowspan="1" colspan="1">8.8</td>
<td align="left" rowspan="1" colspan="1">13796</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Others</td>
<td align="left" rowspan="1" colspan="1">13.4</td>
<td align="left" rowspan="1" colspan="1">2162</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>Wealth quintile</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Poorest</td>
<td align="left" rowspan="1" colspan="1">12.0</td>
<td align="left" rowspan="1" colspan="1">4215</td>
<td align="left" rowspan="1" colspan="1">299.2 (<0.0001)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Poor</td>
<td align="left" rowspan="1" colspan="1">9.8</td>
<td align="left" rowspan="1" colspan="1">4785</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Medium</td>
<td align="left" rowspan="1" colspan="1">8.8</td>
<td align="left" rowspan="1" colspan="1">6126</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Rich</td>
<td align="left" rowspan="1" colspan="1">6.5</td>
<td align="left" rowspan="1" colspan="1">6134</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Richest</td>
<td align="left" rowspan="1" colspan="1">2.9</td>
<td align="left" rowspan="1" colspan="1">4533</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="nt101">
<label></label>
<p>Numbers and percentage missing per variable are provided at the foot of the table.</p>
</fn>
<fn id="nt102">
<label></label>
<p>
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e027.jpg"></inline-graphic>
</inline-formula>
Missing:
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e028.jpg"></inline-graphic>
</inline-formula>
(1.9%).</p>
</fn>
<fn id="nt103">
<label></label>
<p>
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e029.jpg"></inline-graphic>
</inline-formula>
Missing:
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e030.jpg"></inline-graphic>
</inline-formula>
(3.1%).</p>
</fn>
<fn id="nt104">
<label></label>
<p>
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e031.jpg"></inline-graphic>
</inline-formula>
Missing:
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e032.jpg"></inline-graphic>
</inline-formula>
(1.1%).</p>
</fn>
<fn id="nt105">
<label></label>
<p>
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e033.jpg"></inline-graphic>
</inline-formula>
Missing:
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e034.jpg"></inline-graphic>
</inline-formula>
(1.7%).</p>
</fn>
<fn id="nt106">
<label></label>
<p>
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e035.jpg"></inline-graphic>
</inline-formula>
Missing:
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e036.jpg"></inline-graphic>
</inline-formula>
(1.3%).</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Model estimation was carried out using the Bayesian approach and the following prior distributions were specified for all parameters of the model (5). For the intercept, diffuse priors were assumed, that is,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e037.jpg"></inline-graphic>
</inline-formula>
, while for the other fixed effects,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e038.jpg"></inline-graphic>
</inline-formula>
, highly dispersed normal distribution priors were chosen, that is,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e039.jpg"></inline-graphic>
</inline-formula>
. The smooth functions of continuous covariates were modelled using a second-order random walk prior given by
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e040.jpg"></inline-graphic>
</inline-formula>
for
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e041.jpg"></inline-graphic>
</inline-formula>
with noninformative priors for the initials. Again
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e042.jpg"></inline-graphic>
</inline-formula>
controlled the amount of smoothing, with larger values leading to less smoothing. The unstructured spatial effects
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e043.jpg"></inline-graphic>
</inline-formula>
were assumed to follow a multivariate normal distribution, i.e.,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e044.jpg"></inline-graphic>
</inline-formula>
, with covariance matrix
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e045.jpg"></inline-graphic>
</inline-formula>
. The spatial structured effects
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e046.jpg"></inline-graphic>
</inline-formula>
were assigned a multivariate conditional autoregressive (MCAR) prior, i.e.,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e047.jpg"></inline-graphic>
</inline-formula>
, again
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e048.jpg"></inline-graphic>
</inline-formula>
is a covariance matrix
<xref ref-type="bibr" rid="pone.0073500-Neelon1">[18]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Jin1">[22]</xref>
.</p>
<p>The covariance matrices have their diagonal elements equal to the variances and the off-diagonals are correlation components between the outcome processes. Thus, for example
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e049.jpg"></inline-graphic>
</inline-formula>
is variance components, while
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e050.jpg"></inline-graphic>
</inline-formula>
are cross-covariance components between the prevalence of mortality (part 1) and severity of mortality (part 2). Correspondingly,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e051.jpg"></inline-graphic>
</inline-formula>
, for example, gives a measure of spatial correlation between the processes. The variance components were assigned inverse Wishart priors, i.e.,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e052.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e053.jpg"></inline-graphic>
</inline-formula>
where
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e054.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e055.jpg"></inline-graphic>
</inline-formula>
are scalars, while
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e056.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e057.jpg"></inline-graphic>
</inline-formula>
are symmetric and positive definitive matrices. The hyperpriors were assigned
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e058.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e059.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e060.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e061.jpg"></inline-graphic>
</inline-formula>
where
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e062.jpg"></inline-graphic>
</inline-formula>
is an identity matrix.</p>
<p>The regression tool for full Bayesian inference was based on the posterior distribution of all parameters. Markov Chain Monte Carlo techniques were used to draw samples from the full conditionals of all parameter distribution which were then summarized to obtain model estimates in the posterior analysis, i.e.,
<disp-formula id="pone.0073500.e063">
<graphic xlink:href="pone.0073500.e063"></graphic>
<label>(6)</label>
</disp-formula>
where
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e064.jpg"></inline-graphic>
</inline-formula>
is the joint distribution of all parameters in the observation model and the corresponding priors, and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e065.jpg"></inline-graphic>
</inline-formula>
is the likelihood for all observable data (
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e066.jpg"></inline-graphic>
</inline-formula>
). Gibbs sampler was used to draw samples from the full conditionals.</p>
<p>Model implementation was carried out in WinBUGS 1.4
<xref ref-type="bibr" rid="pone.0073500-Spiegelhalter1">[23]</xref>
. The sample code used for our bivariate problem was adapted from Neelon et al.
<xref ref-type="bibr" rid="pone.0073500-Neelon1">[18]</xref>
. See
<xref ref-type="supplementary-material" rid="pone.0073500.s001">Appendix S1</xref>
. We fitted a Poisson and negative binomial hurdle models. For each model, three separate chains were run to help assess convergence, starting from different initial values for the priors. As pointed out by Gelman
<xref ref-type="bibr" rid="pone.0073500-Gelman1">[24]</xref>
, the performance of the model is sensitive to the choice of the hyperparameters. We therefore considered alternative specifications on variance components hyperparameters,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e067.jpg"></inline-graphic>
</inline-formula>
and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e068.jpg"></inline-graphic>
</inline-formula>
, and carried out sensitivity of our model, by assuming
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e069.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e070.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e071.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e072.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e073.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e074.jpg"></inline-graphic>
</inline-formula>
, and
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e075.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e076.jpg"></inline-graphic>
</inline-formula>
 = 
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e077.jpg"></inline-graphic>
</inline-formula>
. The results were similar, therefore, the last choice was maintained. Convergence was monitored by visual examination of time series plots of samples for each chain, and confirmed by plotting the Gelman-Rubin statistic. The first 5,000–10,000 samples were discarded as a “burn-in” and then each chain was run for a further 30,000 iterations or till convergence was achieved. Models were compared using the Deviance Information Criterion (DIC: Spiegelhalter et al.
<xref ref-type="bibr" rid="pone.0073500-Spiegelhalter2">[25]</xref>
), which was simultaneously computed in the estimation process. Smaller values of DIC indicated a better fitting model.</p>
</sec>
</sec>
<sec id="s3">
<title>Results</title>
<p>
<xref ref-type="fig" rid="pone-0073500-g001">Figure 1</xref>
displays the age-specific histogram of adult and old-age deaths. As proposed the number of reported zero deaths was high consisting of 92% of the total observations, with 1 death given in about 6.6% household, while 2, 3 and 6 deaths reported per household represented 0.8%, 0.3% and 0.01% respectively.
<xref ref-type="fig" rid="pone-0073500-g002">Figure 2</xref>
shows the percentage who died by region (panel (a)) and the mean number of persons reported dead per household by region (panel (b)). Table1 presents a descriptive summary of percentage who died by socio-demographic covariates. Relatively small numbers were missing for some variables mainly due to item non-response. There was significant difference across categories. In part, there were more deaths in rural households, in female-headed households and in the poorest of the poor households, and the nearer the health facility the fewer the deaths reported.</p>
<fig id="pone-0073500-g001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.g001</object-id>
<label>Figure 1</label>
<caption>
<title>Histogram of percentage of reported adult and old-age mortality in Namibia, 2006/07 DHS.</title>
</caption>
<graphic xlink:href="pone.0073500.g001"></graphic>
</fig>
<fig id="pone-0073500-g002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.g002</object-id>
<label>Figure 2</label>
<caption>
<title>(a) Percentage of household with at least one member dead within 12</title>
<p>
<bold>months of the survey date by region; (b) Mean number of household members reported dead by region.</bold>
</p>
</caption>
<graphic xlink:href="pone.0073500.g002"></graphic>
</fig>
<p>We fitted four models, two were of random effects only based on the Poisson and Negative binomial, and the other two extended the two random effects models by including fixed and nonlinear effects. The DIC was used to for model selection. The DIC values are given in
<xref ref-type="table" rid="pone-0073500-t002">Table 2</xref>
. The basic Poisson logit hurdle had a
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e078.jpg"></inline-graphic>
</inline-formula>
(
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e079.jpg"></inline-graphic>
</inline-formula>
,
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e080.jpg"></inline-graphic>
</inline-formula>
) which was of poor fit compared to the basic NB logit hurdle (
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e081.jpg"></inline-graphic>
</inline-formula>
). Comparatively, models with fixed and nonlinear effects outperformed the basic models. However, the extended NB logit hurdle model was superior to the extended Poisson logit hurdle (
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e082.jpg"></inline-graphic>
</inline-formula>
 = 8439.92 versus
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e083.jpg"></inline-graphic>
</inline-formula>
 = 9739.95 respectively).</p>
<table-wrap id="pone-0073500-t002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.t002</object-id>
<label>Table 2</label>
<caption>
<title>Model comparison values based on Deviance Information Criterion (DIC) for the Poisson logit hurdle and negative binomial (NB) logit hurdle models.</title>
</caption>
<alternatives>
<graphic id="pone-0073500-t002-2" xlink:href="pone.0073500.t002"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Model</td>
<td align="left" rowspan="1" colspan="1">Description</td>
<td align="left" rowspan="1" colspan="1">
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e084.jpg"></inline-graphic>
</inline-formula>
</td>
<td align="left" rowspan="1" colspan="1">
<italic>p
<sub>D</sub>
</italic>
</td>
<td align="left" rowspan="1" colspan="1">
<italic>DIC</italic>
</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Poisson logit hurdle</td>
<td align="left" rowspan="1" colspan="1">Region (Random effects: RE)+ Cluster (RE)</td>
<td align="left" rowspan="1" colspan="1">9625.65</td>
<td align="left" rowspan="1" colspan="1">347.66</td>
<td align="left" rowspan="1" colspan="1">10320.97</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Poisson logit hurdle</td>
<td align="left" rowspan="1" colspan="1">Fixed + Nonlinear + Region (RE) + Cluster(RE)</td>
<td align="left" rowspan="1" colspan="1">9008.61</td>
<td align="left" rowspan="1" colspan="1">365.67</td>
<td align="left" rowspan="1" colspan="1">9739.95</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">NB logit Hurdle</td>
<td align="left" rowspan="1" colspan="1">Region (RE) + Cluster(RE)</td>
<td align="left" rowspan="1" colspan="1">7799.10</td>
<td align="left" rowspan="1" colspan="1">339.06</td>
<td align="left" rowspan="1" colspan="1">8577.21</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">NB logit Hurdle</td>
<td align="left" rowspan="1" colspan="1">Fixed + Nonlinear+ Region (RE) + Cluster(RE)</td>
<td align="left" rowspan="1" colspan="1">7741.38</td>
<td align="left" rowspan="1" colspan="1">349.27</td>
<td align="left" rowspan="1" colspan="1">8439.92</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>
<xref ref-type="table" rid="pone-0073500-t003">Table 3</xref>
presents posterior summaries from the Negative binomial hurdle model. Given are the fixed and random effects for both the logit and negative binomial components. Odds of death in the family was positively associated with female headed households (log odds of 0.48, 95% CI: 0.38, 0.59), negatively with being married (−0.48, 95% CI: −0.64, −0.32) and never married (−0.21, 95% CI: −0.38, −0.04), positively with being poorest (0.95, 95% CI: 0.60, 1.28), poor (0.81, 95% CI: 0.46, 1.11), medium (0.98, 95% CI: 0.67, 1.26) and rich (0.78, 95% CI: 0.54, 1.04) relative to being the richest, increased when one takes hours to reach a hospital (0.65, 95% CI: 0.60, 1.28), and positively if the family head had primary education (0.18, 95% CI: 0.06, 0.23).</p>
<table-wrap id="pone-0073500-t003" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.t003</object-id>
<label>Table 3</label>
<caption>
<title>Posterior means (post. mean) for Fixed and Random effects estimates with corresponding 95% credible intervals (CI) from the spatial negative binomial hurdle model of adult mortality.</title>
</caption>
<alternatives>
<graphic id="pone-0073500-t003-3" xlink:href="pone.0073500.t003"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Variable</td>
<td colspan="2" align="left" rowspan="1">Bernoulli</td>
<td colspan="2" align="left" rowspan="1">Negative binomial</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Post. Mean</td>
<td align="left" rowspan="1" colspan="1">Post. 95% CI</td>
<td align="left" rowspan="1" colspan="1">Post. Mean</td>
<td align="left" rowspan="1" colspan="1">Post. 95% CI</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>Fixed effects</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Constant</td>
<td align="left" rowspan="1" colspan="1">−4.73</td>
<td align="left" rowspan="1" colspan="1">(−5.44, −4.22)</td>
<td align="left" rowspan="1" colspan="1">−4.45</td>
<td align="left" rowspan="1" colspan="1">(−5.02, −3.91)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Urban</td>
<td align="left" rowspan="1" colspan="1">−0.23</td>
<td align="left" rowspan="1" colspan="1">(−0.54, 0.08)</td>
<td align="left" rowspan="1" colspan="1">−0.35</td>
<td align="left" rowspan="1" colspan="1">(−0.66, −0.07)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Rural</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Hospital</td>
<td align="left" rowspan="1" colspan="1">0.21</td>
<td align="left" rowspan="1" colspan="1">(−0.06, 0.61)</td>
<td align="left" rowspan="1" colspan="1">0.09</td>
<td align="left" rowspan="1" colspan="1">(−0.20, 0.41)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Clinic</td>
<td align="left" rowspan="1" colspan="1">0.29</td>
<td align="left" rowspan="1" colspan="1">(−0.01, 0.61)</td>
<td align="left" rowspan="1" colspan="1">0.15</td>
<td align="left" rowspan="1" colspan="1">(−0.11, 0.42)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Health centre</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Walk</td>
<td align="left" rowspan="1" colspan="1">0.05</td>
<td align="left" rowspan="1" colspan="1">(−0.14, 0.31)</td>
<td align="left" rowspan="1" colspan="1">0.16</td>
<td align="left" rowspan="1" colspan="1">(−0.05, 0.35)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Public transport</td>
<td align="left" rowspan="1" colspan="1">0.17</td>
<td align="left" rowspan="1" colspan="1">(−0.06, 0.38)</td>
<td align="left" rowspan="1" colspan="1">0.23</td>
<td align="left" rowspan="1" colspan="1">(0.02, 0.43)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Car/motorcycle</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Female head</td>
<td align="left" rowspan="1" colspan="1">0.48</td>
<td align="left" rowspan="1" colspan="1">(0.38, 0.59)</td>
<td align="left" rowspan="1" colspan="1">0.33</td>
<td align="left" rowspan="1" colspan="1">(0.23, 0.42)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Male head</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Female member</td>
<td align="left" rowspan="1" colspan="1">0.005</td>
<td align="left" rowspan="1" colspan="1">(−0.08, 0.16)</td>
<td align="left" rowspan="1" colspan="1">0.02</td>
<td align="left" rowspan="1" colspan="1">(−0.07, 0.09)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Male member</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Married</td>
<td align="left" rowspan="1" colspan="1">−0.48</td>
<td align="left" rowspan="1" colspan="1">(−0.64, −0.32)</td>
<td align="left" rowspan="1" colspan="1">−0.37</td>
<td align="left" rowspan="1" colspan="1">(−0.54, −0.21)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Not married</td>
<td align="left" rowspan="1" colspan="1">−0.21</td>
<td align="left" rowspan="1" colspan="1">(−0.38, −0.04)</td>
<td align="left" rowspan="1" colspan="1">−0.14</td>
<td align="left" rowspan="1" colspan="1">(−0.30, 0.02)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Other (married)</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Poorest</td>
<td align="left" rowspan="1" colspan="1">0.95</td>
<td align="left" rowspan="1" colspan="1">(0.60, 1.28)</td>
<td align="left" rowspan="1" colspan="1">0.75</td>
<td align="left" rowspan="1" colspan="1">(0.47, 1.09)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Poor</td>
<td align="left" rowspan="1" colspan="1">0.81</td>
<td align="left" rowspan="1" colspan="1">(0.46, 1.11)</td>
<td align="left" rowspan="1" colspan="1">0.75</td>
<td align="left" rowspan="1" colspan="1">(0.47, 1.05)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Medium</td>
<td align="left" rowspan="1" colspan="1">0.98</td>
<td align="left" rowspan="1" colspan="1">(0.67, 1.26)</td>
<td align="left" rowspan="1" colspan="1">0.83</td>
<td align="left" rowspan="1" colspan="1">(0.58, 1.12)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Rich</td>
<td align="left" rowspan="1" colspan="1">0.78</td>
<td align="left" rowspan="1" colspan="1">(0.54, 1.04)</td>
<td align="left" rowspan="1" colspan="1">0.73</td>
<td align="left" rowspan="1" colspan="1">(0.51, 0.99)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Richest</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Time to facility (Min)</td>
<td align="left" rowspan="1" colspan="1">0.33</td>
<td align="left" rowspan="1" colspan="1">(−0.17, 0.90)</td>
<td align="left" rowspan="1" colspan="1">0.29</td>
<td align="left" rowspan="1" colspan="1">(−0.17, 0.79)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Time to facility (Hr)</td>
<td align="left" rowspan="1" colspan="1">0.65</td>
<td align="left" rowspan="1" colspan="1">(0.14, 1.27)</td>
<td align="left" rowspan="1" colspan="1">0.62</td>
<td align="left" rowspan="1" colspan="1">(0.18, 1.12)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Time to facility (Day)</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">No education</td>
<td align="left" rowspan="1" colspan="1">−0.001</td>
<td align="left" rowspan="1" colspan="1">(−0.19, 0.16)</td>
<td align="left" rowspan="1" colspan="1">−0.03</td>
<td align="left" rowspan="1" colspan="1">(−0.19, 0.08)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Primary education</td>
<td align="left" rowspan="1" colspan="1">0.18</td>
<td align="left" rowspan="1" colspan="1">(0.06, 0.31)</td>
<td align="left" rowspan="1" colspan="1">0.12</td>
<td align="left" rowspan="1" colspan="1">(0.01, 0.23)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Secondary and higher education</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<italic>Random effects</italic>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Spatial structured</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<sub>11</sub>
<sub>22</sub>
)</td>
<td align="left" rowspan="1" colspan="1">0.13</td>
<td align="left" rowspan="1" colspan="1">(0.002, 0.58)</td>
<td align="left" rowspan="1" colspan="1">0.44</td>
<td align="left" rowspan="1" colspan="1">(0.06, 1.58)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<sub>12</sub>
)</td>
<td align="left" rowspan="1" colspan="1">0.88</td>
<td align="left" rowspan="1" colspan="1">(0.59, 1.31)</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Spatial unstructured</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<sub>11</sub>
, Ω
<sub>22</sub>
)</td>
<td align="left" rowspan="1" colspan="1">1.69</td>
<td align="left" rowspan="1" colspan="1">(1.25, 2.04)</td>
<td align="left" rowspan="1" colspan="1">1.77</td>
<td align="left" rowspan="1" colspan="1">(1.40, 2.20)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<sub>33</sub>
)</td>
<td align="left" rowspan="1" colspan="1">0.34</td>
<td align="left" rowspan="1" colspan="1">(0.21, 0.59)</td>
<td align="left" rowspan="1" colspan="1">-</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>With regards to severity of mortality at household level, we observed that risk of deaths decreased in urban areas (−0.35, 95% CI: −0.66, −0.07), increased if public transport was used (0.23, 95% CI: 0.02, 0.43), increased in female headed households (0.33, 95% CI: 0.23, 0.42), was higher in all socio-economic strata other than in the richest stratum, was positively associated with primary education (0.12, 95% CI: 0.01, 0.23), and increase if time to facility was in hours (0.62, 95% CI: 0.18, 1.12).</p>
<p>The variance components for the random effects, in
<xref ref-type="table" rid="pone-0073500-t003">Table 3</xref>
, showed strong spatial correlation in the prevalence estimated as 0.13 (95% CI: 0.002, 0.58), while severity correlation was 0.44 (95% CI: 0.06, 1.58). The covariance was estimated as 0.88 (95% CI: 0.59, 1.31) showing a strong correlation between prevalence and severity. For the unstructured variance component we obtained 1.69 (95% CI: 1.25, 2.04) and 1.77 (95% CI: 1.40, 2.20) for the logit and NB parts respectively. The unstructured covariance for the two components (prevalence and severity) was 0.34 (95% CI: 0.21, 0.59). The calculated correlation coefficients for the structured and unstructured components were 0.04 and 0.09 respectively.</p>
<p>
<xref ref-type="fig" rid="pone-0073500-g003">Figure 3</xref>
displays age curves for the household head and household member for the count component [since the logit and count components yield similar curves]. In both panels we observed a significant departure from zero, as well as from linearity, although this was more pronounced in age of household head (
<xref ref-type="fig" rid="pone-0073500-g003">Figure 3a</xref>
). In the left panel, the risk of mortality increased between 15–20 years and then decreased up to age 30 years, then rose again steadily up to age 65 years, with a little dip at age 50 years. A similar pattern of up and down continued from age of 65–70 years with a final decrease at age of 80 years. From age of 15 years to 55 years, the risk of death lied below zero, suggesting a reduced mortality risk in such households, while at 60 years to the end we observed a risk of above 0, indicating an increased risk of mortality. Overall the dip in risk was at age of 30 years, and a peak in risk was at 65 years. For the age of household member (
<xref ref-type="fig" rid="pone-0073500-g003">Figure 3b</xref>
), there was an overall decreasing risk with increasing age. Up to age 30–35 years the risk was significantly above zero, an age of increased mortality among household members. From age of 40 years, although the risk remained below zero, this remained non-significant based on the confidence bands. Overall a drop in risk was noted at age of 50 years, while a peak in mortality was displayed at age 80 years.</p>
<fig id="pone-0073500-g003" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.g003</object-id>
<label>Figure 3</label>
<caption>
<title>Nonlinear effects of: (a) age of household head; and (b) age of household member, given by the solid centre line as log relative risk (RR) with corresponding 80% confidence band (dotted outer lines).</title>
</caption>
<graphic xlink:href="pone.0073500.g003"></graphic>
</fig>
<p>
<xref ref-type="fig" rid="pone-0073500-g004">Figure 4</xref>
presents residual spatial effects for both model components. As the patterns from the two plots suggest there were similarity in risk between the logit (prevalence) model and negative binomial (severity) model. The increased odds or risk of mortality were observed in the Caprivi and Kavango regions, while decreased odds and risk of mortality were predicted in Erongo and Omaheke regions. The other regions clearly showed a predicted risk not different from zero i.e
<inline-formula>
<inline-graphic xlink:href="pone.0073500.e085.jpg"></inline-graphic>
</inline-formula>
.
<xref ref-type="fig" rid="pone-0073500-g005">Figure 5</xref>
displays the significance map corresponding to the spatial effects given in
<xref ref-type="fig" rid="pone-0073500-g004">Figure 4</xref>
. The significance map displays three colour schemes: black, white and grey. Black colour denotes regions with strictly negative credible intervals, whereas white denotes regions with strictly positive credible intervals and grey denotes regions with no significance association with the outcome. For the logit spatial effects (
<xref ref-type="fig" rid="pone-0073500-g005">Figure 5a</xref>
), only two regions, Caprivi and Kavango, had significant positive effects implying that the odds of adult mortality were significantly higher in these regions than in others. For the NB part (
<xref ref-type="fig" rid="pone-0073500-g005">Figure 5b</xref>
), we obtained both positive and negative significance areas. Positively significant effects were obtained in Caprivi, Kavango, Ohangwena and Omusati regions suggesting that severity of mortality was higher in these regions compared to others. Regions with negatively significant effects were in Erongo, Khomas, Otjozondjupa and Omaheke. This means the risk of adult mortality was relatively lower in these four regions than in others.</p>
<fig id="pone-0073500-g004" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.g004</object-id>
<label>Figure 4</label>
<caption>
<title>Regional model estimates of residual total spatial effects for the logit part (left panel) and the negative binomial count part (right panel).</title>
</caption>
<graphic xlink:href="pone.0073500.g004"></graphic>
</fig>
<fig id="pone-0073500-g005" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0073500.g005</object-id>
<label>Figure 5</label>
<caption>
<title>Significance effects of region for a nominal level of 80%.</title>
<p>(a) Posterior probabilities for logit part; (b) Posterior probabilities for NB part. Black denotes regions with strictly negative credible intervals. White denotes regions with strictly positive credible intervals. Grey denotes regions with no significance association with the outcome.</p>
</caption>
<graphic xlink:href="pone.0073500.g005"></graphic>
</fig>
</sec>
<sec id="s4">
<title>Discussion</title>
<p>We presented a bivariate or two-part model estimating risk factors of adult and old-age mortality in Namibia based on a nationally representative health survey, 2006/07 DHS, which captured mortality and survivorship data of all members in sampled households. This is a rich source of data, which produced consistent mortality estimates when compared to the census or indirect estimates of WHO and UNDP, despite limitations due to sampling errors and selection bias. We refer interested readers to the work by Bendavid et al.
<xref ref-type="bibr" rid="pone.0073500-Bendavid2">[3]</xref>
which presented tables showing comparison of various estimates using WHO, UNDP and this source.</p>
<p>Our approach used the two-part model by assuming two processes governing mortality observed at households. We conjectured that the process of death occurring (extent) would be different from those influencing multiple deaths (intensity) in the households reporting deaths. We therefore postulated that the risk set, although similar, will have different association for the two processes. Indeed, the different significance covariates obtained in the two sets shows that this is the case. For example “urban” and “public transport” were not significant on the logit model, yet these were significant on the count model (
<xref ref-type="table" rid="pone-0073500-t003">Table 3</xref>
). Similarly, variable like “not married” was significant on the logit model but was not on the count model. The significance of “urban” under the severity model indicates that that although death may occur at a particular household, the risk of repeat within a year is significantly lower in urban areas than in rural areas. Moreover, the difference in magnitude of the estimates between the logit and count model (
<xref ref-type="table" rid="pone-0073500-t003">Table 3</xref>
), suggests that the odds of death occurring is relatively higher than the risk of observing multiple deaths in same household.</p>
<p>In estimating risk factors of adult mortality, we have included both individual and household variables as demand-side factors, and clinic factors as supply-side factors within the conceptual framework of health care
<xref ref-type="bibr" rid="pone.0073500-Andersen1">[26]</xref>
. That is, here we assumed that adult mortality is aggravated by availability and access to health care
<xref ref-type="bibr" rid="pone.0073500-Tanser1">[27]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Stock1">[28]</xref>
. However, the epidemiology of adult mortality is not this simplistic. It is an interaction of a myriad of factors, be it health care, socio-economic, demographic and behavioural factors. Despite limited variables in the DHS we have included variables that fall within all these risk categories. In fact, we extended our model to include spatial random effects, purporting that unobserved or unmeasured covariates present other potential source of risk to mortality. The significance of spatial effects does support our proposition that disparities in health are engrossed in Namibian regions
<xref ref-type="bibr" rid="pone.0073500-Namibia1">[9]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Kazembe1">[10]</xref>
. These results should generate further research to unearth potential risk factors of varied adult health in the country.</p>
<p>From a statistical point of view, several issues arise here. First, we used a bivariate conditional autoregressive process with permits correlation between the logit and count model. Without such an approach, a two-part model would give biased estimates
<xref ref-type="bibr" rid="pone.0073500-Neelon1">[18]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Su1">[21]</xref>
. Second, we applied a structured additive regression (STAR) model in an attempt to explain the complex relationship between adult mortality risk and various risk factors. The use of such models is increasingly being applied in epidemiology. See Kazembe
<xref ref-type="bibr" rid="pone.0073500-Kazembe1">[10]</xref>
; Neelon et al.
<xref ref-type="bibr" rid="pone.0073500-Neelon1">[18]</xref>
; Fahrmeir and Lang
<xref ref-type="bibr" rid="pone.0073500-Fahrmeir1">[29]</xref>
and references therein. STAR models simultaneously model spatially structured random effects, unstructured random effects, nonlinear effects of metrical covariates, together with the usual fixed categorical variables. Third, with respect to the likelihood, the zero-augmented models show a better fit to the data, moreover, the observed zero counts are well captured in Hurdle and ZINB as discussed in Zeileis et al.
<xref ref-type="bibr" rid="pone.0073500-Zeileis1">[30]</xref>
.</p>
<p>Our analysis, however, has the following limitation. While the use of structured and unstructured spatial effects provides robust spatial estimates when the locations are many, in this analysis, we only have 13 provinces to estimate spatial effects and most of the provinces do not have more than one neighbour. This may bias the spatial pattern observed here. Moreover, the large regional areas may conceal or over-step variability of risk within that region in which all areas within are depicted as having common risk of mortality. An ideal analysis would be to use small-areas (districts or constituencies) to assess spatial variability in adult mortality, unfortunately these spatial units were not available at the time of this analysis, but would be worthwhile to pursue this further.</p>
<p>Nevertheless, the results of this study has various important implications. First, it might assist epidemiologists with understanding potential risk factors of adult mortality that needs to be factored in when planning interventions. Second, it could help understand the health consequences of social inequality, human behaviour and demographic factors on adult mortality, which in turn is crucial towards comprehending population dynamics
<xref ref-type="bibr" rid="pone.0073500-Rogers2">[8]</xref>
. Third, it might help social planners with understanding where resources should be targeted. Fourth, it might generate hypotheses as to factors explaining spatial variability in adult mortality, be it HIV epidemic which has a strong age-specific impact on adult mortality
<xref ref-type="bibr" rid="pone.0073500-Bendavid1">[2]</xref>
;
<xref ref-type="bibr" rid="pone.0073500-Ngom1">[4]</xref>
.</p>
<p>In conclusion, this paper explored the use of advanced structural additive models to study adult mortality risks. Our study used the most recent data, although six years old, notwithstanding provide a glimpse of multivariate relationship existing between various factors. It provides a natural picture of forces affecting adult mortality, in contrast to the UNDP/WHO indirect estimates. There is need to consider the risk of cause-specific adult mortality. For example, it may be of interest to ascertain what is the spatial distribution of HIV related mortality as this is the major cause of adult mortality in Namibia
<xref ref-type="bibr" rid="pone.0073500-Bendavid1">[2]</xref>
. However, the lack of HIV related factors in data in general, and in the model in particular, is possibly missing a lot of information and makes it difficult to assess the significance of such factors which is crucial toward informing further public health effort. Be as it may, multilevel models can be used to incorporate reliable HIV related data often aggregated at areal level. Furthermore, multivariate models that categorize death into communicable, non-communicable and accidents/injuries will be worthwhile investigating to explicitly display the spatial variability of risk of adult mortality by cause.</p>
</sec>
<sec sec-type="supplementary-material" id="s5">
<title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0073500.s001">
<label>Appendix S1</label>
<caption>
<p>
<bold>WinBugs Code.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0073500.s001.docx">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back>
<ack>
<p>We acknowledge permission granted by Macro International to use the 2006/2007 Namibia DHS data.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0073500-Jamison1">
<label>1</label>
<mixed-citation publication-type="other">Jamison DT, Feachen RG, Makgoba MW, Bos ER, Baigana FK,
<etal>et al</etal>
. (2006) Disease and mortality in sub-Saharan Africa. 2nd Edition. Washington DC: The World Bank.</mixed-citation>
</ref>
<ref id="pone.0073500-Bendavid1">
<label>2</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bendavid</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Holmes</surname>
<given-names>CB</given-names>
</name>
,
<name>
<surname>Bhattacharya</surname>
<given-names>I</given-names>
</name>
,
<name>
<surname>Muller</surname>
<given-names>G</given-names>
</name>
(
<year>2012</year>
)
<article-title>HIV development assistance and adult mortality in Africa</article-title>
.
<source>J Am Medical Assoc</source>
<volume>307</volume>
:
<fpage>2060</fpage>
<lpage>2067</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Bendavid2">
<label>3</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bendavid</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Seligman</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Kubo</surname>
<given-names>J</given-names>
</name>
(
<year>2011</year>
)
<article-title>Comparative analysis of old-age mortality estimations in Africa</article-title>
.
<source>PLoS ONE</source>
<volume>6(10)</volume>
:
<fpage>e26607</fpage>
<pub-id pub-id-type="pmid">22028921</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Ngom1">
<label>4</label>
<mixed-citation publication-type="other">Ngom P, Clark S (2003) Adult mortality in the era of HIV/AIDS: Sub-Saharan Africa. Workshop on HIV/AIDS and Adult mortality in developing countries. Department of Economic and Social Affairs. United Nations Secretariat. New York, 8–13 September, 2003.</mixed-citation>
</ref>
<ref id="pone.0073500-INDEPTH1">
<label>5</label>
<mixed-citation publication-type="other">INDEPTH Network (2004) In depth model life tables for sub-Saharan Africa. Aldershot: Aldershot Ashgate Publications Ltd.</mixed-citation>
</ref>
<ref id="pone.0073500-Ministry1">
<label>6</label>
<mixed-citation publication-type="other">Ministry of Health and Social Social Services (MoHSS) [Namibia] and Macro International Inc. (2008) Namibia Demography and Health Survey 2006–2007.Windhoek, Namibia and Claverton, Maryland: MoHSS and Macro.</mixed-citation>
</ref>
<ref id="pone.0073500-Rogers1">
<label>7</label>
<mixed-citation publication-type="other">Rogers RG, Hummer RA, Nam CB (2000) Living and dying in the USA: Behavioural, health and social differentials of adult mortality. San Diego: Academic Press.</mixed-citation>
</ref>
<ref id="pone.0073500-Rogers2">
<label>8</label>
<mixed-citation publication-type="other">Rogers RG, Hummer RA, Krueger PM (2002) Adult mortality. In: Poston Jr. Handbook of Population. pp: 283–309.</mixed-citation>
</ref>
<ref id="pone.0073500-Namibia1">
<label>9</label>
<mixed-citation publication-type="other">Namibia Statistics Agency (2012) Namibia household income and expenditure survey 2009/2010. Windhoek: NSA.</mixed-citation>
</ref>
<ref id="pone.0073500-Kazembe1">
<label>10</label>
<mixed-citation publication-type="other">Kazembe LN (2012) Analysis of non-food household expenditures using a multivariate structured additive regression. Proceedings of 54th Conference of South African Statistics Association (SASA). Nelson Mandela Metropolitan University, Port Elizabeth, South Africa. 7th-9th Nov, 2012.</mixed-citation>
</ref>
<ref id="pone.0073500-Obermeyer1">
<label>11</label>
<mixed-citation publication-type="journal">
<name>
<surname>Obermeyer</surname>
<given-names>Z</given-names>
</name>
,
<name>
<surname>Rajaratnam</surname>
<given-names>JK</given-names>
</name>
,
<name>
<surname>Park</surname>
<given-names>CH</given-names>
</name>
,
<name>
<surname>Gakidou</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Hogan</surname>
<given-names>MC</given-names>
</name>
,
<etal>et al</etal>
(
<year>2010</year>
)
<article-title>Measuring Adult Mortality Using Sibling Survival: A New Analytical Method and New Results for 44 Countries, 19742006</article-title>
.
<source>PLoS Med</source>
<volume>7(4)</volume>
:
<fpage>e1000260</fpage>
<pub-id pub-id-type="pmid">20405004</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Houle1">
<label>12</label>
<mixed-citation publication-type="other">Houle BC (2011) Bio-social Determinants of Child and Adult Mortality in South Africa [PhD thesis]. University of Washington. pp. 320.</mixed-citation>
</ref>
<ref id="pone.0073500-Ali1">
<label>13</label>
<mixed-citation publication-type="journal">
<name>
<surname>Ali</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Jin</surname>
<given-names>Y</given-names>
</name>
,
<name>
<surname>Kim</surname>
<given-names>DR</given-names>
</name>
,
<name>
<surname>De</surname>
<given-names>ZB</given-names>
</name>
,
<name>
<surname>Park</surname>
<given-names>JK</given-names>
</name>
,
<etal>et al</etal>
(
<year>2007</year>
)
<article-title>Spatial risk for gender specific adult mortality in an area of southern China</article-title>
.
<source>Int J Health Geogr</source>
<volume>6</volume>
:
<fpage>31</fpage>
<pub-id pub-id-type="pmid">17645807</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Turra1">
<label>14</label>
<mixed-citation publication-type="journal">
<name>
<surname>Turra</surname>
<given-names>CM</given-names>
</name>
,
<name>
<surname>Goldman</surname>
<given-names>N</given-names>
</name>
(
<year>2007</year>
)
<article-title>Socio-economic differences in mortality among U.S. adults: insights into the Hispanic paradox</article-title>
.
<source>J Gerontol B Psychol Sci Soc Sci</source>
<volume>62(3)</volume>
:
<fpage>S184</fpage>
<lpage>92</lpage>
<pub-id pub-id-type="pmid">17507594</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Fantahun1">
<label>15</label>
<mixed-citation publication-type="journal">
<name>
<surname>Fantahun</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Berhane</surname>
<given-names>Y</given-names>
</name>
,
<name>
<surname>Högberg</surname>
<given-names>U</given-names>
</name>
,
<name>
<surname>Wall</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Byass</surname>
<given-names>P</given-names>
</name>
(
<year>2008</year>
)
<article-title>Young adult and middle age mortality in Butajira demographic surveillance site, Ethiopia: lifestyle, gender and household economy</article-title>
.
<source>BMC Public Health</source>
<volume>8(7)</volume>
:
<fpage>268</fpage>
<pub-id pub-id-type="pmid">18671854</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-TuckerSeeley1">
<label>16</label>
<mixed-citation publication-type="journal">
<name>
<surname>Tucker-Seeley</surname>
<given-names>RD</given-names>
</name>
,
<name>
<surname>Li</surname>
<given-names>Y</given-names>
</name>
,
<name>
<surname>Sorensen</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Subramanian</surname>
<given-names>SV</given-names>
</name>
(
<year>2011</year>
)
<article-title>Lifecourse socioeconomic circumstances and multimorbidity among older adults</article-title>
.
<source>BMC Public Health</source>
<volume>11</volume>
:
<fpage>313</fpage>
<pub-id pub-id-type="pmid">21569558</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Winkelmann1">
<label>17</label>
<mixed-citation publication-type="other">Winkelmann R (2008) Analysis of count data. Fifth Edition. Hiedelberg: Springer.</mixed-citation>
</ref>
<ref id="pone.0073500-Neelon1">
<label>18</label>
<mixed-citation publication-type="journal">
<name>
<surname>Neelon</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Ghosh</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Loebs</surname>
<given-names>PF</given-names>
</name>
(
<year>2013</year>
)
<article-title>A spatial Poisson hurdle model for exploring geographical variation in emergency department visits</article-title>
.
<source>J R Statist Soc Ser A</source>
<volume>176</volume>
:
<fpage>389</fpage>
<lpage>413</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Winkelmann2">
<label>19</label>
<mixed-citation publication-type="journal">
<name>
<surname>Winkelmann</surname>
<given-names>R</given-names>
</name>
(
<year>2004</year>
)
<article-title>Health care reform and the number of doctor visits – an econometric analysis</article-title>
.
<source>J Appl Econometrics</source>
<volume>19</volume>
:
<fpage>455</fpage>
<lpage>472</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Olsen1">
<label>20</label>
<mixed-citation publication-type="journal">
<name>
<surname>Olsen</surname>
<given-names>MK</given-names>
</name>
,
<name>
<surname>Schafer</surname>
<given-names>J</given-names>
</name>
(
<year>2001</year>
)
<article-title>A two-part random-effects model for semi-continuous longitudinal data</article-title>
.
<source>J Am Statist Assoc</source>
<volume>96</volume>
:
<fpage>247</fpage>
<lpage>281</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Su1">
<label>21</label>
<mixed-citation publication-type="journal">
<name>
<surname>Su</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Tom</surname>
<given-names>BDM</given-names>
</name>
,
<name>
<surname>Farewell</surname>
<given-names>VT</given-names>
</name>
(
<year>2009</year>
)
<article-title>Bias in two-part mixed models for longitudinal semi-continuous data</article-title>
.
<source>Biostatist</source>
<volume>10</volume>
:
<fpage>374</fpage>
<lpage>389</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Jin1">
<label>22</label>
<mixed-citation publication-type="journal">
<name>
<surname>Jin</surname>
<given-names>X</given-names>
</name>
,
<name>
<surname>Carlin</surname>
<given-names>BP</given-names>
</name>
,
<name>
<surname>Banerjee</surname>
<given-names>S</given-names>
</name>
(
<year>2005</year>
)
<article-title>Generalized hierarchical multivariate CAR models for areal data</article-title>
.
<source>Biometrics</source>
<volume>61</volume>
:
<fpage>950</fpage>
<lpage>961</lpage>
<pub-id pub-id-type="pmid">16401268</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Spiegelhalter1">
<label>23</label>
<mixed-citation publication-type="other">Spiegelhalter DJ, Thomas A, Best NG (1999)
<italic>WinBUGS Version 1.2 User Manual</italic>
MRC Biostatistics Unit. Cambridge University.</mixed-citation>
</ref>
<ref id="pone.0073500-Gelman1">
<label>24</label>
<mixed-citation publication-type="journal">
<name>
<surname>Gelman</surname>
<given-names>A</given-names>
</name>
(
<year>2006</year>
)
<article-title>Prior distributions for variance parameters in hierarchical models</article-title>
.
<source>Bayesian Anal</source>
<volume>1</volume>
:
<fpage>515534</fpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Spiegelhalter2">
<label>25</label>
<mixed-citation publication-type="journal">
<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>
,
<name>
<surname>van der Linde</surname>
<given-names>A</given-names>
</name>
(
<year>2002</year>
)
<article-title>Bayesian measures of model complexity and fit</article-title>
.
<source>J R Statist Soc Ser B (with discussion)</source>
<volume>64</volume>
:
<fpage>1</fpage>
<lpage>34</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Andersen1">
<label>26</label>
<mixed-citation publication-type="journal">
<name>
<surname>Andersen</surname>
<given-names>RM</given-names>
</name>
(
<year>1995</year>
)
<article-title>Revisiting the behavioural model and access to medical care: does it matter?</article-title>
<source>J Health Social Behaviour</source>
<volume>36</volume>
:
<fpage>1</fpage>
<lpage>10</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Tanser1">
<label>27</label>
<mixed-citation publication-type="journal">
<name>
<surname>Tanser</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Hosegood</surname>
<given-names>V</given-names>
</name>
,
<name>
<surname>Benzler</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Solarsh</surname>
<given-names>G</given-names>
</name>
(
<year>2001</year>
)
<article-title>New approaches to spatially analyse primary health care usage patterns in rural South Africa</article-title>
.
<source>Trop Med Int Health</source>
<volume>6</volume>
:
<fpage>826</fpage>
<lpage>838</lpage>
<pub-id pub-id-type="pmid">11679131</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0073500-Stock1">
<label>28</label>
<mixed-citation publication-type="journal">
<name>
<surname>Stock</surname>
<given-names>R</given-names>
</name>
(
<year>1983</year>
)
<article-title>Distance and the utilization of health facilities in rural Nigeria</article-title>
.
<source>Social Sci Med</source>
<volume>17</volume>
:
<fpage>563</fpage>
<lpage>570</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Fahrmeir1">
<label>29</label>
<mixed-citation publication-type="journal">
<name>
<surname>Fahrmeir</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Lang</surname>
<given-names>S</given-names>
</name>
(
<year>2001</year>
)
<article-title>Bayesian Inference for Generalized Additive Mixed Models Based on Markov Random Field Priors</article-title>
.
<source>J R Statist Soc Series C (Appl Stat)</source>
<volume>50</volume>
:
<fpage>201</fpage>
<lpage>220</lpage>
</mixed-citation>
</ref>
<ref id="pone.0073500-Zeileis1">
<label>30</label>
<mixed-citation publication-type="other">Zeileis A, Kleiber C, Jackman S (2008) Regression Models for Count Data in R. J Statist Software 27(8). Available:
<ext-link ext-link-type="uri" xlink:href="http://www.jstatsoft.org/v27/i08/">http://www.jstatsoft.org/v27/i08/</ext-link>
</mixed-citation>
</ref>
</ref-list>
</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 002A53 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 002A53 | 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:3774685
   |texte=   A Bayesian Two Part Model Applied to Analyze Risk Factors of Adult Mortality with Application to Data from Namibia
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:24066052" \
       | 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