Determination of the Prevalence of Infection with Mycobacterium tuberculosis among Persons Vaccinated against Bacillus Calmette-Guérin in South Korea
Beat E. Neuenschwander1,2,
Marcel Zwahlen1,
Sang Jae Kim3,
Eun Gyu Lee3 and
Hans L. Rieder4
1 Federal Office of Public Health, Berne, Switzerland.
2 Institute of Social and Preventive Medicine, University of Berne, Berne, Switzerland.
3 Korean Institute of Tuberculosis, Seoul, Republic of Korea.
4 International Union Against Tuberculosis and Lung Disease, Paris, France.
 |
ABSTRACT
|
---|
The prevalence of tuberculous infection was estimated among 12,032 persons with a Bacillus Calmette-Guérin (BCG) vaccination scar and 7,788 persons without such a scar who participated in a nationwide tuberculin skin test survey conducted in the Republic of Korea in 1975. The analysis was built upon mixture models that captured the heterogeneity of indurations arising from tuberculous infection, cross-reactions due to infection with environmental mycobacteria, and BCG vaccination. The three distributions were allowed to vary by age, sex, and BCG vaccination status in the Bayesian manner, according to the prior opinion of the authors. Estimated prevalences of tuberculous infection were similar among persons with a BCG scar and persons without one: 7.5% (95% credibility interval (CI): 3.1, 12.5) and 5.2% (95% CI: 4.2, 6.3), respectively, at age 04 years and 87.3% (95% CI: 84.0, 90.2) and 84.0% (95% CI: 81.9, 85.8), respectively, at age 2529 years. From this analysis it can be concluded that mixture models allow investigators, for the first time, to estimate the prevalence of tuberculous infection not only in unvaccinated persons but also in the BCG-vaccinated population. Mixture models are a versatile tool for analyzing diagnostic test data and more general classification problems of considerable complexity.
Bayes theorem;; BCG vaccine;; mixture model;; Mycobacterium bovis;; prevalence;; tuberculin test;; tuberculosis
Abbreviations:
BCG, Bacillus Calmette-Guérin;; PPD, purified protein derivative;; PPD-S, purified protein derivative standard;; TU, tuberculin units.
 |
