Modelling HIV/AIDS epidemic in Nigeria

Eze, Jude Ikechukwu (2009) Modelling HIV/AIDS epidemic in Nigeria. PhD thesis, University of Glasgow.

Abstract

Nigeria is one of the countries most affected by the HIV/AIDS pandemic, third only to India and South Africa. With about 10% of the global HIV/AIDS cases estimated to be in the country, the public health and socio-economic implications are enormous.

This thesis has two broad aims: the first is to develop statistical models which adequately describe the spatial distribution of the Nigerian HIV/AIDS epidemic and its associated ecological risk factors; the second, to develop models that could reconstruct the HIV incidence curve, obtain an estimate of the hidden HIV/AIDS population and a short term projection for AIDS incidence and a measure of precision of the estimates.

To achieve these objectives, we first examined data from various sources and selected three sets of data based on national coverage and minimal reporting delay. The data sets are the outcome of the National HIV/AIDS Sentinel Surveillance Survey conducted in 1999, 2001, 2003 and 2005 by the Federal Ministry of Health; the outcome of the survey of 1057 health and laboratory facilities conducted by the Nigerian Institute of Medical Research in 2000; and case by case HIV screening data collected from an HIV/AIDS centre of excellence.

A thorough review of methods used by WHO/UNAIDS to produce estimates of the Nigerian HIV/AIDS scenario was carried out. The Estimation and Projection Package (EPP) currently being used for modelling the epidemic partitions the population into at-risk, not-at-risk and infected sub-populations. It also requires some parameter input representing the force of infection and behaviour or high risk adjustment parameter. It may be difficult to precisely ascertain the size of these population groups and parameters in countries as large and diverse as Nigeria. Also, the accuracy of vital rates used in the EPP and Spectrum program is doubtful. Literature on ordinary back-calculation, nonparametric back-calculation, and modified back-calculation methods was reviewed in detail. Also, an indepth review of disease mapping techniques including multilevel models and geostatistical methods was conducted.

The existence of spatial clusters was investigated using cluster analysis and some measure of spatial autocorrelation (Moran I and Geary c coefficients, semivariogram and kriging) applied to the National HIV/AIDS Surveillance data. Results revealed the existence of spatial clusters with significant positive spatial autocorrelation coefficients that tended to get stronger as the epidemic developed through time. GAM and local regression fit on the data revealed spatial trends on the north-south and east - west axis.

Analysis of hierarchical, spatial and ecological factor effects on the geographical variation of HIV prevalence using variance component and spatial multilevel models was performed using restricted maximum likelihood implemented in R and empirical and full Bayesian methods in WinBUGS. Results confirmed significant spatial effects and some ecological factors were significant in explaining the variation. Also, variation due to various levels of aggregation was prominent.

Estimates of cumulative HIV infection in Nigeria were obtained from both parametric and nonparametric back-calculation methods. Step and spline functions were assumed for the HIV infection curve in the parametric case. Parameter estimates obtained using 3-step and 4-step models were similar but the standard errors of these parameters were higher in the 4-step model. Estimates obtained using linear, quadratic, cubic and natural splines differed and also depended on the number and positions of the knots. Cumulative HIV infection estimates obtained using the step function models were comparable with those obtained using nonparametric back-calculation methods. Estimates from nonparametric back-calculation were obtained using the EMS algorithm. The modified nonparametric back-calculation method makes use of HIV data instead of the AIDS incidence data that are used in parametric and ordinary nonparametric back-calculation methods. In this approach, the hazard of undergoing HIV test is different for routine and symptom-related tests. The constant hazard of routine testing and the proportionality coefficient of symptom-related tests were estimated from the data and incorporated into the HIV induction distribution function. Estimates of HIV prevalence differ widely (about three times higher) from those obtained using parametric and ordinary nonparametric back-calculation methods. Nonparametric bootstrap procedure was used to obtain point-wise confidence interval and the uncertainty in estimating or predicting precisely the most recent incidence of AIDS or HIV infection was noticeable in the models but greater when AIDS data was used in the back-projection model.

Analysis of case by case HIV screening data indicate that of 33349 patients who attended the HIV laboratory of a centre of excellence for the treatment of HIV/AIDS between October 2000 and August 2006, 7646 (23%) were HIV positive with females constituting about 61% of the positive cases. The bulk of infection was found in patients aged 15-49 years, about 86 percent of infected females and 78 percent of males were in this age group. Attendance at the laboratory and the proportion of HIV positive tests witnessed a remarkable increase when screening became free of charge. Logistic regression analysis indicated a 3-way interaction between time period, age and sex. Removing the effect of time by stratifying by time period left 2-way interactions between age and sex.
A Correction factor for underreporting was ascertained by studying attendance at the laboratory facility over two time periods defined by the cost of HIV screening. Estimates of HIV prevalence obtained from corrected data using the modified nonparametric back-calculation are comparable with UN estimates obtained by a different method.

The Nigerian HIV/AIDS pandemic is made up of multiple epidemics spatially located in different parts of the country with most of them having the potential of being sustained into the future given information on some risk factors. It is hoped that the findings of this research will be a ready tool in the hands of policy makers in the formulation of policy and design of programs to combat the epidemic in the country. Access to data on HIV/AIDS are highly restricted in the country and this hampers more in-depth modelling of the epidemic. Subject to data availability, we recommend that further work be done on the construction of stratification models based on sex, age and the geopolitical zones in order to estimate the infection intensity in each of the population groups. Uncertainties surrounding assumptions of infection intensity and incubation distribution can be minimized using Bayesian methods in back-projection.

• Item Type: Thesis (PhD)
• Qualification Level: Doctoral
• Keywords: HIV/AIDS,back-calculation, disease mapping,geostatistical methods, multilevel models, cluster analysis, spatial autocorrelation, semivariogram, kriging, ecological factors, Bayesian methods,underreporting
• Subjects: R Medicine > RA Public aspects of medicine > RA0421 Public health. Hygiene. Preventive Medicine; H Social Sciences > HA Statistics; G Geography. Anthropology. Recreation > GE Environmental Sciences
• Colleges/Schools: College of Science and Engineering > School of Mathematics and Statistics > Statistics
• Supervisor's Name: McColl, Prof. John
• Date of Award: February 2009
• Depositing User: Mr Jude I Eze
• Unique ID: glathesis:2009-642
• Copyright: Copyright of this thesis is held by the author.
• Date Deposited: 25 Mar 2009
• Last Modified: 24 Apr 2019 13:56
• URI: https://theses.gla.ac.uk/id/eprint/642