- Research article
- Open Access
- Open Peer Review
Analysis of latent tuberculosis and mycobacterium avium infection data using mixture models
BMC Public Healthvolume 6, Article number: 240 (2006)
Estimation of the frequency of latent tuberculosis infection (LTBI) is difficult in areas with low tuberculosis infection rates and high exposure to non-tuberculous mycobacteria (NTM), including BCG vaccination. The objective was to assess LTBI and M avium infection and to estimate their probability based on skin tests responses in an infant population from a region with the aforementioned characteristics.
A population-based tuberculin skin test (TST) and sensitin (M avium) survey was conducted on seven years old infants in Biscay, a province from The Basque Country (Spain). 2268 schoolchildren received sensitin and 5277 TST. Participation rate was 89%. Commonly used estimation methods were compared with a method based on the fit of mixture models using the Expectation Maximization algorithm. Functions estimating the probabilities of LTBI and M avium infection given the observed skin tests responses were developed for vaccinated and unvaccinated children.
LTBI prevalences varied widely according to the estimation method. The mixture model provided prevalences higher than expected although intermediates between those obtained by currently recommended approaches. Exposure to previous BCG vaccine produces an upward shift of an average of about 3 mm on the induration size to attain the same probability of infection.
Our results confirm the commonplace exposure to NTM which effect should be taken into account when performing and assessing tuberculin surveys. The use of mixture analysis under the empirical Bayes framework allows to better estimate the probability of LTBI in settings with presence of other NTM and high BCG-vaccination coverage. An estimation of the average effect of BCG vaccination on TST induration is also provided. These models maximise information coming from classical tuberculin surveys and could be used together with the newly developed blood tests to improve survey's specificity and cost-effectiveness.
Tuberculosis is a well-recognised worldwide health problem . However, the picture differs dramatically depending on the socioeconomic development of the area in question. While in developing countries tuberculosis is still a first order challenge, it is being kept increasingly under control in developed countries. In the latter, the diagnosis and treatment of Latent Tuberculosis Infection (LTBI) is now necessary in order to eliminate tuberculosis . New blood tests [3, 4] are already changing the approach to this problem, and are to be implemented in the near future in western industrialised countries. However, tuberculin skin tests (TST) remain today the most used tools for LTBI diagnostic  and treatment  decisions at a global scale.
Unfortunately, the interpretation of TST results suffers from several drawbacks when performed on populations where the BCG vaccination rate is high . This situation becomes even more awkward when cross-reactivity due to infection with nontuberculous mycobacteria (NTM) is present [7, 8] (a subject which has not been addressed much so far) and the survey is carried out in a region with low tuberculosis prevalence .
In the region under study, a province in the Basque Country (Spain), all the aforementioned characteristics that make the interpretation of TST results difficult are present. The Basque Country is the only community that has had a programme of systematic BCG vaccination in use for more than two decades , in spite of being one of the communities with the highest socioeconomic level. Despite the current debate regarding whether it is appropriate to maintain this policy in our region due to a substantial decline in tuberculosis prevalence rates, methods that allow a continuous and accurate monitoring of LTBI in settings of this kind are still required.
Several statistical approaches have been attempted to cope with the confounding effect of NTM infection (BCG vaccination included) in the accurate estimation of LTBI. Traditionally, methods that propose fixed general cut-off points have been in use [11–14]. However, problems in determining the prevalence of LTBI using these methods have long been recognised, making the exploration of new statistical methods  necessary. Along these lines, Neuenschwander et al. [15–17] have applied the so-called mixture model analysis  to this problem, first focusing on unvaccinated individuals and then extending this approach to BCG-vaccinated subjects. They have shown mixture models to be a flexible tool for analysing data of this kind.
The purpose of the present study is threefold: a) to estimate the prevalence of LTBI and Mycobacterium avium (M avium) infection in the infant population of The Basque Country (Spain), b) to find a method that could help us to accurately discriminate LTBI from infection with other mycobacteria, including BCG vaccination and c) to help determine the probability of LTBI and M avium infection according to tuberculin or sensitin induration sizes respectively.
Population and sample size
The region under study is called Biscay, a province of The Basque Country with a population of 1.2 million inhabitants. In 1997, 13 smear-positive sputum cases among the total 38 incident tuberculosis cases per 100,000 population were notified in the region . Our study population consisted of all children living in Biscay and turning seven years of age during 1998. We decided to study school-age children based on two grounds: to allow for an effect of the BCG vaccination to still be present and accessibility criteria. This accounts for 8430 children.
To estimate the required sample size, we used results from a previous survey in the same region . A vaccination coverage rate of 89%, a participation rate of 80%, a maximum prevalence of 20% among BCG-vaccinated and 2% among unvaccinated children, respectively, with TST induration size ≥ 10 mm were anticipated. Standard errors of 0.4 and 0.35% were respectively accepted for the prevalence estimations. A required sample size of 6250 children was obtained, a quantity close to the total population size assuming a recruitment rate of 80%, so eventually we decided to recruit as many schoolchildren as possible, since the study could provide the potential benefit of detecting additional children with very large indurations, requiring further examination in order to assess the need for prophylactic treatment.
This study focuses on some aspects of a more general research project devoted to investigate the distribution of tuberculin and sensitin induration sizes as well as the relationship among them and the so-called booster effect. In order to accomplish this endeavour, recruited children were randomly distributed in three groups. Two thirds received tuberculin as the first test and one third received the M avium sensitin, which was the most prevalent MNT in our region at the study time (data provided by the Basque Country Public Health Department). Out of the children receiving tuberculin as the first test, half were given a repeat dose of tuberculin one week later, and half were given the sensitin. Children who had received the sensitin as the first test were given a tuberculin one week later. In this paper, only results regarding administration of the first skin test are shown.
The field work was planned to be conducted in 5 months so as to prevent substantial changes in the children's characteristics and personnel involved in the process that could influence results.
Firstly, an informative campaign among parents and teachers was conducted, presenting the aim and methods of the study, and encouraging them to participate. In addition, an informative leaflet was handed out. Parents had to give written informed consent prior to inclusion of their children in the study. The vaccine information recorded in the Personal Paediatric Record was requested. The protocol was approved by the hospital's Research Ethics Committee.
There were two criteria that led to children not receiving the skin test: absence of written informed permission and the existence of documented evidence of previous positive TST.
Demographic data, past medical history and vaccination status were recorded at study entry. As the latter was deemed a key point in this research, we used three methods to improve its classification: a) certificate of BGC vaccination in the Personal Paediatric Record; b) confirmation of BGC vaccination in the province register; c) BGC scar. To classify a child in the non-BCG group, all three methods had to have negative results and to classify a child as BCG-vaccinated, at least two of three criteria had to be positive. Those who did not fulfil these requirements were considered to have an unknown vaccination status.
Intradermal injections with 0.1 ml of 2 TU PPD-RT23 tuberculin (bioequivalent to 5 TU PPD-S) and 0.1 ml of RS-10 (M. avium) sensitin were used. The PPD RS-10 was chosen because it is the most frequently isolated NTM among the patients in our area. All tuberculin doses belonged to the same lot (L-03) and sensitin doses came from two lots (L-66 and L-68) . Both tuberculin and sensitin were prepared in and obtained from Statens Serum Institut in Copenhagen, Denmark.
The tuberculin skin test was administered intradermally  and the induration diameter was measured transversely to the long axis of the arm 72 hours after application.
Test application and reading were performed by a specifically trained team led by a specialist in respiratory medicine and carried out by eight highly experienced nurses. Four two-nurse groups were assembled, and each group performed double and independent readings.
The distribution of induration diameters was described graphically using histograms and categorised using 5 mm-intervals. The prevalence of LTBI was estimated using three methods proposed in the literature [14, 23, 24]. The first criterion, probably the most widely used, considers all reactions ≥ 10 mm to be a marker of LTBI . The second one multiplies the proportion of reactions of size ≥ 14 mm by 1.22 , and the last one, also known as the mirror method, considers all indurations of 17 mm counted once and indurations of ≥ 18 mm counted twice to obtain the estimated number of infections [23, 24].
Prevalence estimation using mixture models
This statistical approach has been extensively used in situations where a population is composed of several subgroups and the main issue is to find out which subgroup one individual belongs to . Neuenschwander et al  applied this method to the analysis of LTBI prevalence, initially restricted to an unvaccinated population and later on was extended to BCG-vaccinated populations [16, 17].
In our setting there are four possible subgroups according to tuberculin induration response: individuals with null induration after intradermal injection, individuals with reaction due to NTM, individuals with reaction due to previous BCG vaccination and individuals with reaction due to LTBI. Therefore, for BCG-vaccinated children the mixture model will be comprised of four components in the following sequential order :
y j ≈ p 0 + p 1 f EM (y j , θ1) + p 2 f BCG (y j , θ2) + p 3 f TB (y j , θ3) (1).
where y j represents the induration (in millimetres) for the j thindividual, is an indicator function (takes value 1 if y j = 0 and 0 otherwise), f EM , f BCG and f TB represent the distribution for each subgroup respectively: environmental mycobacteria reaction, BCG vaccination reaction and LTBI reaction, being, p 1, p 2 and p 3 the corresponding probabilities and p 0 = 1-(p 1+ p 2+p 3) the probability associated with null induration. For unvaccinated children, the f BCG component is excluded resulting in a three-component mixture model. In this model, the order of the components reflects the assumption of increasing component influences on the induration sizes. A similar model was derived for the sensitin analysis, changing the component order (f BCG , f TB and f EM ).
To analyse the tuberculin and sensitin responses, we tested different parametric distributions for each component, namely, normal, lognormal and Weibull, as suggested by Neuenschwander et al. . Their Bayesian approach to parameter estimation led to convergence problems in our data, and therefore we used an alternative procedure within the empirical-Bayes framework to tackle this problem: the EM estimation algorithm  in the presence of incomplete data. This method need not incorporate prior information from subject-matter knowledge. Parameter confidence intervals have been derived using parametric bootstrap techniques , which have already been used in mixture models applications . These methods were programmed in the freely available R software  [see Additional file 2].
I-Participation and BCG coverage
The survey fieldwork lasted five months as scheduled. Out of the 8,430 children with a registered age of seven, 885 children were excluded for the following reasons: 638 for not having obtained parents' written consent, 146 because the results of the tests could not be read, 21 because their vaccination status was unclassifiable and 95 for miscellaneous reasons (not being at school on the survey days, refusing the test, delay in receiving permission, concurrent minor illness and history of preventive therapy). Fifteen out of the 95 children were nevertheless included because a previous tuberculin skin test reaction size diameter had been recorded. Hence, our analyses are based on 7,545 children (49.6% girls and 50.4% boys). Table 1 displays the distribution of children according to BCG vaccination status The ungrouped induration data is provided in an additional [see Additional file 1].
II-Tuberculin and sensitin induration distributions
Using 5-mm interval categorisation (Table 2), it is clear that the proportion of children having large induration diameters is larger amongst BCG-vaccinated children. Regardless of the vaccination status, more people react to sensitin compared to tuberculin. Analyses of induration distribution by sex and local health area did not show statistically significant differences.
III-Prevalence estimation using cut-off points
LTBI prevalences according to the commonly used cut-off points are displayed in Table 3. Tuberculin results show broad variation both according to vaccination status and the cut-off point method selected. However, sensitin results are quite stable regardless of vaccination status if strict criteria are used. If a less stringent cut-off point is chosen (including doubtful reactors, that is, those with 6–11 mm induration), this stability is lost (unvaccinated prevalence 11.59%, BCG-vaccinated prevalence 26.25%).
IV-Mixture models analysis
Among the theoretical distributions tested, the lognormal and Weibull distributions seemed too flexible to provide an adequate fit to our data, which are characterised by showing a high proportion of zero indurations and not having clearly identifiable modes for the components. In contrast, the normal mixture models provided a proper fit as reflected in comparatively good BIC results and good adjustment to the observed distribution. In Table 4 we provide the parameters' estimates of the component of interest for each mixture model.
The LTBI estimated prevalence in the BCG-vaccinated group, although showing wide confidence intervals, is higher than expected, especially as compared with its counterpart prevalence figure in the unvaccinated group.
The observed induration distributions, the components of the mixture model and the normal mixture model are represented in figures 1a and 1b. On one hand, the fit of the mixture model is better in the BCG-vaccinated group, because it is based on many more people and three components are in use. On the other hand, in the unvaccinated group the components are more clearly distinguished. Finally, figures 1c and 1d show the probability of LTBI as a function of induration among unvaccinated and BCG-vaccinated children respectively. It is inferred from the first figure that indurations of 10 mm carry an estimated probability of infection of 35%, and it sharply arises as the induration increases from 10 to 14 mm, attaining a probability of infection of 97% for indurations of 14 mm (see Figure 1c). In contrast, the induration needed to consider LTBI with a probability of 98% among BCG-vaccinated people is 17 mm, 3 mm higher than among unvaccinated.
M avium estimated prevalence (Table 4) is very similar regardless of the vaccination status, around 7.5%. Similarly to figures 1a and 1b, figures 2a and 2b show a better fit of the mixture model in the BCG-vaccinated group. The relative influence of the last component seems to be bigger for the sensitin test than for the tuberculin test.
Finally, figures 2c and 2d show the estimated probabilities of M avium infection as a function of induration diameters among unvaccinated and BCG-vaccinated children respectively. The first figure shows that indurations of 8 mm carry a probability of infection of 25%, and it sharply arises as the induration increases from 8 to 12 mm, attaining a probability of infection of 98% for indurations of 12 mm. In contrast, the induration needed to consider M avium infection with a probability of 98% among BCG-vaccinated people is 15 mm, 3 mm higher than among BCG-unvaccinated people.
Discussions and conclusions
We present the results of a tuberculin survey conducted in a western European population characterised by a low-intermediate TB incidence rate and high vaccination coverage. The prevalences of LTBI and M avium infection are estimated using different mixture models based on empirical-Bayes estimation methods and compared to traditional estimation approaches. Additionally, a decision-aid tool that capitalises on our own data and provides estimates of the LTBI probability as a function of induration instead of externally drawn cut-off proposals is presented. This approach and the results obtained shed some light on a complex issue with remarkable importance both from clinical and public health perspectives.
In a setting with low-intermediate tuberculosis incidence rates, high BCG coverage and a suspected presence of MNT, data of good quality is essential. For this reason, we designed this study taking special care in controlling the quality of all its procedures. Skilled personnel were enrolled, and a rigorous training process was undertaken in order to ensure homogeneous administration and reading of the derivatives . Double and independent readings are not, as far as we know, a routine procedure in tuberculin and sensitin surveys.
Two features of this procedure deserve to be highlighted: the rapid execution of the fieldwork which ensures stability of conditions and the high participation rate (about 90%). Furthermore, the classification of vaccination statuses was very stringent, requiring at least two out of the three criteria to be positive: certificate of BCG vaccination, register confirmation and BCG scar .
Regarding the analysis of our TST results, it is remarkable how the method used influences the estimated LTBI prevalence (Table 3). Almost 7% out of the 5277 children presented an induration size greater than or equal to 10 mm., and therefore this cut-off point does not seem to be specific enough in the detection of LTBI due to the presence of NTM, BCG included . When the "ø ≥ 14 mm × 1.22" method is utilised, almost 3% of the children are classified as infected by M tuberculosis. Using the mirror method [11, 14, 24], the LTBI prevalence would be 1.46%, which is more in agreement with tuberculosis rates in our region. Similar conclusions are reached by Salaniponi et al.  who also faced similar problems when applying these criteria in a different socio-demographic context.
Concerning environmental mycobacteria, there has not been any study so far in our country with sensitin using a sample as large as ours. Using the Palmer and Edwards criteria  referred to by Bosman , at least 5.2% out of the 2268 children who received sensitin as the first skin test would be infected with M. avium. These children were given TST one week later, and 3% of those who have had an induration size equal to or larger than 5 mm to the sensitin had null reaction to the tuberculin .
When comparing sensitin reactions with tuberculin ones using 5 mm-interval induration categories, the proportion of children with positive responses is clearly higher for the sensitin (Table 2). These findings provide additional support to the idea that exposure to environmental mycobacteria is common in areas with low tuberculosis infection rates and its effect should be accounted for when studying LTBI prevalence [34, 35].
Because of the difficulty of disentangling the origin of the induration, we chose mixture models since they provide an appropriate framework in situations where a measurable variable (diameter of induration) is observed but the subpopulation (representing the main cause of skin test reaction) to which each individual belongs is unknown. From this analysis, we can obtain not only prevalence estimations but also information about the characteristics of each subpopulation and the individual classification probabilities. The estimation procedure of the parameters of a mixture model in general can be difficult when the number of components is large and/or when the component distributions overlap to a great extent. For this reason, it is more difficult to analyse BCG-vaccinated children than unvaccinated children, where convergence problems may arise in the estimation procedure. To minimise this problem, we analysed the data assuming several distributions for the components. Eventually, the normal mixture model provided both a reasonably good fit to the data and clinically interpretable estimates. We fitted the models using an empirical Bayesian analysis, which is a time-saving alternative to the full Bayesian approach used in this context by other authors , with the additional advantage that there is no need to include prior knowledge. To test this alternative estimation method for robustness, it has also been applied to data obtained from a rather different setting (summary data from the 1975 tuberculin survey in Korea analysed in  using full Bayesian techniques and kindly provided by HL Rieder) with very similar results and interpretation provided by both methods.
The interpretation of the mixture model results show that the prevalence of LTBI among unvaccinated children is about 1.5%. On the other hand, the prevalence of LTBI among BCG-vaccinated children is about 4.9%. These figures, although lower than those obtained with the 10 mm cut-off approach, still remain higher than we expected. We think that this fact could be due to the low expected prevalence in this region, which makes it difficult for the mixture model to estimate the relevant last component (attributed to LTBI) with sufficient precision. In spite of the fact that our study enrolled almost 90% of the target population, it may be possible that among non-participants LTBI would be more frequent, contributing to the imprecision of the last component. However, for populations with a higher presence of Mycobacterium tuberculosis, the model should perform better, even in BCG-vaccinated populations. This is reinforced by our findings regarding M. avium infection, with a higher prevalence in our study.
This mixture model approach, hence, in our view provides important information concerning the effect of NTM exposure on skin reactions and also about LTBI and M avium infection probabilities according to induration size (figures 1 and 2).
Focusing on the probability of LTBI (figures 1c and 1d), it seems that on average, exposure to previous BCG vaccine shifts around 3 mm towards greater induration size to attain the same probability of infection in our study population. In particular, to attain a LTBI probability of around 97%, the observed induration is 14 mm. for unvaccinated and 17 mm. for BCG-vaccinated children, results that were similar regardless of the distributional assumptions adopted. This finding agrees with current guidelines for LTBI treatment  and can be considered as a useful, induration size calibrated aid to be used in addition to other epidemiological and clinical information to make individual decisions.
Regarding M avium infection, some 7–8% of the children could be infected according to the mixture model results, without significant differences among BCG-vaccinated and unvaccinated, with both prevalences being higher than those obtained from the TST. BCG vaccine seems to have a similar effect on sensitin and tuberculin induration responses, requiring an upward shift of an average of about 3 mm to attain the same probability of infection.
Comparing the probability of infection with M tuberculosis and M avium as a function of induration, it is clear that the diameters needed to consider infection are lower for sensitin than for tuberculin. According to the results for the sensitin test, 8 mm of induration for unvaccinated and 10 mm for BCG-vaccinated children show a probability of M avium infection of about 25%, attaining a probability of infection above 98% with indurations of 12 and 15 mm for unvaccinated and BCG-vaccinated respectively.
In summary, we can conclude that the survey carried out using both tuberculin and sensitin skin tests supports the fact that the presence of Environmental Mycobacteria in our region cannot be neglected. The approach proposed here to analyse the results of the TST provides an alternative method for LTBI prevalence estimation based on the information obtained from our own local setting data instead of using external criteria and endorses current LTBI treatment decisions in populations similar to ours . Recently developed and currently available blood tests represent the biggest change since introduction of the TST test for the diagnostic and therapeutic management of LTBI in developed countries . CDC guidelines recommend utilising whole blood IFN-γ test QFT-G instead of TST although a combined use of both tests may increase the diagnostic yield insomuch as specificity is improved . This combined strategy may also maximise cost-effectiveness [37, 38], at least in adults as recent research in children suggests in them tests interpretation can be more complex . If this is the case our method may add valuable information and we recommend its use to analyse TBC surveys, especially in areas with intermediate-high TBC prevalences in order to better deal with this still challenging problem.
Maher D, Raviglione M: Global epidemiology of tuberculosis. Clin Chest Med. 2005, 26: 167-182. 10.1016/j.ccm.2005.02.009.
American Thoracic Society, Centers for Disease Control and Prevention: Targeted tuberculin testing and treatment of latent tuberculosis infection. Am J Respir Crit Care Med. 2000, 161: S221-47.
Mori T, Sakatani M, Yamagishi F, Takashima T, Kawabe Y, Nagao K, Shigeto E, Harada N, Mitarai S, Okada M, Suzuki , et al: Specific detection of tuberculosis infection: an interferon-γ-based assay using new antigens. Am J Respir Crit Care Med. 2004, 170: 59-64. 10.1164/rccm.200402-179OC.
Davies PDO, Drobniewski F: The use of interferon-γ-based blood tests for the detection of latent tuberculosis iinfection. Eur Respir J. 2006, 28: 1-3. 10.1183/09031936.06.00047606.
Brodie D, Schluger NW: The diagnosis of tuberculosis. Clin Chest Med. 2005, 26: 247-271. 10.1016/j.ccm.2005.02.012.
Menzies R, Vissandjee B: Effect of Calmette Guérin vaccination on tuberculin reactivity. Am Rev Respir Dis. 1992, 145: 621-625.
Rieder H: Epidemiologic basis of tuberculosis control. 1999, Paris: International Union Against Tuberculosis and lung Disease
Palmer CE, Long MW: Effects of infection with atypical mycobacteria on BCG vaccination and tuberculosis. Am Rev Respir Dis. 1996, 94: 553-568.
Rieder HL: Methodological issues in the estimation of the tuberculosis problem from tuberculin surveys. Tubercle Lung Dis. 1995, 76: 114-121. 10.1016/0962-8479(95)90552-9.
Departamento de Sanidad del Gobierno Vasco. Manual de vacunaciones. (Spanish). 2001, Vitoria-Gasteiz: Servicio Central de Publicaciones del Gobierno Vasco, 1
Bleiker MA, Sutherland I, Styblo K, ten Dam HG, Misljenovic O: Guidelines for estimating the risks of tuberculous infection from tuberculin test results in a representative sample of children. Bulletin of the International Union Against Tuberculosis and Lung Disease. 1989, 64: 7-12.
Rust P, Thomas J: A method for estimating the prevalence of tuberculous Infection. Am J Epidemiol. 1975, 101: 311-322.
Palmer CE, Edwards LB: Epidemiological studies of tuberculin sensitivity. Am J Hyg. 1958, 68: 213-
Bosman MCJ, Swai OB, Kwamanga DO, Agwanda R, Idukitta G, Mislijenovic O: National tuberculin survey of Kenya. Int J Tuber Lung Dis. 1998, 2 (4): 272-280.
Neuenschwander BE, Zwahlen M, Kim SJ, Engel RR, Rieder HL: Trends in the prevalence of infection with Mycobacterium tuberculosis in Korea from 1965 to 1995: an analysis of seven surveys by mixture models. Int J Tuber Lung Dis. 2000, 4: 719-729.
Neuenschwander BE, Zwahlen M, Kim SJ, Lee EG, Rieder HL: Determination of the prevalence of Infection with Mycobacterium tuberculosis among persons vaccinated with Bacillus Calmette-Guérin in South Korea. Am J Epidemiol. 2002, 155: 654-663. 10.1093/aje/155.7.654.
Dubuis M, Fiekert K, Johnston M, Neuenschwander BE, Rieder HL: A tuberculin skin test survey among Afghan children in Kabul. Int J Tuber Lung Dis. 2004, 8 (9): 1065-1072.
Titterington DM, Smith AFM, Makov UE: Statistical Analysis of Finite Mixture Distributions. 1985, New York: John Wiley & Sons
MPTR Working Group: Tuberculosis in Spain: Results of the Multicentre Project for Tuberculosis Research. 2000, Madrid: The Carlos III Institute of Public Health
Sanz MA, Errézola M: La infección tuberculosa en Euskadi. Resultados de la encuesta tuberculínica, Vizcaya 1983–1985. (Spanish). 1986, Vitoria: Departamento de trabajo, Sanidad y Seguridad Social. Gobierno Vasco
Won Reyn CF, Magnusson M: Use and standardisation of the mycobacterium avium sensitin skin-test reagent. CID. 1995, 21: 1527-1528.
American Thoracic Society, Centers for Disease Control: The tuberculin skin test. Am Rev Respir Dis. 1981, 124: 356-363.
Styblo K: The first round of the National Tuberculin Survey in Tanzania, 1983–1987. TSRU Progress Report. 1992, 1: 127-164.
Salaniponi FML, Kwanjana J, Veen J, Misjelnovic O, Borgdoff MW: Risk of infection with Mycobacterium tuberculosis in Malawi: national tuberculin survey 1994. Int J Tuber Lung Dis. 2004, 8 (6): 718-723.
Böhning D, Schalattman P, Lindsay B: Computer – assisted analysis of mixtures (C.A.MAN):statistical algorithms. Biometrics. 1992, 48: 283-303. 10.2307/2532756.
Dempster AP, Laird NM, Rubin DB: Maximum likelihood from incomplete data via the EM algorithms. J Roy Statist Soc. 1977, 39: 1-38.
Davison AC, Hinkley DV: Bootstrap Methods and Their Application. 1997, London: Cambridge University Press
Ugarte MD, Militino AF, Ibáñez B: Confidence intervals for relative risks in disease mapping. Biom J. 2003, 45: 410-425. 10.1002/bimj.200390022.
R Development Core Team: The Comprehensive R Archive Network. October 2005 (cited 2 November 2005). 2005, [http://lib.stat.cmu.edu/R/CRAN]October (cited 2 November 2005)
Arnadottir T, Rieder HL, Trébucq A, Waaler HT: Guidelines for conducting tuberculin skin test surveys in high prevalence countries. Tuber Lung Dis. 1996, 77 (Suppl): 1-20. 10.1016/S0962-8479(96)90127-6.
Young TK, Mirdad S: Determinants of tuberculin sensitivity in a child population covered by mass BCG vaccination. Tuber Lung Dis. 1992, 73: 94-100. 10.1016/0962-8479(92)90062-O.
Wang L, Turner MO, Elwood RK, Schulzer M, FitzGerald JM: A meta-analysis of the effect of Bacille Calmette Guérin vaccination on tuberculin skin test. Thorax. 2002, 57: 804-809. 10.1136/thorax.57.9.804.
Villate JI, Cabriada V, Sanz A, Urcelay MI, Galarza A, Díez I: Infección por Mycobacterium avium en la población infantil de Bizkaia. Influencia de la BCG. (Spanish). Archc Bronconeumol. 1999, 35 (S2): 52-
Larson LO, Magnusson M, Skoogh BE, Lind A: Sensitivity to sensitins and tuberculin in Swedish children. IV. The influence of BCG-vaccination. Eur Respir J. 1992, 5: 584-586.
Kröger L, Katila ML, Korppi M, Pietikäinen M: Skin reactivity to atypical mycobacteria among healthy Finnish preschool children vaccinated with BCG vaccine at birth. Acta Paedriatr. 1992, 81–5: 891-895.
Diel R, Ernst M, Döscher G, Visuri-Karbe L, Greinert U, Niemann S, Nienhaus A, Lange C: Avoiding the effect of BCG vaccination in detecting Mycobacterium tuberculosis infection with a blood test. Eur Respir J. 2006, 28: 16-23. 10.1183/09031936.06.00107005.
Diel R, Nienhaus A, Lange C, Schaberg T: Cost optimization of screening for latent tuberculosis in close contacts. Eur Respir J. 2006, 28 (1): 35-44. 10.1183/09031936.06.00011806.
Wrighton-Smith P, Zellweger J-P: Direct cost of three models for the screenin of latent tuberculosis infection. Eur Respir J. 2006, 28: 45-50. 10.1183/09031936.06.00005906.
Connell TG, Curtis N, Ranganathan SC, Buttery JP: Performance of a whole blood interferon gamma assay for detecting latent infection with Mycobacterium tuberculosis in children. Thorax. 2006, 61: 616-620. 10.1136/thx.2005.048033.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/6/240/prepub
We strongly acknowledge the invaluable help provided by Dr. HL Rieder (International Union Against Tuberculosis and Lung Disease, Paris, France) by kindly providing us with requested datafiles and making useful comments and remarks to draft versions of this paper.
We also wish to express our gratitude to other people for their contribution for the development of this work and we want to emphasize this acknowledgment to J.R. Palenzuela, who was the designer of the database, to M.I. Urcelay who coordinated the skin survey and to M.J. Olazaguirre, who entered the data in the database.
This research project was supported by the Health Department of the Basque Government and Osakidetza-Servicio vasco de Salud and by a grant from the "Fondo de Investigaciones Sanitarias" from the Spanish Ministry of Health, code 96/1171.
The author(s) declare that they have no competing interests.
JIV had the original idea for this research and participated in its design and coordination, BI performed the statistical analysis and helped to draft the manuscript, VC participated in the design of the research and conducted the skin survey, JIP participated in the statistical analysis and helped to draft the manuscript, JT worked in the data magement process and helped to draft the manuscript and AU participated in the statistical analysis and helped to draft the manuscript. All authors read and approved the final manuscript.