Skip to main content

Sample size re-estimation in paired comparative diagnostic accuracy studies with a binary response

Abstract

Background

The sample size required to power a study to a nominal level in a paired comparative diagnostic accuracy study, i.e. studies in which the diagnostic accuracy of two testing procedures is compared relative to a gold standard, depends on the conditional dependence between the two tests - the lower the dependence the greater the sample size required. A priori, we usually do not know the dependence between the two tests and thus cannot determine the exact sample size required. One option is to use the implied sample size for the maximal negative dependence, giving the largest possible sample size. However, this is potentially wasteful of resources and unnecessarily burdensome on study participants as the study is likely to be overpowered. A more accurate estimate of the sample size can be determined at a planned interim analysis point where the sample size is re-estimated.

Methods

This paper discusses a sample size estimation and re-estimation method based on the maximum likelihood estimates, under an implied multinomial model, of the observed values of conditional dependence between the two tests and, if required, prevalence, at a planned interim. The method is illustrated by comparing the accuracy of two procedures for the detection of pancreatic cancer, one procedure using the standard battery of tests, and the other using the standard battery with the addition of a PET/CT scan all relative to the gold standard of a cell biopsy. Simulation of the proposed method illustrates its robustness under various conditions.

Results

The results show that the type I error rate of the overall experiment is stable using our suggested method and that the type II error rate is close to or above nominal. Furthermore, the instances in which the type II error rate is above nominal are in the situations where the lowest sample size is required, meaning a lower impact on the actual number of participants recruited.

Conclusion

We recommend multinomial model maximum likelihood estimation of the conditional dependence between paired diagnostic accuracy tests at an interim to reduce the number of participants required to power the study to at least the nominal level.

Trial registration

ISRCTN ISRCTN73852054. Registered 9th of January 2015. Retrospectively registered.

Peer Review reports

Background

An assessment of diagnostic accuracy is crucial in the development of medical testing procedures [1]. Comparing the accuracy of these procedures in terms of their sensitivities and specificities [2, 3] relative to a gold standard, is essential to ensuring that the most appropriate tests are deployed in the clinical setting [4, 5]. The focus of this paper is sample size re-estimation in the comparison of two candidate tests to a gold standard. However, diagnostic accuracy studies do not necessarily involve comparisons; many such studies report the accuracy of a single test.

At the outset of a study, a sample size is calculated based on assumptions made about the expected changes in sensitivity and specificity and, in a prospective design, the likely prevalence of the condition to be tested for in the sample. However, the initial assumptions about parameters in the study, especially the conditional dependence between the two tests, may be revealed to be inaccurate, resulting in a potentially over- or under-powered study. A planned interim analysis can allow the study’s sample size to be updated based on the data already collected. This involves utilising the information observed at the interim stage to refine the sample size estimate. A resulting increase in sample size allows the time, cost and patient discomfort already invested in the study to yield valid results while a decrease in sample size means that less time and cost will be expended overall and patients will not needlessly undergo unnecessary testing [6].

There are well-established methodologies for interim sample size re-estimation in treatment studies for continuous and normally distributed response variables [7,8,9,10,11], some of which provide mechanisms to maintain blinding in the study [8,9,10]. Methods also exist for the re-estimation with binary response variables [12, 13], and mechanisms to maintain blinding have been proposed in this more complex situation where the variance and mean parameters are not separable [14]. Proschan [15] gives an overview of sample size re-estimation procedures based on a nuisance parameter. Specifically, procedures for determining the difference of means between two samples with a common, unknown, variance and difference in proportions between two groups, with an unknown overall proportion, are considered. In the case of normally distributed data, the independence of the sample variance and sample mean ensures that the validity of estimates is unaffected by the interim sample size re-estimation and this is shown to hold asymptotically in the binary case. However, Proschan does not consider the case of paired data which is the focus of the current paper. Furthermore, the implications of sample size re-estimation in the context of comparative diagnostics studies, inherently different from those in treatment (randomised controlled) studies [16], have not been fully explored in the statistical literature.

A number of salient differences in interim analysis between studies comparing diagnostic tests and those comparing treatments are highlighted in Gerke et al. [5] and Gerke et al. [16]. Firstly, in paired diagnostics accuracy studies, full blinding is often not possible, specifically, certain types of test may not be able to be blinded from the patient, the person administering the test, the person interpreting the test, or the person measuring the outcome. However, as long as the results of the two-tests which are being compared are temporarily blinded from the person measuring the outcome, this is not a major threat to a study’s validity [17]. In fact, it has the advantage that the patients can benefit from their clinicians knowing the results of both diagnostic tests after testing has taken place. Secondly, in diagnostic accuracy studies, early cessation of the study due to futility is not as easy to establish as in treatment studies. The reasons for this are 1) the fact that treatment studies often test a single outcome while diagnostic studies test two outcomes, sensitivity and specificity, and futility must be established for both simultaneously, and 2) patient outcomes may only be seen further downstream from the test results [18]. Thirdly, the sample size required for a hypothesis test in diagnostic studies, powered to a given level, is closely related to the conditional dependence between the two testing procedures which has been shown to present problems in a number of contexts [5, 19,20,21,22,23,24]. More specifically, the lower the conditional dependence between the tests, the greater the sample size will be, with the largest sample size being implied by the maximum negative dependence, given the specified alternative hypotheses. This level of conditional dependence between the tests is one of the primary factors driving the required sample size estimate and it is often difficult to estimate a priori. Gerke et al. [5] assert that for comparative diagnostic studies, as long as an interim sample size re-estimation is planned it bears no threat to the validity of the study. However, Gerke et al. [5] do not provide justification for this assertion and, furthermore, their assertion does not take the inherent uncertainty of the interim data into account. This study aims to present a method and give practical guidelines for its application, for the initial estimation and interim re-estimation of sample size in a paired diagnostic study which will allow utilisation of information on the conditional dependence between tests at the interim to potentially reduce the required sample size while maintaining the approximate nominal statistical power of the experiment as a whole. While we present a method of estimating the size of the conditional dependence to reduce sample size, it should also be noted that there is a body of literature dealing with the problems caused by conditional dependence in other areas [25,26,27].