INTRODUCTION
|
---|
Estimation of the prevalence of infection with Mycobacterium tuberculosis by means of tuberculin skin testing is difficult in situations where sensitization with environmental mycobacteria or Mycobacterium bovis Bacillus Calmette-Guérin (BCG) is prevalent (1
). BCG vaccination is a well-recognized cause of nonspecific reactions to tuberculin, leading to a distribution closely similar to that due to M. tuberculosis infection in the period shortly after vaccination (2
), but the degree of postvaccinal tuberculin sensitivity may vary with the strain of BCG (3
). BCG may also induce cross-reactions for a prolonged time after vaccination (4
). For this reason, persons with a BCG vaccination scar have usually been excluded from the analysis of tuberculin skin test surveys (5


9
). However, BCG vaccination coverage has increased in many regions of the world. Therefore, results of analyses of tuberculin skin test indurations limited to children without prior BCG vaccination (i.e., without a BCG scar) may be increasingly nonrepresentative for the population at large. New methods of analyzing tuberculin skin test survey data in BCG-free individuals have recently been introduced and have been applied to surveys in the Republic of Korea (10
). These methods rely on mixture modeling, a statistical approach designed for situations with inherent structural heterogeneity that is unobserved but of crucial interest for the problem under scrutiny. For tuberculin skin testing in BCG-free individuals, such heteroge-neity consists of the two types of indurations arising from infection with M. tuberculosis and cross-reactions due to infection with environmental mycobacteria, where primary interest is centered around the prevalence of the former, although two-component mixture models also provide valuable insight into the distributions of both reactions.
The analysis presented here goes a step further. To our knowledge, it is the first attempt to estimate prevalence of infection for BCG-vaccinated individuals, a group that has hitherto been considered out of analytical reach because of the problem of cross-reactions arising from prior BCG vaccination. This additional source of heterogeneity necessitates the use of three-component mixture models that are more complex than two-component mixture models used in the population of unvaccinated individuals. Moreover, BCG reactions are particularly disturbing, because they are usually larger and more prevalent shortly after vaccination and fade away later. A further assumption is that BCG vaccination does not prevent infection but protects against the consequences of infection (11
). However, because of the considerable overlap between the distribution of BCG reactions and those resulting from M. tuberculosis infection, it might be unrealistic to assume that an analysis of the vaccinated population alone would provide useful estimates of infection prevalence. In this article, we put forward an alternative and more promising framework that borrows information from the group of unvaccinated individuals. It consists of a simultaneous analysis of both the vaccinated and unvaccinated populations, where in the latter, estimation of the distributions of infections due to M. tuberculosis and cross-reactions due to environmental mycobacteria turns out to be successful (10
); it further assumes that the two distributions are similar in vaccinated and unvaccinated individuals. The latter assumption's being fulfilled seems to be a necessity in order to make an estimation of the third mixture component (distribution of cross-reactions) feasible. However, if this assumption is not valid, or if the size of the unvaccinated population is too small for estimating the distributions of cross-reactions and infections, estimation of the three-component mixture distribution in the vaccinated population is hardly possible, meaning that prevalence and mixture component estimates would be unstable or too imprecise to be of practical value.
For a population in which assignment of BCG vaccination is close to random, the suggested methodology will be regarded as successful if the ratio of the infection prevalence in the BCG-vaccinated population to the prevalence in the unvaccinated population is close to 1. The South Korean survey of 1975 provided us with a good opportunity to validate the chosen methodology, because both groups (persons with and without a BCG scar) were fairly large, lending credibility to the hypothesis of similar infection prevalences in both groups. Furthermore, if the proposed methodology is considered successful, it can also be applied in populations with high BCG coverage wherein the vaccinated and unvaccinated might differ in terms of the prevalence of tuberculous infection. This would also allow much more cost-effective conduct of such surveys.
 |
