Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Geographical variations and contextual effects on age of initiation of sexual intercourse among women in Nigeria: a multilevel and spatial analysis</title>
<author><name sortKey="Uthman, Olalekan A" sort="Uthman, Olalekan A" uniqKey="Uthman O" first="Olalekan A" last="Uthman">Olalekan A. Uthman</name>
<affiliation><nlm:aff id="I1">Center for Evidence-Based Global Health, Save the Youth Initiative, Ilorin, Kwara state, Nigeria</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">18513414</idno>
<idno type="pmc">2414823</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2414823</idno>
<idno type="RBID">PMC:2414823</idno>
<idno type="doi">10.1186/1476-072X-7-27</idno>
<date when="2008">2008</date>
<idno type="wicri:Area/Pmc/Corpus">000D21</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000D21</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Geographical variations and contextual effects on age of initiation of sexual intercourse among women in Nigeria: a multilevel and spatial analysis</title>
<author><name sortKey="Uthman, Olalekan A" sort="Uthman, Olalekan A" uniqKey="Uthman O" first="Olalekan A" last="Uthman">Olalekan A. Uthman</name>
<affiliation><nlm:aff id="I1">Center for Evidence-Based Global Health, Save the Youth Initiative, Ilorin, Kwara state, Nigeria</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">International Journal of Health Geographics</title>
<idno type="eISSN">1476-072X</idno>
<imprint><date when="2008">2008</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>The age of initiation of sexual intercourse is an increasingly important issue to study given that sexually active young women are at risk of multiple outcomes including early pregnancies, vesico-vaginal fistula, and sexually transmitted infections. Much research has focused on the demographic, familial, and social factors associated with sexual initiation and reasons adolescents begin having consensual intercourse. Less is known, however, about the geographical and contextual factors associated with age of initiation of sexual intercourse. Therefore, the purpose of this study was to examine the extent of regional and state disparities in age of initiation of sexual intercourse and to examine individual- and community-level predictors of early sexual debut.</p>
</sec>
<sec sec-type="methods"><title>Methods</title>
<p>Multilevel logistic regression models were applied to data on 5531 ever or currently married women who had participated in 2003 Nigeria Demographic and Health Survey. Coital debut at 15 years or younger was used to define early sexual debut. Exploratory spatial data analysis methods were used to study geographic variation in age at first sexual intercourse.</p>
</sec>
<sec><title>Results</title>
<p>The median age at first sexual intercourse for all women included in the study was 15 years (range; 14 – 19). North West and North East had the highest proportion of women who had reported early sexual debut (61% – 78%). The spatial distribution of age of initiation of sexual intercourse was nonrandom and clustered with a Moran's I = 0.635 (p = .001). There was significant positive spatial relationship between median age of marriage and spatial lag of median age of sexual debut (Bivariate Moran's I = 0.646, (p = .001). After adjusting for both individual-level and contextual factors, the probability of starting sex at an earlier age was associated with respondents' current age, education attainment, ethnicity, region, and community median age of marriage.</p>
</sec>
<sec><title>Conclusion</title>
<p>The study found that individual-level and community contextual characteristics were independently associated with early sexual debut, suggesting that interventions to reduce adolescent high-risk sexual behaviour should focus on high-risk places as well as high-risk groups of people.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></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">Int J Health Geogr</journal-id>
<journal-title>International Journal of Health Geographics</journal-title>
<issn pub-type="epub">1476-072X</issn>
<publisher><publisher-name>BioMed Central</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">18513414</article-id>
<article-id pub-id-type="pmc">2414823</article-id>
<article-id pub-id-type="publisher-id">1476-072X-7-27</article-id>
<article-id pub-id-type="doi">10.1186/1476-072X-7-27</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research</subject>
</subj-group>
</article-categories>
<title-group><article-title>Geographical variations and contextual effects on age of initiation of sexual intercourse among women in Nigeria: a multilevel and spatial analysis</article-title>
</title-group>
<contrib-group><contrib id="A1" corresp="yes" contrib-type="author"><name><surname>Uthman</surname>
<given-names>Olalekan A</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<email>uthlekan@yahoo.com</email>
</contrib>
</contrib-group>
<aff id="I1"><label>1</label>
Center for Evidence-Based Global Health, Save the Youth Initiative, Ilorin, Kwara state, Nigeria</aff>
<pub-date pub-type="collection"><year>2008</year>
</pub-date>
<pub-date pub-type="epub"><day>30</day>
<month>5</month>
<year>2008</year>
</pub-date>
<volume>7</volume>
<fpage>27</fpage>
<lpage>27</lpage>
<ext-link ext-link-type="uri" xlink:href="http://www.ij-healthgeographics.com/content/7/1/27"></ext-link>
<history><date date-type="received"><day>29</day>
<month>1</month>
<year>2008</year>
</date>
<date date-type="accepted"><day>30</day>
<month>5</month>
<year>2008</year>
</date>
</history>
<permissions><copyright-statement>Copyright © 2008 Uthman; licensee BioMed Central Ltd.</copyright-statement>
<copyright-year>2008</copyright-year>
<copyright-holder>Uthman; licensee BioMed Central Ltd.</copyright-holder>
<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"></ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p>
<pmc-comment>
Uthman
A
Olalekan
uthlekan@yahoo.com
Geographical variations and contextual effects on age of initiation of sexual intercourse among women in Nigeria: a multilevel and spatial analysis
2008 International Journal of Health Geographics 7(1): 27-. (2008) 1476-072X(2008)7:1<27> urn:ISSN:1476-072X </pmc-comment>
</license>
</permissions>
<abstract><sec><title>Background</title>
<p>The age of initiation of sexual intercourse is an increasingly important issue to study given that sexually active young women are at risk of multiple outcomes including early pregnancies, vesico-vaginal fistula, and sexually transmitted infections. Much research has focused on the demographic, familial, and social factors associated with sexual initiation and reasons adolescents begin having consensual intercourse. Less is known, however, about the geographical and contextual factors associated with age of initiation of sexual intercourse. Therefore, the purpose of this study was to examine the extent of regional and state disparities in age of initiation of sexual intercourse and to examine individual- and community-level predictors of early sexual debut.</p>
</sec>
<sec sec-type="methods"><title>Methods</title>
<p>Multilevel logistic regression models were applied to data on 5531 ever or currently married women who had participated in 2003 Nigeria Demographic and Health Survey. Coital debut at 15 years or younger was used to define early sexual debut. Exploratory spatial data analysis methods were used to study geographic variation in age at first sexual intercourse.</p>
</sec>
<sec><title>Results</title>
<p>The median age at first sexual intercourse for all women included in the study was 15 years (range; 14 – 19). North West and North East had the highest proportion of women who had reported early sexual debut (61% – 78%). The spatial distribution of age of initiation of sexual intercourse was nonrandom and clustered with a Moran's I = 0.635 (p = .001). There was significant positive spatial relationship between median age of marriage and spatial lag of median age of sexual debut (Bivariate Moran's I = 0.646, (p = .001). After adjusting for both individual-level and contextual factors, the probability of starting sex at an earlier age was associated with respondents' current age, education attainment, ethnicity, region, and community median age of marriage.</p>
</sec>
<sec><title>Conclusion</title>
<p>The study found that individual-level and community contextual characteristics were independently associated with early sexual debut, suggesting that interventions to reduce adolescent high-risk sexual behaviour should focus on high-risk places as well as high-risk groups of people.</p>
</sec>
</abstract>
</article-meta>
</front>
<body><sec><title>Background</title>
<p>The age of initiation of sexual intercourse is an increasingly important issue to study given that sexually active young women are at risk of multiple outcomes including early pregnancies, vesico-vaginal fistula, and sexually transmitted infections. To date, most studies on the effects of early initiation of sexual intercourse among women have focused on those who had their sexual debut before marriage. Bruce et al[<xref ref-type="bibr" rid="B1">1</xref>
] found that 87 out of 100 girls in Nigeria who had sex in the previous week were married. In the northern regions of Nigeria, where child marriage is common, virtually all sexual activity among girls occurs within the context of marriage. In Nigeria, where women are more susceptible to HIV, child marriage and early sexual intercourse may be significant risk factors for adolescent girls [<xref ref-type="bibr" rid="B2">2</xref>
]. Some 3.9 percent of Nigeria's population is infected with HIV [<xref ref-type="bibr" rid="B3">3</xref>
]. Girls aged 15–24 are about twice as likely as boys the same age to be infected [<xref ref-type="bibr" rid="B4">4</xref>
]. Teenage girls are also more susceptible than older women to sexually transmitted infections (STI)[<xref ref-type="bibr" rid="B4">4</xref>
].</p>
<p>Much research has focused on the demographic, familial, and social factors associated with sexual initiation and reasons adolescents begin having consensual intercourse [<xref ref-type="bibr" rid="B5">5</xref>
-<xref ref-type="bibr" rid="B7">7</xref>
]. Less is known, however, about the geographical and contextual factors associated with age of initiation of sexual intercourse. In the present investigation, more appropriate and recent multi-level modelling techniques were used to examine the association between age of initiation of sexual intercourse and contextual factors. The analysis of geographic variation in age of sexual initiation is important. First, such a focus is consistent with the Nigerian national health initiative, which aims to reduce and ultimately eliminate health inequalities among gender, racial/ethnic, socioeconomic, and geographic groups. In addition, a geographic analysis should help identify states and regions of the country that have a relatively high proportion of women with high risk sexual behaviour, which may in turn lead to the development and implementation of more effective geographically differentiated intervention programs. Exploratory spatial data analysis (ESDA)[<xref ref-type="bibr" rid="B8">8</xref>
] methods were used to study geographic variations in age of initiation of sexual intercourse in Nigeria. ESDA is an extension of exploratory data analysis (EDA) that focuses on detecting spatial patterns in data and the generation of hypotheses based on the spatial patterns in the data.</p>
<p>Therefore, the purpose of this study was (1) to estimate prevalence of early sexual debut among ever or currently married women in six geographic regions (37 states) of Nigeria using a large, nationally representative sample, (2) to examine the extent of regional and state disparities in age of initiation of sexual intercourse, (3) to identify individual and contextual predictors of early sexual debut, and (4) to study the association between age of marriage and age of initiation of sexual intercourse.</p>
</sec>
<sec sec-type="methods"><title>Methods</title>
<sec><title>Setting</title>
<p>Nigeria is located in western Africa on the Gulf of Guinea and has a total area of 923,768 kilometer squared (km<sup>2</sup>
), making it the world's 32nd-largest country (after Tanzania). Nigeria is the most populous country in Africa. The United Nations estimates that the population in 2004 was at 131,530,000, with the population distributed as 48.3% urban and 51.7% rural and population density at 139 people per km<sup>2</sup>
. Nigeria has more than 250 ethnic groups, with varying languages and customs, creating a country of rich ethnic diversity. The largest ethnic groups are the Fulani/Hausa, Yoruba, Igbo, accounting for 68% of population, while the Edo, Ijaw, Kanuri, Ibibio, Ebira Nupe and Tiv comprise 27%; other minorities make up the remaining 5%. The middle belt of Nigeria is known for its diversity of ethnic groups, including the Pyem, Goemai, and Kofyar.</p>
<p>One third of Nigeria's total population are youth between the ages of 10 and 24[<xref ref-type="bibr" rid="B9">9</xref>
]. By 2025, the number of Nigerian youth will exceeds 57 million [<xref ref-type="bibr" rid="B9">9</xref>
]. Lack of sexual health information and services places these young people at risk for pregnancy, abortion, sexually transmitted infections (STI), and HIV/AIDS [<xref ref-type="bibr" rid="B9">9</xref>
]. In addition, early marriage and childbearing limit youth's educational and employment opportunities [<xref ref-type="bibr" rid="B9">9</xref>
]. Among women ages 20 to 24, 19.8 percent reported having married by age 15, 39.6 percent by age 18, and 52.7 percent by age 20. In 1999, Nigeria's adolescent fertility rate was 111 births per 1,000 women ages 15 to 19, and Nigerian women averaged more than five births during their lifetime[<xref ref-type="bibr" rid="B9">9</xref>
].</p>
</sec>
<sec><title>Study design</title>
<p>Cross-sectional and population-based study using data from the 2003 Nigeria Demographic and Health survey (NDHS).</p>
</sec>
<sec><title>Sampling technique</title>
<p>Methods used in the NDHS have been published elsewhere [<xref ref-type="bibr" rid="B10">10</xref>
]. Briefly, the survey used a two-stage cluster sampling technique. The country was stratified into 36 states and the Federal Capital Territory (FCT) of Abuja. Each domain is made up of enumeration areas (EAs) established by a general population and housing census in 1991. The sampling frame was a list of all EAs (clusters). Within each domain, a two-stage sample was selected. The first stage involved selecting 466 clusters (primary sampling units) with a probability proportional to the size, the size being the number of households in the cluster. The second stage involved the systematic sampling of households from the selected clusters. A nationally representative probability sample of 7864 households was then selected from the clusters, in which all women aged 15 to 49 years were eligible to be interviewed.</p>
</sec>
<sec><title>Data collection</title>
<p>Data collection procedures have been published elsewhere [<xref ref-type="bibr" rid="B10">10</xref>
]. Briefly, data were collected by visiting households and conducting face-to-face interviews to obtain information on demographic characteristics, wealth, anthropometry, female genital cutting, HIV knowledge, and sexual behaviour between March and August 2003.</p>
</sec>
<sec><title>Ethical consideration</title>
<p>This study is based on an analysis of existing survey data with all identifier information removed. The survey was approved by the Ethics Committee of the ORC Macro at Calverton in the USA and by the National Ethics Committee in the Ministry of Health in Nigeria. All study participants gave informed consent before participation and all information was collected confidentially.</p>
</sec>
<sec><title>Variables</title>
<sec><title>Outcome variable</title>
<p>For the present study, 5531 ever or currently married women, all of whom who have had at least one episode of sexual intercourse in their lifetime, were drawn from the overall sample of 7620 women to examine factors associated with early sexual initiation. Thus, if a woman has never had sexual intercourse in her life (e.g., a "virgin") and never married, she was not included in the analysis. Coital debut at 15 years or younger (Yes or No) was used to define early sexual debut. The cut-off point for age at coital debut was based on the median values of Nigeria women included in the study.</p>
</sec>
<sec><title>Explanatory variables</title>
<p>Four variables characterising individual-level features shown to affect age of initiation at first sexual intercourse were included in this study. The variables were age, educational attainment, ethnicity and religion affiliation. There were five community-level variables. These include: the degree of ethnic diversity in a community, the community urbanization measure (urban residence) and regional dummies (region). The other contextual variables include community median age of marriage and community median spousal age gap. The full definition and coding for each explanatory variable is given in Table <xref ref-type="table" rid="T1">1</xref>
.</p>
<table-wrap position="float" id="T1"><label>Table 1</label>
<caption><p>Variables, measures, and coding</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><td align="left"><bold><italic>Variable</italic>
</bold>
</td>
<td align="left"><bold><italic>Measure</italic>
</bold>
</td>
</tr>
</thead>
<tbody><tr><td align="left"><bold><italic>Dependent variable</italic>
</bold>
</td>
<td></td>
</tr>
<tr><td align="left">Early sexual debut</td>
<td align="left">Coital debut at 15 years or younger (yes = 1, no = 0)</td>
</tr>
<tr><td align="left"><bold><italic>Individual characteristics</italic>
</bold>
</td>
<td></td>
</tr>
<tr><td align="left">Age</td>
<td align="left">Completed years [1) 15–24, 2) 25–34, 3) 35+]</td>
</tr>
<tr><td align="left">Education</td>
<td align="left">Highest level of education [1) No education, 2) primary, 3) secondary or higher]</td>
</tr>
<tr><td align="left">Ethnicity</td>
<td align="left">Ethnicity [1) Hausa or Fulani, 2) Igbo, 3) Yoruba, 4) Others]</td>
</tr>
<tr><td align="left">Religion</td>
<td align="left">Religion affiliation [1) Islam, 2) Catholic, Protestant, Other Christian, 3) traditionalist, others]</td>
</tr>
<tr><td align="left"><bold><italic>Community characteristics</italic>
</bold>
</td>
<td></td>
</tr>
<tr><td align="left">Type of residence</td>
<td align="left">Cluster classification [1) rural or remote, 2) urban]</td>
</tr>
<tr><td align="left">Community ethnic diversity</td>
<td align="left">Number of ethnic group in a community, represents the likelihood that two persons, chosen at random from the same area, belong to different ethnic group.</td>
</tr>
<tr><td align="left">Community median age of marriage</td>
<td align="left">Median age of women at first marriage in community</td>
</tr>
<tr><td align="left">Community median spousal age gap</td>
<td align="left">Median spousal age gap in community</td>
</tr>
<tr><td align="left">Region</td>
<td align="left">The six geopolitical regions [1) north-central, 2) north-east, 3) north-west, 4) south-east, 5) south-south, 6) south-west]</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec><title>Statistical analyses</title>
<sec><title>Descriptive and bivariate association</title>
<p>The descriptive statistics show the distribution of respondents by the key variables. Values were expressed as absolute numbers (percentages) and mean (standard error) for categorical and continuous variables respectively. Contingency tables were analyzed using Pearson's chi-squared test. For continuous variable, a Mann-Whitney test was used to examine the association between early sexual debut and each explanatory variable. All cases in the DHS data were given weights to adjust for differences in probability of selection and to adjust for non-response in order to produce the proper representation. Individual weights were used for descriptive statistics in this study, using Stata 10 for Windows [<xref ref-type="bibr" rid="B11">11</xref>
].</p>
</sec>
<sec><title>Multilevel modelling</title>
<p>Given the hierarchical structure of the sample and the binary outcome, a logistic multilevel modeling approach was adopted [<xref ref-type="bibr" rid="B12">12</xref>
]. A three-level model with a binary response (<italic>y</italic>
, reported as having started sexual intercourse at an early age) for women <italic>i </italic>
living in community <italic>j </italic>
in state <italic>k </italic>
of the form:</p>
<p><disp-formula id="bmcM1"><label>(1)</label>
<italic>π</italic>
<sub><italic>ijk</italic>
</sub>
: <italic>y</italic>
<sub><italic>ijk</italic>
</sub>
~Bemoulli(1, <italic>π</italic>
<sub><italic>ijk</italic>
</sub>
)</disp-formula>
</p>
<p>The probability was related to a set of categorical predictors, <italic>X</italic>
; and a random effect for each level, by a logit-link function as</p>
<p><disp-formula id="bmcM2"><label>(2)</label>
logit(<italic>π</italic>
<sub><italic>ijk</italic>
</sub>
) = log [<italic>π</italic>
<sub><italic>ijk</italic>
</sub>
/(1-<italic>π</italic>
<sub><italic>ijk</italic>
</sub>
)] = <italic>β</italic>
<sub>0 </sub>
+ <italic>βX</italic>
<sub><italic>ijk </italic>
</sub>
+ <italic>u</italic>
<sub>0<italic>jk </italic>
</sub>
+ <italic>v</italic>
<sub>0<italic>k</italic>
</sub>
</disp-formula>
</p>
<p>The linear predictor on the right-hand side of the equation consisted of a fixed part (<italic>β</italic>
<sub>0 </sub>
+ <italic>βX</italic>
<sub><italic>ijk</italic>
</sub>
) estimating the conditional coefficients for the explanatory, and two random intercepts attributable to communities (<italic>u</italic>
<sub>0<italic>jk</italic>
</sub>
) and states (<italic>v</italic>
<sub>0<italic>k</italic>
</sub>
) with each assumed to have an independent and identical distribution and variance estimated at each level.</p>
<p>The analysis was done in three steps. In Model 1 (empty model), no explanatory variable was included. In Model 2, only individual-level factors were included. In Model 3, community contextual factors were added to Model 2. The results of fixed effects (measures of association) were shown as odds ratios (ORs) with 95% confidence intervals (CIs). The results of random effects (measures of variation) were presented as variance partition coefficient and percentage change in variance.</p>
<p>The MLwiN software, version 2.0.2 [<xref ref-type="bibr" rid="B13">13</xref>
], was used for the analyses. Parameters were estimated using the Markov Chain Monte Carlo (MCMC) procedure. The default settings in MLwiN were used for the analyses, i.e., chains of length 5000 after a burn-in of 500. The Deviance Information Criterion (DIC) was used as a measure of how well our different models fitted the data. A lower value on DIC indicates a better fit of the model [<xref ref-type="bibr" rid="B14">14</xref>
,<xref ref-type="bibr" rid="B15">15</xref>
].</p>
</sec>
<sec><title>Exploratory spatial data analysis (ESDA)</title>
<p>The spatial distribution of the rate of early sexual initiation was first visualized and spatially smoothed prevalence was obtained by means of spatial empirical Bayes estimation[<xref ref-type="bibr" rid="B8">8</xref>
]. The spatial pattern of prevalence of early sexual debut was analyzed using Univariate LISA Moran Scatter plot[<xref ref-type="bibr" rid="B16">16</xref>
]. The bivariate spatial relationship between age of sexual initiation and age of marriage was studied using bivariate Moran's I[<xref ref-type="bibr" rid="B17">17</xref>
]. The bivariate LISA (BiLISA) Moran's I "gives an indication of the degree of linear association (positive or negative) between the value for one variable at a given location <italic>I </italic>
and the average of another variable at neighboring locations." Similar to LISA, BiLISA suggests two classes of positive spatial correlation, or spatial clusters (High-High and Low-Low), and two classes of negative spatial correlation, or spatial outliers (High-Low and Low-High). Inference for Moran's I was based on a permutation approach, in which a reference distribution is calculated for spatially random layouts with the same data as observed. The randomization uses an algorithm to generate spatially random simulated data sets as outlined by Anselin[<xref ref-type="bibr" rid="B18">18</xref>
]. 999 random permutations were used in constructing the reference distribution. ESDA is implemented through the <italic>GeoDa </italic>
software[<xref ref-type="bibr" rid="B19">19</xref>
]. GeoDa provides a very user-friendly environment to implement ESDA methods and is freely downloadable.</p>
</sec>
</sec>
</sec>
<sec><title>Results</title>
<sec><title>Sample characteristics and bivariate associations</title>
<p>The weighted descriptive statistics are presented in Table <xref ref-type="table" rid="T2">2</xref>
, 5531 ever or currently married women nested within 177 communities nested within 37 states were analyzed in this study. Among women included in the study, nearly half (48%) reported first sex by age 15. Most of the respondents reside in the rural areas (69%), and were Muslim (60%). The results of bivariate associations with early sexual debut are presented in Table <xref ref-type="table" rid="T3">3</xref>
. Compare to Christian women, Muslim women were more likely to have started sex intercourse at an earlier age (62% versus 38%, p = .001). Counterintuitively, urban women were less likely to have reported early sexual debut (39% versus 54%; p = .001). Compared to Hausa women, Yoruba (14% vs. 71%) and Igbo women (23% vs. 71%) were less likely to have started sex early. In addition, compared to their counterparts in the southern part of the country, women from northern regions were more likely to have started sex at an early age.</p>
<table-wrap position="float" id="T2"><label>Table 2</label>
<caption><p>Socio-demographic characteristics of ever or currently married women, Nigeria 2003*</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><td align="left"><bold><italic>Measure/variable</italic>
</bold>
</td>
<td align="right"><bold><italic>Weighted sample (percentage)</italic>
</bold>
</td>
</tr>
</thead>
<tbody><tr><td align="left"><bold><italic>Dependent variable</italic>
</bold>
</td>
<td></td>
</tr>
<tr><td align="left">Early sexual debut (15 years or less)</td>
<td></td>
</tr>
<tr><td align="left"> Yes</td>
<td align="right">2755 (48.4)</td>
</tr>
<tr><td align="left"> No</td>
<td align="right">2937 (51.6)</td>
</tr>
<tr><td align="left"><bold><italic>Individual characteristics</italic>
</bold>
</td>
<td></td>
</tr>
<tr><td align="left">Age (completed years)</td>
<td></td>
</tr>
<tr><td align="left"> 15 – 24</td>
<td align="right">1523 (26.8)</td>
</tr>
<tr><td align="left"> 25 – 34</td>
<td align="right">2106 (37.0)</td>
</tr>
<tr><td align="left"> 35+</td>
<td align="right">2063 (36.2)</td>
</tr>
<tr><td align="left">Education</td>
<td></td>
</tr>
<tr><td align="left"> No education</td>
<td align="right">3044 (53.5)</td>
</tr>
<tr><td align="left"> Primary</td>
<td align="right">1266 (22.2)</td>
</tr>
<tr><td align="left"> Secondary or higher</td>
<td align="right">1381 (24.3)</td>
</tr>
<tr><td align="left">Ethnicity</td>
<td></td>
</tr>
<tr><td align="left"> Hausa/Fulani</td>
<td align="right">2387 (41.9)</td>
</tr>
<tr><td align="left"> Igbo</td>
<td align="right">582 (10.2)</td>
</tr>
<tr><td align="left"> Yoruba</td>
<td align="right">535 (9.4)</td>
</tr>
<tr><td align="left"> Others</td>
<td align="right">2168 (38.1)</td>
</tr>
<tr><td align="left">Religion</td>
<td></td>
</tr>
<tr><td align="left"> Christian</td>
<td align="right">2156 (37.9)</td>
</tr>
<tr><td align="left"> Islam</td>
<td align="right">3444 (60.5)</td>
</tr>
<tr><td align="left"> Others</td>
<td align="right">89 (1.6)</td>
</tr>
<tr><td align="left"><bold><italic>Community characteristics</italic>
</bold>
</td>
<td></td>
</tr>
<tr><td align="left">Type of residence</td>
<td></td>
</tr>
<tr><td align="left"> Rural</td>
<td align="right">3907 (68.6)</td>
</tr>
<tr><td align="left"> Urban</td>
<td align="right">1785 (31.4)</td>
</tr>
<tr><td align="left">Community ethnic diversity<sup>†</sup>
, Mean (SE)</td>
<td align="right">3.7 (0.09)</td>
</tr>
<tr><td align="left">Community median age of marriage, Mean (SE)</td>
<td align="right">16.2 (0.11)</td>
</tr>
<tr><td align="left">Community median spousal age gap, Mean (SE)</td>
<td align="right">9.7 (0.08)</td>
</tr>
<tr><td align="left">Geopolitical region</td>
<td></td>
</tr>
<tr><td align="left"> North-central</td>
<td align="right">807 (14.2)</td>
</tr>
<tr><td align="left"> North-east</td>
<td align="right">1197 (21.2)</td>
</tr>
<tr><td align="left"> North-west</td>
<td align="right">1949 (34.2)</td>
</tr>
<tr><td align="left"> South-east</td>
<td align="right">396 (7.0)</td>
</tr>
<tr><td align="left"> South-south</td>
<td align="right">764 (13.4)</td>
</tr>
<tr><td align="left"> South-west</td>
<td align="right">579 (10.2)</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>Abbreviation – SE – Standard Error</p>
<p><sup>†</sup>
Number of ethnic group in a community, represents the likelihood that two persons, chosen at random from the same area, belong to different ethnic group</p>
<p>*Data source – 2003 Nigerian Demographic and Survey. Sample size – 5531</p>
</table-wrap-foot>
</table-wrap>
<table-wrap position="float" id="T3"><label>Table 3</label>
<caption><p>Bivariate association between early sexual debut and selected individual compositional and community contextual characteristics of ever or currently married Nigerian women, 2003*</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><td></td>
<td align="center" colspan="3"><bold><italic>Early sexual debut</italic>
</bold>
</td>
</tr>
</thead>
<tbody><tr><td></td>
<td align="left">No [%(n)]</td>
<td align="left">Yes [%(n)]</td>
<td align="left">p-value</td>
</tr>
<tr><td colspan="4"><hr></hr>
</td>
</tr>
<tr><td align="left"><bold><italic>Individual characteristics</italic>
</bold>
</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left">Age (completed years)</td>
<td></td>
<td></td>
<td align="left">.001</td>
</tr>
<tr><td align="left"> 15 – 24</td>
<td align="left">46.6 (640)</td>
<td align="left">53.4 (735)</td>
<td></td>
</tr>
<tr><td align="left"> 25 – 34</td>
<td align="left">56.0 (1159)</td>
<td align="left">44.0 (910)</td>
<td></td>
</tr>
<tr><td align="left"> 35+</td>
<td align="left">50.6 (1057)</td>
<td align="left">49.4 (1030)</td>
<td></td>
</tr>
<tr><td align="left">Education</td>
<td></td>
<td></td>
<td align="left">.001</td>
</tr>
<tr><td align="left"> No education</td>
<td align="left">35.3 (1013)</td>
<td align="left">64.7 (1856)</td>
<td></td>
</tr>
<tr><td align="left"> Primary</td>
<td align="left">57.1 (724)</td>
<td align="left">42.9 (543)</td>
<td></td>
</tr>
<tr><td align="left"> Secondary or higher</td>
<td align="left">80.2 (1119)</td>
<td align="left">19.8 (276)</td>
<td></td>
</tr>
<tr><td align="left">Ethnicity</td>
<td></td>
<td></td>
<td align="left">.001</td>
</tr>
<tr><td align="left"> Hausa/Fulani</td>
<td align="left">29.2 (606)</td>
<td align="left">70.8 (1467)</td>
<td></td>
</tr>
<tr><td align="left"> Igbo</td>
<td align="left">76.7 (578)</td>
<td align="left">23.3 (176)</td>
<td></td>
</tr>
<tr><td align="left"> Yoruba</td>
<td align="left">85.5 (541)</td>
<td align="left">14.5 (92)</td>
<td></td>
</tr>
<tr><td align="left"> Others</td>
<td align="left">54.7 (1122)</td>
<td align="left">45.3 (930)</td>
<td></td>
</tr>
<tr><td align="left">Religion</td>
<td></td>
<td></td>
<td align="left">.001</td>
</tr>
<tr><td align="left"> Christian</td>
<td align="left">69.9 (1593)</td>
<td align="left">38.1 (686)</td>
<td></td>
</tr>
<tr><td align="left"> Islam</td>
<td align="left">38.4 (1202)</td>
<td align="left">61.6 (1932)</td>
<td></td>
</tr>
<tr><td align="left"> Others</td>
<td align="left">52.6 (60)</td>
<td align="left">47.4 (54)</td>
<td></td>
</tr>
<tr><td align="left"><bold><italic>Community characteristics</italic>
</bold>
</td>
<td></td>
<td></td>
<td align="left">.001</td>
</tr>
<tr><td align="left">Type of residence</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> Urban</td>
<td align="left">60.7 (1247)</td>
<td align="left">39.3 (807)</td>
<td></td>
</tr>
<tr><td align="left"> Rural</td>
<td align="left">46.3 (1609)</td>
<td align="left">53.7 (1868)</td>
<td></td>
</tr>
<tr><td align="left">Community ethnic diversity<sup>†</sup>
, Mean (SD)</td>
<td align="left">3.5 (2.6)</td>
<td align="left">3.9 (2.5)</td>
<td align="left">.001</td>
</tr>
<tr><td align="left">Community median age of marriage, Mean (SD)</td>
<td align="left">17.5 (2.5)</td>
<td align="left">15.4 (1.9)</td>
<td align="left">.001</td>
</tr>
<tr><td align="left">Community median spousal age gap, Mean (SD)</td>
<td align="left">9.3 (2.3)</td>
<td align="left">10.2 (1.8)</td>
<td align="left">.001</td>
</tr>
<tr><td align="left">Geopolitical region</td>
<td></td>
<td></td>
<td align="left">.001</td>
</tr>
<tr><td align="left"> North-central</td>
<td align="left">60.4 (553)</td>
<td align="left">39.6 (363)</td>
<td></td>
</tr>
<tr><td align="left"> North-east</td>
<td align="left">38.1 (523)</td>
<td align="left">61.9 (750)</td>
<td></td>
</tr>
<tr><td align="left"> North-west</td>
<td align="left">32.0 (523)</td>
<td align="left">67.8 (1101)</td>
<td></td>
</tr>
<tr><td align="left"> South-east</td>
<td align="left">78.5 (442)</td>
<td align="left">21.5 (121)</td>
<td></td>
</tr>
<tr><td align="left"> South-south</td>
<td align="left">58.8 (314)</td>
<td align="left">41.2 (220)</td>
<td></td>
</tr>
<tr><td align="left"> South-west</td>
<td align="left">82.4 (562)</td>
<td align="left">17.6 (120)</td>
<td></td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>Chi-squared test was used for categorical data, and Mann-Whitney test continuous data for the bivariate analysis. For continuous variables, results are shown as mean and standard deviation.</p>
<p><sup>†</sup>
Number of ethnic group in a community, represents the likelihood that two persons, chosen at random from the same area, belong to different ethnic group</p>
<p>Data source – 2003 Nigeria Demographic and Survey. Sample size – 5531</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec><title>Exploratory spatial data analysis</title>
<p>Figure <xref ref-type="fig" rid="F1">1</xref>
shows geographic variation in state-specific median age of sexual initiation. The median age of sexual debut was lowest in North West (14 years) and highest in Kwara, Osun, and Abia states (19 years). Figure <xref ref-type="fig" rid="F2">2</xref>
presents a state map of the smoothed proportion of women who had reported early sexual debut. The smoothed proportions present clearer patterns and show clearly, where the problem was most severe. The spatial empirical bayes "smooth" estimates technique was able to deal with spatial heterogeneity. The estimation technique guarantees that estimates of neighboring states are more alike than estimates of states that are far away. North West and North East had the highest proportion of women who had reported early sexual debut (61% to 78%). The proportion of women who had reported early sexual intercourse was lowest in South West, South South, South East, and Kwara State (8% to 25%). The prevalence of early sexual debut was checked for spatial patterns by use of a Moran scatter plot. The results of this plot show statistically significant autocorrelation (Moran's I = 0.635, p = .0001) (Figure <xref ref-type="fig" rid="F3">3</xref>
). In addition, the bivariate Local LISA Moran's I statistics for spatial correlation between median age of marriage and spatial lag of median age of sex initiation is 0.646, (p = .001), indicating a significant positive spatial relationship between median age of marriage and median age of initiation of sexual intercourse. As expected, in the northern regions there is significant correlation between low median age of marriage and low median age of initiation of sexual intercourse neighbours (Figure <xref ref-type="fig" rid="F4">4</xref>
).</p>
<fig position="float" id="F1"><label>Figure 1</label>
<caption><p>Map showing median age of sexual initiation by states.</p>
</caption>
<graphic xlink:href="1476-072X-7-27-1"></graphic>
</fig>
<fig position="float" id="F2"><label>Figure 2</label>
<caption><p>Map showing the spatial empirical bayes smoothed proportions of women who reported early sexual initiation at state level.</p>
</caption>
<graphic xlink:href="1476-072X-7-27-2"></graphic>
</fig>
<fig position="float" id="F3"><label>Figure 3</label>
<caption><p>Local Indicator of Spatial Association (LISA) cluster map for Nigeria women age of initiation sexual intercourse.</p>
</caption>
<graphic xlink:href="1476-072X-7-27-3"></graphic>
</fig>
<fig position="float" id="F4"><label>Figure 4</label>
<caption><p>Bivariate LISA (BiLISA) cluster map for community median age of marriage and spatial lag of community median age of sex initiation in Nigeria, 2003.</p>
</caption>
<graphic xlink:href="1476-072X-7-27-4"></graphic>
</fig>
</sec>
<sec><title>Multilevel analyses</title>
<p>The result of the random effects model is shown in Table <xref ref-type="table" rid="T4">4</xref>
(Empty Model 1). There is significant variation in the log odds of early sex debut across the communities (<italic>τ </italic>
= 0.289, <italic>p </italic>
= .001) and states (<italic>τ </italic>
= 1.058, <italic>p </italic>
= .001). According to the intra-state correlation coefficient implied by the estimated intercept component variance, 22.8% and 6.0% of the variance in the age of initiation of sexual intercourse could be attributed to state-level and community-level variables respectively. This variation remained significant, even after controlling for individual-level factors (Model 2). However, the unexplained variations that remained after controlling for both individual and community-level characteristics were not significant (Model 3). As judged by the proportional change in variance, 79% and 95% of the variance in the log odds of early sexual debut variance across states was explained by individual compositional factors (Model 2) and both individual compositional and contextual factors (Model 3) respectively. 63% and 92% of the variance in the log odds of early sexual debut variance across communities was explained by individual compositional factors (Model 2) and both individual compositional and contextual factors (Model 3) respectively. In addition, the deviance information criterion (DIC) – was significant to reveal that the individual compositional factors (Model 2) and individual compositional and contextual variables (Model 3) increased the multivariable multilevel model's ability to explain variation in the log odds of early sexual debut, as indicated by lower DIC.</p>
<table-wrap position="float" id="T4"><label>Table 4</label>
<caption><p>Individual compositional and community contextual factors associated with early sexual debut among ever or currently marriage Nigeria women identified by multivariable multilevel logistic regressions<sup>‡</sup>
</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><td></td>
<td align="left"><bold><italic>Empty model</italic>
</bold>
</td>
<td align="left"><bold><italic>Individual variables</italic>
</bold>
</td>
<td align="left"><bold><italic>Individual and community variables</italic>
</bold>
</td>
</tr>
<tr><td></td>
<td colspan="3"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left"><bold>Model 1</bold>
</td>
<td align="left"><bold>Model 2</bold>
</td>
<td align="left"><bold>Model 3</bold>
</td>
</tr>
<tr><td></td>
<td colspan="3"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left"><bold>OR (95% CI)</bold>
</td>
<td align="left"><bold>OR (95% CI)</bold>
</td>
<td align="left"><bold>OR (95% CI)</bold>
</td>
</tr>
</thead>
<tbody><tr><td align="left"><bold>Fixed effects</bold>
</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"><bold><italic>Individual characteristics</italic>
</bold>
</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left">Age (completed years)</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> 15 – 24</td>
<td></td>
<td align="left">1</td>
<td align="left">1</td>
</tr>
<tr><td align="left"> 25 – 34</td>
<td></td>
<td align="left">1.02 (0.87, 1.21)</td>
<td align="left">1.05 (0.89, 1.24)</td>
</tr>
<tr><td align="left"> 35+</td>
<td></td>
<td align="left">1.19 (1.02, 1.41)*</td>
<td align="left">1.24 (1.05, 1.47)**</td>
</tr>
<tr><td align="left">Education</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> No education</td>
<td></td>
<td align="left">1</td>
<td align="left">1</td>
</tr>
<tr><td align="left"> Primary</td>
<td></td>
<td align="left">0.69 (0.59, 0.80)***</td>
<td align="left">0.74 (0.62, 0.88)***</td>
</tr>
<tr><td align="left"> Secondary or higher</td>
<td></td>
<td align="left">0.26 (0.22, 0.70)***</td>
<td align="left">0.30 (0.25, 0.36)***</td>
</tr>
<tr><td align="left">Ethnicity</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> Hausa/Fulani</td>
<td></td>
<td align="left">1</td>
<td align="left">1</td>
</tr>
<tr><td align="left"> Igbo</td>
<td></td>
<td align="left">0.35 (0.23, 0.52)***</td>
<td align="left">0.53 (0.34, 0.80)**</td>
</tr>
<tr><td align="left"> Yoruba</td>
<td></td>
<td align="left">0.22 (0.14, 0.33)***</td>
<td align="left">0.32 (0.21, 0.48)***</td>
</tr>
<tr><td align="left"> Others</td>
<td></td>
<td align="left">0.56 (0.45, 0.70)***</td>
<td align="left">0.63 (0.51, 0.78)***</td>
</tr>
<tr><td align="left">Religion</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> Islam</td>
<td></td>
<td align="left">1</td>
<td align="left">1</td>
</tr>
<tr><td align="left"> Christian</td>
<td></td>
<td align="left">0.77 (0.61, 0.96)*</td>
<td align="left">0.81 (0.64, 1.02)</td>
</tr>
<tr><td align="left"> Others</td>
<td></td>
<td align="left">1.02 (0.65, 1.59)</td>
<td align="left">0.92 (0.60, 1.45)</td>
</tr>
<tr><td align="left"><bold><italic>Community characteristics</italic>
</bold>
</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left">Type of residence</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> Urban</td>
<td></td>
<td></td>
<td align="left">1.10 (0.95, 1.25)</td>
</tr>
<tr><td align="left"> Rural</td>
<td></td>
<td></td>
<td align="left">1</td>
</tr>
<tr><td align="left">Community ethnic diversity<sup>†</sup>
</td>
<td></td>
<td></td>
<td align="left">0.98 (0.95, 1.01)</td>
</tr>
<tr><td align="left">Community median age of marriage</td>
<td></td>
<td></td>
<td align="left">0.77 (0.77, 0.81)***</td>
</tr>
<tr><td align="left">Community median spousal age gap</td>
<td></td>
<td></td>
<td align="left">1.01 (0.98, 1.05)</td>
</tr>
<tr><td align="left">Geopolitical region</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> North-central</td>
<td></td>
<td></td>
<td align="left">0.44 (0.31, 0.64)***</td>
</tr>
<tr><td align="left"> North-east</td>
<td></td>
<td></td>
<td align="left">0.40 (0.24, 0.66)***</td>
</tr>
<tr><td align="left"> North-west</td>
<td></td>
<td></td>
<td align="left">0.48 (0.30, 0.77)**</td>
</tr>
<tr><td align="left"> South-east</td>
<td></td>
<td></td>
<td align="left">0.43 (0.24, 0.75)**</td>
</tr>
<tr><td align="left"> South-south</td>
<td></td>
<td></td>
<td align="left">1</td>
</tr>
<tr><td align="left"> South-west</td>
<td></td>
<td></td>
<td align="left">0.57 (0.34, 0.96)*</td>
</tr>
<tr><td align="left">Constant</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"><bold>Random effects</bold>
</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left">State-level</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> Variance (SE)</td>
<td align="left">1.058 (0.289)***</td>
<td align="left">0.220 (0.084)**</td>
<td align="left">0.057 (0.03)</td>
</tr>
<tr><td align="left"> VPC (%)</td>
<td align="left">22.8</td>
<td align="left">6.1</td>
<td align="left">1.7</td>
</tr>
<tr><td align="left"> Explained variation (%)</td>
<td align="left">Reference</td>
<td align="left">79.2</td>
<td align="left">94.6</td>
</tr>
<tr><td align="left">Community-level</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> Variance (SE)</td>
<td align="left">0.280 (0.057)***</td>
<td align="left">0.104 (0.039)**</td>
<td align="left">0.021 (0.017)</td>
</tr>
<tr><td align="left"> VPC (%)</td>
<td align="left">6.0</td>
<td align="left">2.9</td>
<td align="left">0.6</td>
</tr>
<tr><td align="left"> Explained variation (%)</td>
<td align="left">Reference</td>
<td align="left">62.8</td>
<td align="left">92.5</td>
</tr>
<tr><td align="left"><bold>Model fit statistics</bold>
</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> DIC</td>
<td align="left">6543</td>
<td align="left">6255</td>
<td align="left">6219</td>
</tr>
<tr><td align="left"> Sample size</td>
<td></td>
<td></td>
<td></td>
</tr>
<tr><td align="left"> State-level</td>
<td align="left">36</td>
<td align="left">36</td>
<td align="left">36</td>
</tr>
<tr><td align="left"> Community-level</td>
<td align="left">177</td>
<td align="left">177</td>
<td align="left">177</td>
</tr>
<tr><td align="left"> Individual-level</td>
<td align="left">5531</td>
<td align="left">5508</td>
<td align="left">5508</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>Abbreviations: OR – Odds ratio, CI – Confidence interval, VPC – Variance Partition Coefficient, SE – Standard error, DIC – Deviance information criterion</p>
<p>†Number of ethnic group in a community, represents the likelihood that two persons, chosen at random from the same area, belong to different ethnic group</p>
<p>*p < .05, **p < .01, and ***p < .001</p>
<p><sup>‡</sup>
Data source: 2003 Nigerian Demographic and Health Survey</p>
</table-wrap-foot>
</table-wrap>
<p>The results of fitting the model including individual-level variables appear in Table <xref ref-type="table" rid="T4">4</xref>
(Model 2). The log odds of early sexual debut increased with increasing age and decreased with increasing level of education. Christian women were less likely to have started sex at early age. Compared to their counterparts that are Hausa or Fulani, Igbo and Yoruba women were 65% and 88% less likely to have reported early initiation of sexual intercourse respectively. The effect of the inclusion of contextual factors is shown in Table <xref ref-type="table" rid="T4">4</xref>
(Model 3). Inclusion of the community-level variables had minimal effect on the contribution of individual-level variables to the likelihood of starting sex at an early age. After individual compositional and contextual factors were added, regional dummies and community median age of marriage were statistically significant with early sexual debut. Independent of other factors, higher community median age of marriage was significantly associated with a decreased likelihood of started sex at early age, with 26% lower likelihood with every one-year increase in the community median age of marriage.</p>
</sec>
</sec>
<sec><title>Discussion</title>
<sec><title>Main findings</title>
<p>Using an explicit multilevel analytic framework and exploratory spatial data analysis, the study has shown that both individual composition and community contextual characteristics are important predictors of women's age at first sexual intercourse, and demonstrates geographic variation in rates of early initiation of sex in Nigeria. After adjusting for both individual-level and contextual factors, the probability of starting sex at an earlier age was associated with respondents' current age, educational attainment, ethnicity, region, and community median age of marriage. This paper has offered an alternative to more traditional ways of thinking about the distribution of sexual behaviour and sexual risk at the population level. The potential of multilevel modelling and exploratory spatial data analysis for exploring population contextual influences on individual sexual behaviour is clearly established.</p>
<p>Nigeria is made up of six major geopolitical regions. It is ethnically and religiously diverse and economic development and education levels vary widely across the country[<xref ref-type="bibr" rid="B20">20</xref>
]. Not unexpectedly, this study found that Nigerian women's age of initiation of sexual intercourse varies widely by region. Specifically, traditional and cultural beliefs that seem to encourage early marriage and lead to early initiation of sexual intercourse may be principal reasons for the observed geographical disparities. The study found that women living in neighborhoods with low median age of marriage had increase odds of reporting early initiation of sexual intercourse than those in neighborhoods with higher median age of marriage. This is consistent with the preponderance of evidence suggesting that Northern regions of Nigeria have some of the highest rates of early marriage in the world.</p>
</sec>
<sec><title>Study strengths and Limitations</title>
<p>There are a number of caveats to be considered when interpreting these results. The cross-sectional nature of the data limits ability to draw casual inferences. For this study, regions were administratively defined. The main drawback of this is that these boundaries may not match residents' perceptions and lived experiences of their "regional" boundaries. Another important limitation is validity constraints of self-reported sexual activity. By definition, self-reports of initiation of sexual activity cannot be externally validated and studies have revealed considerable inconsistencies in individual's self-reported sexual activity [<xref ref-type="bibr" rid="B21">21</xref>
,<xref ref-type="bibr" rid="B22">22</xref>
]. According to a study of eight African Demographic and Health Surveys by Meekers [<xref ref-type="bibr" rid="B22">22</xref>
], recall errors and misinterpretation of questions about first sexual intercourse may be responsible for misreporting of sexual behaviour.</p>
<p>Instead of reporting the age at which they first ever had sexual intercourse, some women report the age at which they first had sex with either their first or current husband [<xref ref-type="bibr" rid="B22">22</xref>
]. Lauritsen et al[<xref ref-type="bibr" rid="B21">21</xref>
] found that adolescents may under report sexual activity, depending on how recently that activity took place relative to the timing of the interview. Recalling age at first intercourse has been found to be significantly related to race, gender and family structure [<xref ref-type="bibr" rid="B21">21</xref>
]. In addition to standard limitations such as self-reported data and recall bias, the study was limited in that an existing data set used was not designed specifically to answer the research questions. Finally, spatial clustering methods are exploratory tools that help researchers and policy makers make sense of complex geographic patterns. Once a cluster is identified, only detailed epidemiological investigation can determine if the cluster is a "random" event. Therefore, it is essential that statistical significant clusters be examined in more detail. Despite these limitations, the study strengths are significant. It relies on a large, population-based study with national coverage. Individual compositional and community contextual characteristic accounted for 95% and 92% of the state and community disparities in early sexual debut, illustrating that the models were effective in predicting the risk of early sexual debut, and inclusion of contextual factors eliminated the significance of the state- and community-level variance (<bold>Model 3</bold>
).</p>
</sec>
<sec><title>Policy implications</title>
<p>The study findings have some important and relevant policy messages. The findings demonstrate that both individual compositional and community contextual factors are key factors that predict adolescent's risk of becoming sexually active at an earlier age. Better identifying and understanding the interplay between these individual-level characteristics and a variety of social and cultural contexts can provide valuable information for planners who determine the allocation of scarce resources. In addition, understanding of regional differences may aid in the identification of regions that may need to be particularly targeted with education and prevention programs. Multifaceted geographically differentiated intervention may represent a potentially effective approach for addressing issues related to adolescent sexuality. Program components could include: youth development, media campaigns, educational interventions, male involvement initiatives and legislative reform efforts.</p>
</sec>
</sec>
<sec><title>Conclusion</title>
<p>This first known national multilevel and spatial study of age at first sexual intercourse found strong evidence that women's odds of starting sexual intercourse early was significantly associated with respondents' current age, education attainment, religion affiliation, ethnicity, geopolitical region, and community median age of marriage. The study also shows that substantial geographical variations in rates of early initiation of sex exist within Nigeria. The findings suggest that multifaceted geographically differentiated interventions may represent a potentially effective approach for addressing issues related to adolescent sexuality. Scholars trying to understand variation in individual high risk sexual behaviour should pay attention to the characteristics of both individuals and places of residence.</p>
</sec>
<sec><title>Competing interests</title>
<p>The author declares that they have no competing interests.</p>
</sec>
<sec><title>Authors' contributions</title>
<p>OAU conceived the study, extracted the data, did the analyses and interpretation, and wrote the first and final draft of the manuscript.</p>
</sec>
</body>
<back><ack><sec><title>Acknowledgements</title>
<p>Without wishing to implicate them in any way, the author is grateful to the editor and two anonymous referees for their helpful comments on an earlier version of the paper. The author is grateful to Measure DHS for providing them with the 2003 Nigeria DHS data.</p>
</sec>
</ack>
<ref-list><ref id="B1"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Bruce</surname>
<given-names>J</given-names>
</name>
<name><surname>Shelley</surname>
<given-names>C</given-names>
</name>
</person-group>
<source>Including married adolescents in adolescent reproductive health and HIV/AIDS policy</source>
<year>2003</year>
<publisher-name>Geneva , World Health Organization</publisher-name>
</citation>
</ref>
<ref id="B2"><citation citation-type="book"><person-group person-group-type="author"><collab>Population Council</collab>
</person-group>
<source>Nigeria: child Marriage briefing</source>
<year>2005</year>
<publisher-name>New York , Population council</publisher-name>
</citation>
</ref>
<ref id="B3"><citation citation-type="book"><person-group person-group-type="author"><collab>UNAIDS</collab>
</person-group>
<source>Report on the Global AIDS Epidemic</source>
<year>2006</year>
<publisher-name>Geneva , UNAIDS</publisher-name>
</citation>
</ref>
<ref id="B4"><citation citation-type="book"><person-group person-group-type="author"><collab>United Nations Children's Fund</collab>
</person-group>
<source>Early marriage: child spouses</source>
<year>2001</year>
<publisher-name>Florence, Italy , UNICEF</publisher-name>
</citation>
</ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Andersson-Ellstrom</surname>
<given-names>A</given-names>
</name>
<name><surname>Forssman</surname>
<given-names>L</given-names>
</name>
<name><surname>Milsom</surname>
<given-names>I</given-names>
</name>
</person-group>
<article-title>Age of sexual debut related to life-style and reproductive health factors in a group of Swedish teenage girls</article-title>
<source>Acta Obstet Gynecol Scand</source>
<year>1996</year>
<volume>75</volume>
<fpage>484</fpage>
<lpage>489</lpage>
<pub-id pub-id-type="pmid">8677776</pub-id>
<pub-id pub-id-type="doi">10.3109/00016349609033359</pub-id>
</citation>
</ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rink</surname>
<given-names>E</given-names>
</name>
<name><surname>Tricker</surname>
<given-names>R</given-names>
</name>
<name><surname>Harvey</surname>
<given-names>SM</given-names>
</name>
</person-group>
<article-title>Onset of sexual intercourse among female adolescents: the influence of perceptions, depression, and ecological factors</article-title>
<source>J Adolesc Health</source>
<year>2007</year>
<volume>41</volume>
<fpage>398</fpage>
<lpage>406</lpage>
<pub-id pub-id-type="pmid">17875466</pub-id>
<pub-id pub-id-type="doi">10.1016/j.jadohealth.2007.04.017</pub-id>
</citation>
</ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kinsman</surname>
<given-names>SB</given-names>
</name>
<name><surname>Romer</surname>
<given-names>D</given-names>
</name>
<name><surname>Furstenberg</surname>
<given-names>FF</given-names>
</name>
<name><surname>Schwarz</surname>
<given-names>DF</given-names>
</name>
</person-group>
<article-title>Early sexual initiation: the role of peer norms</article-title>
<source>Pediatrics</source>
<year>1998</year>
<volume>102</volume>
<fpage>1185</fpage>
<lpage>1192</lpage>
<pub-id pub-id-type="pmid">9794952</pub-id>
<pub-id pub-id-type="doi">10.1542/peds.102.5.1185</pub-id>
</citation>
</ref>
<ref id="B8"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Anselin</surname>
<given-names>L</given-names>
</name>
</person-group>
<person-group person-group-type="editor"><name><surname>Longley PA, Goodchild MF, Maguire DJ, Rhind DW</surname>
</name>
</person-group>
<article-title>Interactive techniques and exploratory spatial data analysis</article-title>
<source>Geographical Information Systems: Principles, Techniques, Management and Applications</source>
<year>1999</year>
<publisher-name>New York , John Wiley</publisher-name>
<fpage>251</fpage>
<lpage>264</lpage>
</citation>
</ref>
<ref id="B9"><citation citation-type="other"><person-group person-group-type="author"><collab>Advocacy for youths</collab>
</person-group>
<article-title>Adolescent Reproductive Health in Nigeria: The facts</article-title>
<ext-link ext-link-type="uri" xlink:href="http://www.advocatesforyouth.org/PUBLICATIONS/factsheet/fsnigeria.htm"></ext-link>
</citation>
</ref>
<ref id="B10"><citation citation-type="other"><article-title>Nigeria: DHS, 2003 - Final Report</article-title>
<ext-link ext-link-type="uri" xlink:href="http://www.measuredhs.com/pubs/pdf/FR148/00FrontMatter.pdf"></ext-link>
</citation>
</ref>
<ref id="B11"><citation citation-type="book"><person-group person-group-type="author"><collab>StataCorp</collab>
</person-group>
<source>Stata Statistical Software: Release 10</source>
<year>2007</year>
<publisher-name>College Station, TX , StataCorp LP</publisher-name>
</citation>
</ref>
<ref id="B12"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Snijders</surname>
<given-names>TAB</given-names>
</name>
<name><surname>Bosker</surname>
<given-names>RJ</given-names>
</name>
</person-group>
<source>Multilevel analysis - an introduction to basic and advanced multilevel modeling</source>
<year>1999</year>
<publisher-name>Thousand Oaks, California , SAGE Publications</publisher-name>
</citation>
</ref>
<ref id="B13"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Rasbash</surname>
<given-names>J</given-names>
</name>
<name><surname>Steele</surname>
<given-names>F</given-names>
</name>
<name><surname>Browne</surname>
<given-names>W</given-names>
</name>
<name><surname>Prosser</surname>
<given-names>B</given-names>
</name>
</person-group>
<source>A user's guide to MLwiN. Version 2.0.</source>
<year>2004</year>
<publisher-name>London , Centre of Multilevel Modelling, Institute of Education, University of London</publisher-name>
</citation>
</ref>
<ref id="B14"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Browne</surname>
<given-names>W</given-names>
</name>
</person-group>
<source>MCMC estimation in MLwiN. Version 2.0</source>
<year>2003</year>
<publisher-name>London , University of London, Centre of Multilevel Modelling, Institute of Education</publisher-name>
</citation>
</ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Spiegelhalter</surname>
<given-names>DJ</given-names>
</name>
<name><surname>Best</surname>
<given-names>NG</given-names>
</name>
<name><surname>Carlin</surname>
<given-names>BP</given-names>
</name>
<name><surname>van der Linde</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Bayesian measures of model complexity and fit</article-title>
<source>J R Stat Soc B</source>
<year>2002</year>
<volume>64</volume>
<fpage>583</fpage>
<lpage>639</lpage>
<pub-id pub-id-type="doi">10.1111/1467-9868.00353</pub-id>
</citation>
</ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anselin</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Local indicators of spatial association – LISA</article-title>
<source>Geographical Analysis</source>
<year>1995</year>
<volume>27</volume>
<fpage>93</fpage>
<lpage>115</lpage>
</citation>
</ref>
<ref id="B17"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Anselin</surname>
<given-names>L</given-names>
</name>
<name><surname>Syabri</surname>
<given-names>I</given-names>
</name>
<name><surname>Smirnov</surname>
<given-names>O</given-names>
</name>
</person-group>
<source>Visualizing Multivariate Spatial
Correlation with Dynamically Linked Windows</source>
<year>2002</year>
<publisher-name>Urbana, IL , University of Illinois, Urbana-Campaign</publisher-name>
</citation>
</ref>
<ref id="B18"><citation citation-type="book"><person-group person-group-type="author"><name><surname>Anselin</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>MicroQAP: a microcomputer implementation of generalized measures of spatial association</article-title>
<source>In Working Paper</source>
<year>1986</year>
<publisher-name> University of California, Santa Barbara, Department of Geography</publisher-name>
</citation>
</ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anselin</surname>
<given-names>L</given-names>
</name>
<name><surname>Syabri</surname>
<given-names>I</given-names>
</name>
<name><surname>Kho</surname>
<given-names>Y</given-names>
</name>
</person-group>
<article-title>GeoDa: An introduction to spatial data analysis</article-title>
<source>Geographical Analysis</source>
<year>2006</year>
<volume>38</volume>
<fpage>5</fpage>
<lpage>22</lpage>
<pub-id pub-id-type="doi">10.1111/j.0016-7363.2005.00671.x</pub-id>
</citation>
</ref>
<ref id="B20"><citation citation-type="other"><article-title>Early Childbearing in Nigeria: A Continuing Challenge</article-title>
<ext-link ext-link-type="uri" xlink:href="http://www.guttmacher.org/pubs/rib/2004/12/10/rib2-04.pdf"></ext-link>
</citation>
</ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lauritsen</surname>
<given-names>JL</given-names>
</name>
<name><surname>Swicegood</surname>
<given-names>CG</given-names>
</name>
</person-group>
<article-title>The consistency of self-reported initiation of sexual activity</article-title>
<source>Fam Plann Perspect</source>
<year>1997</year>
<volume>29</volume>
<fpage>215</fpage>
<lpage>221</lpage>
<pub-id pub-id-type="pmid">9323498</pub-id>
<pub-id pub-id-type="doi">10.2307/2953398</pub-id>
</citation>
</ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Meekers</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Immaculate conceptions in sub-saharan Africa: exploratory analysis of inconsistencies in the timing of first sexual intercourse and first birth</article-title>
<source>Social biology</source>
<year>1995</year>
<volume>42</volume>
<fpage>151</fpage>
<lpage>161</lpage>
<pub-id pub-id-type="pmid">8738542</pub-id>
</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 000D21 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000D21 | 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é= |texte= }}
![]() | This area was generated with Dilib version V0.6.32. | ![]() |