The remainder of the article is organised as follows. The methods section outlines sample size estimation methods for paired diagnostic test studies, introduces a motivating example application, and then goes on to propose a new method for re-estimation based on a multinomial likelihood. The results section first provides extensive simulations of the method under various real world conditions and then moves back to apply the sample size re-estimation method proposed in this paper to the motivating example. The article then continues with a brief discussion of the place of this study in the literature and the optimal interim sample size to choose. Finally, the conclusion, summarises and restates the major outcomes of this study.

Methods

A representation of data from a paired comparative diagnostic accuracy study is given in Table 1. The subjects are initially divided according to whether they are discovered, via the gold standard test, to be diseased or non-diseased. They are then further subdivided as to whether they test positive or negative on tests A and B. For example, the cell n A represents subjects that were found to have the disease via the gold standard test and also tested positive on both test A and B, while cell n F denotes subjects who tested negative on the gold standard and test B but positive on test A.

Table 1 Paired study design

A possible initial sample size calculation, using a normal approximation of the logarithm of the ratio of sensitivities and specificities, and assuming a comparison between a new test, test A, and an existing test, test B, follows from Alonzo et al. [21] and a full derivation can be found therein. The experiment, as a whole tests jointly both sensitivity and specificity improvement to pre-specified levels, the sample size is calculated for each and the largest sample size is chosen to power the study. Note that this paper concentrates on the situation in which superiority is tested for both sensitivity and specificity. However, the method elaborated below should be extendable to situations where we are interested in testing non-inferiority in either or both of sensitivity and specificity. For details on the construction of the confidence intervals and hypothesis tests in these situations see Alonzo et al. [21]. In the case of the estimation of a sample size for superiority, the initial sample size calculation for sensitivity is given by:

$$ {n}_{p1}={\left(\frac{Z^{\left(1-\beta \right)}+{Z}^{\left(1-\alpha /2\right)}}{\mathit{\log}{\gamma}_1}\right)}^2\left(\frac{\left({\gamma}_1+1\right){TPR}_B-2 TPPR}{\gamma_1{TPR}_B^2}\right)/\pi $$
(1)

where, α is the type I error rate of the study and β is the power of the study. The main quantity of interest, γ 1, is the ratio of true positive rates=TPR A /TPR B , TPR B is the true positive rate (sensitivity) on test B, i.e. TPR B  = (n A  + n C ) / (n A  + n B  + n C  +  n D ), TPR A is the true positive rate (sensitivity) on test A, i.e. TPR A = (n A  + n B ) / (n A  + n B  + n C  +  n D ), TPPR is the proportion of diseased patients who test positive on both tests, i.e. TPPR = n A /( n A  + n B  + n C  +  n D ) and π is the prevalence of disease. The null hypothesis is that γ 1 = 1, the alternative hypothesis is that γ 1≠1.

For testing superiority of specificity we are interested in the true negative rates so the formula is instead:

$$ {n}_{n1}={\left(\frac{Z^{\left(1-\beta \right)}+{Z}^{\left(1-\alpha /2\right)}}{\mathit{\log}{\gamma}_2}\right)}^2\left(\frac{\left({\gamma}_2+1\right){TNR}_B-2 TNNR}{\gamma_2{TNR}_B^2}\right)/\left(1-\pi \right) $$
(2)

where, γ 2, the main quantity of interest is the ratio of true negative rates =TNR A /TNR B , TNR A is the true negative rate (specificity) on test A = (n G  + n H ) / (n E  + n F  + n G  +  n H ), TNR B is the true negative rate (specificity) on test B = (n F  + n H ) / (n E  + n F  + n G  +  n H ), and TNNR is the proportion of non-diseased patients who test negative on both tests = n H /( n E  + n F  + n G  +  n H ).

It is interesting to note that, following the notation of Vacek [25] and considering the population 2 × 2 table (in Table 1), the conditional dependence of the two tests can be denoted by e b and e a. , the conditional covariance when the gold standard disease status is positive or negative, respectively [25]. Therefore, the probability of both tests being positive can be expressed as TPPR = TPR A  ∙ TPR B  + e b and the probability of both tests being negative TNNR = (1 − TNR A ) ∙ (1 − TNR B ) + e a . When e a and e b  = 0 the tests are conditionally independent, when e a and/or e b  ≠ 0 the response on one test changes the probability of that response on the other test. For example, when e b  > 0 an individual who responds positively on test A is more likely to respond positively on test B.

