Small-area estimation and analysis of HIV/AIDS indicators for precise geographical targeting of health interventions in Nigeria. a spatial microsimulation approach
International Journal of Health Geographics volume 22, Article number: 23 (2023)
Precise geographical targeting is well recognised as an indispensable intervention strategy for achieving many Sustainable Development Goals (SDGs). This is more cogent for health-related goals such as the reduction of the HIV/AIDS pandemic, which exhibits substantial spatial heterogeneity at various spatial scales (including at microscale levels). Despite the dire data limitations in Low and Middle Income Countries (LMICs), it is essential to produce fine-scale estimates of health-related indicators such as HIV/AIDS. Existing small-area estimates (SAEs) incorporate limited synthesis of the spatial and socio-behavioural aspects of the HIV/AIDS pandemic and/or are not adequately grounded in international indicator frameworks for sustainable development initiatives. They are, therefore, of limited policy-relevance, not least because of their inability to provide necessary fine-scale socio-spatial disaggregation of relevant indicators.
The current study attempts to overcome these challenges through innovative utilisation of gridded demographic datasets for SAEs as well as the mapping of standard HIV/AIDS indicators in LMICs using spatial microsimulation (SMS).
The result is a spatially enriched synthetic individual-level population of the study area as well as microscale estimates of four standard HIV/AIDS and sexual behaviour indicators. The analysis of these indicators follows similar studies with the added advantage of mapping fine-grained spatial patterns to facilitate precise geographical targeting of relevant interventions. In doing so, the need to explicate socio-spatial variations through proper socioeconomic disaggregation of data is reiterated.
In addition to creating SAEs of standard health-related indicators from disparate multivariate data, the outputs make it possible to establish more robust links (even at individual levels) with other mesoscale models, thereby enabling spatial analytics to be more responsive to evidence-based policymaking in LMICs. It is hoped that international organisations concerned with producing SDG-related indicators for LMICs move towards SAEs of such metrics using methods like SMS.
The intensity of the HIV/AIDS pandemic, its key drivers, and the outcomes of public interventions all exhibit significant spatial heterogeneity, not only at large regional scales, but also at microscale levels like towns, communities and neighbourhoods . This motivates the overwhelming emphasis on precise geographical targeting by several authors and organisations, including the Joint United Nations Programme on HIV/AIDS (UNAIDS) [2,3,4]. The lack of small-area estimates (SAEs) of standard HIV-related indicators in LMICs greatly undermines the possibility of tracking trends and patterns of HIV-related metrics at fine spatial scales, as is typically done for large regions like countries . This is particularly problematic in LMICs like Nigeria, which suffer from a scarcity of geocoded data amidst the high prevalence rate of HIV/AIDS. A clear implication is that regardless of overall records of positive outcomes for large regions, the localities therein (particularly disadvantaged ones) not only may be suffering a much slower pace of improvement, but also could be experiencing deteriorating outcomes. In addition to greatly hampering the spatial allocative efficiency of public resources for the control of the HIV/AIDS pandemic, the dearth of fine-scale metrics also limits the possibility of determining contextual drivers of these indicators. Consequently, the differential effects of intervention programmes upon various localities remain obscure. In practical terms, the ability to obtain granular estimates of standard HIV-related indicators in this paper also means that a wide variety of robust spatial analytics can be implemented in LMICs using a realistic spatially enriched synthetic population.
Even though there has been a recent upsurge in the number of studies concerned with fine-scale mapping of HIV-related indicators , the vast majority of these works do not adequately synthesise locational and social (population at risk) aspects in a unified framework. This owes largely to limitations regarding the methods employed, which are seldom as robust as the spatial microsimulation (SMS) employed in the present study. Furthermore, SAE of HIV/AIDS indicators are yet to benefit from the immense analytical power of SMS, which typically requires rich data that are scarce in LMICs [7, 8]. Nevertheless, this has been employed in the analysis of other health phenomena, including smoking, obesity, mental illness, alcohol consumption, diabetes, and healthcare access, especially in developed countries, which do not experience the same data limitations as in the present study context [10,11,12,13,14,15]. Similar indicators and data have been used by other studies, except that these are not disaggregated at small-area scales like in the present study [16,17,18,19,20]. Mweemba et al.  is a recent similar attempt at estimating the HIV prevalence in Africa (Zambia); however, the SAE method employed is not as robust as SMS. Overall, many of these studies are not based on standard SDG-related estimates such as the Multiple Indicator Cluster Survey (MICS), Millennium Development Goals (MDG) or Sustainable Development Goals (SDG) Indicators framework [21,22,23]. This is a major limitation to their reproducibility and comparison in different international contexts as well as their utility for the monitoring of sustainable development initiatives.
To advance ongoing research in this area, there are two (2) main objectives of this paper. The first objective is to simulate a spatially enriched synthetic population by linking survey microdata with small-area zoning systems using SMS. The second is to utilise the derived spatially enriched synthetic individual-level population to estimate and analyse four SDG-related multivariate indicators of HIV/AIDS at small-area levels. These indicators are (1) sex with non-regular partners (MICS Indicator 9.14), (2) condom use with non-regular partners (MICS Indicator 9.15 or MDG Indicator 6.2), (3) knowledge about HIV prevention among young persons (MICS Indicator 9.1 or MDG Indicator 6.3) and (4) sexually active young persons who have been tested for HIV and know the results (MICS Indicator 9.6). In the MICS indicators framework, these are amongst the 17 sub-indicators which capture traits of HIV/AIDS and sexual behaviour. These standard multivariate MDG or SDG indicators of HIV/AIDS and sexual behaviour are estimated and analysed for small areas with a view to facilitating precise geographical targeting of HIV-related interventions [1, 24,25,26]. Not only is the estimation of HIV/AIDS and sexual behaviour indicators in this study premised on the framework of social determinants of sexual behaviour and sexual health [25, 75], these are also amongst the internationally adopted indicators for monitoring HIV/AIDS the world over.
The MICS indicators framework which is congruent with relevant SDG indicators, is operationalised in this study using the fifth wave of MICS microdata collected by UNICEF in 2016/2017. This also has the benefit of temporality and comparability for reproducing similar research both with historical MICS data (since the 1990s) and for the 118 countries currently being covered. About 130 indicators are collected by MICS to analyse the situation of children, women and men in the areas of health, education and child protection. With roughly 33 SDG indicators, MICS is a rich data source for analysing aspects of the SDG including HIV/AIDS of this study. At the time of writing, the fifth wave was the latest MICS microdata of the study area available for public download on UNICEF’s website (https://mics.unicef.org/surveys). There are indications that MICS 5 is representative at the subnational levels like the study area, Kogi State (see https://mics.unicef.org/faq#sampling). This explains why the official MICS 5 report of Nigeria disaggregates several statistics at the state level . If subnational units of the MICS data were not representative, we suppose they would not be appropriate grouping/clustering levels for use in multilevel spatial analysis  or other types of sub-national spatial analysis [29, 30]. Aside from these, a cardinal objective of SMS is to mitigate the problem of data sparsity (of sample surveys) through data enrichment by synthesizing them with spatial constraints that contain the entire study population in aggregate form. Thus, sample size limitations of surveys are of minimal concern for this type of analysis since this is the challenge that SMS is meant to overcome.
As illustrated in Fig. 1, this study follows the standard SMS procedure for developing a spatially enriched synthetic population and deriving SAEs of relevant multivariate outcomes such as health-relevant indicators . ArcGIS Desktop 10.5 and IBM SPSS v26 were used for preparing the microdata and spatial constraints, while Flexible Modelling Framework (FMF) v1.3 [32, 33] was used to simulate a generalised synthetic micro-population of the study area, which was then used for SAEs of relevant MICS 5 indicators. FMF software is mature SMS software  as well as the most user-friendly of the available options. This implements the combinatorial optimization approach which has been adjudged to outperform many competing techniques, including Iterative Proportional Fitting (IPF) and Generalised Regression Weighting (GREGWT) [35, 36]. Other specialist software for implementing SMS include LIAM2 , the SMS package of R  and JAS-mine .
Selecting and preparing spatial constraints
The selection of appropriate spatial constraints can be determined via regression analysis, analysis of intra-area homogeneity, literature reviews, and consultation with end users [36, 40, 41]. Where the outcomes of an SMS operation are multivariate, such as in the present study, it is essential to include key demographic attributes as spatial constraints, especially in equity-sensitive projects . This can be informed by the framework of PROGRESS stratifiers, which represents the internationally adopted standard for statistical data disaggregation in SDG-related studies [43, 44]. However, the availability and quality of relevant datasets for the study area also limits the selection of variables for this study.
All of the spatial constraints and sub-constraints in this study were derived from the following WorldPop grids: age/sex structure, education (per sex) (as illiteracy ratio) and poverty. From the age/sex structure grid, three (3) spatial constraints were produced, namely: age [5 year groupings], sex (female or male), and 5 year age groupings by sex. From the education grid, two (2) spatial constraints were derived, namely: education (educated or uneducated) and education by sex (i.e., educated female, uneducated female, educated male or uneducated male). The poverty grid produced one spatial constraint, poverty (i.e. poor or not poor). As required, these were aggregated to the small-area analytical zoning system of the study area using ArcGIS Desktop 10.5 software. The configuration of the variables used for linking the microdata with the spatial constraints of the SMS operation is outlined in Table 1.
Most processes of data preparation for SMS involve the categorisation of continuous variables or the recoding of variables into fewer classes. Though necessary for model execution, error mitigation, and model validity, data preparation often results in a loss of nuance and a consequent reduction in the robustness of the derived synthetic population.
Preparing the microdata/population
Only respondents from Kogi State are included in the microdata, although SMS allows the use of microdata from a wider geographical area (and collected at different times) than that which is being modelled [45, 46]. Relevant variables of the microdata were recoded accordingly to match the configuration of the available spatial constraints. For instance, since only two categories of socioeconomic status are provided by spatial constraints (namely Poor and Not Poor), the wealth index quintile of the microdata was regrouped into two corresponding categories. ‘Poorest’ (coded 1) and ‘Second’ (coded 2) were regrouped as ‘Poor’, while ‘Middle’ (coded 3), ‘Fourth’ (coded 4) and ‘Richest’ (coded 5) were regrouped as ‘Not Poor’. Similarly, the multiple levels of education present in the MICS 5 microdata were regrouped into two categories, namely ‘Educated’ and ‘Uneducated’, to match available spatial constraints. The selected variables and structure of recoded microdata used as linkage variables are outlined in Table 1. Only respondents aged 15–49 years with all of the relevant linking variables are included in the SMS operation in order to minimise modelling errors. Thus, the input MICS 5 microdata contain 1305 persons, comprising 912 females and 393 males.
Implementing the SMS of Kogi state population
SMS requires two (2) types of appropriately structured data input, namely: spatial constraint files and a microdata/population file. With these, the SMS operation in this paper was performed using FMF v1.3 software, following the detailed user guide of Harland . The same process can also be implemented in R following the comprehensive guidelines by Lovelace and Dumont .
The execution of the SMS within the FMF software environment is a computing-intensive process which produced a spatially enriched synthetic population of the study area, as expected. The synthetic small-area analytical zone to which each synthetic person belongs is also included as a column of the resulting table. These enabled the computation and mapping of SAEs of four (4) MICS 5 indicators of HIV/AIDS for young people aged 15–24 years in the study area. These are (1) comprehensive knowledge about HIV prevention (MICS 5 Indicator 9.1 or MDG Indicator 6.3), (2) sexually active young persons who have been tested for HIV and know the results (MICS 5 Indicator 9.6), (3) sex with non-regular partners in the last 12 months (MICS 5 Indicator 9.14), and (4) condom use with non-regular partners (MICS 5 Indicator 9.15).
The computation of these indicator estimates followed the standard MICS 5 indicator computation algorithms but was adapted so as to account for the new small-area geographies ascribed to the synthetic individual-level population. Furthermore, other relevant disaggregation variables not considered in the original algorithms were also included, particularly the ethnicity and religion of household heads. This entailed making requisite modifications to the standard MICS 5 indicator computation algorithms provided as SPSS syntax files (sps) on the MICS 5 website: http://mics.unicef.org/tools?round=mics5#data-processing. The goodness of fit of the derived synthetic population is examined in the next section, after which the SAEs derived thereof are discussed.
The spatially enriched synthetic micro-population
The result of the SMS operation in this study is a spatially enriched synthetic individual-level population of the study area comprising 2,249,170 microunits (i.e. persons of age 15–49 years). These are composed of 1,115,283 females and 1,133,887 males, with about 425 MICS 5 attributes. The synthetic small-area analytical zone to which each synthetic person belongs is also included as a column of the resulting table. These data are available at . With the derived spatially enriched synthetic population, it became possible to compute and map SAEs (of relevant MICS 5 indicators) for the optimised analytical zoning system of the study area.
Goodness of fit of the synthetic population: internal validation
It is crucial for a derived synthetic population to be validated before being utilised for further analytical steps . Not only does this ensure that a given SMS is of acceptable quality, it also ensures that potential analytical errors are not propagated unto further analytical operations. This explains why Step 4 of Fig. 1, which is the validation of a synthetic population, precedes ‘the creation of estimates at small-area levels’ (i.e. Step 5) as well as subsequent analytical steps in SMS. Internal validation is performed in this study by comparing the outputs of the SMS (i.e. the synthetic population) with the original inputs (i.e. the sample survey micro-population and the spatial constraints) to estimate the extent to which the original inputs are emulated by the resulting synthetic population [49, 50]. The degree of modelling accuracy is assessed using all twelve  goodness-of-fit statistics available in FMF software, namely: standard root mean square error (SRMSE), absolute entropy difference (AED), R-squared (R2), entropy (standard entropy measure), chi-squared (X2), total absolute error (TAE), standard absolute error (SAE), percentage error (PE), total error (TE), cell percentage error (CPE), Z2, and Z.
These goodness-of-fit statistics are produced for each spatial constraint and small-area zone, constituting a massive table of 342 rows by 73 columns. Presented in are the goodness-of-fit statistics for each of the spatial constraints considered, while Fig. 2 presents maps of two of these goodness-of-fit statistics (i.e. CPE and Z2) for three spatial constraints, namely: ‘age’, ‘age by sex’, and ‘sex by education’. Z-squared is chosen for mapping because of its popularity in cognate literature [51, 52], while CPE is relatively easy to explain. With higher values of SRMSE, AED, SAE, PE and CPE come more error in a model fit. The three spatial constraints mapped tend to have the lowest goodness-of-fit values, as shown in Table 2 (also highlighted in orange colour); therefore, they deserve further consideration. Nevertheless, these error metrics are extremely low, being less than 0.08. This implies very good internal validity of the derived synthetic micro-population of Kogi State.
Overall, from the computed goodness-of-fit statistics, the SMS of Kogi State shows an excellent level of accuracy across all of the spatial constraints utilised. For instance, all SAE (and TAE and TE) values are extremely low, with ‘sex’ exhibiting the best model fit and ‘sex by age’ having the lowest (albeit a negligible) model fit. It is unsurprising that ‘sex’ has the best model fit because it was specified as the reference spatial constraint used by FMF software for determining zone totals; thus, it is the reference benchmark by which the accuracy of the synthetic population in this study is adjudged. Any constraint can be chosen as the reference variable for internal validation; however, sex was chosen because it is bivariate and tends to constitute relatively accurate input data. This near-perfect SMS validity is exhibited by all of the other computed goodness-of-fit statistics. For instance, R2 values (i.e. coefficient of determination) typically range from 0 to 1, with 1 signifying a perfect model fit. Notice that the R2 for all of the constraint variables is approximately 1. Furthermore, a synthetic population is deemed to be unfit if the Z-squared statistic is greater than the critical value (i.e. |Z|> 1.96) . Clearly, no Z2 is above the critical value for the current SMS operation, thus confirming the validity of the synthetic population developed in this study.
Figure 2 shows the spatial patterns of CPE and Z-squared for three spatial constraints, namely: ‘age’, ‘age by sex’ and ‘age by education’. These three spatial constraints tend to have the lowest goodness-of-fit values; therefore, they are likely to reveal problem small areas as well as to show substantial spatial variations in mapped data of SMS errors. Like Table 2, Fig. 2 shows that ‘Age by Sex’ has the highest error of any small-area zone, with 16 zones having CPEs of between 5.42% and 10%. Age has the second-highest error, with 25 zones having CPEs of between 3.1% and 7.58%. Nevertheless, not only are these CPE values relatively low, only very few zones have these levels of error, that is, 16, 25 and 19 zones (out of the 341 analytical small-area zones in the study area), for the ‘sex by age’, ‘age’, and ‘sex by education’ spatial constraints, respectively. On the flipside, 45, 144 and 193 zones are 100% accurate in terms of CPE for the ‘Sex by Age’, ‘Age’, and ‘Sex by Education’ spatial constraints, respectively. These corroborate Table 2 in demonstrating the very high internal validity of the SMS in the current study. The resulting synthetic population is, therefore, fit for use in subsequent sections of this paper.
Small-area estimates of select HIV/AIDS indicators
From the derived synthetic micro-population of the study area, the four (4) SAEs of standard MICS 5 indicators are presented in this section as a prelude to the associated discussion in the subsequent section.
The spatial pattern of comprehensive HIV/AIDS knowledge among young persons aged 15–24 years (MICS 5 Indicator 9.1 and MDG Indicator 6.3) shown in Fig. 3A suggests that there are a relatively greater number of young persons within the middle horizontal strip of the study area who have comprehensive knowledge about HIV than of people in the upper and lower flanks. The small-area spatial pattern of HIV testing amongst sexually active young persons (MICS 5 Indicator 9.6) is shown in Fig. 3B. This suggests that there is more HIV testing amongst sexually active young persons on the western side of the River Niger (and more variations in the same region). The percentage of young persons who have tested for HIV in the last 12 months and know the results ranges from 3.6% to 12.4% in the study area. MICS 5 Indicators 9.1 and 9.6 are further disaggregated socioeconomically for the study area, as presented in Fig. 4. This shows that more young people in urban areas have comprehensive knowledge of HIV (about 12.6%) than have those in rural areas (about 5.9%). Similarly, in the last 12 months, there has been much more HIV testing amongst sexually active young persons in urban areas (about 25.4%) in comparison to rural areas (about 5.8%).
The percentage of multiple sexual partnerships among young people in Kogi State (MICS 5 Indicator 9.14) ranges from 31.4% to 48.4%, whereas 31.9% to 49.9% of this population had used a condom during their last sexual intercourse with a non-regular partner (i.e. MICS 5 Indicator 9.15). This is shown in Fig. 5A, B respectively. These are further disaggregated according to relevant socioeconomic markers, as shown inFig. 6Footnote 1 There is an indication that with higher levels of education come increased proportions of young people reporting having had sex with multiple sexual partners in the last 12 months. For instance, 64.9% of young people with higher education reported having had sex with a non-regular partner in the last 12 months, compared to 29.9% of young people with primary education reporting a similar practice. The disaggregation of MICS 5 Indicators 9.14 and 9.15 by marital status, sex, and urban/rural dwelling reveals other interesting patterns which are discussed in the next section.
HIV/AIDS knowledge and testing amongst young people (MICS 5 Indicators 9.1 and 9.6)
It is not certain as to why there are a relatively greater number of young persons within the middle horizontal strip of the study area who have comprehensive knowledge about HIV (MICS 5 Indicator 9.1 and MDG Indicator 6.3) than of people in the northern and southern flanks. However, better knowledge about HIV is likely to be associated with both urban dwelling and more education [54, 55]. More educational opportunities as well as urban dwelling tend to be concentrated in the same middle horizontal strip of the study area. Nevertheless, the percentage of young persons with comprehensive HIV knowledge in the study area is generally low, ranging from 4 to 9%. This implies a need to improve HIV-related education in the study area, particularly in the northern and southern flanks.
From the mapping of SAEs of MICS 5 Indicator 9.6 shown in Fig. 3B, not only is the percentage of young persons who have tested for HIV in the last 12 months and know the results very low (3.6% to 12.4%), places with exceptionally poor testing for HIV are also revealed. This indicates that while the study area will benefit from universal improvement in the coverage of HIV testing for sexually active young persons, localities with exceptionally low records of HIV testing (3.6% to 6.3%) should be targeted for prioritised provision of HIV testing services as well as awareness creation .
The suggested urban advantage (shown in Fig. 3) is plausible for a number of reasons, including relatively higher health literacy and better accessibility of HIV services in comparison to rural areas [56, 57]. This indicates that rural areas should be given prioritised attention for HIV/AIDS-related interventions. Furthermore, higher levels of both education and wealth quintiles are associated with remarkable increases in both Indicators 9.1 and 9.6. Together with urban residence, wealth and education are key markers of socioeconomic advantage. As these have been shown to be greatly associated with better knowledge of, as well as more testing for, HIV amongst young people, poor knowledge of HIV and a lack of HIV testing amongst sexually active young people can be considered matters of disadvantage in respect of socioeconomic status. This supports studies which show that even though highly educated and/or well-off people are more likely to engage in risky sexual behaviours (such as keeping multiple sexual partners), they also tend to both be more health-literate and to practise safe sex, thereby being less prone to HIV infection than are poor and/or less educated people [54, 55, 58]. Theoretically, there is a complex relationship between risky sexual behaviours and both schooling and wealth, which also varies by gender . Nevertheless, the spatial and social disaggregation of these indicators suggests socio-spatial inequalities that are worth investigating further in subsequent studies.
Multiple sexual partnerships and condom use amongst young people (MICS 5 Indicators 9.14 and 9.15)
Small-area estimates of MICS 5 Indicators 9.14 and 9.15 (shown in Fig. 5A, B respectively) indicate that even though a large proportion of the young population in the study area had had multiple sexual partners in the last 12 months, a good percentage of these had used a condom during their last sexual intercourse with a non-regular partner. However, the observed spatial patterns in Fig. 5 suggest an inverse relationship between multiple sexual partnerships and condom use amongst young persons in the study area. For instance, relative to other senatorial districts, a smaller proportion of young persons in Kogi Central have multiple sexual partners, while a greater percentage of them reported having used a condom during their last sexual intercourse with a non-regular partner. Conversely, Kogi East, which records a high proportion of multiple sexual partnerships among young people, is also the senatorial district with a relatively low proportion of condom use during their last sexual intercourse with a non-regular partner. This shows that although there is a need to improve sex education and HIV/AIDS-related healthcare/screening services in the study area, this is particularly more crucial in Kogi East to ensure that the higher proportion of multiple sexual partnerships are matched with safe-sex practices.
Despite the increased sexual partnerships associated with higher levels of education depicted in Fig. 5, these are matched with more reports of regular usage of condoms [56, 58]. In fact, 84.8% of this population had used a condom during their last sexual affair with a non-regular partner, compared to 12.6% of young people with primary education who reported having used a condom during a similar sexual encounter. This shows that although more educated young people tend to be more exploratory regarding having sex with multiple sexual partners, they are also more cautious in practising safe sex (by using a condom with non-regular partners) than are their less educated peers [57, 59]. This suggests a need for more sex education for young people with secondary education or lower. For young people aged 15 – 24 years, a key reason for increased multiple sexual partnerships with higher education is that time in higher education and beyond is often associated with more liberty/autonomy, being away from the restrictions of parents or guardians [60, 61]. Apparently, the educational level of young people is likely to be directly correlated with age. Results show that with increased age comes more tendency towards sex with multiple sexual partners amongst young people; however, a greater proportion of non-teenage young people (42.1% and 43.8% for ages 20–22 and 23–24 respectively) reported having used a condom during their last sexual affair with a non-regular partner in comparison to the proportion of teenage young people (23.5% and 30.6% for ages 15–17 and 18–19 respectively) reporting the usage of condoms. This is consistent with extant empirical literature, which suggests that awareness of safe-sex practices is positively related to age, education, and socioeconomic status . As with having a higher level of education, post-teenage young people are likely to be more autonomous than teenage young people, since they may no longer be subject to as many parental restrictions as those imposed on teenagers [60, 61]. Despite this, they tend to be twice as likely to practise safe sex as teenage young people. Thus, teenage young people should be prioritised for HIV-related interventions, especially on the need for safer sexual behaviour.
Other interesting patterns are also observable when Indicators 9.14 and 9.15 are disaggregated according to marital status, sex, and urban/rural dwelling. Being in non-committal relationships, many more single young people (44.7%) have had sex with multiple sexual partners in the last 12 months in comparison to their married peers (7.9%), as expected [63, 64]. Furthermore, amongst married young people who have had an extramarital sexual partner in the last 12 months, the large majority (72.2%) reported having used a condom, compared to the proportion of single young people who reported having used a condom with their non-regular sexual partner in the same period (36.8%). For cultural and religious reasons, married young people in the study area are much less likely to engage in extramarital sex than are their single peers . Moreover, having a married regular sexual partner means that most married young people are not likely to be keenly searching for new sexual partners. Whenever such extramarital affairs happen, condoms are often used to prevent both pregnancy of the non-regular (female) partner as well as the transmission of an STI/STD to their married spouse.
It is unsurprising that a greater proportion of young males have had sex with a non-regular sexual partner (48.0%) in comparison to the proportion of females with multiple sexual partners in the last 12 months (31.9%) [66,67,68]. This is because Nigeria, like many sub-Saharan African countries, is notoriously patriarchal, thus being culpable of hegemonic masculinity [69,70,71]. This concept develops on ideas of patriarchy in explaining entrenched patterns of social practices (including actions and expectations) that perpetuate male dominance of females, often facilitated by culture, institutions, and political influence . Consequently, while females are highly discouraged from keeping multiple sexual partners for cultural reasons, males do not experience the same levels of restrictions, even when married [73, 74]. Indeed, in many localities (including the study area), while it is taboo for married females to engage in extramarital sex, this is not the case for married males. It is ironic, however, that a higher proportion of married young males who reported having had sex with a non-regular partner in the last 12 months (48.4%) indicated having used a condom in such affairs, compared to the percentage of their female peers who reported having used a condom during their last sexual affair with a non-regular partner (21.0%) [75,76,77]. This may be because in the study area there is a tendency for (married) females who have sex with non-regular partners to do so for a variety of transactional reasons (including in exchange for gifts or other favours from men), in which case they are less able to negotiate for safe sex [78,79,80]. Consequently, young males are twice as likely to practise safe sex with a non-regular partner as females, thereby partly explaining why young females (aged 15–24 years) in sub-Saharan Africa are twice as likely as young males to have HIV . This suggests a need to promote safer-sex practices amongst young people, especially females, as well as to target empowerment interventions at women to make them less vulnerable .
In rural areas, more young people reported having had sex with a non-regular partner (41.1%) in the last 12 months, compared to their peers in urban areas who reported a similar practice (33.9%) . This suggests that young people in rural areas enjoy more autonomy than do their peers in urban centres, probably because of the very informal and communal nature of rural areas in the study area. With this comes an increased risk of HIV infection because people who become sexually active at a younger age are more prone to having multiple sexual partners in their lifetime, which is associated with a higher tendency towards indulging in other risky sexual behaviours [58, 63]. This is, however, slightly at odds with the patterns recorded based on disaggregation by educational level, as a greater proportion of young people in urban areas are expected to possess higher education in comparison to their peers in rural areas. Despite more reports of having had sex with non-regular partners amongst young people in rural areas, a much lower proportion of these people use condoms in comparison to their peers in urban areas . While 67.5% of young people with multiple sexual partners in urban areas had used a condom during their last sexual intercourse with a non-regular partner, only 33.3% of young people in rural areas reported the usage of condoms. Thus, young rural dwellers are twice as unlikely to practise safe sex with a non-regular partner than their urban peers. In addition to poorer access to condoms in rural areas, this could indicate lower levels of sex education and, by extension, less awareness of safe-sex practices than in urban areas [54, 57]. On the one hand, this suggests a need for increased sex education, condom accessibility, and the promotion of safe-sex practices in rural areas. On the other hand, being more vulnerable, rural areas should be targeted for increased accessibility to medical services related to the sexual health of young people, such as HIV/AIDS-related services as well as services for other STDs/STIs, such as relevant screenings/tests [1, 83].
Unlike purely statistical methods, SMS is arguably the only explicitly geographical method for SAEs of multivariate indicators, especially when there is a necessity to simultaneously incorporate multiple variables . From MICS 5 microdata, four (4) standard indicators that relate to HIV/AIDS and sexual behaviour are estimated (and mapped) at a small-area scale and comprehensively analysed. In tandem with a need to illuminate inequalities using the PROGRESS framework , these indicators are also disaggregated according to relevant socioeconomic dimensions in order to explicate pertinent socio-spatial variations . In addition to creating SAEs of standard health-related indicators from disparate multivariate data, the outputs make it possible to establish more robust links (even at individual levels) with other mesoscale models, thereby enabling spatial analytics to be more responsive to evidence-based policymaking in LMICs [53, 86,87,88]. With SMS, it is also possible to further link and co-analyse various other types of contextual information at multiple spatial scales with the rich socioeconomic data of sample surveys. What is more, this study provides a robust analytical framework for estimating and analysing other health—and/or SDG-related indicators at very granular scales in various data-sparce contexts through the use of a variety of publicly available sample survey microdata and gridded demographic datasets.
From anecdotal evidence, spatial patterns with regard to both ethnicity and religion are likely to be unreliable because none of the spatial constraints are known to have a significant correlation with these two in the study area. In other words, sex, age, wealth, and educational level are not known to be significantly correlated with ethnicity or religion in the study area. Future works would benefit from having either ethnicity or religion as a spatial constraint because in the study area any of these would serve as a good predictor of the other. This would go a long way in enhancing the fidelity of the resulting synthetic population in these dimensions. For some of the spatial constraints, such as wealth/poverty and literacy, only a binary categorisation was used instead of the full 4–5 categories in the MICS 5 sample survey microdata. This is a limitation of the gridded demographic data used in this project. New data sources offer opportunities to provide more detailed spatial constraints, which will enhance the model fidelity. For instance, comprehensive data on rural/urban classification are now available from the Global Human Settlement—Settlement Model Grid (GHS-SMOD) . Furthermore, future works should consider using a variety of innovative data sources to implement dynamic SMSs of standard SDG-related indicators beyond the health-related exemplars in the present study. This would benefit from the use of more robust SMS tools in geocomputation environments like Python and R. It is hoped that in the years ahead, the simulation of small-area metrics of SDG-related indicators will become the standard practice of international organisations concerned with data provision in LMICs.
Availability of data and materials
The datasets generated and/or analysed in this manuscript are available in the Harvard Dataverse repository, https://dataverse.harvard.edu/privateurl.xhtml?token=b47f9178-ce21-428b-94f0-757b2096c2aa
Indicator 9.15 is depicted relative to Indicator 9.14 such that the golden bars have percentages that are independent of the y-axis of the chart, but relative to the purple bars (of indicator 9.14) in order to be more meaningful and aid proper comparisons.
Anderson S-J, Ghys PD, Ombam R, Hallett TB. HIV prevention where it is needed most: comparison of strategies for the geographical allocation of interventions. J Int. 2017;20(4):e25020.
Kedziora DJ, Stuart RM, Pearson J, Latypov A, Dierst-Davies R, Duda M, et al. Optimal allocation of HIV resources among geographical regions. BMC Public Health. 2019;19(1):1509.
Grantham KL, Kerr CC, Wilson DP. Local responses to local epidemics for national impact need advanced spatially explicit tools. AIDS. 2016;30(9):1981.
Meyer-Rath G, McGillen JB, Cuadros DF, Hallett TB, Bhatt S, Wabiri N, et al. Targeting the right interventions to the right people and places: the role of geospatial analysis in HIV program planning. AIDS. 2018;32(8):957.
Mweemba C, Hangoma P, Fwemba I, Mutale W, Masiye F. Estimating district HIV prevalence in Zambia using small-area estimation methods (SAE). Popul Health Metrics. 2022;20(1):8.
Prabhu S, Mehta SH, McFall AM, Srikrishnan AK, Vasudevan CK, Lucas GM, et al. Role of geospatial mapping in the planning of HIV programs a case study from Southern India. Medicine. 2021;100(34):e2709.
Mohammed JI, Comber A, Brunsdon C, (eds). Population estimation in small areas: combining dasymetric mapping with pycnophylactic interpolation. GIS Research UK (GISRUK) Conference; 2012 11th-13th April. Lancaster University.
Tomintz MN, García-Barrios VM. Location-allocation planning. In: Cockerham W, Dingwall R, Quah SR, editors. The Wiley Blackwell encyclopedia of health, illness, behavior and society. Hoboken: Wiley-Blackwell; 2014.
Entwisle B, Walsh R, Rindfuss R, Chamratritharong A. 1998. Land-use/land-cover and population dynamics National Research Council, Division of Behavioral and Social Sciences and Education Board on Environmental Change and Society Committee on the Human Dimensions of Global Change People and Pixels Linking Remote Sensing and Social Science. National Academies Press. Washington.
Tomintz MN, Barnett R. Geosimulation approach for filling the gap of non-response smoking data from the census 2013: A spatial analysis of census area unit geographies. New Zealand Geographer. 2018.
Eberth JM, McLain AC, Hong Y, Sercy E, Diedhiou A, Kilpatrick DJ. Estimating county-level tobacco use and exposure in South Carolina: a spatial model-based small area estimation approach. Ann Epidemiol. 2018;28(7):481–8.
Koh K, Grady SC, Darden JT, Vojnovic I. Adult obesity prevalence at the county level in the United States, 2000–2010: downscaling public health survey data using a spatial microsimulation approach. Spatial Spatio-Temp Epidemiol. 2018;26:153–64.
Morrissey K, Clarke G, Williamson P, Daly A, O’Donoghue C. Mental illness in Ireland: simulating its geographical prevalence and the role of access to services. Environ Plann B Plann Des. 2015;42(2):338–53.
Riva M, Smith DM. Generating small-area prevalence of psychological distress and alcohol consumption: validation of a spatial microsimulation method. Soc Psychiatry Psychiatr Epidemiol. 2012;47(5):745–55.
Zmölnig J, Tomintz MN, Clarke GP. 2016. Mapping the distribution of diabetes for public health support in Austria. gisScience–Die Zeitschrift fur Geoinformatik 3: 86–97.
Lucas AM, Wilson NL. Schooling, wealth, risky sexual behaviour, and HIV/AIDS in Sub-Saharan Africa. J Develop Stud. 2019;55(10):2177–92.
Akinyemi YC. Spatio-temporal pattern and gender disparity of HIV/AIDS prevalence in Nigeria. Int Soc Sci J. 2016;66(221–222):303–22.
Hajizadeh M, Sia D, Heymann SJ, Nandi A. Socioeconomic inequalities in HIV/AIDS prevalence in sub-Saharan African countries: evidence from the demographic health surveys. Int J Equity Health. 2014;13(1):18.
Doyle AM, Mavedzenge SN, Plummer ML, Ross DA. The sexual behaviour of adolescents in sub-Saharan Africa: patterns and trends from national surveys. Trop Med Int Health. 2012;17(7):796–807.
Ghys PD, Garnett GP. Trends in hiv prevalence and sexual behaviour among young people aged 15–24 years in countries most affected by hiv. Sex Trans Infect. 2010;86(Suppl. 2):72–83.
UN Statistics Division. Official list of MDG indicators: department of economic and social affairs UN 2008. https://mdgs.un.org/unsd/mdg/Host.aspx?Content=Indicators/OfficialList.htm.
UNICEF. MICS5 indicator list: UNICEF. 2022. https://mics.unicef.org/tools?round=mics5.
UN Statistics Division. Global indicator framework for the Sustainable Development Goals and targets of the 2030 Agenda for Sustainable Development. In: Affairs DoEaS, editor.: United Nations; 2022.
Slaymaker E. A critique of international indicators of sexual risk behaviour. Sex Trans Infect. 2004;80(suppl 2):ii13.
Mmari K, Sabherwal S. A review of risk and protective factors for adolescent sexual and reproductive health in developing countries: an update. J Adolesc Health. 2013;53(5):562–72.
UNAIDS, United Nations General Assembly Special Session on HIV/AIDS. Monitoring the declaration of commitment on HIV/AIDS guidelines on construction of core indicators 2008 reporting. Geneva: World Health Organization; 2007.
National Bureau of Statistics and United Nations Children's Fund. Nigeria multiple indicator cluster survey 2016–2017. 2018. https://nigerianstat.gov.ng/download/670.
Mahmood T, Abbas F, Kumar R, Somrongthong R. Why under five children are stunted in Pakistan? A multilevel analysis of Punjab multiple indicator cluster survey (MICS-2014). BMC Public Health. 2020;20(1):952.
Pezzulo C, Nilsen K, Carioli A, Tejedor-Garavito N, Hanspal SE, Hilber T, et al. Geographical distribution of fertility rates in 70 low-income, lower-middle-income, and upper-middle-income countries, 2010–16: a subnational analysis of cross-sectional surveys. Lancet Glob Health. 2021;9(6):e802–12.
Danert K, Healy A. Monitoring groundwater use as a domestic water source by urban households: analysis of data from lagos state, nigeria and sub-saharan africa with implications for policy and practice. Water. 2021;13(4):568.
Clarke G, Harding A. Conclusions and future research directions Spatial microsimulation A reference guide for users. Berlin: Springer; 2013.
Harland K. microsimulation model user guide: flexible modelling framework national centre for research methods working Paper 06/13. 2013.
Evans AJ. Flexible modelling framework (FMF). v1.3 ed2014. p. A tool for the microsimulation component to link synthetic populations with the original microdata.
Lovelace R, Dumont M. Spatial microsimulation with R. Boca Raton: CRC Press; 2016.
Edwards KL, Clarke G. SimObesity: combinatorial optimisation (deterministic) model a reference guide for users spatial microsimulation. Berlin: Springer; 2013.
Harland K, Heppenstall A, Smith D, Birkin MH. Creating realistic synthetic populations at varying spatial scales: a comparative critique of population synthesis techniques. J Artif Soc Soc Simul. 2012. https://doi.org/10.1564/jasss.1909.
De Menten G, Dekkers G, Bryon G, Liégeois P, O’Donoghue C. Liam2: a new open source development tool for discrete-time dynamic microsimulation models. J Artif Soc Soc Simul. 2014;17(3):1–9.
Kavroudakis D. sms: microdata for geographical analysis in R. J Stat Softw. 2015;68(2):1–23.
Richiardi M, Richardson R. JAS-mine: a new platform for microsimulation and agent-based modelling. Int J Microsimul. 2017;10(1):106–34.
Burden S, Steel D. Constraint choice for spatial microsimulation. Popul Space Place. 2016;22(6):568–83.
Birkin M, Clarke G. The enhancement of spatial microsimulation models using geodemographics. Ann Reg Sci. 2012;49(2):515–32.
Birkin M, Clarke M. Spatial microsimulation models: a review and a glimpse into the future population dynamics and projection methods. Berlin: Springer; 2011.
Freistein K, Mahlert B. The potential for tackling inequality in the sustainable development goals. Third World Quarterly. 2016;37(12):2139–55.
O’Neill J, Tabish H, Welch V, Petticrew M, Pottie K, Clarke M, et al. Applying an equity lens to interventions: using PROGRESS ensures consideration of socially stratifying factors to illuminate inequities in health. J Clin Epidemiol. 2014;67(1):56–64.
Moon G, Quarendon G, Barnard S, Twigg L, Blyth B. Fat nation: deciphering the distinctive geographies of obesity in England. Soc Sci Med. 2007;65(1):20–31.
Smith DM, Clarke GP, Harland K. Improving the synthetic data generation process in spatial microsimulation models. Environ Plan A. 2009;41(5):1251–68.
Abubakar EO, Cunningham N. 2022. A spatially enriched synthetic population developed by spatial microsimulation of 2016/2017 Multiple Indicator Cluster Survey microdata of Kogi State Nigeria. In: Abubakar EO. Draft version (eds). Harvard Dataverse. Kogi
Edwards KL, Clarke GP, Thomas J, Forman D. Internal and external validation of spatial microsimulation models: small area estimates of adult obesity. Appl Spat Anal Policy. 2011;4(4):281–300.
Edwards KL, Tanton R. Validation of spatial microsimulation models spatial microsimulation: a reference guide for users. Berlin: Springer; 2013.
Lovelace R, Ballas D, Watson M. A spatial microsimulation approach for the analysis of commuter patterns: from individual to regional levels. J Transp Geogr. 2014;34:282–96.
Morrissey K, O’Donoghue C. Validation issues and the spatial pattern of household income spatial microsimulation for rural policy analysis. Berlin: Springer; 2013. p. 87–102.
Edwards KL, Clarke GP. The design and validation of a spatial microsimulation model of obesogenic environments for children in Leeds, UK: SimObesity. Soc Sci Med. 2009;69(7):1127–34.
Morrissey K, Clarke G, O’Donoghue C. Linking static spatial microsimulation modelling to meso-scale models: the relationship between access to GP services and long-term illness spatial microsimulation: a reference guide for users. Berlin: Springer; 2013.
de Walque D. Does education affect HIV status? evidence from five African countries. World Bank Econo Rev. 2009;23(2):209–33.
Bärnighausen T, Tanser F, Gqwede Z, Mbizana C, Herbst K, Newell M-L. High HIV incidence in a community with high HIV prevalence in rural South Africa: findings from a prospective population-based study. AIDS. 2008;22(1):139.
Kasirye I. HIV/AIDS Sero-prevalence and Socio-economic status: evidence from Uganda. Afr Dev Rev. 2016;28(3):304–18.
Kayeyi N, Fylkesnes K, Wiium N, Sandøy IF. Decline in sexual risk behaviours among Young People in Zambia (2000–2009): do neighbourhood contextual effects play a role? PLoS ONE. 2013;8(5): e64881.
de Walque D, Nakiyingi-Miiro JS, Busingye J, Whitworth JA. Changing association between schooling levels and HIV-1 infection over 11 years in a rural population cohort in south-west Uganda. Trop Med Int Health. 2005;10(10):993–1001.
Uchudi J, Magadi M, Mostazir M. A multilevel analysis of the determinants of high-risk sexual behaviour in sub-Saharan Africa. J Biosoc Sci. 2012;44(3):289–311.
Somefun OD, Odimegwu C. The protective role of family structure for adolescent development in sub-Saharan Africa. PLoS ONE. 2018;13(10): e0206197.
Pilgrim NA, Ahmed S, Gray RH, Sekasanvu J, Lutalo T, Nalugoda F, et al. Family structure effects on early sexual debut among adolescent girls in Rakai. Uganda Vulnerable Child Youth Stud. 2014;9(3):193–205.
Uthman OA, Kayode GA, Adekanmbi VT. Individual and contextual socioeconomic determinants of knowledge of the ABC approach of preventing the sexual transmission of HIV in Nigeria: a multilevel analysis. Sexual Health. 2013;10(6):522–9.
Marston M, Slaymaker E, Cremin I, Floyd S, McGrath N, Kasamba I, et al. Trends in marriage and time spent single in sub-Saharan Africa: a comparative analysis of six population-based cohort studies and nine demographic and health surveys. Sex Trans Infect. 2009;85(Suppl 1): i64.
Żaba B, Isingo R, Wringe A, Marston M, Slaymaker E, Urassa M. Influence of timing of sexual debut and first marriage on sexual behaviour in later life: findings from four survey rounds in the Kisesa cohort in northern Tanzania. Sex Trans Infect. 2009;85(Suppl 1):i20–6.
Olivier J, Wodon Q. Religion, reproductive health, and sexual behavior in ghana: why statistics from large surveys don’t tell the whole story. Rev Faith Int Affairs. 2015;13(2):64–73.
Romero-Estudillo E, González-Jiménez E, Mesa-Franco MC, García-García I. Gender-based differences in the high-risk sexual behaviours of young people aged 15–29 in Melilla (Spain): a cross-sectional study. BMC Public Health. 2014;14(1):745.
Bamidele JO, Asekun-Olarinmoye EO, Odu OO, Amusan OA, Egbewale BE. Sociodemographic characteristics and health risk behaviours among students of a tertiary institution in south western Nigeria. Afr J Med Med Sci. 2007;36(2):129–36.
Somefun OD. Religiosity and sexual abstinence among Nigerian youths: does parent religion matter? BMC Public Health. 2019;19(1):416.
Jewkes R, Morrell R, Hearn J, Lundqvist E, Blackbeard D, Lindegger G, et al. Hegemonic masculinity: combining theory and practice in gender interventions. Cult Health Sex. 2015;17(sup2):112–27.
Hearn J, Morrell R. Reviewing hegemonic masculinities and men in Sweden and South Africa. Men Mascul. 2012;15(1):3–10.
Messerschmidt JW. Hegemonic masculinity: formulation, reformulation, and amplification. Lanham: Rowman Littlefield Publishers; 2018.
Connell RW, Messerschmidt JW. Hegemonic masculinity: rethinking the concept. Gend Soc. 2005;19(6):829–59.
MacPhail C. Challenging dominant norms of masculinity for HIV prevention. Afr J AIDS Res. 2003;2(2):141–9.
Bowleg L, Teti M, Massie JS, Patel A, Malebranche DJ, Tschann JM. ‘What does it take to be a man? What is a real man?’: ideologies of masculinity and HIV sexual risk among black heterosexual men. Cult Health Sex. 2011;13(5):545–59.
Wellings K, Collumbien M, Slaymaker E, Singh S, Hodges Z, Patel D, et al. Sexual behaviour in context: a global perspective. The Lancet. 2006;368(9548):1706–28.
Ahmed S, Lutalo T, Wawer M, Serwadda D, Sewankambo NK, Nalugoda F, et al. HIV incidence and sexually transmitted disease prevalence associated with condom use: a population study in Rakai. Uganda. 2001;15(16):2171–9.
de Walque D, Kline R. Variations in condom use by type of partner in 13 Sub-Saharan African Countries. Stud Fam Plann. 2011;42(1):1–10.
Soler H, Quadagno D, Sly DF, Riehman KS, Eberstein IW, Harrison DF. Relationship dynamics, ethnicity and condom use among low-income women. Fam Plann Perspect. 2000;32(2):82–101.
Ramjee G, Daniels B. Women and HIV in Sub-Saharan Africa. AIDS Res Ther. 2013;10(1):30.
Dunkle KL, Jewkes RK, Brown HC, Gray GE, McIntryre JA, Harlow SD. Transactional sex among women in Soweto, South Africa: prevalence, risk factors and association with HIV infection. Soc Sci Med. 2004;59(8):1581–92.
UNAIDS. The gap report. Geneva: United Nations Programme on HIV/AIDS; 2014.
Oyediran KA, Feyisetan OI, Akpan T. Predictors of condom-use among young never-married males in Nigeria. J Health Popul Nutr. 2011;29(3):273–85.
Nompumelelo ZD, Shissana O, Udjo E, Mosala T, Seager J. An Audit of HIV/AIDS policies in Botswana, Lesotho, Mozambique, South Africa. Cape Town: Republic of South Africa; 2004.
Rahman A, Harding A. Small area estimation and microsimulation modeling. Boca Raton: Chapman and Hall CRC; 2017.
World Health Organization. State of inequality reproductive maternal newborn and child health interactive visualization of health data. Geneva: World Health Organization; 2015.
Foley R, Charlton M, Fotheringham S. GIS in Health and Social Care Planning. In: Bavaud F, Mager C, editors. Handbook of Theoretical and Quantitative Geography. UNIL-FGSE-Workshop series. Lausanne: Univ. de Lausanne-Faculté des géosciences et de l'environnement. 2009. p. 73–115.
Eide AH, Dyrstad K, Munthali A, Van Rooy G, Braathen SH, Halvorsen T, et al. Combining survey data, GIS and qualitative interviews in the analysis of health service access for persons with disabilities. BMC Int Health Hum Rights. 2018;18(1):26.
United Nations Development Programme (UNDP). What does it mean to leave no one behind? A UNDP discussion paper and framework for implementation. New York: UNDP; 2018.
Pesaresi M, Florczyk A, Schiavina M, Melchiorri M, Maffenini L. GHS settlement grid, updated and refined REGIO model 2014 in application to GHS-BUILT R2018A and GHS-POP R2019A, multitemporal (1975–1990–2000–2015), R2019A. European Commission, Joint Research Centre (JRC). 2019.
Ethics approval and consent to participate
Consent for publication
Both authors consent to the publication of this manuscript in the International Journal of Health Geographics should this be accepted by the journal’s editor.
The authors have no competing interests as defined by BMC, or other interests that might be perceived to influence the results and/or discussion reported in this manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Abubakar, E.O., Cunningham, N. Small-area estimation and analysis of HIV/AIDS indicators for precise geographical targeting of health interventions in Nigeria. a spatial microsimulation approach. Int J Health Geogr 22, 23 (2023). https://doi.org/10.1186/s12942-023-00341-8