Sample size and power considerations in network meta-analysis
Systematic Reviews volume 1, Article number: 41 (2012)
Network meta-analysis is becoming increasingly popular for establishing comparative effectiveness among multiple interventions for the same disease. Network meta-analysis inherits all methodological challenges of standard pairwise meta-analysis, but with increased complexity due to the multitude of intervention comparisons. One issue that is now widely recognized in pairwise meta-analysis is the issue of sample size and statistical power. This issue, however, has so far only received little attention in network meta-analysis. To date, no approaches have been proposed for evaluating the adequacy of the sample size, and thus power, in a treatment network.
In this article, we develop easy-to-use flexible methods for estimating the ‘effective sample size’ in indirect comparison meta-analysis and network meta-analysis. The effective sample size for a particular treatment comparison can be interpreted as the number of patients in a pairwise meta-analysis that would provide the same degree and strength of evidence as that which is provided in the indirect comparison or network meta-analysis. We further develop methods for retrospectively estimating the statistical power for each comparison in a network meta-analysis. We illustrate the performance of the proposed methods for estimating effective sample size and statistical power using data from a network meta-analysis on interventions for smoking cessation including over 100 trials.
The proposed methods are easy to use and will be of high value to regulatory agencies and decision makers who must assess the strength of the evidence supporting comparative effectiveness estimates.
Over the past 2 decades, meta-analysis has become increasingly accepted by clinicians, decision-makers and the public as providing high-quality assessments of evidence. Network meta-analysis, a new expansion of meta-analysis that allows for simultaneous comparison of several treatments, is similarly becoming increasingly accepted in the clinical research community[2–11]. Having been available for more than 3 decades, meta-analysis has been studied extensively, and several hundred articles published in this period have identified and resolved a vast array of basic and advanced methodological issues[1, 12]. Network meta-analysis inherits all the challenges present in a standard meta-analysis (e.g., issues of bias, heterogeneity and precision), but with increased complexity due to the multitude of comparisons involved[5, 13]. Since network meta-analysis is still a relatively new technique, the number of publications addressing methodological challenges is still relatively sparse.
One important issue that has received much attention in individual trials and meta-analysis is the issue of sample size and statistical power[14–26]. Several studies have demonstrated the importance of interpreting pooled meta-analysis estimates and confidence intervals according to the statistical level of evidence (i.e., precision)[14, 15, 18–24, 26, 27], and sound recommendations have been provided[16, 22, 28]. So far, however, only a small number of studies have addressed the issue of power and precision in network meta-analysis[8, 13, 29], and no comprehensive guidance exists on the topic. Network meta-analyses typically include many more trials than standard meta-analyses because of the multitude of comparisons involved and for this reason may artificially appear to provide a stronger evidence base. Likewise, the accompanying graphical representation of a treatment network can provide a similar compelling but potentially false impression of a strong evidence base.
Network meta-analysis utilizes evidence from direct (head-to-head) comparisons (i.e., trials directly comparing treatment A and B) and indirect comparisons (e.g., the combination of trials comparing A with C and trials comparing B with C)[4, 6, 30]. The major challenge in interpreting the power and precision of a network meta-analysis stems from the fact that there are (typically) varying levels of power and precision across all comparisons. In addition, the power and precision of indirect evidence are more complex to assess than for direct evidence, and thus, without proper guidance, it will be difficult for most authors to evaluate the precision gain from use of indirect evidence as well as the strength of evidence in a treatment network.
In this article, we provide guidance on quantifying the power and precision in network meta-analysis using simple sample size considerations. We first describe how to quantify the precision in indirect comparison meta-analysis and subsequently in network meta-analysis with combinations of direct and indirect evidence. We then outline the concept of sample size requirements and power calculations in pairwise meta-analysis. Finally, we show how to combine these measures in order to quantify the power and strength of evidence available for all treatment comparisons in a network meta-analysis. We illustrate the described methods using data from a recent network meta-analysis on interventions for smoking cessation.
Basic methodological framework
Indirect effect estimates are obtained with the effect estimates from two comparisons sharing a common comparator. For example, when two treatments A and B have both been compared to some common comparator C (e.g., placebo) in a number of randomized clinical trials, an indirect effect estimate of treatment A versus B can be obtained using the meta-analysis effect estimate of A versus C and the meta-analysis effect estimate of B versus C. In particular, the indirect effect estimate of A versus B (d AB ) is calculated as the estimated effect of A versus C (d AC ) minus the estimated effect of B versus C (d BC ). Mathematically, this corresponds to the equation
(Note, when dealing with ratio effect measures, such as relative risks and odds ratios, all calculations are done on the log scale to preserve linearity and approximate normality). To produce confidence intervals for the indirect estimate, we first need to estimate its variance. The variance of in the indirect estimate of A versus B (V AB ) is simply equal to the sum of the variance of the effect estimate of A versus C (V AC ) and the variance of the effect estimate of A versus B (V AB ). Mathematically this corresponds to the equation V AB = V AC + V BC . It is therefore clear that the variance of a (direct) meta-analysis effect estimate based on some number of trials, say k, will always be smaller than the variance of an indirect meta-analysis based on the same number of trials, k (all trial sample sizes being equal). In other words, direct estimates come with higher precision and power (trial count and trial sample sizes being equal). In many situations, however, using indirect estimates can add considerable power and precision.
Combining direct and indirect evidence
When both direct and indirect evidence is available, it may often be advantageous to combine the two statistically[2, 4–7, 30, 33]. For example, if only two small trials have investigated two active interventions A and B head to head, but 20 trials have compared A or B with placebo, the indirect evidence will be able to add much power and precision to the comparative estimate of A and B. The combination of indirect and direct evidence requires advanced statistical regression techniques (i.e., network meta-analysis) that are beyond the scope of this article[4, 6, 30]. However, in the context of sample size and power considerations, it suffices to understand that indirect evidence, when combined with direct evidence, increases the power and precision of treatment effect estimates[4, 6, 7, 9, 30, 33]. The extent to which is does so can be evaluated readily by using the methods we describe below.
Sample size in indirect comparisons
In this section we introduce three methods for gauging how much statistical precision an indirect estimate provides when no direct evidence is available. In particular, we describe how to approximate the amount of information required in a direct (head-to-head) meta-analysis to produce the same precision as that in the available indirect evidence. Simply put, what direct meta-analysis sample size would provide a similar degree of information? We dub this the effective sample size of the indirect evidence or, interchangeably, the effective indirect sample size. We describe three different methods for approximating the effective indirect sample size. Each of these methods differs with respect to simplicity and validity (the simpler one being the least valid), so we outline the simplicity-validity trade-offs at the end of the section.
Method 1: the effective number of trials
A simple approach to gauging the degree of power and precision available in indirect evidence is to approximate how many trials are required in an indirect comparison to produce a matching degree of power and precision from a single head-to-head trial. This type of approximation is possible under the simple assumptions that the variances (of the mean) are equal for each trial and that no heterogeneity is present. Glenny et al. showed that when the number of trials is the same in both of two comparisons informing the indirect evidence (e.g., two trials of A vs. C and two trials of B vs. C), it takes four trials in the indirect evidence to produce the same precision as one direct head-to-head trial. In indirect comparisons, however, it is common that one comparison will include more trials than the other. When this happens, the above 1:4 precision ratio no longer holds true. For example, if the number of trials is twice as high in one comparison (i.e., a 1:2 trial count ratio), the indirect comparison will need exactly 4.5 trials to produce the same precision as one head-to-head trial (see mathematical derivation in Appendix 1.a). In reality, however, to maintain a ratio of 1:2 in the trial count, one would need six trials (2:4) in the indirect comparison to produce the same precision as one head-to-head trial. To produce the same precision as two head-to-head trials, one would need 2 × 4.5 = 9 trials, which allows maintaining the 1:2 ratio with three trials in one comparison and six in the other (i.e., 3:6 as the trial count ratio). Table 1 presents the approximate number of trials required in an indirect comparison under different scenarios where the number of trials in the two comparisons is unbalanced. The mathematical derivations for all exact precision ratios are presented in Appendix 1.a.
The cells in underlined italics indicate where the indirect evidence produces the exact precision of the corresponding number of trials. The remaining cells indicate where the indirect evidence produces precision slightly above that of the corresponding number of head-to-head trials.
In some cases the required number of trials in an indirect comparison for a specific trial count ratio produces a precision corresponding to more than that of a stated number of single head-to-head trials. For example, with a trial count ratio of 1:3, one would require 2 × 5.33 = 10.66 indirect comparison trials to produce the precision of two head-to-head trials. However, since we cannot have fractions of trials, we take the closest integer above 10.66 where the trial count ratio is maintained: 12 trials with a trial count ratio of 3:9.
Table 1 can readily be used for quickly and easily checking how many head-to-head trials the indirect evidence ‘effectively’ corresponds to. That is, if the indirect evidence produces the precision of, say, three trials, we can think of the evidence as being as strong as a meta-analysis of three head-to-head trials. For example, if one has an indirect comparison with 4 trials comparing A with C, and 12 trials comparing B with C, the precision of the indirect comparison corresponds to a meta-analysis of 3 trials directly comparing A with B (Table 1). It should be noted that Table 1 is only valid to the extent that trial sample sizes and trial population variances are similar across trials, as well as the extent to which heterogeneity is absent or ignorable.
Method 2: the effective sample size
Another relatively simple approach to gauging the degree of power and precision from an indirect comparison is to consider the collection of trials included in each comparison as one (large) clinical trial. From a sample size perspective, following similar mathematic derivations as the above trial count perspective, the relationship between the precision of an indirect comparison and the precision of a direct meta-analysis turns out the same (see Appendix 1.b for mathematical derivations). For example, to produce the same precision as a head-to-head meta-analysis including 1,000 patients, one would need a total of 4,000 (4 × 1,000) patients in the indirect comparison, provided the number of patients is the same for the two comparisons (2,000:2,000). Taking the sample size perspective comes with the flexibility of a possible reversal of the calculation. For example, an indirect comparison with 500 patients in the A vs. C comparison and 500 patients in the B vs. C comparisons would produce the same precision as a direct comparison with 250 patients [(500 + 500)/4]. Likewise, in a scenario with 1,000 patients in comparison A vs. C, and 10,000 patients in comparison B vs. C, the exact precision ratio is 12.1 (see Table 1), and so the effective direct meta-analysis sample size would be (1,000 + 10,000)/12.1 = 909.
Often, the sample sizes in the two comparisons do not line up to produce the exact precision ratio presented in Table 1 and Table 2. Letting n AC and n BC denote the sample sizes for the comparisons of A vs. C and B vs. C, respectively, a more general formula for the effective indirect sample size is (see Appendix 1.b).
For the above example, the effective sample size using this formula is therefore
The above simple approach to calculate the effective indirect sample size does not consider the possibility that statistical heterogeneity exists across trials. When heterogeneity is present, the effect estimates that go into the indirect estimate incur a higher degree of variation, and so the effective indirect sample size corresponding to a head-to-head meta-analysis will be smaller than with the above simple approach. In line with already-established heterogeneity corrections for meta-analysis required sample sizes[23, 24, 28], we put forward that the actual number of patients in each of the comparisons informing the indirect estimate can be penalized by the additional variation explained by heterogeneity. In line with previous proposals, we penalize for the ‘lack of homogeneity’[23, 24] using the popular measure of heterogeneity, I 2, as a basis for the penalization.
Consider the example where a meta-analysis of A vs. C includes 6,000 patients and a meta-analysis of B vs. C includes 8,000 patients, and assume the estimated degree of heterogeneity for A vs. C is 50% (I AC 2 = 50%) and 25% for B vs. C (I BC 2 = 25%). Then the lack of homogeneity is 100%-50% = 50% for A vs. C and 100%-25% = 75% for B vs. C. We penalize the actual sample size by multiplying the actual sample size by the lack of homogeneity, so that the penalized sample size of A vs. C is 50% × 6,000 = 3,000, and the penalized sample size for B vs. C is 75% × 8,000 = 6,000. The total penalized number of patients in the indirect comparison is then 3,000 + 6,000 = 9,000, the patient count ratio is 1:2 (same as 3,000:6,000), the precision ratio is 4.5 (see Table 1), and so the effective heterogeneity-corrected sample size in this indirect comparison is 9,000/4.5 = 2,000.
Following the above example, the general formula for a heterogeneity-corrected effective sample size for indirect evidence is
where n AC and n BC are the actual sample sizes (before correction) in the meta-analyses of A vs. C and B vs. C, respectively, and where the precision ratio is based on the heterogeneity-corrected sample sizes (see Table 1 and Table 2). In the appendix we provide a general formula for the precision ratio.
As with the above example of non-penalized sample sizes, the penalized sample sizes may not always line up to match the precision ratio given in Table 1 and Table 2. The more general formula for the heterogeneity-corrected effective sample size is (see Appendix 1.b).
One immediate limitation of the above-proposed sample size heterogeneity correction is the fact that I 2s are typically unreliable and unstable in meta-analyses including a limited number of trials and will depend on the effect metric used[34–36]. In most cases, it will therefore be preferable to simply assume some plausible degree (percentage) of heterogeneity derived from a mix of clinical considerations and the I2 estimate at hand. Typically, an assumption of 25% or 50% heterogeneity will be reasonable in the context of sample size considerations. Table 2 illustrates the effective sample size for various indirect comparison scenarios including a total of 10,000 patients, under different combinations of heterogeneity corrections (one including no correction).
Another limitation of the sample size approach is the inherent assumption that sample size is a good proxy for precision. This may not be true if there are some important differences in event rates (for binary data) or counts (count data), or if there are population differences in trials that result in notably different standard deviations, but not necessarily different effect estimates. To curb this limitation for binary data, one may for example choose to focus on the effective number of events. A more universally applicable approach focuses on a measure called statistical information. We describe this below.
The effective statistical information
Statistical information, also sometimes referred to as Fisher information, is a more complex statistical measure for gauging the degree of precision present in a data set. For pairwise meta-analyses, the statistical information is equal to the inverse of the pooled variance (i.e., one divided by the variance), which is also the measure of precision. For indirect comparisons, the statistical information (precision) is equal the inverse of the pooled indirect variance. That is, with variances V AC and V BC for comparisons A vs. C and B vs. C, the indirect variance is V AC + V BC , so indirect statistical information is 1/(V AC + V BC ). Because the variance incorporates heterogeneity and is dependent on the number of trials and the sample size, no further calculations or adjustments are needed.
A major disadvantage of statistical information in meta-analysis is that it operates on a scale that no statistically non-sophisticated reader could easily grasp. The statistical information may however be useful in the context of sample size requirements since it is possible to calculate the required statistical information (analogous to the required sample size) and compare the actual statistical information present in an indirect comparison with such a yardstick for sufficient power.
Strength and limitations of the approaches
Each of the above approaches comes with strengths and limitations. These are outlined in Table 3.
So far we have only described situations where the indirect estimate is obtained through one common comparator. However, it is worth noting that indirect evidence may often come from scenarios where the link between two treatments of interest must go through two or more comparators. For example, if we wish to compare A with D and the following three comparisons are available—A vs. B, B vs. C and C vs. D—then the link between A and D goes through both B and C. The indirect variance is now a sum of three (rather than two) direct variances, . We would therefore expect comparably smaller precision. Calculating the effective number of trials in the two common comparator example, we find that with an equal number of trials, say 3:3:3, the precision of the indirect comparison corresponds to that of only one direct trial, that is, an exact precision ratio of 1:9. With a more unbalanced number of trials, say 8:1:8, we get an exact precision ratio of 1:21.
Such consistently large direct to indirect precision ratios indicate that indirect comparisons with two (or more) comparators in the link will typically add very little precision. For this reason we dub them ‘weak links.’ In the context of combining direct and indirect evidence as we address below, it seems that weak links will typically only add an ignorable small amount of precision to the final estimate and thus, for simplicity, may be ignored for sample size and power considerations in indirect comparisons and network meta-analysis.
Effective sample size in treatment networks
The three-treatment loop
The simplest example of a combination of direct and indirect evidence is the three-treatment loop where precision is added to the comparison of A and B by borrowing strength from an indirect comparison based on some common comparator C. Whether our measure of precision (information) is the number of trials, the sample size or the statistical information, the total amount of precision available for a particular comparison in a three-treatment loop is conceptually the sum of information in the direct evidence and in the indirect evidence.
To calculate the effective number of trials informing a particular comparison (e.g., A vs. B) in the combined evidence, we simply add the effective number of trials in the indirect evidence to the number of trials in the head-to-head evidence.
To calculate the effective number of patients informing a particular comparison, we simply add the effective number of patients in the indirect evidence to the number of patients in the head-to-head evidence. If heterogeneity adjustments have been applied for the indirect evidence, a similar adjustment should be applied to the direct evidence (we illustrate this in our worked example below).
To calculate the statistical information informing a particular comparison, we simply take the sum of the inverse indirect variance and the inverse direct variance. Alternatively, the statistical information may be extracted directly from the statistical software package used to run the network meta-analysis. (Note that if multi-arm trials are included in the analysis, some adjustments for correlations are needed).
Multiple sources of evidence
In many treatment networks two or more sources of indirect evidence informing some (or all) of the comparisons may exist. For example, two active treatments A and B may both have been compared to standard-of-care and placebo. In this case, indirect evidence exists from two sources. Extending the three-treatment example above, estimating the total effective amount of information is simply a task of summing all indirect and direct evidence.
In a similar situation, multiple sources of indirect evidence may exist where no direct evidence exists. In this case, the effective number of trials, sample size or statistical information is simply obtained by summing all indirect information.
We previously discussed ‘weak links.’ Since these add a relatively small amount of information, one may chose to ignore them without notable loss of information (but with ease in the calculations). This goes for both situations where direct evidence and indirect evidence are being combined, and where there are multiple sources of indirect evidence.
Power and sample size requirements in network meta-analysis
Before getting into sample size and power considerations for indirect comparison meta-analysis and network meta-analysis, we first outline the already well-established framework for pairwise meta-analysis. We then extend the concept to indirect comparison meta-analysis and lastly to network meta-analysis.
Sample size requirements for direct meta-analysis
Several methodological studies have explored sample size and power considerations for direct (head-to-head) meta-analysis[14, 15, 19–24, 26]. By now, it has been well established that the required sample size (i.e., the required number of patients) for a meta-analysis, should be at least that of a large well-designed clinical trial. Sample size calculations are derived from an a priori estimate of a treatment effect, d, that investigators wish to demonstrate; the associated variance around that treatment effect, V 2; and a maximum risk of type I error, α, (i.e., maximum false-positive risk) and type II error, β, (i.e., maximum false-negative risk). As a basis, one can use the required sample size corresponding to a large multicenter clinical trial as the required sample size for the head-to-head meta-analysis[16, 19, 20, 23, 24].
Here z 1-α/2 and z 1-β are the (1-α/2)th and (1-β)th percentiles of a standard normal distribution, and C is a constant depending on the randomization ratio and number of treatment arms (C = 4 with a randomization ratio of 1:1 and two treatment arms).
If statistical heterogeneity exists across the included trials in a meta-analysis, one can adjust the calculated sample size to account for the additional variation (i.e., increased uncertainty)[21–24, 28]. This is achieved by multiplying the required sample size, N, by a heterogeneity correction factor 1/(1-H), where H has a similar interpretation as the well-known measure I 2 (the percentage of variation in a meta-analysis explained by heterogeneity) and is the a priori or maximum acceptable degree of heterogeneity[21–24, 28]. Empirical evidence and simulations have demonstrated that such adjustments perform well in maintaining the desired statistical power[21, 22, 24].
An alternative approach to dealing with heterogeneity is to calculate the required statistical information (also known as Fisher information)[17, 25]. In pairwise meta-analysis the statistical information is simply the inverse of the pooled variance, that is, the pooled precision. The required statistical information resembles the required sample size.
A simulation study has demonstrated adequate performance of the required statistical information when the heterogeneity is modeled with certain Bayesian priors.
Information fractions and power in direct meta-analysis
At any point in a direct meta-analysis before the cumulative amount of evidence has surpassed the required sample size (or required statistical information), we can calculate two useful measures to gauge the strength of evidence. The first measure, the information fraction (IF), is the accrued number of patients, n (or statistical information), divided by the required sample size (or required statistical information)[19, 23].
This measure gives us an idea of how far we have come and how much farther of a distance there is to the yardstick—our required sample size.
The second measure is a retrospective power calculation. Re-arranging the expression of the required sample size, we can retrospectively estimate the power (1-β) of the current data
where Φ is the cumulative standard normal distribution function.
Information fractions and power in indirect comparisons
So far we have described method for estimating the effective sample size in indirect and combined evidence, as well as methods for estimating the required sample size and gauging the strength of evidence in pairwise meta-analyses. We now combine these measures to evaluate the strength of indirect comparisons and network meta-analyses. For the remainder of this article we concentrate on the number of patients, but our example and calculations can easily be carried out using the number of trials or statistical information.
We have introduced the effective sample size for indirect evidence. In line with its definition, we can derive the effective indirect information fraction and the effective power in an indirect comparison. The steps required to do this are as follows. First, calculate the effective indirect sample size. Second, calculate the required sample size for a direct meta-analysis. Third, to get the effective indirect information fraction, simply divide the effective number of patients in the indirect comparison by the required sample size for a direct comparison. Fourth, to calculate the power of the available indirect evidence, simply insert the effective indirect sample size as the n in the above formula.
Information fractions and power in network meta-analysis
The steps required to calculate the information fraction and power of treatment comparisons that are both informed by direct and indirect evidence or by multiple sources of indirect evidence are similar to the step required for indirect comparisons. First, calculate the effective sample size for the comparison of interest by summing up the evidence from the available sources (we described how to do this above). Second, as before, calculate the required sample size for a direct meta-analysis. Third, as before, calculate the effective information fraction by dividing the effective number of patients by the required sample size. Fourth, calculate the power of the available evidence by inserting the effective sample size as the n in the above formula.
Results and discussion
Worked example – interventions for smoking cessation
A recent MTC explored the efficacy of five different interventions for smoking cessation, low-dose NRT (<22 mg nicotine patches), which is an over-the-counter intervention, and four newer and more expensive interventions, combination NRT (i.e., patch plus gum), high-dose NRT (>22 mg patches), buproprion and varenicline.The published MTC provides a useful data-set to illustrate the issues we have raised in this article. Low-dose NRT is already well established and well known to yield about a 50% higher success rate than inert control interventions (i.e., relative risk of approximately 1.50). To consider any of the four newer treatments worthwhile, one could argue for the need to demonstrate at least an additional 20% efficacy. These considerations make up the basis for the required sample size considerations. Keep in mind that the published MTC has much higher rates of effective sample size than the assumptions we make in this example.
The median smoking cessation success rate at 6 months for patients in the inert control group is 15.0% across the 79 included trials reporting at this time point. Assuming a 1.5 relative risk, we would expect a 22.5% success rate with low-dose NRT. An additional 20% relative efficacy would suggest a 26.0% success rate (with any of the newer treatments). Accepting a maximum type I error of α = 5% and a power of (1-β) = 90%, the required sample size to demonstrate that any of the newer interventions are at least 20% better than low-dose NRT is 6,303. To assess the strength of the evidence contained in the treatment network, we calculate the effective sample size (number of patients) for each of the newer treatments compared with low-dose NRT. We do this with and without taking heterogeneity into account. We subsequently calculate the information fraction and power based on the effective sample size. Figure 1 presents the number of trials, number of patients and degree of heterogeneity in the comparisons informed by head-to-head evidence. Figure 2 presents the sources of direct and indirect evidence for the four comparisons of interest (newer treatments vs. low-dose NRT).
For low-dose NRT vs. combination NRT, the direct evidence includes 1,664 patients (and no heterogeneity). Indirect evidence exists with inert control as the common comparator. The comparison of low-dose NRT and inert control includes 19,929 patients, but with 63% heterogeneity, so the heterogeneity penalized sample size is 19,929 × (1–0.63) = 7,374. The comparison of combination NRT vs. inert control includes 1,848 patients (and no heterogeneity). The effective indirect sample size without heterogeneity penalization (n indirect ) and the effective indirect sample size with heterogeneity penalization (n indirect-Pen ) are therefore
Adding this to the direct evidence sample size, we get an effective total sample size of
These two total effective sample sizes correspond to an information fraction of 53% and 50% (half of the require sample size has been accumulated) and a statistical power of 66% and 63%. Table 4 presents the direct effective indirect and total sample sizes and the corresponding information fractions and statistical power for all four comparisons of the newer treatments vs. low-dose NRT. The calculations for the remaining three comparisons are presented in Appendix 2.
In the network meta-analysis of interventions for smoking cessation, high-dose NRT and varenicline were both significantly better than low-dose NRT and demonstrated effect estimates larger than the a priori considered minimally important difference. For these comparisons, high-dose NRT is supported by 88% power, and varenicline is supported by 76% power. Considering that the true effects of high dose NRT and varenicline over low dose NRT are much higher than the 20% increase that was assumed for these calculations, the true power of these comparisons are also much higher. For example, if the MTC actual effects for varenicline (i.e., a 38% increase in smoking cessation compared with low dose NRT) are used, the statistical power to detect this difference then exceeds 99%. Combination NRT and buproprion both yielded effect estimates very similar to low-dose NRT. Considering that the two are supported by 66% and 95% power, respectively, and that none of the two effect estimates appear superior, it should be reasonable to infer that the two interventions do not offer any noteworthy benefit over low-dose NRT.
In this article we have outlined available methods for gauging the strength of the evidence in a network meta-analysis using sample size and power considerations. We recommend sample size considerations in the context of the number of patients, as the required calculations are relatively straightforward and will resonate well with most clinicians and decision-makers. The methods we have outlined are of high value to regulatory agencies and decision makers who must assess the strength of the evidence supporting comparative effectiveness estimates.
1.a Calculating the effective number of trials
Consider the situation where three treatments, A, B and C, have been compared head to head in randomized clinical trials. For any one trial, assume that the estimated treatment effect has variance v. For a meta-analysis of 2k trials, using the inverse variance approach would produce an estimated variance of the pooled treatment effect of σ 2 /2k. By the expected variance of an indirect comparison, if we have two comparison including k trials, we would expect an indirect variance estimate of σ 2 /k + σ 2 /k = 2σ 2 /k. Now letting R denote a ratio describing the relationship between the precision of indirect and direct evidence; we can derive R as follows
That is, in the scenario where the number of trials are equal in the two comparisons informing the indirect comparison (and the other above assumptions are met), it would require four trials in the indirect evidence to produce the same precision as that corresponding to a single head-to-head trial. We can generalize this ratio to the situation where the number of trials is not equal in the two comparisons informing the indirect evidence. Let k AC and k BC be the number of trials informing the comparison of A vs. C and B vs. C, respectively. For a single meta-analysis, with k AC + k BC trials we would expect a variance of the pooled effect of σ 2 /(k AC + k BC ). Moreover, we would expect a variance from the indirect comparison of σ 2 /k AC + v/k BC . Proceeding as above we then have
This formula creates the basis for the results presented in Table 1.
1.b Calculating the effective number of patients
Consider the situation where three treatments, A, B and C, have been compared head to head in randomized clinical trials. Assume that the population variance of comparative treatment effects is the same for A vs. B, A vs. C and B vs. C, and assume the population variance produced by a fixed-effect pairwise meta-analysis can be regarded as a large well-designed clinical trial. Let n AB , n AC and n BC denote the meta-analysis sample size (total number of patients) for the three comparisons A vs. B, A vs. C and B vs. C, respectively.
We are interested in finding the ratio between the variance of the direct meta-analysis pooled treatment effect estimate and the variance of the indirect meta-analysis pooled treatment estimate. Let R denote this ratio, and let σ AB 2, σ AC 2 and σ BC 2 denote the population variances for the three comparisons (where we assume σ AB 2 = σ AC 2 = σ BC 2 = σ 2). Then we have
Thus, by multiplying this ratio with the total indirect sample size (nAC + nBC) we have that the formula for the effective indirect sample size is
When heterogeneity exists for one or both of the comparisons in the indirect evidence, one can penalize the sample size by multiplying by the ‘lack of homogeneity,’ much similar to what is done for a heterogeneity correction of a required meta-analysis sample size. With estimates of the percentage of variation in the meta-analysis due to between-trial heterogeneity for A vs. C, I AC 2 , and for B vs. C, I BC 2 , we can derive penalized sample sizes within each comparison
and subsequently use these penalized sample sizes in the formula for the effective indirect sample size.
2. Information fraction and power calculations – worked example
For low-dose NRT vs. high-dose NRT, the direct evidence includes 3,605 patients (and no heterogeneity). Indirect evidence exists with inert control as the common comparator. The comparison of low-dose NRT and inert control includes 19,929 patients, but with 63% heterogeneity, so the heterogeneity penalized sample size is 19,929×(1–0.63) = 7,373. The comparison of high-dose NRT vs. inert control includes 2,487 patients, but with 60% heterogeneity, so the heterogeneity penalized sample size is 2,487×(1–0.60) = 1,492. The effective sample size from this indirect comparison is therefore
A second indirect comparison with varenicline as the common comparator only includes 32 patients in one of the two involved comparisons. The effective sample size of this indirect comparison (n indirect = 31) is so comparably small that we choose to ignore it. Adding the above calculated indirect sample sizes to the direct evidence sample size, we get effective total sample sizes of
This total effective sample sizes correspond to information fractions of 92% and 60% and statistical power estimates of 88% and 72%.
For low-dose NRT vs. buproprion, no direct evidence exists. Indirect evidence exists through inert control as the common comparator. As above, the sample size for low-dose NRT vs. inert control is 19,929, or 7,373 if heterogeneity penalized. The sample size for buproprion vs. inert control is 12,567, or 12,567×(1–0.39) = 7,666 when heterogeneity is penalized. Therefore, the total effective sample sizes (which are equal to the effective indirect sample sizes) are
This total effective sample sizes correspond to information fractions of >100% and 60% and statistical power estimates of 95% and 71%.
For low-dose NRT vs. varenicline, the direct evidence includes 740 patients (and no heterogeneity). As above, the sample size for low-dose NRT vs. inert control is 19,929, or 7,373 if heterogeneity is penalized. The sample size for varenicline vs. inert control is 4,331, or 4,331 × (1–0.69) = 1,343 if heterogeneity is penalized. Therefore, the total indirect sample sizes are
and so the total effective sample are
All power, and information fraction calculations above are geared to detect an assumed relative improvement in smoking cessation of 20%. All calculations are highly sensitive to the assumed relative improvement. In particular, assuming larger improvements would result in substantially larger power and information fraction estimates.
Sutton AJ, Higgins JP: Recent developments in meta-analysis. Stat Med. 2008, 27: 625-650. 10.1002/sim.2934.
Ioannidis JP: Integration of evidence from multiple meta-analyses: a primer on umbrella reviews, treatment networks and multiple treatments meta-analyses. Cmaj. 2009, 181: 488-493. 10.1503/cmaj.081086.
Jansen JP, Crawford B, Bergman G, Stam W: Bayesian meta-analysis of multiple treatment comparisons: an introduction to mixed treatment comparisons. Value Health. 2008, 11: 956-964. 10.1111/j.1524-4733.2008.00347.x.
Lu G, Ades AE: Combination of direct and indirect evidence in mixed treatment comparisons. Stat Med. 2004, 23: 3105-3124. 10.1002/sim.1875.
Mills EJ, Bansback N, Ghement I, Thorlund K, Kelly S, Puhan MA, Wright J: Multiple treatment comparison meta-analyses: a step forward into complexity. Clin Epidemiol. 2011, 3: 193-202.
Salanti G, Higgins JP, Ades AE, Ioannidis JP: Evaluation of networks of randomized trials. Stat Methods Med Res. 2008, 17: 279-301.
Sutton A, Ades AE, Cooper N, Abrams K: Use of indirect and mixed treatment comparisons for technology assessment. Pharmacoeconomics. 2008, 26: 753-767. 10.2165/00019053-200826090-00006.
Glenny AM, Altman DG, Song F, Sakarovitch C, Deeks JJ, D'Amico R, Bradburn M, Eastwood AJ: Indirect comparisons of competing interventions. Health Technol Assess. 2005, 9: 1-134. iii-iv
Song F, Altman DG, Glenny AM, Deeks JJ: Validity of indirect comparison for estimating efficacy of competing interventions: empirical evidence from published meta-analyses. BMJ. 2003, 326: 472-10.1136/bmj.326.7387.472.
Song F, Loke YK, Walsh T, Glenny AM, Eastwood AJ, Altman DG: Methodological problems in the use of indirect comparisons for evaluating healthcare interventions: survey of published systematic reviews. BMJ. 2009, 338: b1147-10.1136/bmj.b1147.
Song F, Xiong T, Parekh-Bhurke S, Loke YK, Sutton AJ, Eastwood AJ, Holland R, Chen YF, Glenny AM, Deeks JJ, Altman DG: Inconsistency between direct and indirect comparisons of competing interventions: meta-epidemiological study. BMJ. 2011, 343: d4909-10.1136/bmj.d4909.
Higgins JPT, Green S (editors): Cochrane Handbook for Systematic Reviews of Interventions Version 5.0.2 [updated September 2009]. The Cochrane Collaboration. 2009
Mills EJ, Ioannidis JPA, Thorlund K, Schünemann HJ, Puhan MA, Guyatt GH: How to use an article reporting a multiple treatment comparison meta-analysis. JAMA. 2012, 308: 1246-53. 10.1001/2012.jama.11228.
Brok J, Thorlund K, Gluud C, Wetterslev J: Trial sequential analysis reveals insufficient information size and potentially false positive results in many meta-analyses. J Clin Epidemiol. 2008, 61: 763-769. 10.1016/j.jclinepi.2007.10.007.
Brok J, Thorlund K, Wetterslev J, Gluud C: Apparently conclusive meta-analyses may be inconclusive–Trial sequential analysis adjustment of random error risk due to repetitive testing of accumulating data in apparently conclusive neonatal meta-analyses. Int J Epidemiol. 2009, 38: 287-298.
Guyatt GH, Oxman AD, Kunz R, Brozek J, Alonso-Coello P, Rind D, Devereaux PJ, Montori VM, Freyschuss B, Vist G: GRADE guidelines 6. Rating the quality of evidence--imprecision. J Clin Epidemiol. 2011, 64: 1283-1293. 10.1016/j.jclinepi.2011.01.012.
Higgins JP, Whitehead A, Simmonds M: Sequential methods for random-effects meta-analysis. Stat Med. 2011, 30: 903-921. 10.1002/sim.4088.
Ioannidis J, Lau J: Evolution of treatment effects over time: empirical insight from recursive cumulative metaanalyses. Proc Natl Acad Sci USA. 2001, 98: 831-836. 10.1073/pnas.98.3.831.
Pogue J, Yusuf S: Overcoming the limitations of current meta-analysis of randomised controlled trials. Lancet. 1998, 351: 47-52. 10.1016/S0140-6736(97)08461-4.
Pogue JM, Yusuf S: Cumulating evidence from randomized trials: utilizing sequential monitoring boundaries for cumulative meta-analysis. Control Clin Trials. 1997, 18: 580-593. 10.1016/S0197-2456(97)00051-2. discussion 661–586
Thorlund K, Devereaux PJ, Wetterslev J, Guyatt G, Ioannidis JP, Thabane L, Gluud LL, Als-Nielsen B, Gluud C: Can trial sequential monitoring boundaries reduce spurious inferences from meta-analyses?. Int J Epidemiol. 2009, 38: 276-286.
Thorlund K, Imberger G, Walsh M, Chu R, Gluud C, Wetterslev J, Guyatt G, Devereaux PJ, Thabane L: The number of patients and events required to limit the risk of overestimation of intervention effects in meta-analysis–a simulation study. PLoS One. 2011, 6: e25491-10.1371/journal.pone.0025491.
Wetterslev J, Thorlund K, Brok J, Gluud C: Trial sequential analysis may establish when firm evidence is reached in cumulative meta-analysis. J Clin Epidemiol. 2008, 61: 64-75. 10.1016/j.jclinepi.2007.03.013.
Wetterslev J, Thorlund K, Brok J, Gluud C: Estimating required information size by quantifying diversity in random-effects model meta-analyses. BMC Med Res Methodol. 2009, 9: 86-10.1186/1471-2288-9-86.
van der Tweel I, Bollen C: Sequential meta-analysis: an efficient decision-making tool. Clin Trials. 2010, 7: 136-146. 10.1177/1740774509360994.
Thorlund K, Anema A, Mills E: Interpreting meta-analysis according to the adequacy of sample size. An example using isoniazid chemoprophylaxis for tuberculosis in purified protein derivative negative HIV-infected individuals. Clin Epidemiol. 2010, 2: 57-66.
Pereira TV, Ioannidis JP: Statistically significant meta-analyses of clinical trials have modest credibility and inflated effects. J Clin Epidemiol. 2011, 64: 1060-1069. 10.1016/j.jclinepi.2010.12.012.
Thorlund K, Engstrom J, Wetterslev J, Brok J, Imberger G, Gluud C: User manual for trial sequential analysis (TSA). In Book User manual for trial sequential analysis (TSA). 2011, Copenhagen Trial Unit, City
Mills EJ, Ghement I, O'Regan C, Thorlund K: Estimating the power of indirect comparisons: a simulation study. PLoS One. 2011, 6: e16237-10.1371/journal.pone.0016237.
Lumley T: Network meta-analysis for indirect treatment comparisons. Stat Med. 2002, 21: 2313-2324. 10.1002/sim.1201.
Mills EJ, Wu P, Lockhart I, Thorlund K, Puhan MA, Ebbert JO: Comparison of high-dose and combination nicotine replacement therapy, varenicline, and buproprion for smoking cessation: A systematic review and multiple treatment meta-analysis. Ann Med. 2012, 44: 588-97. 10.3109/07853890.2012.705016.
Bucher HC, Guyatt GH, Griffith LE, Walter SD: The results of direct and indirect treatment comparisons in meta-analysis of randomized controlled trials. J Clin Epidemiol. 1997, 50: 683-691. 10.1016/S0895-4356(97)00049-8.
Higgins JP, Whitehead A: Borrowing strength from external trials in a meta-analysis. Stat Med. 1996, 15: 2733-2749. 10.1002/(SICI)1097-0258(19961230)15:24<2733::AID-SIM562>3.0.CO;2-0.
Higgins JP, Thompson SG: Quantifying heterogeneity in a meta-analysis. Stat Med. 2002, 21: 1539-1558. 10.1002/sim.1186.
Borenstein M, Hedges L, Higgins JP: Introduction to meta-analysis. 2009, John Wiley and Sons, Chichester, 25
Thorlund K, Imberger G, Johnston B, Walsh M, Awad T, Thabane L, Gluud C, Devereaux PJ, Wetterslev J: Evolution of heterogeneity (I2) estimates and their 95% confidence intervals in large meta-analyses. PLoS One. 2012, 7: e39471-10.1371/journal.pone.0039471.
Whitehead A, Whitehead J: A general parametric approach to the meta-analysis of randomized clinical trials. Stat Med. 1991, 10: 1665-1677. 10.1002/sim.4780101105.
Thorlund K, Imberger G, Wetterslev J, Brok J, Gluud C: Comments on 'Sequential meta-analysis: an efficient decision-making tool' by I van der Tweel and C Bollen. Clin Trials. 2010, 7: 752-753. 10.1177/1740774510386255. author reply 754
The authors declare that they have no competing interests.
KT conceived the idea, mathematically derived the proposed formulas, drafted the first manuscript and performed all statistical analyses. EM contributed to the methodological development, writing of the manuscript and interpretation of results. Both authors read and approved the final manuscript.
About this article
Cite this article
Thorlund, K., Mills, E.J. Sample size and power considerations in network meta-analysis. Syst Rev 1, 41 (2012). https://doi.org/10.1186/2046-4053-1-41