For initial estimates of TPPR and TNNR, from Alonzo et al. [21] we can use the fact that TPPR ≥ (1 + γ 1)TPR B  − 1 and TNNR ≥ (1 + γ 2)TNR B  − 1 to estimate the lower bounds of the possible values of TPPR and TNNR, under the specified hypotheses. The required sample size is largest when TPPR = (1 + γ 1)TPR B  − 1 and TNNR = (1 + γ 2)TNR B  − 1, thus, these estimates represent the “worst case scenarios” of maximal negative conditional dependence between the tests, conditional on the fixed values of TPR A and TPR B . The sample size implied by using these levels of TPPR and TNNR would very likely overpower the study, i.e. more participants will be recruited than is strictly necessary to achieve the power specified by β. The required sample size is smallest when the conditional dependence between tests A and B are maximal, conditional on the fixed values of TPR A and TPR B , i.e. when TPPR = TPR B and TNNR = TNR B . The implied sample size in this case would likely underpower the study, i.e. too few participants recruited to reach the power specified by β. The sample size in this “best case scenario” can be substantially lower than that in the worst case scenario. Conservatively, it might be thought a good idea to always use the “worst case scenario” implied sample size estimate which will always power the study sufficiently. However, in cases where the recruitment and testing of participants comes at a premium, both financially and in terms of discomfort to the patients, it might be preferable to apply a more nuanced strategy. Furthermore, the sample size implied by the “worst case scenario” implies the highly unlikely condition of a maximal negative conditional dependence between two tests, which are performed on the same patients to detect the same disease. The implied sample size based on this condition is not recommended [28]. One possibility, to enable a more accurate evaluation of the conditional dependence between the two tests, and thus the required sample size, is to perform a planned interim sample size re-estimation using this information to refine the sample size estimate.

At a planned interim, where a proportion of the overall sample size has been collected, we would have some information about the true values of TPPR, TNNR, π, TPR B and TNR B , however, these values would only come from a limited sample size. The crucial parameters to use in re-estimation are those related to the conditional dependence between the tests, i.e., TPPR and TNNR, as these values are difficult to estimate and, for these parameters, it is unlikely that research exists which can provide an approximate value. Conversely, the values of, TPR B and TNR B , the sensitivities and specificities of an established test, may have known values in the literature and these should preferably be used over those from the relatively small interim sample. For the value of π,the prevalence, a judgement must be made as to whether the researcher feels that any pre-existing estimate of prevalence would be a more accurate reflection of the true prevalence in the specific study population than any interim estimate. In the example given below, we use values for TPPR, TNNR and π at the interim in the sample size calculation.

Naively, it might appear that interim sample size re-estimation would entail a straightforward replication of eqs. (1) and (2) with π, and in the case of (1), TPPR or in the case of (2), TNNR, replaced with the estimates at the interim point. However, this approach does not effectively take into account the inherent uncertainty in the interim parameter estimates of TPPR, TNNR and π, nor the fact that only a specific range of values for TPPR and TNNR are actually possible under the alternative hypothesis. An approach which does take these factors into account is re-estimation of the sample size based on maximum likelihood estimation, at the interim, of the parameters in question under a multinomial model. This model is constrained by the hypothesised values of TPR A  , TPR B , TNR A , and TNR B , i.e. the marginals in Table 1.

Application

The numerical example we use involves an interim sample size recalculation of a study comparing the incremental benefits to sensitivity and specificity of augmenting current methods for diagnosing pancreatic cancer with Positron Emission Tomography (PET) and computed tomography (CT) technologies. The alternative hypotheses were that sensitivity would rise from 81% to 90%, and specificity would rise from 66% to 80%, additionally, the expected prevalence of pancreatic cancer from the literature was 47%.

To calculate the sample size for sensitivity equation 1 was used, taking \( \alpha =0.05,\kern0.5em \beta =0.2,\kern0.5em {\widehat{\gamma}}_1=\frac{0.9}{0.81}, \) \( \widehat{TPR_B}=0.81, \) \( \widehat{TPPR}=0.71 \), and \( \widehat{\pi}=0.47 \) gives a sample size of 598. To calculate the sample size for specificity equation 2 was used taking \( \alpha =0.05,\kern0.5em \beta =0.2,{\ \widehat{\gamma}}_2=\frac{0.8}{0.66}, \) \( \widehat{\ {TNR}_B}=0.66, \) \( \widehat{TNNR}=0.46 \), and \( \widehat{\pi}=0.47 \) gives a sample size of 409. The minimum sample sizes for sensitivity and specificity, given \( \widehat{TPPR}=0.81 \) and \( \widehat{TNNR}=0.66 \), are 186 and 106, respectively. Given the disparity between the minimum and maximum sample size estimates it was decided to re-assess the sample size at a planned interim.