MATERIALS AND METHODS
|
---|
Survey methodology
Between 1965 and 1995, the Korean National Tuberculosis Association conducted seven tuberculosis and tuberculin skin test surveys at 5-year intervals. The survey methodology has been described elsewhere in detail (12
). Briefly, all surveys were based on multistage stratified cluster sampling. A total of 190 clusters were sampled, using as a sample unit a modified enumerating district, defined according to the census. Two thirds of the clusters were in an urban stratum and one third were in a rural stratum, all clusters having the same sampling fraction. Enumeration districts with special mass accommodations (such as military bases, dormitories, or hospitals) were excluded.
For this analysis, we selected the survey conducted in 1975, because our interest was in the analysis of the estimated prevalence of M. tuberculosis infection at a time when BCG coverage was sufficiently high to allow meaningful analysis of BCG-vaccinated individuals.
Tuberculin skin testing
The dose of international standard tuberculin was calibrated. It was established to be 5 tuberculin units (TU) of purified protein derivative (PPD) standard (PPD-S) (13
, 14
). Subsequently, the introduction of a stabilizing agent (a polyoxyethylene derivative of sorbitan mono-oleate) demonstrated that the number of units of tuberculin PPD RT-23 (Statens Serum Institut, Copenhagen, Denmark) could be reduced to 2 TU to achieve a response that was bioequivalent to 5 TU of PPD-S (15
). International recommendations for the dosage using tuberculin PPD RT-23 have changed over time (16
), which complicates the issue further. For this reason, the Korean Institute of Tuberculosis opted to use the dosage of 1 TU of PPD RT-23 originally recommended by the World Health Organization throughout all of the surveys to ensure comparability of the data. Thus, in addition to clinical, radiographic, and bacteriologic examinations, individuals under 30 years of age were tuberculin-tested by trained personnel with 1 TU of PPD RT-23 in all seven surveys. Approximately 72 hours after application of the test, the diameter of induration of the reaction was measured transversely to the long axis of the arm. The number of reactors by exact millimeter of induration was recorded up to a reaction size of 29 mm. Persons with indurations of 30 mm or more were recorded as having indurations of 30 mm.
BCG vaccination policy in South Korea
Before 1952, there was no systematic BCG vaccination program in South Korea. Subsequently, systematic vaccination targeting primary school, junior high school, and high school students was implemented. Coverage was fairly low until 1963, when the scheme was expanded to include preschool children. From 1966 onward, vaccination policy was altered, and preschool children were vaccinated without prior tuberculin testing, while schoolchildren were only vaccinated (or revaccinated) if they had a negative tuberculin skin test result. From 1967 onward, vaccination without prior tuberculin testing was extended to first graders if they had no visible BCG scar, but for sixth graders, BCG vaccine was given in the presence of a negative tuberculin skin test irrespective of the presence or absence of a BCG scar. This scheme remained unaltered until 1997. In 1975, when the population of South Korea was 35 million, 2.5 million doses of BCG vaccine were given (Korean Institute of Tuberculosis, unpublished data).
Data set for the current analysis
The population surveyed in 1975 was stratified by 5-year age group (04, 59, 1014, 1519, 2024, and 2529 years), by sex, and by the presence or absence of a BCG scar upon physical inspection (table 1). A total of 19,820 individuals with recorded tuberculin skin test results were included in the 1975 tuberculin skin test survey. Of these, 12,032 (60.7 percent) had a BCG scar and 7,788 (39.3 percent) did not (17
).
View this table:
[in this window]
[in a new window]
|
TABLE 1. Distribution of individuals in a tuberculin skin test survey, by age, sex, and the presence or absence of a BCG* scar, South Korea, 1975
|
|
Analysis by mixture models
The statistical analysis was based on mixture models (18
) that captured the heterogeneity of indurations arising from infection with M. tuberculosis, cross-reactions due to infection with environmental mycobacteria, and reactions due to BCG vaccination. The distributions of these three types of reactions, together with their relative weights, were the building blocks of the mixture model (table 2). It was assumed that the distribution of indurations in the group of individuals with a BCG scar was a mixture of three components: a distribution of infections with M. tuberculosis (f), a distribution of cross-reactions due to environmental mycobacteria (g), and a distribution of BCG reactions (h) due to prior vaccination. In the group without a BCG scar, only the distributions of M. tuberculosis infections (f') and cross-reactions (g') needed to be considered. The prevalences of M. tuberculosis infection were q1 in the group with a BCG scar and p1 in the group without a BCG scar. Because of the possibility of multiple reactions, it was further assumed that an individual's reaction to M. tuberculosis infection always dominates the other two reactions (in the sense of being larger in size than the other two). As a consequence, q2 (p2) and q3 cannot be interpreted as the prevalences of cross-reactions and BCG reactions, respectively, since a reaction due to M. tuberculosis infection never excludes the possibility of a smaller (unobserved) cross-reaction or BCG reaction. The mixture component distributions were assumed to be Weibull, lognormal (both parameterized by their 25 percent and 75 percent quantiles), or gamma distributions (parameterized by mean and standard deviation). Hence, the two mixture distributions m(x) and m'(x) comprise 15 unknown model parameters (five mixing weights and 10 quantile parameters).
View this table:
[in this window]
[in a new window]
|
TABLE 2. Model structure for mixture components and mixture distributions in a tuberculin skin test survey, South Korea, 1975
|
|
It was assumed that sex and age not only affect the prevalences (mixture weights) but might have an influence on the mixture component distributions (f, g, h, f', g') as well. When taking these covariates into account (12 strata), the full model would consist of 180 free parameters. However, this model would unrealistically assume that the mixture component distributions are completely unrelated to each other with regard to the different covariate strata. For example, the distribution of infections for males aged 1014 years with a BCG scarfM3(x), for exampleis likely to be close to 1) the distribution of infections for males aged 59 years, fM2(x), and males aged 1519 years, fM4(x); 2) the distribution of infections for females aged 1014 years, fF4(x); and 3) the distribution of infections for males aged 1014 years without a BCG scar, f'M3(x). On the other hand, assuming that the component distributions are not affected at all by sex, age, and BCG scar status might also be questionable. Therefore, a compromise driven by knowledge of the subject matter and the observed data seems reasonable, where the former has to be included in the statistical model as external knowledge and the latter enters the analysis automatically during the process of statistical estimation.
Since the analysis was conducted in a Bayesian way (19
, 20
), subject-matter knowledge was included in the form of prior opinion. Noninformative (uniform) prior distributions were used for all prevalence parameters (mixture weights), whereas prior opinion about the relations between the param-eters (25 (75) percent quantiles) of the component distri-butions was specified in the form of constraints (intervals) for the ratios of related parameters. Uniform distributions were assumed within these intervals. Elicitation of prior opinion among the authors led to the following constraints for the ratios of the 25 (75) percent quantiles (table 3):
- Differences of quantiles with regard to sex were assumed to be, at most, 10 percent for all distributions.
- Differences of quantiles between subsequent age classes were assumed to be, at most, 10 percent for the distribution of M. tuberculosis infections and, at most, 15 percent for the distribution of cross-reactions. For the distribution of BCG reactions, a decrease of 1025 percent between subsequent age classes was assumed.
- Differences of quantiles between the distributions of M. tuberculosis infections in the groups with and without a BCG scar (f vs. f') were assumed to be, at most, 10 percent. The same assumption was made for the distributions of cross-reactions (g vs. g').
View this table:
[in this window]
[in a new window]
|
TABLE 3. Prior specifications (uniform distributions) for ratios of 25 (75) percent quantiles of related mixture component distributions in a tuberculin skin test survey, South Korea, 1975
|
|
The posterior distribution of all model parameters was obtained by combining the likelihood function (given by the mixture distributions m(x) and m'(x) for all covariate strata) and the prior distribution via Bayes' theorem. Because of the complexity of the model, computation of the posterior distribution was analytically intractable, but Markov chain Monte Carlo methods (21
, 22
) offered a feasible alternative. The simulations were done by a random-walk Metropolis algorithm, with a burn-in period of 100,000 and a sample size of 500,000 iterations. All posterior quantities were extracted from a thinned sample comprising 2,000 values. Convergence of the algorithm was good; we checked for it by comparing runs starting from different initial values.
 |
RESULTS
|
---|
Analysis of the class of mixture models based on Weibull, lognormal, and gamma distributions suggested a best model with Weibull distributions for M. tuberculosis infections and BCG reactions and lognormal distributions for cross-reactions (table 4). The log-likelihood values showed that the second-best model (with lognormal distributions instead of Weibull distributions for BCG reactions) might be a reasonable choice as well. However, since the posterior odds in favor of the first model were overwhelming (independent of the model selection criterion used), the results of the final analysis were based on the Weibull-lognormal-Weibull model. Posterior predictive model checks were done by comparing observed induration frequencies for seven intervals (0, 12, 37, 812, 1317, 1822, and >22 mm) and 24 covariate strata with their 95 percent prediction intervals. An observed frequency falling outside of its prediction interval was considered a predictive failure. The model showed a predictive failure rate of approximately 5 percent, and we thus concluded that calibration of the model was reasonably good. Examples of observed induration frequencies and model estimates are displayed in figures 1 and 2 by age, sex, and vaccination status.

View larger version (22K):
[in this window]
[in a new window]
|
FIGURE 1. Observed distribution of tuberculin skin test induration sizes (histograms), mixture distribution (solid line on left panel), distribution due to infection with M. tuberculosis (solid thick line on right panel), distribution due to Bacillus Calmette-Guérin (BCG) vaccination (solid thin line on right panel), and distribution due to infection with environmental mycobacteria (dashed line on right panel) among children aged 04 years, by sex and presence or absence of a BCG scar, tuberculin skin test survey, South Korea, 1975.
|
|

View larger version (24K):
[in this window]
[in a new window]
|
FIGURE 2. Observed distribution of tuberculin skin test induration sizes (histograms), mixture distribution (solid line on left panel), distribution due to infection with M. tuberculosis (solid thick line on right panel), distribution due to Bacillus Calmette-Guérin (BCG) vaccination (solid thin line on right panel), and distribution due to infection with environmental mycobacteria (dashed line on right panel) among persons aged 2529 years, by sex and presence or absence of a BCG scar, tuberculin skin test survey, South Korea, 1975.
|
|
For the 1975 tuberculin skin test survey in South Korea, the model supported the assumption that the mean induration size of reactions due to M. tuberculosis infection varied by age, sex, and the presence of a BCG scar (table 5). Mean values were consistently higher (by 0.51 mm) among females than among males across all age groups, and they were somewhat lower among persons with a vaccination scar as compared with those without one. On the other hand, mean induration sizes of cross-reactions due to environmental mycobacteria seemed fairly constant with regard to sex and BCG vaccination status and showed only a slight increase over age (22.5 mm). The mean values for reactions due to BCG vaccination were smaller than those for M. tuberculosis infection and decreased with increasing age (and time elapsed since vaccination), from approximately 11 mm to 5 mm (see figures 1 and 2).
View this table:
[in this window]
[in a new window]
|
TABLE 5. Mean (expected) values (in mm) for distributions of Mycobacterium tuberculosis infections, cross-reactions due to environmental mycobacteria, and BCG* reactions due to prior vaccination in a tuberculin skin test survey, South Korea, 1975
|
|
The prevalence of M. tuberculosis infection increased rapidly with age, from an average of 5 percent in the youngest age group (04 years) to almost 90 percent in the oldest age group (2529 years) (table 6). The credibility intervals were wide in the youngest age group of BCG-vaccinated individuals. This finding was due to the very low prevalence of M. tuberculosis infection and the fact that the highest mean value for reactions due to BCG vaccination appeared in this age group. The total prevalence of tuberculous infection (BCG-vaccinated and unvaccinated individuals combined) is shown in figure 3.
View this table:
[in this window]
[in a new window]
|
TABLE 6. Prevalence (%) of infection with Mycobacterium tuberculosis, by age and sex, in a tuberculin skin test survey, South Korea, 1975
|
|

View larger version (9K):
[in this window]
[in a new window]
|
FIGURE 3. Estimated prevalence of infection with M. tuberculosis, combined for individuals with and without a Bacillus Calmette-Guérin scar, by age and sex, tuberculin skin test survey, South Korea, 1975. Symbols, posterior medians; lines, 95% credibility intervals.
|
|
When comparing the prevalences of tuberculous infection in the BCG-vaccinated and unvaccinated populations, the ratio of the two prevalences was close to 1, with the exception of the youngest age group, which had a highly uncertain ratio (figure 4). Moreover, the results for the unvaccinated population were very similar to those obtained in a previous analysis of unvaccinated individuals only (10
). Both findings support the conclusion that the chosen methodology based on mixture models allows one to estimate prevalences of tuberculous infection in the population of BCG-vaccinated and unvaccinated individuals.

View larger version (11K):
[in this window]
[in a new window]
|
FIGURE 4. Estimated ratio of the prevalence of tuberculous infection among those with and those without the presence of a Bacillus Calmette-Guérin scar, by age and sex, tuberculin skin test survey, South Korea, 1975. Symbols, posterior medians; lines, 95% credibility intervals.
|
|
 |
DISCUSSION
|
---|
The extent of transmission of M. tuberculosis in the community is obviously important, since it influences the extent to which future cases of tuberculosis will emerge. Styblo et al. (23
) have proposed utilizing a determination of the average annual risk of M. tuberculosis infection in the community as an index of transmission. This average annual risk of infection can be derived algebraically, under suitable conditions, from the prevalence of M. tuberculosis infection (23
, 24
). The approach is hampered by the limitations of the operating characteristics of the tuberculin skin test (1
) and the high BCG coverage in most societies, precluding a representative ascertainment of the prevalence of M. tuberculosis infection. Thus, limiting the estimation of tuberculous infection to persons without a BCG scar (as a proxy for vaccination) becomes increasingly questionable with increasing BCG coverage; furthermore, it would require ever-increasing sample sizes in tuberculin skin test surveys. It would therefore be desirable to identify a method that would allow estimation of the prevalence of tuberculous infection irrespective of prior BCG vaccination.
The opportunity to investigate the new method arose in studying the surveys conducted systematically and regularly in South Korea (17
, 25



30
). While all surveys provided information on the existence or absence of a scar resulting from BCG vaccination, the 1975 survey data allowed us to estimate the prevalence of tuberculous infection among both persons with a scar and persons without a scar at a time when BCG coverage was high but not yet great enough to make either the vaccinated group or the unvaccinated group obviously unrepresentative of the Korean population. The 1975 survey fell into a period when approximately 60 percent of the surveyed population had been vaccinated. Thus, we expected that the prevalences of M. tuberculosis infection would be similar in the vaccinated and nonvaccinated populations. The estimated ratio of infection prevalences (BCG-vaccinated:unvaccinated) was highly uncertain in the youngest age group (with documented recent vaccination) but close to 1 for the remaining age groups. This supports the conclusion that the chosen methodology was successful in estimating the prevalence of tuberculous infection in both BCG-vaccinated and unvaccinated individuals.
The statistical inference of tuberculous infection prevalences was based on mixture modeling. Mixture models provide the appropriate inferential framework in situations where heterogeneity (subpopulations) in an outcome variable is present and where the subpopulation to which an individual belongs is unknown. Quantities of interest can be the distributions of the subpopulations (mixture components), the prevalences of the subpopulations (mixture weights), or quantities related to them (e.g., individual classification probabilities). Mixture models offer great flexibility (number of components, types of component distributions) in describing epidemiologic applications. In the context of tuberculin skin test results, the population of persons with a BCG scar can be described reasonably well by a three-component mixture representing infection with M. tuberculosis, cross-reactions due to environmental myco-bacteria, and BCG reactions due to prior vaccination.
Even though the mixture models merely capture the indispensable characteristics of tuberculin skin testing, their complexity is, relative to the information provided by the data, already high. For example, when assuming that prevalences and mixture components are affected by sex and age (12 covariate strata), induration data for a specific sex/age stratum may hardly be informative enough to estimate the nine parameters of the three-component mixture (group with a BCG scar) or the six parameters of the two-component mixture (group without a BCG scar).
In general, estimation of parameters in mixture models can be difficult if the number of mixture components is large, and even more difficult if the component distributions overlap to a great extent. In the extreme case of complete overlap of component distributions, mixture weights (prevalences) are no longer identifiable. In practical situations, identifiability problems can arise to an extent that depends both on the complexity of the model and on the available data, and that is reflected in imprecise parameter estimates. When applying mixture analysis to tuberculin skin test data, we experienced problems of identifiability mostly in situations with aggregated data. Improved estimates can generally be obtained whenever covariate information is available and 1) mixture component distributions are similar or even identical (in location and spread) in the different covariate strata and 2) some covariate strata are particularly informative with regard to one of the mixture component distributions. The latter is typically the case if prevalences in some of the subpopulations are large. For example, in the Korean data on nonvaccinated persons, stratification with regard to the important covariate "age" turned out to be advantageous: It showed large prevalences of cross-reactions in the very young and large prevalences of infection in the elderly, revealing valuable information on the distribution of cross-reactions and infections, respectively. If we had ignored age and pooled the data from the subpopulations, information on the component distributions would have been obscured and estimation of mixture parameters made more difficult (less precise).
For the Korean data, a mixture model was chosen (among several competing models) that was flexible enough to incorporate covariate information (age, sex) while imposing reasonable constraints on the similarity of the mixture component distributions over the covariate strata. Moreover, the distributions of cross-reactions and infections were assumed to be similar for vaccinated and unvaccinated persons as well. The degree of similarity was introduced in the form of prior opinion that was agreed upon by the authors. No formal assessment of prior opinion was done, however.
For the statistical analysis of the mixture models (including prior opinion), we took the Bayesian stance, a statistical paradigm that has been shown to be not only theoretically sound (a long-recognized feature) but also versatile and generally practicable (31
). Bayesianism (full probability modeling) puts forward a general framework of statistical thinking that is based on a probability model for the data, the likelihood function, and a prior distribution for the unknown parameters of this model. Once these parts are specified, the approach is fully automatic, since there is no ambiguity in the estimation process (computation of the posterior distri-bution), which relies exclusively on the laws of probability (an axiomatic system). However, specification of the model and assessment of prior knowledge are not as straightforward because of their being heavily dependent on the epidemiologic problem (32
), and both are open to discussion, but primarily by epidemiologists, not statisticians. Since the focus of the inferential process is thereby shifted away from statistics to epidemiology, more realistic and better calibrated inference statements are likely to result.
The methodology adopted here demonstrates that the prevalence of infection with tubercle bacilli can be estimated successfully in persons with and without prior BCG vaccination. The distribution of reaction sizes is importantly affected by the amount of time elapsed since vaccination (2
), which is reflected by the diminishing size of BCG-attributable indurations with increasing age. The vaccine strain used may also affect the rapidity with which reactions due to BCG vaccination wane (3
). Nevertheless, this analysisto our knowledge, the first of its kindsuggests that tuberculin skin-test survey data from both vaccinated and unvaccinated populations might be utilized. This would be important in the most commonly prevailing situation, where BCG coverage is high and an analysis restricted to the nonvaccinated might seriously hamper an unbiased assessment of the epidemiologic situation in a community.
 |
