 Research article
 Open Access
 Published:
A twostage Bayesian method for estimating accuracy and disease prevalence for two dependent dichotomous screening tests when the status of individuals who are negative on both tests is unverified
BMC Medical Research Methodology volume 14, Article number: 110 (2014)
Abstract
Background
Estimating the disease prevalence and test accuracy (sensitivity and specificity) for two dependent screening tests when the status of individuals who are negative on both tests is unverified represents a considerable challenge, as the disease rates for individuals negative on both tests are not identifiable without additional assumptions.
Methods
This article presents a unified framework for handling this nonidentifiability problem using twostep hierarchical informative prior on the sensitivities by twostage Bayesian modeling with the characterized by joint testing strategies based on the inherent attribute of screening/diagnostic tests. We assign a diffuse and less risky twostep hierarchical informative uniform prior to the sensitivities while assigning a uniform (0,1) prior distribution to the specificities and prevalence. Strategies for model evaluation, general global evaluations, and individual cell checking are presented. Simulations are conducted under various scenarios to evaluate the performance of the proposed method. Applications to real data are also presented to illustrate the potential impact and benefit of the proposed method.
Results
Our results indicate that when the priors of sensitivities are assigned as appropriate twostep hierarchical informative priors, or even in the absence of the priors for the specificities and prevalence, the parameters involved in this study can still be estimated well. The advantages and limitations of this method in solving such problems are discussed and compared with other twostage methods.
Conclusions
We developed a twostage Bayesian method for two dependent dichotomous screening tests with unverified individuals who are negative on both tests, and addressed the ad hoc model evaluation and checking procedures. The method can be understood easily and used conveniently by nonstatisticians.
Background
Screening programs for a specific disease or condition of interest are typically divided into two stages. In the first stage, a population of known size n is screened by simple and rapid but imperfect tests to classify people as likely or unlikely to have the disease or condition of interest. In the second stage, individuals who appear likely to have the disease or condition are examined further using a gold standard test to confirm the disease or condition. To improve the sensitivity of these screenings, two simple and rapid tests are often used in the first stage [1]. If either of the tests is positive, then a full evaluation of the correct disease status using a gold standard classification is undertaken in the second stage [2]. This strategy is widely used in screening for chronic diseases, infectious diseases, and animal diseases [2–4]. One characteristic of this strategy is that negatives in both tests are not verified by a gold standard test because the disease probability in the subpopulation is so low that further investigations are costly, unacceptable, and possibly unethical. Another characteristic is that the two tests in the first stage are often conditionally dependent on disease status and/or nondisease status [5]. If the two screening tests have a similar biological basis, as is often the case, the conditional dependence assumption is suitable. For example, when we use two fecal occult blood tests (occult blood tests and immunohistochemical methods) for colon cancer screening, a positive occult blood test is often accompanied by a positive immunochemical result; thus, the results of the two tests are dependent [2].
Estimating the disease prevalence and test accuracy (sensitivity and specificity) when the two tests are dependent and when individuals negative on both tests are unverified is a considerable challenge, as the disease rate for individuals who are negative on both tests is not identifiable without additional assumptions [5–9]. Li et al. reviewed several methods for resolving the nonidentifiability problem in the framework of classical statistics [8]. They classified these methods into two types. The methods of the first type are intended to solve the nonidentifiability problem by introducing additional assumptions about the association of the two tests. Walter et al. assumed that the two test errors were independent [9], while van der Merwe also assumed that the two tests were independent under the disease status but dependent under nondisease status [10]. The methods in the second case resolve the nonidentifiability problem directly by adding an assumption that is more likely to hold than the independence assumption [7]. Bohning and Patilea proposed a capturerecapture approach for screening using two diagnostic tests with information on disease status available for the test positives only [11]. Chu et al. proposed latent class models for screening studies using two screening tests, with a categorical disease status verified in positives test only [7]. However, these studies did not provide sufficient insight and details to solve the problem [8]. Li et al. presented a unified framework for solving the nonidentifiability problem by saturating the model using certain homogeneous association assumptions [8]. In their paper, although the five distinct models provided different estimators, they were all saturated models and provided the same maximum likelihood. Therefore, they are not differentiable in terms of goodnessoffit [8].
Recently, Bayesian analysis has been used increasingly to solve the nonidentifiability problem using prior information [12–17]. Briefly, the use of this analysis is appropriate for three scenarios. The first involves the use of Bayesian model averaging (BMA) to address this challenging estimation problem, as a different homogeneousdependent model can provide the same goodnessoffit for the data but with substantively different estimates [7]. The second scenario involves the application of BMA to inferences over conditional independence and dependence models for the estimation of disease prevalence in a situation involving two imperfect tests in the absence of a gold standard [18]. The third scenario involves the administration of a moderately effective but relatively cheap screening test to all subjects in the first stage, after which a gold standard test is performed on a subset of the high screen and low screen groups in the second stage [19].
This article presents a unified framework to handle the nonidentifiability problem by a twostage Bayesian model using twostep hierarchical informative uniform prior on the sensitivities. We provide detailed Bayesian modeling for stages one and two, the specification of the prior distribution, and the calculation of the posteriors of the parameter distributions. The strategies for model evaluation and checking are proposed. We illustrate our approach using an example and simulation and conclude this article with a discussion and some remarks.
Methods
Twostage Bayesian modeling
Data structure and estimated parameters
Table 1 presents the outcomes of two screening tests (T _{ 1 } , T _{ 2 }) from a total of n subjects, with disease statuses (D _{ + } , D _{  }) verified only for individuals with at least one positive test. D _{ + } represents a true diseased condition, whereas D _{  } represents the nondiseased, as determined by a gold standard test. T _{ i(+,)} (i = 1,2) represents the results of the test i. Let vector x = (x _{ 11 }, x _{ 10 }, x _{ 01 }, x _{ 00 }) denote the observed number for each combination of two test outcomes (1 or 0 for test positive or test negative). Let SI _{(+,)} (abbreviation of Simultaneous testing) represent the outcomes of the simultaneous testing in stage one, and let a _{ 11, } a _{ 10, } a _{ 01, } a _{ 00 } denote D _{ + } SI _{ + } , D _{  } SI _{ + }, D _{ + } SI, D SI, respectively. Because the full evaluation of the disease status is not performed if neither of the two tests is positive in the first stage, the frequencies of [a _{ 01 }] and [a _{ 00 }] in the bracket are unobserved, whereas the total frequency of both tests being negative a _{ 0 . } = [a _{ 01 }] + [a _{ 00 }] is known.
The prevalence is defined as π = P(D+). The sensitivity and specificity for the ith test are defined as Se _{ i } = P(T _{ i+ } D _{ + } ) and Sp _{ i } = P(T _{ i } D _{  } ), respectively. Furthermore, let CovDp and CovDn be the covariances between the two tests for diseases (CovDp) and for nondiseases (CovDn). Let ρ _{ D+ } and ρ _{ D } be the correlation coefficients under the conditions of disease and nondisease, and let PPV _{ JE } be the positive predictive value of the simultaneous testing. These parameters are defined as follows [12, 20, 21]:
where Se _{ JE } and Sp _{ JE } denote the joint sensitivity and specificity of the simultaneous testing (described in equation (6))
where Se _{11} = P(T _{1 +}, T _{2 +}D _{+}), Sp _{22} = P(T _{1 }, T _{2 }D _{})
Bayesian modeling for stage one
According to the screening programs, in the first stage (see Table 1, left column), two simple and rapid but potentially imperfect screening tests, denoted as T _{ 1 } and T _{ 2 }, are applied simultaneously to each of the units sampled. If either individual test is positive, this condition indicates a positive result of the joint test. The strategy is often called simultaneous testing and seeks to obtain a joint testing strategy with increased sensitivity [22]. Let vector Px = (Px _{ 11 } ,Px _{ 10 } ,Px _{ 01 } ,Px _{ 00 }) be the probabilities, respectively, from the vector x = (x _{ 11 }, x _{ 10 }, x _{ 01 }, x _{ 00 }). Supposing that n subjects are sampled randomly with prevalence π, the vector x has the following multinomial sampling distribution:
Considering that the two tests are conditionally dependent under disease status and nondisease status, the multinomial cell probabilities are given by
In stage one, the individuals are classified as likely to have the disease or condition of interest only if either of the two test outcomes is positive. The joint sensitivity and specificity of the simultaneous testing (denoted by the subscript JE, for ‘Joint Either’) are as follows:
Bayesian modeling for stage two
In the second stage, a gold standard test is used only for positive samples from the simultaneous testing of the first stage, whereas the negative samples are not verified by the gold standard test (see the righthand column in Table 1). Obviously, in this case, the simultaneous testing of the first stage and the gold standard testing of the second stage constitute sequential testing [22]. This sequential testing is used mainly to increase feasibility because applying the gold standard test to the negatives is often costly, unacceptable, and unethical. Given that the gold standard test and the simultaneous testing are independent and that the sensitivity and specificity of the gold standard test are equal to 100%, the joint sensitivity and specificity of the sequential testing (denoted by the subscript JB, for ‘Joint Both’) are as follows:
where Se _{ GS } and Sp _{ GS } are the sensitivity and specificity of the gold standard (denoted by the subscript GS, for ‘Gold Standard’), respectively.
Let vector a = (a _{ 11 } , a _{ 10 } , a _{ 0. }) denote the observed number for each combination in stage two, where a _{ 0. } is the sum of the test negatives from the first stage. Let Pa = (Pa _{ 11 } ,Pa _{ 10 } ,Pa _{ 0. }) represent the respective probabilities from vector a. Consider that elements from the proportion vector of multinomial a are required the probabilities to sum to 1. The multinomial sampling distribution is given by
where the multinomial cell probabilities are given by
Prior distributions
For Bayesian statistics, we wish to use all available information at the design stage but might prefer a more vague, less risky prior at the data analysis stage [23]. To reduce the influence of subjective opinions on the data analysis as much as possible, we placed an informative prior distribution on a minimum number of parameters [24].
In principle, a uniform prior distribution or Beta prior distribution can be used over the set of unknown parameters. We chose the uniform nonBeta prior distribution because a uniform prior distribution is more easily accessible for nonstatisticians and less risky in this study. For example, to determine the value of α and β for a Beta prior distribution, we need to know the corresponding mean and standard deviation or the 2.5 and 97.5 percentiles [25], whereas for a uniform prior distribution, we need only to determine the range of estimated parameters. Especially for the main parameters (Se _{ 1 } , Se _{ 2 } , Sp _{ 1 } , Sp _{ 2 } ,) in this study, lower bound values for the uniform distribution are often estimated securely based on expert opinion, published papers, or even test kit instructions. For example, an epidemiologist can easily be assured that the sensitivity of the enzymelinked immunosorbent assay (ELISA) for HIV antibody screening is higher than 80% [4]; i.e., the lower bound of the uniform distribution is 0.8. However, a reasonable choice for the upper bound of a uniform prior distribution is sometimes difficult to determine because the accuracy of the test kit might differ under various practical conditions. We use the method of twostep hierarchical priors to set the upper bound of the uniform distribution. For an estimated parameter, such as sensitivity (Se _{ i }), we first provide a uniform prior distribution (a _{ Sei } , b _{ Sei }) and then give b _{ Sei } another uniform distribution (b _{ 1Sei } , b _{ 2Sei }). The method of twostep hierarchical priors is a good strategy if it is difficult to determine the prior distribution using a onestep method; even though an incorrect prior is set in step two, the risk of this action resulting in a mistaken result is smaller than in a onestep prior [26].
For sensitivity, we set two steps for the hierarchical uniform priors:
For specificities and prevalence, we might set a uniform prior distribution based on the characteristics of this type of study (see the Discussion section):
For covariances CovDp and CovDn, the feasible range is determined by the sensitivities among the diseased subjects and the specificities among the nondiseased subjects, where 0 ≤ CovDp ≤ min (Se _{ 1 } ,Se _{ 2 })Se _{ 1 } Se _{ 2 } for the diseased subjects and 0 ≤ CovDn ≤ min(Sp _{ 1 } ,Sp _{ 2 })Sp _{ 1 } Sp _{ 2 } for the nondiseased subjects [20]. Because prior information regarding the two covariances is typically unavailable, uniform prior distributions over these ranges can be used for CovDp and CovDn:
Calculation of the posteriors of parameter distributions
The posteriors of parameter distributions were calculated using Markovchain Monte Carlo techniques, in particular the Gibbs sampler in WinBUGS (MRC Biostatistics Unit, Cambridge, UK) [27]. For the analyses presented in this paper, inferences were based on 105,000 iterations after discarding an initial burnin of 5,000 iterations, with convergence assessed by running multiple chains from various starting values [28]. The WinBUGS code used in this paper is available in the Additional files 1 and 2 and can be altered easily for use with different data.
Model evaluation and checking
Model evaluation and checking are highly active areas of Bayesian statistics research. Researchers can use various statistics to determine the plausibility of an assumption of interest in light of the observed data [29]. In this study, we divide the model evaluation and the checking of the assumptions into individual and overall diagnostics. Individual checks are based on the cell in Table 1, and overall diagnostics aim to check the more general assumptions of the model by DIC[21], pD[5], and the Chisquared goodnessoffit test. The technical details were given in Additional file 3. Given the special circumstance that the negatives for both tests are unavailable, none of the methods could be used alone in the preceding model checks because they diagnosed the models from a different perspective. For example, a local Chisquared goodnessoffit test is based only on known cells. For this reason, it is recommended that the above methods be combined for model checking.
Below, we present an ad hoc preliminary, exploratory model that serves to check the criteria, based on empirical information for this type of study. If (1) the values of DIC and pD are reasonable [5]; (2) the pvalue of the local χ ^{2} test is close to 0.5 or far from zero or one [29]; (3) the 95% BCI of the estimated PPV (model key point) includes its actual value; and (4) most frequencies of the cells in Table 1 fall within their corresponding Bayesian credible intervals, we may conclude that the model fits the data well based on the current information and at a specified probability level.
Results and discussion
Results
Screening study for colorectal cancers
The data in Table 2 are based on the data of Castiglione et al. [30]. That study compared rehydrated guaiac testing (Hemoccult) on three consecutive bowel movements with immunochemical testing using reversed passive haemagglutination (RPHA—Hemselect) on the first bowel movement only to detect colorectal cancers or adenomas ≥ 10 mm in 5727 individuals aged 40–59 years. Subjects with a positive hemoccult and/or a positive/borderline hemeselect test were invited to undergo pancolonoscopy. A doublecontrast barium enema was performed when pancolonoscopy was not possible. Further details of the methods used in the study were provided by Castiglione et al. [2], and the data were analyzed by Geoffrey Berry et al [6].
According to the description in of Prior distribution section of this paper, we placed hierarchical uniform priors on Se _{ 1 } and on Se _{ 2 }. These hierarchical priors were elicited based on the opinions of a coauthor (Dr. Hao Yu) and on published papers [2, 30]. The expert opinion and previous information showed that the sensitivity of test 1 (rehydrated hemoccult) was less than the sensitivity of test 2 (RPHAHemeselect) and that the values of their lower bounds values were at least 0.5 and 0.6, respectively. The values of the upper bounds for Se _{ 1 } and Se _{ 2 } were assigned as much diffuse, less risky twostep hierarchical priors, as follows:
For the specificities and prevalence, we set the following uniform prior:
For the covariances CovDp and CovDn, the uniform priors were assigned according to equation (12).
It is clear from Table 3 that the results from the twostage Bayesian model with twostep hierarchical prior on Se _{ 1 } and Se _{ 2 } fit the data well in accordance with the criteria presented in the Model evaluation and checking section. All of the estimated values were close to their true values, and all fell within 95% Bayesian credible intervals over their corresponding true values. The effective number of parameters, a pD of approximately 3.9, indicates that the twostage Bayesian strategy substantially improves the identifiability of the models. The χ ^{2} = 0.002 and P = 0.61 values indicate that the distributions of the replicated and actual data are similar.
Table 3 shows the sensitivity and specificity of joint testing to be approximately 0.81 and 0.94, respectively, with the strategy of jointsimultaneous testing adding approximately 0.13 (Se _{ JE } (0.8103) minus Max(Se _{ 1 }, Se _{ 2 }) (0.6825)) to the sensitivity at the expense of only an approximate loss of 0.04 (Sp _{ JE } (0.9405) minus Max(Sp _{ 1 }, Sp _{ 2 }) (0.9812)) in specificity for the colorectal cancer data. This result conformed to the theory of jointsimultaneous testing. The correlation coefficients under disease and nondisease status, which were approximately 0.4 and 0.1, respectively, suggested moderate and weak dependence under the conditions of disease and nondisease, respectively.
To analyze the influence of the prior distribution of sensitivities (Se _{ i }, (i = 1,2)) on posterior distributions, the noninformative prior, pessimistic and optimistic twostep hierarchical informative priors were assigned, respectively, as below:
For the specificities and prevalence, the following uniform priors were assigned:
Table 4 shows that prior information of the sensitivities produced substantive influence on the posterior distributions. The hierarchical informative priors on the sensitivities must be assigned reasonably in these types of studies.
Berry et al. analyzed the data using maximum likelihood by fitting the four versions of the models: (1) the independence model, (2) independence in the nondiseased group, (3) independence in the diseased group, and (4) the dependence model [6]. Based on the χ ^{2} test statistics, our models provide a better fit to the data than those in the Berry et al. paper. In addition, our method can estimate parameter intervals that were not reported in the Berry et al. paper. In fact, it is difficult for the method of Berry et al. to calculate the confidence intervals of certain complex statistics involved in joint simultaneous testing, such as the joint positive predictive value (PPV _{ JE }), joint sensitivity (Se _{ JE }), and joint specificity (Sp _{ JE }), which are important in evaluating model fitting in this type of study.
Simulation studies
To further illustrate the performance of the twostage Bayesian models, we applied them to a series of simulated data sets with a conditional correlation coefficient of ρ _{ D+ } = 0.5 for the diseased and ρ _{ D } = 0.4 for the nondiseased, as would typically be found in practice. Because our models contain a large number of parameters and because the study design has been associated with a wide range of prevalence values, sample sizes, and test properties, it is impossible to investigate the performance of these models across all possible scenarios. We therefore selected a range related to the actual screening and diagnostic tests, based on the following parameters:

(1)
low prevalence and large sample size, such as cancer screening for community population [30] and HIVantibody screening for blood donors [31].
(1.1) high test accuracy, (π, n, (Se _{1}, Sp _{1}), (Se _{2}, Sp _{2})) = (0.01, 20000, (0.90, 0.95), (0.95, 0.90));
(1.2) low test accuracy, (π, n, (Se _{1}, Sp _{1}), (Se _{2}, Sp _{2})) = (0.01, 20000, (0.60, 0.70), (0.70, 0.60)).

(2)
high prevalence and small sample size, such as the screening of suspicious patients in the hospital or clinic.
(2.1) high test accuracy, (π, n, (Se _{1}, Sp _{1}), (Se _{2}, Sp _{2})) = (0.40, 200, (0.90, 0.95), (0.95, 0.90));
(2.2) low test accuracy, (π, n, (Se _{1}, Sp _{1}), (Se _{2}, Sp _{2})) = (0.40, 200, (0.60, 0.70), (0.70, 0.60))
We also examined the consequences of specifying different prior information, which were divided into the four scenarios (Table 5) given below:

(1)
informative priors for model: the prior density of the estimated model parameters (Se _{ 1 } , Se _{ 2 } , Sp _{ 1 } , Sp _{ 2 } , and π) was centered at their true values.