Table 2 gives the results after data from 187 participants had been collected. The observed values at the interim are: \( \widehat{TPPR}=0.80 \), \( \widehat{TNNR}=0.66 \) and \( \widehat{\pi}=0.44 \). Taking a naive approach and plugging these values directly into equations 1 and 2 the implied sample sizes for sensitivity become 242 and for specificity 100, giving a total sample size for the study of 242 (or 342 and 145, respectively, had we also used the interim values of TPR B and TNR B ). However, this method does not take into account the fact that \( \widehat{TPPR} \) and\( \widehat{\ TNNR} \) are random variables and we are actually interested in the true value of the probability of TPPR and TNNR under the specified alternative hypothesis. In fact, had the observed value for TPPR been equal to 0.86, the sample size given via the naive method would have been −22, given the fact that \( \widehat{TPPR} \) would have been larger than both TPR A and TPR B . Clearly, the naive method, which uses the random value of a single cell, is inappropriate and a method that uses information about the value of TPPR from all of the observed cells and the specified marginals is required.

Table 2 Interim PET diagnostic study results

Sample size re-estimation via maximum likelihood estimation of TPPR

For illustration purposes, we will discuss the re-estimation of the sample size for sensitivity, the estimation procedure for specificity is analogous. Taking TPR A as the test with the highest expected diagnostic utility, i.e. the “new” test whose performance we are comparing to the “standard”, the probabilities corresponding to the cells in Table 1, given the situation of the maximally negative conditional dependence between the tests are: p 1 = TPR B  − (1 − TPR A ), p 2 = 1 − TPR B , p 3 = 1 − TPR A , p 4 = 0. The probabilities of the cells when the conditional dependence between TPR A and TPR B is at its maximally positive are given by: p 1 = TPR B , p 2 = TPR A  − TPR B , p 3 = 0, p 4 = 1 − TPR A . We could alternatively specify these cell probabilities according to the covariance between the two tests. Specifically, Vacek [25] gives the maximum value of the covariance as TPR B  (1 − TPR A ) and the minimum value as −(1 − TPR A )(1 − TPR B ). Thus, the maximum and minimum values for the cells can be ascertained by finding the product of the marginal probabilities associated with a cell and adding the minimum or maximum value of covariance, for cells p 1 and p 4, or subtracting the values of covariance for cells p 2 and p 3. For example, the minimum value for p 1 = TPR A  ∙ TPR B  − (1 − TPR A )(1 − TPR B ). Between the minimum and maximum values lies every permissible joint configuration. Let these possible joint configurations be expressed as vector, p, with p 1 = TPPR,where \( {\sum}_{i=1}^4{\mathbf{p}}_i=1,{\ p}_1+{p}_2 = {TPR}_A \) and p 1 + p 3 = TPR B .

When the conditional dependence is maximally positive the sample size required is the smallest, when it is maximally negative the sample size required is at its largest. At the beginning of the experiment we do not know which of these possible levels of conditional dependence our data were generated under and thus we use the, usually overly conservative, largest possible sample size estimate.

However, at the interim we can use our observed data to infer a likelihood of that data having been generated under each of the permissible joint configurations of cell probabilities given the implied range of probabilities under a multinomial model. A simple method of extracting an estimate of TPPR is to maximise the likelihood function of the interim data given the values of p implied by the marginal probabilities:

$$ \mathcal{L}\left(\boldsymbol{p}| x\right)=\kern0.5em \prod_{i=1}^4{\boldsymbol{p}}_i^{x_i} $$
(3)

where p is the vector of joint probabilities defined above and x are the observed cell frequencies. The constraints imposed on the above multinomial likelihood make the parameter space one dimensional, thus, substituting the constraints in order to express the likelihood in terms of p 1, gives:

$$ \mathcal{L}\left({p}_1| x\right)={p}_1^{x_1}{\left(\ {TPR}_A - {p}_1\right)}^{x_2}{\left(\ {TPR}_B - {p}_1\right)}^{x_3}{\left(1 - {TPR}_A - {TPR}_B + {p}_1\right)}^{x_4} $$
(4)
$$ {\ p}_1\in \left[\ {TPR}_B-\left(1 - {TPR}_A\right),{\ TPR}_B\right] $$

Code to estimate this in R, via optimisation of the negative log-likelihood, is in the Appendix. In effect, this method bounds the value for the conditional dependence between the minimum and maximum values under the specified marginals and then uses information from the frequency values of the four cells of the table to infer the most probable value of p 1. We can use this estimate of\( {\ \widehat{p}}_1 \) as our value of \( \widehat{TPPR} \) and use the observed value of the prevalence (if required) as our measure of \( \widehat{\pi} \) in equation 1 to re-estimate the sample size at the interim.

Results

Simulation studies

In order to verify the integrity of the method for sample size re-estimation described and applied above a series of simulation studies were carried out. The objectives of these studies were to assess the implications of re-estimating a sample size based on data already collected on the type I and II error rates under various permutations of parameters. The type II error rate should be as close to nominal as possible (i.e. 0.8 in the example above), and the type I error rate should be minimally affected by the re-estimation.

