Improved benchmark-multiplier method to estimate the prevalence of ever-injecting drug use in Belgium, 2000–10
© Bollaerts et al.; licensee BioMed Central Ltd. 2013
Received: 15 January 2013
Accepted: 23 March 2013
Published: 3 May 2013
Accurate estimates of the size of the drug-using populations are essential for evidence-based policy making. However, drug users form a ‘hidden’ population, necessitating the use of indirect methods to estimate population sizes.
The benchmark-multiplier method was applied to estimate the population size of ever injecting drug users (ever-IDUs), aged 18–64 years, in Belgium using data from the national HIV/AIDS register and from a sero-behavioral study among injecting drug users. However, missing risk factor information and absence of follow-up of the HIV+/AIDS– cases, limits the usefulness of the Belgian HIV/AIDS register as benchmark. To overcome these limitations, statistical corrections were required. In particular, Imputation by Chained Equations was used to correct for the missing risk factor information whereas stochastic mortality modelling was applied to account for the mortality among the HIV+/AIDS– cases. Monte Carlo simulation was used to obtain confidence intervals, properly reflecting the uncertainty due to random error as well as the uncertainty associated with the two statistical corrections mentioned above.
In 2010, the prevalence (/1000) of ever-IDUs was estimated to be 3.5 with 95% confidence interval [2.5;4.8]. No significant time trends were observed for the period 2000–2010.
To be able to estimate the ever-IDU population size using the Belgian HIV/AIDS register as benchmark, statistical corrections were required without which seriously biased estimates would result. By developing the improved methodology, Belgium is again able to provide ever-IDU population estimates, which are essential to assess the coverage of treatment and to forecast health care needs and costs.
Accurate estimates of the size of the drug-using population are indispensable to govern and evaluate drug policy. They are particularly important to forecast health care needs and costs and to assess the coverage of treatment and harm reduction measures e.g. [1, 2]. However, drug users form a `hidden’ population, hiding their membership because it involves illegal and stigmatized behaviour. As a result, `hidden’ populations lack sampling frames and classical epidemiological methods (e.g. population-based surveys) fail. Therefore, indirect estimation methods are to be used. Generally, three types of indirect estimation methods are distinguished to estimate the size of drug using populations: capture-recapture (CR), benchmark-multiplier (BM) and multivariate indicator (MI) methods [3–5]. An application of the BM method is presented in this paper.
The BM method is a very intuitive method, relatively easy to apply and hence, frequently used . The BM method combines (an estimate of) the size of the known part of the target population, the benchmark, with an estimate of the proportion of individuals from the target population belonging to the benchmark. The reciprocal of the proportion is called the multiplier. The size of the target population is then estimated as the product of the benchmark and the multiplier. Common sources used as benchmark are mortality registers, HIV/AIDS registers, drug treatment or police data e.g. [7, 8]. During the beginning of the 1970s, the mortality multiplier method was already used . Multiplier methods e.g. [10, 11] have produced estimates consistent with estimates produced by CR methods, which have been put forward as a preferred methodology . The key assumptions underlying the BM method are that the benchmark should be exhaustive, the sample used to estimate the multiplier should be representative of the target population and the case definition used for the benchmark should exactly match the one used for the multiplier, including the time window .
In this paper, the benchmark-multiplier (BM) method was applied to estimate the prevalence of ever injecting drug users (ever-IDUs) within the population aged 15–64 years in Belgium, 2000–10, using the HIV/AIDS register as benchmark. In Belgium, HIV-screening is widely used and all sera of which the screening test result was positive, are submitted for confirmation to one of the seven AIDS Reference Laboratories (ARLs) in Belgium. The registration results of all ARLs are then send to the national HIV/AIDS register. As the seven ARLs are the only laboratories subsidized for performing HIV confirmation tests, the Belgian HIV/AIDS register is deemed to be exhaustive and hence, constitutes a potential suitable benchmark.
For Belgium, the latest estimate of the ever-IDU prevalence date from 1995 and was obtained adopting the BM approach using the HIV/AIDS register as benchmark . However, the Belgium HIV/AIDS register suffers from missing risk factor information. Particularly, of all registered cases in 2000–10, 28.6% were reported without information on probably mode of HIV transmission (e.g. sexual contact, blood transfusion, injecting drug use). Simply discarding these cases from analysis would lead to an underestimation of the IDU prevalence. Simply extrapolating the risk factor fractions of the subpopulation for which the risk factors are known to the subpopulation for which they are unknown, would lead to overconfident results and would only lead to unbiased results under very strict assumptions (i.e. the Missing Completely ad Random assumption). In addition, the national HIV/AIDS register lacks follow-up of the non-AIDS cases, implying absence of information on the vital status of the HIV+/AIDS– cases. Not accounting for the mortality among these cases would result in an overestimation of the number of alive seropositive IDUs, with the bias increasing as the time since the onset of the HIV-epidemic (mid-eighties) increases. The methodology proposed in  used a simple extrapolation method to account for the missing risk factor information and failed to account for the mortality among HIV+/AIDS– cases. Therefore, the current paper presents an improved application of HIV/AIDS BM method, using up-to-date statistical methods to correct for the missing risk factor information and lacking follow-up of the HIV+/AIDS– cases.
Benchmark: national HIV/AIDS register
During 2000–10, an average of 56 screening tests per 1000 inhabitants per year were performed in Belgium, excluding tests related to blood donations (National Institute for Sickness and Invalidity Insurance). The confirmation tests are performed by one of the seven ARLs in Belgium, being the only laboratories subsidized for performing these tests. The registration results of the seven ARLs are validated for duplicate recording and are all included in the national HIV/AIDS register, being hosted by the Scientific Institute of Public Health, Brussels (WIV-ISP). Since its foundation in 1985–86, the Belgian HIV/AIDS register is deemed to be exhaustive.
For each confirmed HIV-positive test, a standardized form is sent to the patient’s clinician to collect additional information on nationality, residence, sexual orientation, probable mode of HIV transmission and CD4 count at time of HIV diagnosis. The response categories for the self-reported probable mode of HIV transmission are homo- and heterosexual transmission, transmission through blood transfusion, transmission through injecting drug use and mother-to-child transmission. A person can indicate multiple modes of transmission. Unfortunately, the standardized forms are not always fully completed returned to the WIV-ISP, resulting in missing risk factor information.
Follow-up is conducted for patients who developed AIDS with each year data being collected on last consultation and possible death. The variables available for the AIDS cases are year at AIDS-diagnosis, year at death and year at lost to follow-up, which are only completed upon occurrence of the event and hence, are time dependent. The HIV+/AIDS– cases are not subject to follow-up.
Multiplier: sero-behavioral prevalence study
In Belgium, a sero-behavioral study among drug users in contact with drug treatment facilities or imprisoned was carried out in 2004–05 . In total, 1005 drug users in treatment and 117 incarcerated drug users (15–40 years) enrolled at 65 different drug treatment facilities and 15 different prisons geographically dispersed over Belgium, respectively, participated to the study. Of the drug users in treatment and in prison, 57% (n = 573) and 68% (n = 80) declared to have injected drugs at least once during their life. Intravenous blood samples were taken to determine the HIV- as well as the Hepatitis B and -C status of the participants. The HIV-seroprevalence among IDUs in treatment and in prison was estimated to be 2.8% (95% Wald CI: [1.5;4.2]) and 5% (95% Wald CI:[0.2;9.9]), respectively. These prevalences were not significantly different (p-value = 0.30), yielding an overall estimated prevalence of 3.1% (95% Wald CI: [1.8;4.8]).
In addition to serological studies, the HIV prevalence among IDUs in Belgium can be obtained from routine diagnostic testing, of which the results are yearly available. This allows the investigation of time trends. However, a concern regarding the (geographical) representativeness of the Belgian data exists. In line with (Western) European trends , no significant time trends in HIV prevalences among IDUs were observed during the last 10 years in Belgium based on the results from routine diagnostic testing . Therefore, the HIV prevalence from the sero-behavioral study was assumed to apply for the entire period 2000–10.
Schematic overview of the 3-step Monte Carlo simulation model to estimate ever IDU prevalences Ɨ
Imputation by Chained Equations: missing risk factor information
Missing risk factor information was imputed using Imputation by Chained Equations. The imputation model contains the variables: injecting drug use, sex, nationality, year at registration and age at registration. The imputation results in one complete dataset X k , containing original and imputed values.
Stochastic Mortality Modeling: lacking follow-up of the HIV + /AIDS – cases
For a complete dataset k, the number of registered HIV-cases for whom injecting drug use was the most probable route of transmission and who were alive at time t is calculated as
where I i indicates the ‘vital’ status with I i = 1 if person i is still alive and living in Belgium and I i = 0 otherwise, where r i is the number of years since HIV registration or r i = t − t hi and where p d is the annual non-AIDS mortality rate among seropositive IDUs with p d ~ betapert*(0.58%, 1.08%, 1.58%).
Benchmark-multiplier method: population size estimation
The number of ever-injecting drug users being alive at time t is given by with obtained from step 2, and n = 639.
Step 1: To correct for the missing risk factor information, Imputation by Chained Equations (ICE) was used. ICE is an iterative technique that starts by filling in the missing values in a simple way, e.g. by using mean imputation [17–19]. Then, each variable included in the imputation model X i , i = 1, 2,…q, is imputed in turn by the predictions of the regression model regressing the observed values of X i on the observed and imputed values of the remaining variables within the imputation model. The iteration through all q variables constitutes one cycle. After a sufficient number of cycles (typically 10 cycles), the final imputations are retained, resulting in one ‘complete’ dataset. There are several software packages available providing (M)ICE procedures, i.e. STATA , S-PLUS or R . For the current application, ICE as implemented in STATA was used to account for the missing risk factor information. In particular, the imputation model included all (auxiliary) variables that had a statistically significant association with the variable indicating whether injecting drug use was the most probable route of transmission; being sex, nationality, year at registration and age at registration. For every run k (k = 1,2,…K = 1000), one ‘complete’ dataset was generated.
Step 2: Binomial stochastic mortality modelling was applied to correct for the lacking information on the vital status of the HIV+/AIDS– cases. More precisely, the vital status of a HIV+/AIDS– case i was simulated using a binomial process having survival probability with being the number of years since HIV registration, implying lower survival probabilities when the time since HIV-registration increases. The non-AIDS crude mortality rate among IDUs (p d ) was assumed to follow a betapert distribution  with minimum, most likely and maximum values of 0.58%, 1.08% and 1.58% per annum, respectively. These non-AIDS crude mortality rates were calculated by Degenhardt et al.  based on 44 cohort-studies. The vital status of the AIDS cases at time t was deduced from the AIDS-specific variables year at AIDS-diagnosis (Ta), year at death (Td), and year at lost to follow-up (Tl). In particular, an AIDS case was indicated to be alive at time t if any information was recorded after time t (i.e. if Ta ≥ t or Td ≥ t or Tl ≥ t). On the other hand, an AIDS case was indicated to be death at time t if the person died before time t or was lost to follow-up. In the latter case, the person is most likely to be death or to have left the country.
with n being the size of sample used to estimate π . Observe that the sample size n is the only additional piece of information needed to construct the bias-corrected estimator as compared to the traditional estimator. Then, the estimate of the benchmark was obtained from step 1 and 2 whereas the estimate and the corresponding sample size n were derived from the sero-behavioral prevalence study . As the study indicated that x = 20 out of the n = 639 ever-IDUs were seropositive, was assumed to follow a beta distribution with or .
Results of the Monte Carlo (MC) simulation study estimating the prevalence of ever Injecting Drug Use (IDU), Belgium, 2000-10
Number of seropositive IDUs
Number of alive seropositive IDUs
Number of alive IDUs
95% MI CIƗ
95% MI CIƗ
95% MI CIƗ
95% MI CIƗ
To obtain recent estimates of the ever-IDU prevalence in Belgium, an improved HIV multiplier method was applied. The improved methodology allows using the Belgian HIV/AIDS register as benchmark while accounting for missing risk factor information (through Imputation by Chained Equations - ICE) and lacking follow-up of the HIV+/AIDS– cases (through stochastic mortality modelling). The statistical corrections were implemented in the best possible way. Indeed, ICE is a very flexible imputation approach, that can deal with several incomplete variables, does not require specific missing data patterns and operates under the less stringent assumption that the probability of the missing values depends only on the observed values and not on unobserved values (the so-called Missing at Random assumption, ) [17–19]. The stochastic mortality model accounts for the uncertainty in the mortality rate and the Monte Carlo simulation model properly acknowledges the uncertainty associated with the statistical corrections mentioned above.
In 2010, the prevalence (per 1000 inhabitants, 15–64 years) of ever injecting drug use in Belgium was estimated to be 3.5 (95% CI: 2.5-4.8), with no significant time trends being observed for the period 2000–10. Similar findings were reported by the other European countries providing time trends of the national prevalence of ever injecting drug use, i.e. Norway, Cyprus and the Czech republic . This absence of a time trend is seemingly in contrast with recent findings from the EMCDDA treatment demand indicator, suggesting a decline in recent-onset heroin users and heroin injectors in the Western European countries . However, it need to be stressed that the prevalence of ever injecting drug use is less sensitive to recent trends in drug use.
The above mentioned results should also be regarded in the light of the assumptions underlying the BM method. Particularly, the BM estimates rely on the conditions that the benchmark is complete and that the multiplier is estimated based on a representative sample of the target population. As benchmark of the current analysis, the national HIV/AIDS register was used. Since its foundation, the Belgian HIV/AIDS register is deemed to be exhaustive as it includes the results of all Belgian laboratories subsidised for performing HIV confirmation tests. As not only AIDS cases but also HIV cases are registered, the Belgian HIV/AIDS register has the additional advantage that back-calculation methods are not needed to obtain the HIV prevalence. The multiplier was obtained from a national sero-prevalence study including 653 IDUs . Although this is the largest and most recent Belgian sero-prevalence study, it did not include IDUs outside treatment/prison, for whom currently no information is available. Despite the fact that low threshold drugs services are common in Belgium, the possibility of bias resulting from a multiplier derived from a sample that is not fully representative of the target population cannot be ruled out.
To summarize, the current paper provides a thorough application of the BM method. Nevertheless, the use of such indirect methods inherently relies on empirically non-verifiable (but reasonable) assumptions, such as the representativeness of the multiplier. Therefore, the current estimate of the ever IDU prevalence should be complemented with estimates based on other data sources (e.g. substitution treatment register) and/or using other methods (e.g.capture-recapture, multivariate indicator method) (e.g. ). For now, it is (reasonably) assumed that the HIV prevalence among IDUs remained stable during the last ten years. However, new prevalence estimates are needed for future applications of the BM method based on the Belgian HIV/AIDS register. Finally, this paper presents estimates of the ever IDU prevalence, which are important to assess future health care needs regarding e.g. Hepatitis C treatment. However, when interest is in assessing the coverage adequacy of needle exchange programmes (NEPs) or opioid substitution treatment (OST), national estimates of the current IDU prevalence are required. The latter prevalence (per 1000 inhabitants, 15–64 years) can be roughly estimated as 1.4 (95% CI: 1.02-1.97) by multiplying the ever-IDU prevalence with the ratio between current injectors and ever injectors (r = 41%) among all treatment demands as provided by the treatment demand indicator .
As the latest estimates of the prevalence of ever-IDUs for Belgium date from 1995 , obtaining new estimates was prioritized. To be able to estimate the ever-IDU population size adopting the BM approach using the Belgian HIV/AIDS register as benchmark, statistical corrections were required. In particular, not accounting for the mortality among HIV+/AIDS– cases would nowadays results in seriously biased estimates. By developing the improved methodology, recent ever-IDU population estimates for Belgium could be obtained. These estimates are essential to monitor trends over time and to forecast health care needs and costs (e.g. for Hepatitis C treatment).
KB holds a master degree in psychology and a PhD in statistics on quantitative epidemiology. She is currently working as senior biostatistician within the research group Substance Use and Related Disorders. AS is a medical doctor and head of the HIV/AIDS register. MA is professor in biostatistics with strong interest in quantitative epidemiology.
We thank the AIDS Reference Laboratories and the AIDS Reference Centres for their invaluable contribution to the HIV/AIDS/STI surveillance programmes.
- Hickman M, Higgins V, Hope V, Bellis M, Tilling K: Injecting drug use in Brighton, Liverpool, and London: best estimates of prevalence and coverage of public health indicators. J Epidemiol Community Health. 2004, 58: 471-71.View ArticleGoogle Scholar
- Kimber J, Hickman M, Degenhardt L, Coulson T, van Beek I: Estimating the size and dynamics of an injecting drug user population and implications for healt service coverage: comparison of indirect prevalence estimation methods. Addiction. 2008, 103 (10): 1604-1613. 10.1111/j.1360-0443.2008.02276.x.View ArticlePubMedGoogle Scholar
- Kraus L, Augustin R, Frischer M, Kümmler P, Uhl A, Wiessing L: Estimating prevalence of problem drug use at national level in countries of the European Union and Norway. Addiction. 2003, 98: 471-485. 10.1046/j.1360-0443.2003.00326.x.View ArticlePubMedGoogle Scholar
- Hickman M, Hope V, Platt L, Higgins V, Bellis M, Rhodes T: Estimating prevalence of injecting drug use: a comparison of multiplier and capture-recapture methods in cities in England and Russia. Drug Alcohol Rev. 2006, 25 (2): 131-140. 10.1080/09595230500537274.View ArticlePubMedGoogle Scholar
- Wickens TD: Quantitative methods for estimating the size of the drug using population. J Drug Issues. 1993, 23: 185-216.View ArticleGoogle Scholar
- Scientific Monograph No 1. Estimating the Prevalence of Problem Drug Use in Europe. Edited by: Stimson GV, Hickman M, Quirk A, Frisher M, Taylor C. 1997, Lisbon: EMCDDA, 113-123.Google Scholar
- Barrio Anta G, Oliva JS, Bravo MAJ, Mateo SD, Domingo-Salvany A: Estimating the prevalence of drug injection using a multiplier method based on a register of new HIV diagnoses. Eur J Public Health. 2011, 21 (5): 646-648. 10.1093/eurpub/ckq076.View ArticlePubMedGoogle Scholar
- Augustin R, Kraus L: Changes in prevalence of problem opiate use in Germany between 1990 and 2000. Eur Addict Res. 2004, 10 (2): 61-67. 10.1159/000076115.View ArticlePubMedGoogle Scholar
- Andima H, Krug D, Bergner L, Patrick S, Whitman S: A prevalence estimation model of narcotics addiction in New York City. Am J Epidemiol. 1973, 98: 56-62.PubMedGoogle Scholar
- Hartnoll R, Lewis R, Mitcheson M, Bryer S: Estimating the prevalence of opioid dependence. Lancet. 1885, 1997: 203-205.Google Scholar
- Hall W, Ross J, Lynskey M, Law M, Degenhardt L: How many heroin users are there in Australia?. Med J Aust. 2000, 73: 528-531.Google Scholar
- Hickman M, Taylor C: Indirect methods to estimate prevalence. Epidemiology of Drug Abuse. Edited by: Sloboda Z. 2005, New York: Springer US, 113-131.View ArticleGoogle Scholar
- Sartor F, Walckiers D, Sasse A, Bils L: Estimate of the prevalence of injecting drug use in Belgium. Arch. Public Health. 2001, 59: 63-75.Google Scholar
- Plasschaert S, Ameye L, De Clercq T, Walckiers D, Sartor F, Micalessi I: Study on HCV, HBV and HIV seroprevalence in a sample of drug users in contact with treatment centres or in prisons in Belgium, 2004–2005. 2005, Brussels: Scientific Institute of Public Health - Unit of Epidemiology, Report No.: IPH/EPI Reports N° 2005–029Google Scholar
- Trends in injecting drug use in Europe. Edited by: EMCDDA. 2010, LuxembourgGoogle Scholar
- Deprez N, Antoine J, Asueta-Lorente J-F, Bollaerts K, Van der Linden T, Van Bussel JCH: Belgian national report on drugs 2011. 2012, Brussels: Scientific Institute Public HealthGoogle Scholar
- Van Buuren S, Boshuizen HC, Knook DL: Multiple imputation of missing blood pressure covariates in survival analysis. Stat Med. 1999, 18 (6): 681-694. 10.1002/(SICI)1097-0258(19990330)18:6<681::AID-SIM71>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Van Buuren S: Multiple imputation of discrete and continuous data by fully conditional specification. Stat Methods Med Res. 2007, 16 (3): 219-242. 10.1177/0962280206074463.View ArticlePubMedGoogle Scholar
- Azur MJ, Stuart EA, Frangakis C, Leaf PJ: Multiple imputation by chained equations: what is it and how does it work?. Int J Methods Psychiatr Res. 2011, 20 (1): 40-49. 10.1002/mpr.329.View ArticlePubMedPubMed CentralGoogle Scholar
- Renschaw E: Stochastic Population Processes: Analysis. 2010, Approximations and Simulations: Oxford University PressGoogle Scholar
- Royston P, White IR: Multiple Imputation by Chained Equation (MICE): Implementation in Stata. J Stat Softw. 2012, 45 (4): 1-20.Google Scholar
- Van Buuren S, Groothuis-Oudshoorn K: MICE: Multivariate Imputation by Chained Equations in R. J Stat Softw. 2011, 45 (3): 1-67.Google Scholar
- Vose D: Risk Analysis: A Quantitative Guide. 2000, Chichester: John Wiley & Sons, 2Google Scholar
- Degenhardt L, Hall W, Warner-Smith M: Using cohort studies to estimate mortality among injecting drug users that is not attributable to AIDS. Sex Transm Infect. 2006, 82 (Suppl 3): iii56-iii63.PubMedPubMed CentralGoogle Scholar
- Gelman A, Carlin JB, Stern HS, Rubin DB: Bayesian Data Analysis. 2003, London: Chapman & Hall/CRCGoogle Scholar
- Molenberghs G, Kenward M: Missing Data in Clinical Studies. 2012, Chichester: John Wiley & SonsGoogle Scholar
- Barrio G, Montanari L, Bravo M, Guarita B, De la Fuente L, Pulido J, Vicente J: Trends of heroin use and heroin injection epidemics in Europe: Findings from the EMCDDA treatment demand indicator (TDI). J Substance Abuse Treatment. 2013, 10.1016/j.jsat.2012.11.002.Google Scholar
- Origer A: Prevalence of problem drug use and injecting drug use in Luxembourg: A longitudinal and methodological perspective. Eur Addiction Res. 2012, 18: 288-296. 10.1159/000337211.View ArticleGoogle Scholar
- Scientific Institute of Public Health (WIV-ISP): Belgian Treatment Demand Indicator Register. 2011, BrusselsGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.