(2)
informative priors for sensitivities alone: the prior density of the sensitivities alone was centered at their true values, whereas the specificities and prevalence were assigned as noninformative priors.

(3)
twostep informative priors for sensitivities alone: the priors of the sensitivities were assigned according to a twostep hierarchical prior with their upper bounds assigned as much vaguer and less risky priors, as described in the Prior distribution section.

(4)
informative priors for specificities alone: the prior density of the specificities alone was centered at their true values, whereas the sensitivities and prevalence were assigned as noninformative priors.
The results of applying our models to the simulated data sets are given in Tables 6 and 7.
Here, we would like to make the following general observations:

(1)
It is clear that, as long as the priors of the sensitivities were correctly set, the models with informative priors for the model, informative priors for the sensitivities alone, and two step informative priors for the sensitivities alone produced similar results, all of which were centered over the true values. This finding suggests that the prior distributions of the prevalence and specificities could be reasonably assigned as uniform (0, 1) in these types of studies.

(2)
For models with informative priors for the specificities alone, the estimates of the prevalence and sensitivities were biased, and the corresponding credible intervals were very wide, suggesting that the precision of the estimated parameters was low. The negative pDvalues meant that informative priors for the specificities alone are insufficient to estimate the parameters. Thus, the Bayesian models are nonidentifiable from the perspective of probabilistic constraints [5]. The DIC was decreased because of the negative pD.