It should be noted that the statistical power provided by the sample size implied by the Alonzo et al. [21] method (when no re-estimation is undertaken) is related to the level of conditional dependence between the tests, Fig. 1 illustrates this relationship. In total 100,000 replications were generated under the specified true alternative hypothesis (i.e. γ 1 = 0.9/0.81 = 1.11), for the example situation above, at various levels of conditional dependence between the two tests. The number of replications 100,000, is more than required, however as the computing time to calculate these was trivial, there was little cost in simulating to this level of accuracy. This number of simulations was used throughout this paper. In all cases in Fig. 1 the simulated power was higher than nominal but where the conditional dependence was highest the power was greatly over specified. As the conditional dependence tends towards becoming maximally positive, i.e. as TPPR tends towards its maximal value, the cell n C tends towards 0. This means that the asymptotic assumptions underlying formulae 1 and 2 and those underlying the significance test no longer hold. However, this should not be of too great a concern, with regards to balancing the minimisation of the required sample size estimate with the statistical power of the experiment, as the instances where the power is over specified are when the sample size is lowest. Additional conservatism at positive levels of conditional dependence has a significantly lesser impact on the overall sample size than it would have at the end of the continuum where the conditional dependence is negative. Whatever the case may be, it should be noted that the results of re-estimation will follow a similar pattern.

Fig. 1
figure 1

Simulated power of sample size specified by the true TPPR in equation 1 when TPRA=0.9, TPRB=0.81 and π=0.45

In the first set of simulations, which aim to assess the stability of the type II error rate, data are generated under the conditions TPR A  = 0.9, TPR B  = 0.81, π = 0.45, while the sample sizes at the interim are varied between 50 and 200 and the values for TPPR are varied between 0.71 and 0.81. The null hypothesis is: TPR A /TPR B  = 1, and our data were simulated under the alternative hypothesis TPR A  = 0.9 and TPR A  = 0.81, with varying levels of conditional dependence within the implied limits. Figure 2 shows how the power of the experiment overall (i.e. using the data from both before and after sample size re-estimation) varies as a function of the interim sample size and the true value of TPPR. As expected the values follow the same pattern as that in Fig. 1. The minimum of the nominal power, or very close to it, was achieved at all levels of conditional dependence and at all interim sample sizes.

Fig. 2
figure 2

Simulated power of re-estimation method across various interim sample sizes and levels of true TPPR when TPRA = 0.9, TPRB = 0.81 and π = 0.45

Table 3 provides information about the mean sample size, bias, coverage and Root Mean Squared Error (RMSE) (from the value specified by equation 1 using the true value of TPPR for the simulated data) under the combinations of conditional dependences and interim sample size. The sample sizes implied by Equation 1 for maximal and minimal levels of conditional dependence are 194 and 625, respectively. The interim sample sizes of 50, 100, 150 and 200 were chosen to illustrate the effects of choosing various interim sample sizes that were smaller than the total sample size of 242 calculated by the Alonzo method described above for our application.

Table 3 Mean sample size (S.D.), bias, coverage and RMSE of simulated sample sizes with varying interim sample size estimates and true levels of TPPR when TPRA= 0.9, TPRB = 0.81 and Prevalence = 0.45. (N = interim sample size)

An increasing interim sample size does not have that great an impact on the average estimated sample size. However, it does have a large impact on the RMSE. Thus, choosing a larger interim sample size at which to re-estimate will ensure a more accurate sample size re-estimate in individual cases, meaning that the experiment will be more likely to be powered to the appropriate level while recruiting as few participants as possible. Of course, if the interim sample size is chosen to be too large then there is a risk of having already recruited too many participants at the interim. Therefore, some sensible trade-off is required. The bias and coverage seem to be at acceptable levels although the coverage does dip when the conditional dependence between the tests is high.

A second set of simulations was run to assess the performance of the method under the null hypothesis where \( {\gamma}_1=\frac{{\ TPR}_A}{{\ TPR}_B}=1 \). Table 4 shows the cell probabilities for these simulations. Rather than report across the entire range only the minimum, 50% (i.e., median) and maximum levels of TPPR are reported.

Table 4 Simulation settings to estimate Type I error

Table 5 shows the type I error rate, mean sample size, bias, coverage and RMSE of simulated sample sizes under various simulation settings. At all levels of conditional dependence and at all interim sample sizes the type I error rate is close to the specified levels. Again, the inference to be made from the RMSE value is that a larger sample size provides a more accurate estimate of the full sample size required, reducing the extent to which an experiment will be over or underpowered in individual cases. The bias and coverage also appear to be at acceptable levels.

Table 5 Type I error rate, Mean sample size (S.D.), bias, coverage and RMSE of simulated sample sizes under various simulation settings

Table 6 gives the results of a range of simulations undertaken at various values of TPR A TPR B and π in both true alternative and null cases. Regarding the best sample size to specify at the interim, a possible balance to be struck between a suitably large interim sample, which would increase the precision of the measure of conditional dependence, and minimising the overall experimental sample size would be to take the minimal possible sample size for the experiment as a whole at the interim. In this way, the interim sample could never be larger than the overall required sample size, which means that it is impossible to collect more data than actually needed. Yet, the minimum possible overall sample size represents a significant proportion of the total experimental sample size. Thus, for the values in Table 6, the sample size re-estimate was conducted at the number implied by equation 2, when TPPR is at maximal value given the marginals. The maximum positive, mid-range and maximum negative levels of TPPR were reported to show a range of values across different levels of TPPR. The mean sample size is provided in parentheses in order to allow intuition about the reduction in the sample size this method brings. In all cases, where data were generated under the true alternative hypothesis, the simulated power is above or very close to the nominal value. Furthermore, in all cases where data were generated under the true null hypothesis the size is close to the nominal value. Comparing the mean sample sizes given for the maximal and mid-point TPPRs against the fixed values that would be used if Alonzo et al. [21] had been followed we can see that the sample size re-estimation method outlined above can dramatically reduce the required sample size to power an experiment to the minimum of a nominal level.