ACKNOWLEDGMENTS
|
---|
Dr. Beat E. Neuenschwander was supported by Swiss Federal Office of Public Health contracts 5919 and 6895.
The authors thank J. I. Bai, who contributed to the preparation of tuberculin data sheets.
 |
NOTES
|
---|
Correspondence to Dr. Hans L. Rieder, Jetzikofenstrasse 12, 3038 Kirchlindach, Switzerland (e-mail: TBRieder{at}cs.com).
 |
REFERENCES
|
---|
-
Rieder HL. Methodological issues in the estimation of the tuberculosis problem from tuberculin surveys. Tuber Lung Dis 1995;76:11421.[ISI][Medline]
-
Baily GV, Narain R, Mayurnath S, et al. Trial of BCG vaccines in south India for tuberculosis prevention. Tuberculosis Prevention Trial, Madras. Indian J Med Res 1980;72(suppl):174.[Medline]
-
Comstock GW. Identification of an effective vaccine against tuberculosis. Am Rev Respir Dis 1988;138:47980.[ISI][Medline]
-
Comstock GW, Livesay VT, Woolpert SF. Evaluation of BCG vaccination among Puerto Rican children. Am J Public Health 1974;64:28391.[ISI][Medline]
-
Styblo K. Recent advances in epidemiological research in tuberculosis. Adv Tuberc Res 1980;20:163.[Medline]
-
Cauthen GM, Pio A, ten Dam HG. Annual risk of infection. Geneva, Switzerland: World Health Organization, 1988. (Document no. WHO/TB/88.154:134).
-
El Ibiary S, de Coster EJ, Tolba FM, et al. Trend in the annual risk of tuberculous infection in Egypt, 19501996. Int J Tuberc Lung Dis 1999;3:2949.[ISI][Medline]
-
Hong YP, Kim SJ, Lew WJ, et al. The seventh nationwide tuberculosis prevalence survey in Korea, 1995. Int J Tuberc Lung Dis 1998;2:2736.[ISI][Medline]
-
Odhiambo JA, Borgdorff MW, Kiambih FM, et al. Tuberculosis and the HIV epidemic: increasing annual risk of tuberculous infection in Kenya, 19861996. Am J Public Health 1999;89:107882.[Abstract]
-
Neuenschwander BE, Zwahlen M, Kim SJ, et al. 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 Tuberc Lung Dis 2000;4:71929.[ISI][Medline]
-
Sutherland I, Lindgren I. The protective effect of BCG vaccination as indicated by autopsy studies. Tubercle 1979;60:22531.[ISI][Medline]
-
Hong YP, Kim SJ, Kwon DW, et al. The sixth nationwide tuberculosis prevalence survey in Korea, 1990. Tuber Lung Dis 1993;74:32331.[ISI][Medline]
-
Furcolow ML, Hewell B, Nelson WE, et al. Quantitative studies of the tuberculin reaction. I. Titration of tuberculin sensitivity and its relation to tuberculous infection. Public Health Rep 1941;56:1082100.
-
World Health Organization. Comité d'Experts pour la Standardisation Biologique. Cinquième rapport. 7. Tuberculine. WHO Tech Rep Ser 1952;56:67.
-
Comstock GW, Edwards LB, Philip RN, et al. A comparison in the United States of America of two tuberculins, PPD-S and RT-23. Bull World Health Organ 1964;31:16170.[ISI]
-
Arnadottir T, Rieder HL, Trébucq A, et al. Guidelines for conducting tuberculin skin test surveys in high prevalence countries. Tuber Lung Dis 1996;77(suppl):120.
-
Ministry of Health and Social Affairs, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the 3rd tuberculosis prevalence survey in Korea1975. Seoul, South Korea: Korean Institute of Tuberculosis, 1976.
-
Titterington DM, Smith AF, Makov UE. Statistical analysis of finite mixture distributions. Chichester, United Kingdom: John Wiley and Sons, 1985.
-
Malakoff D. Bayes offers a "new" way to make sense of numbers. Science 1999;286:14604.[Free Full Text]
-
Gelman A, Carlin JB, Stern HS, et al. Bayesian data analysis. London, United Kingdom: Chapman and Hall Ltd, 1995.
-
Gilks WR, Richardson S, Spiegelhalter DJ. Markov chain Monte Carlo in practice. London, United Kingdom: Chapman and Hall Ltd, 1996.
-
Tanner MA. Tools for statistical inference: methods for the exploration of posterior distributions and likelihood functions. 3rd ed. New York, NY: Springer-Verlag, Inc, 1996.
-
Styblo K, Meijer J, Sutherland I. The transmission of tubercle bacilliits trend in a human population. Tuberculosis Surveillance Research Unit report no. 1. Bull Int Union Tuberc 1969;42:1104.[Medline]
-
Sutherland I. Recent studies in the epidemiology of tuberculosis, based on the risk of being infected with tubercle bacilli. Adv Tuberc Res 1976;19:163.[Medline]
-
Ministry of Health and Social Affairs, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the first tuberculosis prevalence survey in Korea1965. Seoul, South Korea: Korean Institute of Tuberculosis, 1966.
-
Ministry of Health and Social Affairs, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the 2nd tuberculosis prevalence survey in Korea1970. Seoul, South Korea: Korean Institute of Tuberculosis, 1971.
-
Ministry of Health and Social Affairs, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the 4th tuberculosis prevalence survey in Korea1980. Seoul, South Korea: Korean Institute of Tuberculosis, 1981.
-
Ministry of Health and Social Affairs, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the 5th tuberculosis prevalence survey in Korea1985. Seoul, South Korea: Korean Institute of Tuberculosis, 1986.
-
Ministry of Health and Social Affairs, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the 6th tuberculosis prevalence survey in Korea1990. Seoul, South Korea: Korean Institute of Tuberculosis, 1991.
-
Ministry of Health and Welfare, Korean Institute of Tuberculosis, Korean National Tuberculosis Association. Report on the 7th tuberculosis prevalence survey in Korea1995. Seoul, South Korea: Korean Institute of Tuberculosis, 1996.
-
Lindley DV. The philosophy of statistics. Statistician 2000;49:293337.
-
Greenland S. Probability logic and probabilistic induction. Epidemiology 1998;9:32232.[ISI][Medline]
Received for publication May 14, 2001.
Accepted for publication October 30, 2001.