(3)
The parameter specificities were always estimated with greater precision than the sensitivities when the prevalence was low. Because the specificities relate to the negative subjects, when the prevalence is low, there are more truly negative subjects who provide more experimental information for specificities than for sensitivities; thus, the more precise posteriori estimates of specificities are obtained.
All of the above results indicate that prior information on the sensitivities plays a key role in these types of studies. We consider one of the reasons might be that the prevalence was low, because in absence of enough positive individuals, the prior distributions of sensitivities produced substantive influence on the posterior distribution. The results from both practical case (Tables 3 and 4) and simulation study (Tables 6 and 7) supported the consideration.
As described in the Background section, one of characteristics of two stage screening tests is that the two tests in the first stage are often conditionally dependent on disease status and/or nondisease status [5]. To further study how the prevalence and test accuracy estimates vary with the model assumptions (dependent or independent) and to evaluate the impact of the misspecification of different models on the estimation of the prevalence and test accuracy, we performed four sets of simulations under the assumptions of the independent model for independent data, the independent model for dependent data, the dependent model for independent data, and the dependent model for dependent data. For ease of presentation and interpretation, we considered moderate parameter values for the prevalence, sample size, and test accuracy, as shown below:
Table 8 shows that the independent models might generate incorrect statistical inferences for the conditional dependent data and larger DIC (323.08) values, indicating that the model was not appropriate for the data. For conditional independent data, the independent models and dependent models had similar results closer to their true values. Thus, we suggest that, when the researcher is uncertain as to whether the data are dependent or independent, the less risky dependent model should be used despite the slight loss of precision among the estimated parameters (widened credible interval).
Discussion
This article presents additional methods for the estimation problem for screening studies using two screening tests with the disease statuses being verified for test positives only. The core of the methods is that Bayesian modeling is divided into two stages characterized by a joint testing strategy based on the inherent attributes of this type of screening. The first stage is a joint simultaneous testing consisting of two dependent binary tests. The second stage is a special joint sequential testing that consists of simultaneous testing in the first stage and a gold standard test in the second.
Another method presented in this study involves assigning a diffuse and less risky twostep hierarchical prior structure for sensitivities. Because reasonably informative prior information is necessary for resolving the nonidentifiability problem in this type of study, we emphasized the setting of the sensitivity prior distribution based on the following three points. First, the priors of the sensitivity are more important than the priors of the specificity or prevalence when the prevalence is lower than 50%, as is often found in this type of study. Second, we placed the prior information of the sensitivity on the least possible quantity, i.e., diffuse prior or less risky prior information, such as twostep hierarchical prior structure, to reduce the influence of subjective opinion as much as possible on the results. Third, we set the prior structure to be as easy as possible for use by nonstatisticians so that uniform prior distributions could be used in this study.
We found that when the priors of the sensitivities were assigned as appropriate priors, such as twostep hierarchical uniform priors, in the absence of priors for the specificities and prevalence (e.g., assigning a uniform distribution on (0,1)), the parameters involved in this study could still be estimated well. However, changing the prior structure of the sensitivities has a substantive effect on the estimators (Tables 4, 6 and 7). The reasons for these results are presented in the Simulation studies section. Here, we would like to state that this finding has important practical significance because the sensitivities (related to the diseased population) can often, but not always, be obtained more easily than the specificities (related to the healthy population) in medical practice. This characteristic is partially based on the increased feasibility of applying the gold standard test to the diseased population compared with the healthy population, owing to the frequent invasiveness, expense, and sophistication of many gold standard tests, particularly ones with further ethical implications. Therefore, we suggest that the prior structure of the sensitivities should be thoroughly investigated for scenarios in lowprevalence populations.
To our knowledge, this study is the first to estimate test accuracy and disease prevalence for two dependent dichotomous screening tests with unverified negative individuals on both tests using twostage Bayesian modeling. Recently, Bohning and Patilea, Chu et al., and Li et al. proposed twostage methods to address the nonidentifiability problem [7, 8, 11]. Their twostage estimation methods differ from ours. In the first stage, Bohning and Patilea [11], Chu et al., and Li et al. [7, 8] estimated the parameters of known cells; in the second stage, they estimated the parameters of unknown cells by making assumptions regarding the cell probabilities of various homogeneous association models. They circumvented the nonidentifiability problem using a capture–recapture approach or constrained maximum likelihood estimation. Because all of the models are saturated, homogeneous association assumptions are not testable. Therefore, the dependence structure modeling considered can only be viewed as a sensitivity analysis [7]. They did suggest that Bayesian methods that incorporate prior information could be a reasonable alternative to obtain improved estimates, although further research is needed [8]. Their suggestion was an important factor in initiating our study.
Our twostage approach has certain advantages over the previous methods. First, our method can be understood more easily and used with greater convenience by nonstatisticians because the procedure of constructing the model fully agrees with the practical screening procedure. Additionally, the modeling often only involves the prior distribution for the sensitivities, because the prior distributions of the specificities and prevalence are assigned a noninformative uniform distribution of (0, 1), respectively. Moreover, the method of model checking is simple, intuitive, and convenient, as pD and DIC can be obtained directly using WinBUGS. Second, this method improves nonidentifiability more reasonably, as probabilistic constraints not deterministic constraints (simplifying the model) are used. Third, the method could obtain more information to evaluate the screening strategy than the capture–recapture approach or constrained maximum likelihood estimation [7, 8, 11]. For example, the method could conveniently calculate the joint sensitivities and joint specificities and make corresponding statistical inferences, such as 95% credible intervals (Table 3). In the framework of frequency statistics, such as capture–recapture approaches or constrained maximum likelihood estimations, it is often difficult to obtain statistical inferences about the joint sensitivities and joint specificities in these types of studies.
Conclusions
We developed a twostage Bayesian method for two dependent dichotomous screening tests with unverified individuals who are negative on both tests. We also addressed the ad hoc model evaluation and checking procedures based on empirical information for these types of studies. Although the practical example and simulation studies considering various practical situations showed that these models fit well, three points should be clarified. (1) Due to the lack of full data, it is impossible to check the model completely without further information. (2) Although our Bayesian twostage modeling improves nonidentifiability based on probabilistic constraints, the informative prior information on the sensitivities must be assigned reasonably, even though the prior could be diffuse, as in the twostep hierarchical uniform prior. (3) Using our real example and simulated scenarios related to reallife situations, the priors of the sensitivities are more important than the priors of the specificities. However, when the prevalence is large, especially if it is greater than 50% (rare in medical screening tests), the same attention should be paid to the priors of the specificities.
Abbreviations
 BCI:

Bayesian credible interval
 DIC:

Deviance Information Criterion
 pD:

Effective number of parameters
 PPV:

Positive predictive value
 SI:

Simultaneous testing.
References
Macaskill P, Walter SD, Irwig L, Franco EL: Assessing the gain in diagnostic performance when combining two diagnostic tests. Stat Med. 2002, 21 (17): 25272546. 10.1002/sim.1227.
Castiglione G, Zappa M, Grazzini G, Mazzotta A, Biagini M, Salvadori P, Ciatto S: Immunochemical vs guaiac faecal occult blood tests in a populationbased screening programme for colorectal cancer. Brit J Cancer. 1996, 74 (1): 141144. 10.1038/bjc.1996.329.
Branscum AJ, Gardner IA, Johnson WO: Bayesian modeling of animal and herdlevel prevalences. Prev Vet Med. 2004, 66 (1–4): 101112.
Liu P, Xiao S, Shi ZX, Bi XX, Yang HT, Jin H: Bayesian evaluation of the human immunodeficiency virus antibody screening strategy of duplicate enzymelinked immunosorbent assay in Xuzhou Blood Center, China. Transfusion. 2011, 51 (4): 793798. 10.1111/j.15372995.2010.02890.x.
Berkvens D, Speybroeck N, Praet N, Adel A, Lesaffre E: Estimating disease prevalence in a Bayesian framework using probabilistic constraints. Epidemiology. 2006, 17 (2): 145153. 10.1097/01.ede.0000198422.64801.8d.
Berry G, Smith CL, Macaskill P, Irwig L: Analytic methods for comparing two dichotomous screening or diagnostic tests applied to two populations of differing disease prevalence when individuals negative on both tests are unverified. Stat Med. 2002, 21 (6): 853862. 10.1002/sim.1066.
Chu H, Zhou Y, Cole SR, Ibrahim JG: On the estimation of disease prevalence by latent class models for screening studies using two screening tests with categorical disease status verified in test positives only. Stat Med. 2010, 29 (11): 12061218.
Li F, Chu H, Nie L: A twostage estimation for screening studies using two diagnostic tests with binary disease status verified in test positives only. Stat Met Med Res. 2011, Epub ahead of print
Walter SD: Estimation of test sensitivity and specificity when disease confirmation is limited to positive results. Epidemiology. 1999, 10 (1): 6772. 10.1097/0000164819990100000012.
van der Merwe L, Maritz JS: Estimating the conditional falsepositive rate for semilatent data. Epidemiology. 2002, 13 (4): 424430. 10.1097/0000164820020700000010.
Bohning D, Patilea V: A capturerecapture approach for screening using two diagnostic tests with availability of disease status for the test positives only. J Am Stat Assoc. 2008, 103 (481): 212221. 10.1198/016214507000000383.
Branscum AJ, Gardner IA, Johnson WO: Estimation of diagnostic test sensitivity and specificity through Bayesian modeling. Prev Vet Med. 2005, 68 (2–4): 145163.
Cheng D, Stamey JD, Branscum AJ: Bayesian approach to average power calculations for binary regression models with misclassified outcomes. Stat Med. 2009, 28 (5): 848863. 10.1002/sim.3505.
Hanson TE, Johnson WO, Gastwirth JL: Bayesian inference for prevalence and diagnostic test accuracy based on dualpooled screening. Biostatistics. 2006, 7 (1): 4157.
Stamey JD, Young DM, Seaman JW: A Bayesian approach to adjust for diagnostic misclassification between two mortality causes in Poisson regression. Stat Med. 2008, 27 (13): 24402452. 10.1002/sim.3134.
Suess EA, Gardner IA, Johnson WO: Hierarchical Bayesian model for prevalence inferences and determination of a country’s status for an animal pathogen. Prev Vet Med. 2002, 55 (3): 155171. 10.1016/S01675877(02)000922.
Xia M, Gustafson P: A Bayesian method for estimating prevalence in the presence of a hidden subpopulation. Stat Med. 2012, 31 (21): 23862398. 10.1002/sim.5374.
Black MA, Craig BA: Estimating disease prevalence in the absence of a gold standard. Stat Med. 2002, 21 (18): 26532669. 10.1002/sim.1178.
Erkanli A, Soyer R, Stangl D: Bayesian inference in twophase prevalence studies. Stat Med. 1997, 16 (10): 11211133. 10.1002/(SICI)10970258(19970530)16:10<1121::AIDSIM526>3.0.CO;2L.
Dendukuri N, Joseph L: Bayesian approaches to modeling the conditional dependence between multiple diagnostic tests. Biometrics. 2001, 57 (1): 158167. 10.1111/j.0006341X.2001.00158.x.
Spiegelhalter DJ, Best NG, Carlin BP: Bayesian measures of model complexity and fit (with discussion). J R Stat Soc Ser B. 2002, 64: 583640. 10.1111/14679868.00353.
Su CL, Gardner IA, Johnson WO: Diagnostic test accuracy and prevalence inferences based on joint and sequential testing with finite population sampling. Stat Med. 2004, 23 (14): 22372255. 10.1002/sim.1809.
Berry SM, Carlin BP, Lee JJ, Muller P: Bayesian adaptive methods for clinical trials. CRC press; 2010: 72
Lu Y, Dendukuri N, Schiller I, Joseph L: A Bayesian approach to simultaneously adjusting for verification and reference standard bias in diagnostic test studies. Stat Med. 2010, 29 (24): 25322543. 10.1002/sim.4018.
Joseph L, Gyorkos TW, Coupal L: Bayesian estimation of disease prevalence and the parameters of diagnostic tests in the absence of a gold standard. Am J Epidemiol. 1995, 141 (3): 263272.
Mao SS: Bayesian Statistics. 2005, Beijing: China Statistics Press, 105110.
Spiegelhalter DJ, Thomas A, Best NG, Lunn D: WinBUGS Version 1.4.3 user manual. [http://www.mrcbsu.cam.ac.uk/bugs]
Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences. Stat Sci. 1992, 7: 457511. 10.1214/ss/1177011136.
Ntzoufras I: Bayesian Modeling Using WinBUGS. 2009, John Wiley & Sons, Inc: Hoboken, New Jersey
Castiglione G, Zappa M, Grazzini G, Sani C, Mazzotta A, Mantellini P, Ciatto S: Cost analysis in a population based screening programme for colorectal cancer: comparison of immunochemical and guaiac faecal occult blood testing. J Med Screen. 1997, 4 (3): 142146.
Liu P, Yang HT, Qiang LY, Xiao S, Shi ZX: Estimation of the sensitivity and specificity of assays for screening antibodies to HIV: a comparison between the frequentist and Bayesian approaches. J Virol Met. 2012, 186: 8993. 10.1016/j.jviromet.2012.08.018.
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712288/14/110/prepub
Acknowledgments
This work was supported by the National Natural Science Foundation of China [grant number 81273184].
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
LJ and CF conceived and conducted the study, and drafted the manuscript. YH participated in the design, data generation and interpretation of the analyses. ZP and LL participated in the simulation and interpretation of the analyses. All the authors read and approved the final manuscript.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Liu, J., Chen, F., Yu, H. et al. A twostage Bayesian method for estimating accuracy and disease prevalence for two dependent dichotomous screening tests when the status of individuals who are negative on both tests is unverified. BMC Med Res Methodol 14, 110 (2014). https://doi.org/10.1186/1471228814110
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1471228814110
Keywords
 Bayesian modeling
 Screening test
 Test accuracy
 Disease prevalence
 Model checking