Table 6 Simulated type I and II error rates and fixed maximal sample size values under various true values of TPRA, TPRB and prevalence across various levels of conditional dependence (average sample size given in brackets)

Application revisited

Given the robustness of the proposed method of sample size recalculation described and validated in simulation above, we return to apply it to the application presented earlier in this paper. The cell probability values at maximum positive conditional dependence for diseased patients under the specified values of TPR A and TPR B are \( {\widehat{p}}_1=0.81,{\widehat{p}}_2=0.09 \), \( {\widehat{p}}_3=0 \), \( {\widehat{p}}_4=0.1 \). The cell probability values at maximum negative conditional dependence for diseased patients under the specified values of TPR A and TPR B are \( {\widehat{p}}_1=0.71,{\widehat{p}}_2=0.19 \), \( {\widehat{p}}_3=0.10 \), \( {\widehat{p}}_4=0 \). Table 7 shows an example range of the permissible values under the specified values of TPR A and TPR B . Given this, we can create a likelihood of our observed interim data having come from each possible configuration of the alternative hypothesis using equation 3.

Table 7 Example range of cell probabilities based on: TPRA= 0.9 and TPRB = 0.81

Applying the method outlined in section 3, we take; \( \widehat{{\ TPR}_A}=0.9 \), \( \widehat{{\ TPR}_B}=0.81 \), observed \( {\widehat{n}}_A=66 \), \( {\widehat{n}}_B=3,\kern0.5em {\widehat{n}}_C=3 \), \( {\widehat{n}}_D=10 \) and \( {\widehat{n}}_E + {\widehat{n}}_F+{\widehat{n}}_G+\kern0.5em {\widehat{n}}_H=105 \), implying \( \widehat{\pi} = 0.439 \). Using equation 4 the maximum likelihood value of \( \widehat{TPPR} \) is 0.793. Given the fact that π is binomially distributed, the maximum likelihood estimate for the prevalence is equal to the observed prevalence, \( \widehat{\pi\ } \). Taking these values and inserting them into equation 1 we get the value for the sample size required for sensitivity as 275. Taking \( \widehat{\ {TNR}_A}=0.8 \) and \( \widehat{TNR_B}=0.66 \), with the observed values \( {\widehat{n}}_E=21 \), \( {\widehat{n}}_F=4,\kern0.5em {\widehat{n}}_G=11 \), \( {\widehat{n}}_G=69 \) and \( {\widehat{n}}_A + {\widehat{n}}_B+{\widehat{n}}_C+\kern0.5em {\widehat{n}}_D=82 \), implying \( 1-\widehat{\pi} = 0.561 \). Using equation 3 to derive the maximum likelihood of the cell probabilities for specificity we estimate that \( \widehat{TNNR}=0.635 \). Inserting these values into equation 2 gives us a sample size estimate of 136. Thus, the updated sample size, in order to use the interim information about the conditional dependence between the tests and to preserve a minimal nominal power of 0.8 should be 275.

Discussion

This paper has presented a robust method of sample size re-estimation for use in paired diagnostic accuracy studies where the conditional independence between the two tests may be unknown or inaccurately estimated at the start of the study. In terms of the recommendation of sample size estimation for the experiment as a whole a specific protocol is suggested given the results. Rather than basing the estimate for the experiment as a whole on the case where there is the maximal negative conditional dependence between tests – thus the largest possible sample size - as suggested in Alonzo et al. [21], we would suggest an alternative strategy, the robustness of which is highlighted in Table 6. Specifically, initially estimating the sample size at the maximal positive conditional dependence between tests, i.e. using TPPR = TPR B - giving the smallest possible sample size - then, re-estimating the final sample size using the method simulated in Table 6. As long as the initial estimate for prevalence is close to accurate, this protocol is deemed appropriate as it balances the risk of collecting more participants than might actually be needed with collecting the most information about the true conditional dependence at the interim. Table 6 provides strong evidence for the integrity of this method in providing at minimum the nominal power while reducing the sample size when we have a higher than maximally negative true conditional dependence. Should the interim sample size be some other value, the maximum likelihood method will still be appropriate, although it should be kept in mind that the larger the interim sample size, as a proportion of the total possible sample size, the more accurate the interim sample size estimates will be, for individual cases.

Interestingly, the sample size values in the table seem to be somewhat greater, even when using our method than those typically seen in the literature in diagnostic test accuracy studies, see for example van Enst et al. [29] Although it is difficult to know the specifics of the 859 studies mentioned in the van Enst collection of meta-analyses, e.g. clinically significant differences, sample size estimation and hypothesis testing procedures, it is striking that the 50% covariance sample size is only 87 (IQR 45–185) participants. Very few of our sample sizes in Table 6 are this low for the size of effect (ratios) we are considering, even using our method of sample size reduction. It may be that many diagnostic accuracy studies commissioned do not carefully consider their sample sizes. While the method discussed here of estimating the conditional dependence between the tests via maximum likelihood, given constraints imposed by the specified marginals and under a multinomial model, is pertinent to paired diagnostic accuracy tests, there is little reason why similar processes could not be extended to similar problems. The kernel of the method, maximum likelihood estimation of the parameter related to the conditional dependence using a constrained multinomial model, is equally valid in other applications involving sample size re-estimation for paired binary 2 × 2 tables.

Conclusions

In this paper we have described a sample size re-estimation procedure that can be applied in an interim analysis for a diagnostic test study that is comparing two methods of testing on patients that are being followed up over a period of time. The procedure uses information on the levels of conditional dependence between the two tests at the interim in order to refine the required sample size for a paired diagnostic accuracy study with a binary response. Evidence from simulations has been provided to demonstrate its functionality under various parameter values thought to reflect a range of commonly occurring situations. The procedure can be applied in the case of paired comparative diagnostic accuracy studies in order to more accurately gauge the sample size required for a given power thereby reducing both the costs associated with this kind of study and also the burden on patients.

Abbreviations

CT:

computed tomography

PET:

Positron Emission Tomography

RMSE:

Root Mean Squared Error

TBR A :

True negative rate for test A

TBR B :

True negative rate for test B

TNNR :

True negative negative rate

TPPR :

True positive positive rate

TPR A :

True positive rate for test A

TPR B :

True positive rate for test B

References

  1. Knottnerus JA, van Weel C, Muris JWM. Evaluation of diagnostic procedures. Br Med J. 2002;324:477–80.

    Article  Google Scholar 

  2. Swets JA, Pickett RM. Evaluation of Diagnostic Systems. New York: Academic Press; 1982.

    Google Scholar 

  3. Zhou XH, Obuchowski NA, DK MC. Statistical Methods in Diagnostic Medicine. New York: Wiley; 2002.

    Book  Google Scholar 

  4. Freedman LS. Evaluating and comparing imaging techniques: a review and classification of study designs. Br J Radiol. 1987;60:1071–81.

    Article  CAS  PubMed  Google Scholar 

  5. Gerke O, Vach W, Høilund-Carlsen PF. PET/CT in Cancer - Methodological Considerations for Comparative Diagnostic Phase II Studies with Paired Binary Data. Methods Inf. Med. [Internet]. 2008 [cited 2014 Oct 15];470–9. Available from: http://www.schattauer.de/index.php?id=1214&doi=10.3414/ME0540.

  6. Gould AL. Sample size re-estimation: recent developments and practical considerations. Stat Med. 2001;20:2625–43.

    Article  CAS  PubMed  Google Scholar 

  7. Wittes J, Brittain E. The role of internal pilot studies in increasing the efficiency of clinical trials. Stat Med. 1990;9:65–72.

    Article  CAS  PubMed  Google Scholar 

  8. Shih WJ. Sample size reestimation in clinical trials. In: Peace K, editor. Biopharm. Seq. Stat. Appl. New York: Marcel Dekker; 1992. p. 285–301.

    Google Scholar 

  9. Shih WJ. Sample size reestimation for triple blind clinical trials. Drug Inf J. 1993;27:761–4.

    Article  Google Scholar 

  10. Gould AL, Shih WJ. Sample size re-estimation without unblinding for normally distributed outcomes with unknown variance. Commun Stat Methods. 1992;21:2833–53.

    Article  Google Scholar 

  11. Birkett MA, Day SJ. Internal Pilot Studies for Estimating Sample Size. Stat Med. 1994;13:2455–63.

    Article  CAS  PubMed  Google Scholar 

  12. Gould AL. Interim Analyses for Monitoring Clinical Trails that do not Affect Type I Error Rates. Stat Med. 1992;11:55–66.

    Article  CAS  PubMed  Google Scholar 

  13. Herson J, Wittes J. The Use of Interim Analysis for Sample Size Adjustment. Drug Inf J. 1993;27:761–4.

    Article  Google Scholar 

  14. Shih WJ, Zhao P. Design for Sample Size Re-estimation with Interim Data for Double-Blind Clinical Trails. Stat Med. 1997;16:1913–23.

    Article  CAS  PubMed  Google Scholar 

  15. Proschan MA. Two-stage sample size re-estimation based on a nuisance parameter: a review. J Biopharm Stat. 2005;15:559–74.

    Article  PubMed  Google Scholar 

  16. Gerke O, Høilund-carlsen PF, Poulsen MH, Vach W. Interim analyses in diagnostic versus treatment studies : differences and similarities. Am J Nucl Med Mol Imaging. 2012;2:344–52.

    PubMed  PubMed Central  Google Scholar 

  17. Lijmer JG, Jeroen G. Empirical Evidence of Design-Related Bias. JAMA. 1999;282:1061–6.

    Article  CAS  PubMed  Google Scholar 

  18. Lord SJ, Staub LP, Bossuyt PMM, Irwig LM. Target practice : choosing target conditions for test accuracy studies that are relevant to clinical practice. Br Med J. 2011;343:1–5.

    Article  Google Scholar 

  19. Newcombe RG. Improved Confidence Intervals for the Difference between Binomial Proportions Based on Paired Data. Stat Med. 1998;17:2635–50.

    Article  CAS  PubMed  Google Scholar 

  20. Tango T. Equivalence Test and Confidence Interval for the Difference in the Proportions Based on Paired Data. Stat Med. 1998;17:891–908.

    Article  CAS  PubMed  Google Scholar 

  21. Alonzo TA, Pepe MS, Moskowitz CS. Sample Size Calculations for Comparative Studies of Medical Tests for Detecting Presence of Disease. Stat Med. 2002;21:835–52.

    Article  PubMed  Google Scholar 

  22. Lu Y, Jin H, Genant HK. On the Non-Inferiority of a Diagnostic Test Based on Paired Observations. Stat Med. 2006;3:227–79.

    Google Scholar 

  23. Moskowitz CS, Pepe MS. Comparing the Predictive Values of Diagnostic Tests: Sample Size and Analysis for Paired Study Designs. Clin Trials. 2006;3:272–9.

    Article  PubMed  Google Scholar 

  24. Bonett DG, Price RM. Confidence Intervals for a Ratio of Binomial Proportions Based on Paired Data. Stat Med. 2006;25:3039–47.

    Article  PubMed  Google Scholar 

  25. Vacek P. The effect of conditional dependence on the evaluation of diagnostic tests. Biometrics. 1985;41:959–68.

    Article  CAS  PubMed  Google Scholar 

  26. van Smeden M, Naaktgeboren CA, Reitsma JB, Moons KGM, de Groot JA. Latent Class Models in Diagnostic Studies When There is No Reference Standard — A Systematic Review. Am J Epidemiol. 2014;179:423–31.

    Article  PubMed  Google Scholar 

  27. Schiller I, van Smeden M, Hadgu A, Libman M, Reitsma B, Dendukuri N. Bias due to composite reference standards in diagnostic accuracy studies. Stat Med. 2016;35:1454–70.

    Article  PubMed  Google Scholar 

  28. Royston P. Exact conditional and unconditional sample size for pair-matched studies with binary outcome: a practical guide. Stat Med. 1993;12:699–712.

    Article  CAS  PubMed  Google Scholar 

  29. van Enst WA, Naaktgeboren CA, Ochodo EA, de Groot JA, Leeflang MM, Reitsma JB, et al. Small-study effects and time trends in diagnostic test accuracy meta-analyses : a meta-epidemiological study. Syst Rev. 2015;4:66.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This project was funded by the NIHR Health Technology Assessment Programme, REF: 08/29/02.

Availability of data and materials

All data used to illustrate the method can be found in Table 2 of this paper.

Author information

Authors and Affiliations

Authors

Contributions

GL conceived the study. The method was devised by GM, overseen by GL and AT. The manuscript was written by GM. All authors contributed to the reviewing and revising of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Gareth P. J. McCray.

Ethics declarations

Ethics approval and consent to participate

REC reference: 10 H1017 8. The original application was made to North West 1 Research Ethics Committee - Cheshire. This committee has subsequently been superseded by the North West - Greater Manchester East Research Ethics Committee. Consent to participate was provided by all participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Appendix R code for maximum likelihood sample size re-estimation

Appendix R code for maximum likelihood sample size re-estimation

ss.est.mle <- function(obs.a, obs.b, obs.c, obs.d, obs.x, tpra, tprb, alpha, beta){ mle.tppr <- function(theta.1, obs.a, obs.b, obs.c, obs.d, tpra, tprb) {-((obs.a*log(theta.1)) + obs.b*log(tpra-theta.1) + obs.c*log(tprb-theta.1) + obs.d*log(-tpra-tprb+theta.1+1))} tppr <- optim(par=(tpra+(tprb-(1-tpra)))/2 ,fn=mle.tppr, obs.a=obs.a, obs.b=obs.b, obs.c=obs.c, obs.d=obs.d, tpra=tpra, tprb=tprb, method = "Brent", lower=(tprb-(1-tpra)), upper=tprb)$par obs.prev <- (obs.a+obs.b+obs.c+obs.d)/(obs.a+obs.b+obs.c+obs.d+obs.x) alonzo <- function(lambda, prev ,beta, alpha, tprb, gam1){(((qnorm(1-beta) + qnorm(1-alpha))/log(gam1))^2 * (((gam1+1) * tprb)-(2 * lambda))/(gam1*tprb^2))/prev} gam1 <- tpra/tprb ss.est <- alonzo(tppr, obs.prev, beta = beta, alpha = alpha, tprb = tprb, gam1 = gam1 ) return(ss.est)} ### Example sensitivity ss.est.mle(obs.a=66, obs.b=3, obs.c=3, obs.d=10, obs.x=105, tpra=0.9, tprb=0.81, alpha=0.025, beta=0.2) ### Example specificity ss.est.mle(obs.a=69, obs.b=11, obs.c=4, obs.d=21, obs.x=82, tpra=0.8, tprb=0.66, alpha=0.025, beta=0.2)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

McCray, G.P.J., Titman, A.C., Ghaneh, P. et al. Sample size re-estimation in paired comparative diagnostic accuracy studies with a binary response. BMC Med Res Methodol 17, 102 (2017). https://doi.org/10.1186/s12874-017-0386-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12874-017-0386-5

Keywords