 Research article
 Open access
 Published:
Trial Sequential Analysis in systematic reviews with metaanalysis
BMC Medical Research Methodology volume 17, Article number: 39 (2017)
Abstract
Background
Most metaanalyses in systematic reviews, including Cochrane ones, do not have sufficient statistical power to detect or refute even large intervention effects. This is why a metaanalysis ought to be regarded as an interim analysis on its way towards a required information size. The results of the metaanalyses should relate the total number of randomised participants to the estimated required metaanalytic information size accounting for statistical diversity. When the number of participants and the corresponding number of trials in a metaanalysis are insufficient, the use of the traditional 95% confidence interval or the 5% statistical significance threshold will lead to too many false positive conclusions (type I errors) and too many false negative conclusions (type II errors).
Methods
We developed a methodology for interpreting metaanalysis results, using generally accepted, valid evidence on how to adjust thresholds for significance in randomised clinical trials when the required sample size has not been reached.
Results
The LanDeMets trial sequential monitoring boundaries in Trial Sequential Analysis offer adjusted confidence intervals and restricted thresholds for statistical significance when the diversityadjusted required information size and the corresponding number of required trials for the metaanalysis have not been reached. Trial Sequential Analysis provides a frequentistic approach to control both type I and type II errors. We define the required information size and the corresponding number of required trials in a metaanalysis and the diversity (D^{2}) measure of heterogeneity. We explain the reasons for using Trial Sequential Analysis of metaanalysis when the actual information size fails to reach the required information size. We present examples drawn from traditional metaanalyses using unadjusted naïve 95% confidence intervals and 5% thresholds for statistical significance. Spurious conclusions in systematic reviews with traditional metaanalyses can be reduced using Trial Sequential Analysis. Several empirical studies have demonstrated that the Trial Sequential Analysis provides better control of type I errors and of type II errors than the traditional naïve metaanalysis.
Conclusions
Trial Sequential Analysis represents analysis of metaanalytic data, with transparent assumptions, and better control of type I and type II errors than the traditional metaanalysis using naïve unadjusted confidence intervals.
Background
Most metaanalyses include too few randomised participants, to obtain sufficient statistical power that allow reliable assessment of even large anticipated intervention effects [1]. The credibility of statistical significant metaanalyses with too few participants is poor, and intervention effects are often spuriously overestimated (type I errors) or spuriously underestimated (type II errors) [2]. Metaanalyses of, e.g., cardiovascular, anaesthesiologic, and neonatal interventions have many false positive and false negative results, due to low statistical power in a metaanalysis when the required number of randomised participants or trials have not been reached [3–6]. Trial Sequential Analysis (TSA) of a metaanalysis may amend these problems [4, 7]. In this article, we aim to describe the origin, history, adaptation, and criticism of TSA.
Using TSA, we can handle a metaanalysis of several randomised clinical trials in an analogous manner to interim analysis of a single randomised clinical trial. If the accrued cumulative information fails to achieve the required number of randomised participants in order to detect or reject a specific assumed effect, the uncertainty of the estimate of the intervention effect will increase. The uncertainty will decrease the higher the fraction of the required information size the metaanalysis obtain. To statistically solve the problem with uncertainty, we expand the confidence interval, i.e., adjusting the threshold for statistical significance when the required information size has not been reached. The farther from the required number of randomised participants, the wider the confidence interval and the lower the statistical significance level needs to be in order to reliably assess the uncertainty of the point estimate.
In TSA of a metaanalysis, we include the trials in chronological order and we handle the analysis of these trials as an interim analysis relative to the required number of randomised participants. TSA calculates the required number of participants, based on our predefined anticipated intervention effect, i.e., our alternative hypothesis [7–9]. The result of a trial sequential metaanalysis is displayed on a TSA diagram (e.g., Fig. 1a and b) with a TSAadjusted confidence interval and an adjusted level of statistical significance, i.e, a lower threshold for statistical significance compared to the usual nominal of 0.05, if the required information size has not been reached [10].
In the abovementioned adjustments, we take into consideration if the required number of randomised participants and corresponding trials, to show or reject a specific intervention effect, were reached or not. The required information size is defined as the number of participants and events necessary to detect or reject an a priori assumed intervention effect in a metaanalysis [11]. The required information size is not a single sample size, but a summation of sample sizes from a given number of included trials. Therefore, the calculation is performed considering the variability (heterogeneity variance) between the estimates of the intervention effects of the included trials.
In TSA, the sample size, required for a single randomised clinical trial to be conclusive for a specific intervention effect, is adjusted upward by an appropriate measure of the statistical heterogeneity in the metaanalysis in order to become the required information size. This is equivalent to using the variance in the randomeffects model to calculate the required information size (the model variance based calculation of the required information size). In the TSA, we hereafter adjust the confidence interval of the point estimate and the threshold for statistical significance relative to the fraction of the required information size which has been accrued in the actual metaanalysis [11].
First, we will present a motivating example of a metaanalysis on hypothermia versus no hypothermia in comatose patients having survived cardiac arrest. Second, we present an updated metaanalysis with the results of a new trial, and we describe how this update has changed the conclusion of the preceding traditional metaanalysis. We also show how the use of TSA would appropriately have reduced the risk of a wrong conclusion in the first metaanalysis failing to achieve the required information size. Third, we shortly describe the historical development of sequential analyses in a single trial with interim analyses and in a cumulative metaanalysis of several trials. We explain how sequential metaanalysis can be performed with TSA [12]. Finally, we discuss the criticism that has been raised about TSA and we briefly describe the possibility for Bayesian metaanalysis as an alternative to both traditional naïve metaanalysis and TSA of a metaanalysis.
A motivating example: how the Target Temperature ManagementTrial changed the conclusion of the metaanalysis of trials with cooling of patients after out of hospital cardiac arrest
In TSA, we consider each interimanalysis result, produced after the addition of a new trial, a sequential metaanalysis. The possibility to include groups of several new trials at a time is, of course, also possible. This latter approach will decrease the number of interimanalyses in the cumulative metaanalysis [10]. However, updating the metaanalysis in a systematic review each time a new trial is published is a rational decision, and to update a systematic review before a new trial is initiated ought to become mandatory [13–15]. Previous trial results ought to be considered whenever we evaluate the cons and pros of designing new trials, as the evidence on a given intervention may already be sufficient [13–15]. It is surprising to see how little the TSA, conducted after each new trial has been interimanalysed, differs from the last TSA on groups of trials (e.g., TSA only updated every second year).
Figure 1 shows the result of a TSA of metaanalysis of four trials comparing a target temperature of 33°–34 °C versus no cooling, conducted before the initiation of the Target Temperature Management (TTM) Trial (Fig. 1a) [16–18]. The TSA shows that the four trials did not even reach half of the required information size to confirm or reject a 17% relative risk reduction which was the intervention effect indicated in a conventional metaanalysis of the trials [16]. The conventional confidence interval for the relative risk ratio of allcause mortality in a traditional metaanalysis is 0.70 to 1.00 (P = 0.05), suggesting a reduction of mortality. The confidence interval and the Pvalue would not have been sufficient to claim a conclusive interim analysis stopping for benefit in a single randomised trial if analysed with LanDeMets’ group sequential monitoring boundaries [19]. For demonstrating a 17% relative risk reduction, the TSAadjusted confidence interval of the relative risk is 0.63 to 1.12. This confidence interval shows that i) a target temperature of 33°–34 °C versus no cooling can either decrease or increase mortality, and ii) that definitive evidence has not yet been reached. The cumulative Zcurve in the figure does not pass through the trial sequential monitoring boundary for benefit; only the conventional and naïve P = 0.05 (Z = 1.96) level for a beneficial effect has been reached. Therefore, there is not sufficient information to document the effect, or there may not be a beneficial effect at all. Nevertheless, based on this evidence, international guidelines had recommended for ten years that the temperature of comatose cardiac arrest patients should be targeted to 33°–34 °C, calling the intervention »mild therapeutic hypothermia« [20]. No further randomised clinical trials of induced hypothermia versus no temperature control (or normothermia) in comatose cardiac arrest patients after resuscitation and admittance to intensive care units were conducted during this 10year period. This may indicate that a Pvalue of 0.05 in the conventional metaanalysis was used as an unofficial »stopping boundary« for further trials within this same period.
In the TTM Trial, we compared the effect of cooling to target temperature 33 °C versus 36 °C on mortality of cardiac arrest patients [17, 18]. The updated TSA including the TTM Trial showed no statistically significant effect at the conventional level, as the Zcurve returned to the area with P > 0.05 (Z < 1.96) (Fig. 1b). Figure 1b shows that the cumulative Zcurve touches the futilityboundaries in the TSA diagram (see section ‘False Negative Metaanalyses’ below). Therefore, the updated TSA indicates that a 17% relative risk reduction, or an even greater reduction, most likely can be rejected, although the preestimated required information size of 2040 patients has not yet been reached. It is not likely that a metaanalysis will ever show a 17% statistical significant relative risk reduction of mortality, even though the continued conduct of trials until a cumulated number of patients, corresponding to the required metaanalytic information size of 2040 patients, was reached (Fig. 1b). The conclusion is that hypothermia to 33°–34 °C does not seem to have a clinical important effect on mortality compared with no cooling or targeted normothermia (36 °C), as the 17% relative risk reduction only corresponds to a median of 3 weeks’ longer survival [17, 18]. Moreover, the original conventional metaanalysis before inclusion of the TTM Trial had a false positive result; the null hypothesis was falsely rejected. Whether the avoidance of fever is actually beneficial compared with no cooling at all, remains to be tested, as the TTM trial used cooling in both the intervention (target 33 °C) and the control group (target 36 °C).
Interimanalyses during a randomised clinical trial with an accumulating number of participants
If a trial is stopped after an interimanalysis because of a P < 0.05 or the trial is continued if P ≥ 0.05, the real risk of committing a type I error will increase to more than 0.05 with the number of interimanalyses. Introducing an interimanalysis halfway in a randomised clinical trial, using a stopping Pvalue equal to 0.05 in both the halfway analysis and the final analysis, will increase the real maximal type I error risk to 8% [21, 22] (Table 1). If the procedure of interim analysis is performed as four interim analyses and one final analysis, with a constant level of statistical significance of 5%, the real type I error risk will be 14% [21]. A simulation study using repetitive testing on an accumulating number of participants in a single trial, has shown that the Pvalue will inevitably become less than 0.05, despite the true intervention effect being zero [23].
A Bonferroni adjustment of the level of statistical significance, being 5% divided with the number of tests on accumulating data, assumes that all tests are conducted on independent data. As the tests on the accumulating trial population are not statistically independent, the Bonferroniadjusted levels of statistical significance are most often too conservative [24]. The trial participants in an early sequential analysis are also included in the subsequent later sequential analyses. Therefore, there is an increasing overlap of trial participants included in the latest sequential analysis compared to participants included in the previous sequential analyses. The closer we come to the a priori calculated sample size, the Bonferroni adjustment becomes more and more unjustified (too conservative).
Historical development of sequential analyses in a single trial with interim analyses
Methods to avoid an increased risk of a type I error due to repetitive testing on an increasing number of observations was described by Abraham Wald in 1945 in Contributions to the theory of statistical estimation and testing hypotheses [25]. Wald proposed »the sequential probability ratio test« in which the sequential testing continues until a definitive wanted or unwanted effect can be proved [26, 27]. According to this procedure, the trial continues as long as the results of the sequential tests fall within the socalled ‘zone of indifference’ amidst the two alternative hypotheses. This procedure, used as a quality assurance measure of production during the Second World War, has never achieved wide implementation in randomised clinical trials; possibly because the procedure is bound to continue infinitely as long as the true intervention effect lies between the two alternative hypotheses. Consequently, a decision to stop the trial may never become possible [28].
After the Second World War, Peter Armitage suggested more restrictive levels of statistical significance than 5% to stop a trial before the a priori calculated sample size was reached [21]. This procedure was applied in a number of interim analyses of large trials [29]. Furthermore, Stuart Pocock proposed a procedure in which the overall risk of type I error is limited to 5% by setting the statistical significance level to 0.05 divided by k, using k1 interim analyses and a final analysis [22]. This procedure is identical to the Bonferroni procedure for interim analyses and a final analysis of a single trial [30]. Researchers might find it peculiar to only declare statistical significance if P < (0.05/k), despite the estimated sample size has been reached and the required criterion for statistical independence was not fulfilled.
In 1977, Richard Peto suggested the use of a maximal type I error risk (αspending) in each of four interim analyses of 0.001 (1 promille) and 0.05 in the final analysis. As this would produce a summary additional type I error risk of 0.004 to the final 0.05, the total type I error risk would maximally be 5.4% [31] (Fig. 2). However, by a modest increase of the a priori estimated sample size, the summary maximal used type I error risk would remain within the usual 5%. As shown above, the rationale of statistical independence required for this procedure still lacks underlying reason as to why the trial participants in an early sequential analysis are also included in the subsequent sequential analysis.
In 1979, Peter O’Brien and Thomas Fleming proposed the group sequential design of trials with interim analyses, using exponential decreasing levels of statistical significance with the increasing number of patients in the sequentially analysed groups (Fig. 2) [32]. The recommendations of the International Conference on Harmonization – Good Clinical Practice, the U.S.A. Food and Drug Administration, and the European Medicines Agency on the design and analysis of randomised trials with interim analyses are mainly based on works from 1980s, primarily prepared by Gordon Lan, Kuyung Man Kim, and David DeMets (Fig. 2) [18, 33, 34]. Their works allow proper sequential testing at any time during the trial period, without unduly increasing the overall risk of a preset nominal type I error risk [34–36].
Methods
Avoiding the increased risk of random errors in randomised clinical trials with interim analyses
It is and should be mandatory to perform interim analyses in large randomised clinical trials addressing patientcentred outcomes. Even though the preplanned sample size has not been reached, thousands of patients might already have been randomised in a trial. Before we allow the trial to continue, there is a need to secure that no valid evidence showing superiority of one of the compared interventions exists. If one of the interventions (could also be placebo) with a sufficiently small uncertainty is superior to the other one in an interim analysis, it may be unethical to continue the trial. The explanation for this is that the superiority can be so large that it cannot be reversed even though we continue to randomise patients until the total, originally preplanned sample size is obtained. If the trial is continued despite the superiority of the intervention in one of the intervention groups, the patients in the other group will be exposed to an inferior (harmful) intervention and the trial must be stopped [37]. The use of interim analyses in a single randomised trial has to be planned at the design stage of the trial and protocolised upfront as group sequential analyses in the charter for interim analyses [33]. For the conduct of group sequential analyses, a sample size is calculated already at the design stage, based on the anticipation of a minimal important and realistic intervention effect of the primary outcome of the trial [36, 38] (see Appendix).
The sample size calculation considers the level of statistical significance at which we want to test a dichotomous or a continuous outcome when the full sample size has been reached. It is when the precalculated sample size has been reached, and only then, a twosided Pvalue of less than 0.05, corresponding to a teststatistic Zvalue of ±1.96, can be accepted as the statistical significance level when α has been set to 5% in the sample size calculation.
Interim analyses, with the potential to stop a randomised trial before the estimated (or fixed) sample size has been reached due to a positive, negative, or lack of the addressed effect, can be conducted for dichotomous and continuous outcomes by calculating the cumulative Z_{ i }value at the ith analysis (see Appendix). The calculated Z_{ i }value is then related to the more restrictive level of statistical significance, the critical Zvalue being the discrete group sequential boundary according to the actual accrued number of participants.
There is international consensus that the increase of type I error risk with sequential testing, including the risk of overestimating the intervention effect or underestimating the variance, at an interim analysis, should be outweighed by more restrictive levels of statistical significance before the a priori estimated (fixed) sample size has been reached [29, 31–37]. This is why ‘monitoring boundaries’, with significance levels much smaller than a nominal Pvalue of 0.05 (corresponding to much larger Zvalues than ±1.96) are applied as criteria to stop a trial before achieving the estimated sample size [33].
Numerical integration is used to calculate the monitoring boundaries, being the critical levels of statistical significance for the Z_{ i }values (and Pvalues) of the interim analyses [39]. Most often, the O’BrienFleming’s αspending–function is applied and converted to sequential boundaries (critical values) for the Z_{ i }values called LanDeMets’ sequential monitoring boundaries (Fig. 2) [18, 19]. The αspending function allows only a small part of the total nominal type I error risk to be used initially in the sequential analyses, and with a modest increase of the estimated final (fixed) sample size, there is a full 5% type I error risk available for the final analysis when the a priori estimated sample size is reached. LanDeMets’ sequential boundaries allow testing whenever you want during the trial [34, 35]. If we plan, e.g., a halfway analysis in a randomised trial, we can monitor the Pvalue at this time point according to LanDeMets’ monitoring boundaries and suggest that the trial is stopped if the Pvalue is less than 0.003 which corresponds to a 99.7% confidence interval excluding 1.00 for a relative risk or 0.00 for a mean difference [34–36]. Therefore, sequential analyses become a theoretical decision tool to decide whether a trial should be stopped before the estimated (fixed) sample size is achieved, considering the sparse data and the repetitive testing during the trial [37].
Avoiding the increased risk of random errors in cumulative metaanalyses with sparse data and multiple metaanalytic updates
The majority of metaanalyses include less than the required number of randomised participants and trials in order to become conclusive [1, 3, 5, 7]. There are two reasons for this. First, most randomised trials are underpowered [1, 3, 5, 7]. Second, the estimation of the required information size in a randomeffects metaanalysis ought to incorporate the heterogeneity variance (between trial variance) [1, 7, 11]. Only 22% of the metaanalyses in The Cochrane Library have 80% power to conclude whether there is an intervention effect of 30% or not when the usual maximal risks of type I error (α) of 5% and type II error (β) of 20% are applied [1]. This lack of power is primarily caused by small trials and a considerable heterogeneity variance between the estimates of the intervention effect in the included trials [1].
Metaanalyses can be conducted with a fixedeffect model or a randomeffectsmodel [40, 41]. In the fixedeffect model, we assume one true underlying effect in all the included trials. In the randomeffects model, we assume that the true underlying effects vary from trial to trial according to a normal or log normal distribution. Often, the fixedeffect assumption is unrealistic as the possible underlying effect may depend on, e.g., doses of a pharmacological intervention, duration of the interventions, timing of the outcome assessment, and differences between the trial populations. These differences between the included trials are called clinical heterogeneity. Due to these factors and possibly random variation, the included effect estimates often show considerable variation defined as statistical heterogeneity and measured as large inconsistency (I^{2}) [42] and large diversity (D^{2}) [11]. Considerable statistical heterogeneity leads to increased uncertainty, expressed as a wider confidence interval of the intervention effect when the metaanalytic estimate is calculated in a randomeffects model. Early metaanalyses conducted before the required information size and the corresponding number of trials are achieved [43], often wrongly show unrealistic large intervention effects as well as statistical significance which cannot be reproduced when the amount of required information is adequately considered [44, 45]. The reliability in early metaanalyses is lower compared to their updated counterparts years later [2]; the estimated intervention effects, when further trials are included in the metaanalysis update, become considerably lower than previously estimated [2].
A large simulation study of randomeffects metaanalyses shows that there is a considerable risk of overestimating the intervention effect when the required information size has not been reached [6]. These results were based on the assumption that the ‘true’ intervention effect was zero while the frequencies of events in the control groups and the heterogeneity variance were assumed similar to those in large cardiologic metaanalyses [6]. It has been shown empirically that approximately 25% of cardiologic metaanalyses are inconclusive because of lack of power [3]. Turner and colleagues showed that the trials and the metaanalyses of Cochrane systematic reviews have limited power [1]. In Cochrane metaanalyses, each total number of analysed participants provide only 22% of the metaanalyses with an 80% power to detect or refute a 30% relative risk reduction (which is a large intervention effect) [1] (Fig. 3). Recently, Imberger and colleaques confirmed these results in metaanalyses of anaesthesiological interventions [46]. Accordingly, four out of five metaanalyses did not have the statistical power to address even substantial intervention effects. The number of metaanalyses with sufficient power to address smaller and clinically more plausible intervention effects are, undoubtedly, even smaller.
If we test with a constant level of statistical significance (e.g., 5%) on the way towards the required information size, the risk of type I error is increased to more than 5%. The problem for cumulative metaanalyses, due to repeated updating and consecutive calculation of 95% confidence intervals, with inclusion of results from new randomised trials is, therefore, analogous to interim analyses of a single trial [8, 9]. Thus, we, as well as others, recommend that the interpretation of metaanalyses in systematic reviews is done alongside with a sequential analysis, e.g., Trial Sequential Analysis (TSA) [46, 47]. The purpose of using TSA is to avoid the risk of type I and type II errors due to sequential testing on a constant statistical significance level and with inclusion of fewer participants than the required number in order to detect or reject specified effects [7, 10, 11]. It is possible to accommodate Gordon Lan and David DeMets’ group sequential analysis for interim analysis in a single randomised trial to the updating of cumulative metaanalysis as it progresses with the addition of trials. This is done with an appropriate continuous use of type I error risk and an αspending function of the allowed total nominal type I error risk, so that when the required information size and the required number of trials have been reached and beyond, the risk is kept below 5%. The trial sequential monitoring boundaries generated this way make it possible to test if significance is reached and to adjust the confidence intervals every time new trials are added to the metaanalysis. The latter is a prerequisite for using sequential boundaries in cumulative metaanalyses of trials with varying sample sizes [10, 12].
Besides applying the observed estimate of statistical heterogeneity—the observed statistical diversity (D^{2}) [11, 41] in the most recently conducted metaanalysis—it may be reasonable to apply an expected heterogeneity in the calculation of the required information size, especially when the observed heterogeneity is zero [48]. As it is unlikely that diversity will stay zero when larger trials are added, an expected heterogeneity may be used in a sensitivity analysis (e.g., a diversity of 25% or the upper confidence interval of the I^{2} (provided by the TSA program)) when the required information size is calculated [48, 49]. It may also be wise in a post hoc calculation of the required information size to apply the least likely intervention effect, i.e., the confidence limit of the summary estimate in the metaanalysis confidence interval closest to the null effect. The latter is a conservative approach facilitating the evaluation of whether a metaanalysis may show an effect of the least likely magnitude in a TSA. If a TSA with such an approach shows a statistical significant intervention effect, judged by the TSAadjusted confidence interval, there is a very high probability that the intervention has an effect, provided that the included trials are at low risk of bias. In contrast, there will only be very low evidence of effect if the TSAadjusted confidence interval does not exclude the null effect for an intervention effect of a magnitude indicated by the point estimate.
Results
False positive metaanalyses
It is necessary to assume or address a specific magnitude of the intervention effect, different from zero, in order to calculate the sample size in a single trial. Therefore, when a sample size is estimated, we relate not only to the null hypothesis but also to a specific alternative hypothesis. The alternative hypothesis is the assumption or the anticipation of a specific magnitude of the intervention effect different from zero. Most often randomeffects metaanalysis will be the preferred appropriate method to estimate the precision weighted average effect as it does not ignore the statistical heterogeneity variance. If statistical heterogeneity is anticipated, the information size in the conclusive metaanalysis ought to be an upward adjusted sample size of a corresponding adequately powered single trial. The upward adjustment is done with the variance expansion shifting from a ‘fixedeffect’ model to a ‘randomeffects’ model, see Appendix [11].
The described example from cooling of patients after out of hospital cardiac arrest is far from being unique (Fig. 1). Among metaanalyses of interventions for neonatal patients, there were approximately 25% to 30% false positive results [5, 50]. In 2009, we showed empirically that the use of LanDeMets’ trial sequential monitoring boundaries eliminated 25% of the false positive traditional interimmetaanalyses. This analysis included 33 final metaanalyses with sufficient information size to detect or reject a 15% relative risk reduction [44]. In 2013, we showed that 17% of cardiovascular metaanalyses with P < 0.05 were most likely false positive [3]. In 2015, we showed that less than 12% of metaanalyses of anaesthesiological interventions had 80% power to show a 20% relative risk reduction [46].
There may be other important reasons for a traditional metaanalysis to yield a false positive result than only the increased risk of random errors. A risk of systematic error (bias) in the included trials is a frequent cause of overestimation of benefit and underestimation of harm – sequential metaanalyses do not in any way solve problems with bias [51–58]. Therefore, it is recommended that every single trial included in a systematic review with metaanalysis be evaluated for risks of bias. This evaluation should encompass the following domains: generation of the allocation sequence, allocation concealment, blinding of patients and caregivers, blinding of outcome assessment, report on attrition during the trial, report on outcomes, and industry funding. Other types of bias may also need to be considered [51–58].
False negative metaanalyses
Lack of a statistical significant intervention effect in a traditional metaanalysis is not necessarily evidence of no effect of the intervention. »Absence of evidence is not evidence of absence of effect« [59]. Nevertheless, sequential metaanalyses with the TSA software may show that the metaanalysis has sufficient statistical power to reject an intervention effect of a specific magnitude even though the estimated required information size has not yet been reached (Fig. 4).
This can be done by calculating the nonsuperior and noninferior trial sequential monitoring boundaries, the socalled ‘futility boundaries’. Futility boundaries indicate when the assumed effect could be considered unachievable. Futilityboundaries are calculated using a power function analogous to the αspending function for constructing superiority and inferiorityboundaries with application of numerical integration [36]. The example with cooling of comatose patients after cardiac arrest shows a situation where the assumed intervention effect of 17% relative risk reduction can be rejected because the Zcurve crosses the futilityboundary (Fig. 1b). However, this is not always what happens. We found that in 25 of 56 (45%) published cardiovascular systematic reviews in The Cochrane Library, the actual accrued information size failed to reach what was required to refute a 25% relative risk reduction [3]. Only 12 of these reviews (48%) were recognised as inconclusive by the authors. Of the 33 metaanalyses not showing statistical significance, only 12 (36%) were truly negative in the sense that they were able to reject a 25% relative risk reduction [3]. This illustrates that the statistical power is also low in many cardiovascular metaanalyses, and false conclusions are imminent. Within other medical specialities, the problems are likely to be even bigger as trials and metaanalyses usually include less patients. Nevertheless, sequential metaanalyses with calculated futilityboundaries may, in some instances, contribute to adequately declare the a priori anticipated intervention effect to be unachievable, though the required information size was not reached [10].
Analogous to the false positive metaanalyses, a metaanalysis may result in a false negative result due to bias. Bias is a frequent cause for underestimation of harmful intervention effects [51–57], and therefore, the preliminary defined bias risk domains should also be evaluated for all included trials when it comes to serious and nonserious adverse events [51–58].
Discussion
We have explained and shown how the use of TSA may assist the metaanalyst in controlling risks of type I and type II errors when conducting metaanalyses. The use of TSA has now increasingly been advocated by authors, both inside and outside The Cochrane Collaboration [47, 60, 61]. However, the use of TSA is not easy, may be misused, and has been critisised [62].
If TSA is designed and conducted after data were collected, there is a danger that the analysis becomes data driven and that it may not be sufficiently stringent to address a predefined alternative hypothesis [63–65]. However, using datadriven hypotheses and analyses is a critique that could potentially be directed against all metaanalyses. This is why, for each TSA, the anticipated intervention effect, the anticipated between trial heterogeneity, and the proportion of the outcome in the control group, should be part of a peer reviewed protocol, published prior to the conduct of the systematic review and the TSA [49, 64, 65]. These considerations should also impact the choice of the metaanalytic model, e.g., whether to give most credibility to the fixedeffect or the randomeffects model and how to calculate the required information size [11, 65].
TSA has been criticised for transferring a method from a decision theoretical universe in a single randomised clinical trial into a universe where the result does not directly impact the subsequent decisions [63–66]. The postulate seems to be that no matter that a TSA shows benefit, harm, or lack of a relevant effect, it will not impact any part of the already finalised trials, and possibly, not decisions to stop or continue ongoing trials, or to initiate trials. This point of view seems to unduly emphasise the difference between the consequences of an interimanalysis in a single trial and the consequences of a sequential metaanalysis of several trials. First, the systematic review is placed at the top of the generally recognised hierarchy of evidence, meaning that the systematic review is considered the most likely reliable source of evidence, implicating whether an intervention should be implemented in clinical practice or further trials should be launched [52, 53]. Interventions are often recommended in clinical guidelines and implemented in clinical practice when a metaanalysis shows statistical significance on the traditional naïve level (P < 0.05) [16, 18, 67–69]. Furthermore, the chance that a metaanalysis is updated in The Cochrane Library is apparently 57% higher when P ≥ 0.05 than when P < 0.05 [4, 45]. This indicates that metaanalyses with P < 0.05 contribute to the decision to stop doing further trials or to decide if metaanalyses should be updated or not.
Critics of sequential metaanalysis have stressed that the method emphasises too heavily the result of the statistical significance test instead of the 95% confidence interval [70]. However, the fundamental problem is not whether the result is presented as a Pvalue or as a confidence interval, but it is foremost because a (1α)% confidence interval is based upon the choice of the maximally allowed type I error risk (α). If we use naïve unadjusted confidence intervals when the required information size is still not reached, we will be led to make hasty and false declarations of statistical significant effects, likely to be refuted if further trials are added. With TSA we adjust the confidence interval for the incomplete metaanalytic information size and for multiple testing [4]. It has been claimed that a traditional 95% confidence interval is sufficient to evaluate whether the intervention works or not [70], but the traditional 95% confidence interval exclusively relates to the null hypothesis and not to a relevant alternative hypothesis [68, 71]. Thereby, the supporters of the traditional confidence interval forget that the rejection of the null hypothesis (the conventional 95% confidence interval excluding the null effect), does not in itself lead to the acceptance of a relevant alternative hypothesis [71]. Premature rejection of the null hypothesis, in the case of sparse data, may be dismissed if these data become sufficient to conclude on a specific alternative intervention effect that is different from the null hypothesis.
A traditional unadjusted 95% confidence interval excluding the null effect and accepting an effect indicated by, e.g., the point estimate, is sufficient as a criterion for statistical significance only when the required information size has been reached. If the number of randomised participants in the metaanalysed trials is less than the required, the confidence interval needs to be adjusted [34, 36]. By exclusively applying a 95% confidence interval in a metaanalysis, one does not automatically account for the lack of required power in the metaanalysis to conclude on an effect size indicated by, e.g., the point estimate [71]. Therefore, in relation to a relevant and realistic alternative hypothesis, the traditional unadjusted confidence interval will represent a too narrow confidence interval which by chance does not include the null effect, and accordingly, the observed effect of the intervention may be misleading [71, 72]. The credibility of the traditional confidence interval relies on the fact that the required information size for a specific effect has been achieved, and thereby, the ability to conclude on an alternative hypothesis [59, 63–65].
TSA has also been criticised for being a too conservative approach as one may decide to use a too sceptical a priori intervention effect and use the total variance in the randomeffects metaanalysis to calculate the required information size. The use of an a priori intervention effect does not consider the intervention effect estimated from the data already accrued; however, applying such an approach may in fact lead to even larger required information sizes [73]. Moreover, to think of the total variance in the randomeffects model as a result of random variation alone, could be seen as a ‘worstcase scenario’ of risk of random error [73]. However, we may rarely know when a variation is caused by systematic differences or by random variations [52]. Therefore, it seems mandatory to perform an analysis, assuming that all the variance encountered in the randomeffects metaanalysis is arising from ‘play of chance’ [46, 47].
Elena Kulinskaya and John Wood [43] argued, in their important article from 2013, that when estimating the information size in randomeffects model metaanalyses, it is too simplistic to just increase the required information size with the variance increase, going from a fixedeffect to a randomeffects model. Kulinskaya and Wood [43] persuasively showed that the necessary number of future trials to be included should be given with a lower limit (i.e., minimal number), regardless of the sample sizes of the trials, before the power of the randomeffects model metaanalysis becomes sufficient to detect or reject a prespecified clinically relevant intervention effect. Kulinskaya and Wood also showed that increasing the number of future trials in a randomeffects model metaanalysis might decrease the required information size estimated for additional future trials to render sufficient power of the randomeffects metaanalysis [43]. We welcome the proposals for modifying the plan on number of subsequently included trials and their sample size. These considerations are in line with the findings of Joanna in’t Hout et al. [74], Alexander Sutton et al. [73], Jeffrey Valentine et al. [75], and Michael Borenstein et al. [76]. However, we would still argue that the difference between the required information size and the accrued information, may contribute importantly to the estimation of the necessary sample size in future trials, especially if coupled with the considerations proposed by Kulinskaya and Wood [43]. If we use the weighted estimate of the variance in previous trials as being the best estimate of the variance for the future trials, we may need 50% (Appendix) more trials than the minimal number required to cover the information gap of the required minus the acquired information size (RISAIS) (Appendix). Following an example given by Kulinskaya and Wood [43], we will be able to cover the information gap suggested by RISAIS with 12 trials instead of the minimal required number of eight trials. As outlined by Kulinskaya and Wood, we would be able to further decrease the total number of future randomised patients by increasing the number of future planned trials even more. However, this will be at the expense of dramatically decreasing the power of each new trial to detect the difference, indicated so far by the point estimate of the metaanalysis (or even the minimal important difference). Certainly, we could choose to increase the number of future trials with only one or two. However, the corresponding information size will still be huge. The minimal required number of trials calculated as the first integer greater than c ⋅ τ ^{2} (where c is a figure relating to the information already gained and τ ^{2} is the between trial variance, Appendix), and the corresponding metaanalytic information size, may be optimal because it provides each of the new, equally sized, trials with the same power as the ‘planned’ randomeffects metaanalysis aimed to detect or reject a similar intervention effect. However, for most interventions, these huge trials will be unrealistically large to conduct. Alternatively, increasing the number of trials corresponding to a required extra information size of RISAIS will still provide such trials with a power of 80% to detect or reject an intervention effect of 2.5 times the effect indicated in the metaanalysis. Increasing the number of trials even further than the number corresponding to RISAIS will decrease the power of these trials with approximately 10% per additional trial (or increase the detectable alternative to three times or more the effect indicated in the metaanalysis). Such trials will subsequently be substantially underpowered to detect or reject even much larger intervention effects than the realistic difference, or even the minimal important difference. This will obviously destroy the integrity of such small future trials and they will generally, and rightfully so, be disregarded as heavily influenced by random error (‘play of chance’). Therefore, the RIS and thereby the RISAIS seem to be a fair tradeoff between the number of required additional randomised participants and the number of required additional trials. In two examples given by Kulinskaya and Wood, the number of additional randomised participants is reduced from 4700 to 720 and from 11,200,000 to 300,000 when using RISAIS at the expense of four more trials than the minimal number of trials required. However, we agree, that a reasonable strategy for resolving the question of the presence or absence of a specific intervention effect with an adequately powered randomeffects model may include a first trial with a sample size equal to the sample size indicated by formula 1 in the Appendix. This is a sample size corresponding to the minimal number of required trials. Such a trial may very well be substantially larger than the total acquired information size in the metaanalysis conducted before the trial. When the result from such a trial becomes available, the updated cumulative metaanalysis using the a priori anticipated intervention effect and a new estimate of the between trial variance may be used in a fixedeffect or a randomeffects model to evaluate how far we will be from a conclusion of whether the intervention effect exists or not. The fixedeffect model may then turn out to be the most appropriate model to evaluate the pooled intervention effect when one or a few trials heavily dominate the entire accumulated evidence [77].
Nevertheless, we must be aware that including new trials in a cumulative metaanalysis may change the estimate of the ‘between trials variance’ as well as the proportion of events in the control group which are both essential for estimating the required information size and the corresponding number of required future trials. If diversity and the proportion of events in the control group change substantially, the magnitude of the required information size and the corresponding number of required future trials may change accordingly. This is the phenomenon of the ‘moving target’ which critics hold against TSA. However, a moving target seems better than having no target at all. Recently, we documented that in prospective application of TSA in very large cumulative metaanalyses, TSA prevented false positive conclusions in 13 out of 14 metaanalyses when RIS was not reached [45].
Trial Sequential Analysis: a position between frequentist and Bayesian thinking
TSA of metaanalysis like the sequential analysis of a single randomised trial, originates from frequentist statistics [29]. The frequentist way of thinking was initially based on testing of the null hypothesis. This applies to both the Pvalue and its relation to an a priori accepted maximal type I error risk (α) and the possibility of including a null effect in the corresponding (1α)% confidence interval [29]. The anticipation of an intervention effect of a specific magnitude, the alternative hypothesis, and subsequently the calculation of a required information size enabling the conclusion whether such an effect could be accepted or rejected, is, however, intimately related to the Bayesian prior.
TSA contains an element of Bayesian thinking by relating the result of a metaanalysis to the a priori point estimate of the intervention effect addressed in the analysis [77]. Bayes’ factor (BF) for a trial result is the ratio between the probability that the trial data originates under the null hypothesis, and the probability that the trial data originates under the alternative hypothesis or even several alternative hypotheses [72, 78, 79]. The posterior odds ratio for the estimate of the intervention effect after a new trial is added is calculated given the prior odds ratio for the intervention effect before the trial as: posterior odds ratio = BF x prior odds ratio [79]. In a Bayesian analysis, the prior takes form of an anticipated probability distribution of one or more possible alternative hypotheses or intervention effects which multiplied with the likelihood of the trial, results in a posterior distribution [79].
A methodological position between the frequentist and the Bayesian thinking can be perceived both in sequential interimanalyses of a single trial and in TSA of several trials [29]. Both have a decisive anticipation of a realistic intervention effect, although a full Bayesian analysis should incorporate multiple prior distributions with different anticipated distributions of intervention effects: e.g., a sceptical, a realistic, and an optimistic prior [79]. The TSA prioritise one or a few specific alternative hypotheses, specified by point estimates of the anticipated effect in the calculation of the required information size just as in the sample size estimation of a single trial [11].
The incentive to use sequential analyses arise because the true effect is not known and the observed intervention effect may be larger than the effect addressed in the sample size estimation of a single trial as well as in the estimation of the required information size for a metaanalysis of several trials. The need to discover an early, but greater effect than the one anticipated in the sample or information size calculation, or to discard it, thereby originates. If the intervention effect, in relation to its variance, happens to be much larger during the trial or the cumulative metaanalysis, this will be discovered through the breakthrough of the sequential boundary. However, this may also be problematic as too small sample sizes (in relation to the true effect), as mentioned, increase the risk of overestimation of the intervention effect or the risk of underestimation of the variance. In other words, due to a factitious too small sample size, we may erroneously confirm an unrealistic large anticipated intervention effect due to the play of chance.
There seems to be an ancestry between the sceptical prior in a Bayesian analysis and the use of a realistic intervention effect in a sequential analysis when the sample size in a single trial or the information size in a metaanalysis should be calculated [77, 78]. The smaller the effect, the greater the demand for quantity of information, and the sequential statistical significance boundaries become more restrictive. In other words, it becomes more difficult to declare an intervention effective or ineffective, in case the required information size is not achieved.
Christopher Jennison and Bruce Turnbull, however, have shown that on average, when a small, but realistic and important intervention effect is anticipated, a group sequential design requires fewer patients than an adaptive design, e.g., reestimating the (fixed) sample size after the first interim analysis [80]. The group sequential design seems more efficient than the adaptive design. In line with mathematical theory [72], simulation studies [6], and empirical considerations [44, 45, 81, 82], there is evidence that small trials and small metaanalyses by chance tend to overestimate the intervention effect or underestimate the variance. Early indicated large intervention effects are often contradicted in later published large trials or large metaanalyses [6, 45, 81, 82]. The reason might be that statistical confidence intervals and significance tests, relating exclusively to the null hypothesis, ignore the necessity of a sufficiently large number of observations to assess realistic or minimally important intervention effects. The early statistical significance, at the 5% level, may be a result of an early overestimation of the intervention effect or an underestimation of the variance, or both, when the required information size for a realistic effect is not achieved. In general, it is easier to reject the null hypothesis than to reject a small, but realistic and still important, alternative hypothesis [64]. The null hypothesis can never be proven, and in practice, this means that it can never be completely discarded, as this would require an infinitely large number of observations.
The reason for early spurious significant findings may be quite simple, although not selfevident. Even adequate randomisation in a small trial lacks ability to ensure the balance between all the involved, known or unknown, prognostic factors in the intervention groups [81]. When we find a statistically significant intervention effect in a small trial or in a small metaanalysis, it is often due to insufficient balance of important prognostic factors, known or unknown, between the intervention groups. Therefore, it is not necessarily intervention effects that we observe, but rather an uneven distribution of important prognostic factors between groups. In addition to the described risks of random error, the overall risk of bias which includes the risk of publication bias makes it understandable why published trials and metaanalyses often result in unreliable estimates of intervention effects [2, 83].
The power of frequentist inference in a single trial and in a metaanalysis of several trials lies in two basic assumptions. First, the only decisive difference between the intervention groups during the trial is the difference between the interventions. We conclude that ‘despite everything else’, the measured difference in the outcome is due to different properties of the interventions because everything else seems equal in the groups. In a small trial and a small metaanalysis, the assumption, that all other risk factors are equally distributed in the two intervention groups, may not be fulfilled as described above, even though adequate bias control has been exercised. Second, the power of frequentist inference depends on the correctness of applying the ‘reverse law of implication’ from mathematical logic (see Appendix): that a sufficiently small Pvalue, calculated as the probability that we got a specific trial result when the null hypothesis is in fact true, leads us to discard the null hypothesis itself. This assumption, which never totally excludes the possibility that the result of a trial may agree with or be a result of the null hypothesis, demands a specific a priori chosen threshold for statistical significance. That is, a sufficiently small Pvalue leads us to regard the trial result as virtually impossible under the null hypothesis, and, therefore, we regard the opposite to be true and discard the null hypothesis. This automatically raises the question: how small a Pvalue should be before we can apply the ‘reverse law of implication’. Or alternatively expressed, does a Pvalue less than an a priori chosen threshold of statistical significance reject the null hypothesis? Ronald A. Fisher, already in 1956, warned against using a statistical significance level of 5% in all situations [84]. Nevertheless, ever since, it seems to have broadly been implemented as a criterion for conclusion in medical research [83], and this is likely wrong [85].
Sequential interimanalyses in a single trial and TSA of metaanalyses of several trials deal systematically and rationally with the misunderstood application of a constant level of statistical significance (P < 0.05), unrelated to the accrued fraction of the precalculated required (fixed) sample or information size and number of trials.
Conclusions
Most systematic reviews with metaanalyses, including Cochrane systematic reviews, do not have sufficient statistical power to detect or reject even large intervention effects. Metaanalyses are updated continuously, and, therefore, ought to be regarded as interimanalyses on the way towards a required information size. The evaluation of metaanalyses ought to relate the total number of randomised participants to the required metaanalytic information size and the corresponding number of required trials considering statistical diversity. When the number of participants in a metaanalysis is less than the required, based on a realistic and minimally important intervention effect, the constant application of a traditional naïve 95% confidence interval or a naïve 5% statistical significance threshold will lead to too many false positive and false negative conclusions. The LanDeMets’ sequential monitoring boundaries in TSA offer adjusted, expanded confidence intervals and adjusted, restrictive thresholds for statistical significance when the diversityadjusted required information size and the required number of trials for the metaanalysis has not been reached. A Bayesian metaanalysis, using prior distributions for both the intervention effect and the statistical heterogeneity, may even be more reliable for deciding whether an intervention effect is present or not. However, the Bayesian metaanalysis also poses difficulties with interpretation. Until easytouse software programs for full Bayesian metaanalysis become accessible, TSA represents a better assumptiontransparent analysis than the use of traditional metaanalysis with unadjusted confidence intervals and unadjusted thresholds for statistical significance.
Abbreviations
 ϴ :

Intervention effect
 α :

Maximally allowed type 1 error
 AIS :

Acquired information size in a metaanalysis
 β :

Maximally allowed type 2 error
 BF :

Bayes factor
 D ^{2} :

Statistical diversity (Dsquare)
 H _{ 0 } :

Null hypothesis
 H _{ A } :

Alternative hypothesis
 I ^{2} :

Statistical inconsistency (Isquare)
 Ln :

Natural logarithm
 P(DH _{ 0 }):

Probability of getting dataset if the null hypothesis is true
 RIS :

Required information size in a metaanalysis
 RR :

Relative risk
 SD :

Standard deviation
 SE :

Standard error
 SS :

Sample size
 TSA:

Trial sequential analysis
 TTM Trial:

Target Temperature Management Trial
 Z _{ i } :

Cumulative z_{ i }statistics
References
Turner RM, Bird SM, Higgins JP. The impact of study size on metaanalyses: examination of underpowered studies in Cochrane reviews. PLoS One. 2013;8:e59202.
Pereira TV, Ioannidis JP. Statistically significant metaanalyses of clinical trials have modest credibility and inflated effects. J Clin Epidemiol. 2011;64:1060–9.
AlBalawi Z, McAlister FA, Thorlund K, Wong M, Wetterslev J. Random error in cardiovascular metaanalyses: how common are false positive and false negative results? Int J Cardiol. 2013;168:1102–7.
Imberger G. Multiplicity and sparse data in systematic reviews of anaesthesiological interventions: a cause of increased risk of random error and lack of reliability of conclusions? Ph.D. Thesis. Copenhagen: Copenhagen University, Faculty of Health and Medical Sciences; 2014.
Brok J, Thorlund K, Wetterslev J, Gluud C. Apparently conclusive metaanalyses may be inconclusive—trial sequential analysis adjustment of random error risk due to repetitive testing of accumulating data in apparently conclusive neonatal metaanalyses. Int J Epidemiol. 2009;38:287–98.
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 metaanalysis—a simulation study. PLoS One. 2011;6:e25491.
Wetterslev J, Thorlund K, Brok J, Gluud C. Trial sequential analysis may establish when firm evidence is reached in cumulative metaanalysis. J Clin Epidemiol. 2008;61:64–75.
Pogue J, Yusuf S. Cumulating evidence from randomised trials: utilizing sequential monitoring boundaries for cumulative metaanalysis. Control Clin Trials. 1997;18:580–93.
Pogue J, Yusuf S. Overcoming the limitations of current metaanalysis of randomised controlled trials. Lancet. 1998;351:47–52.
Thorlund K, Engstrøm J, Wetterslev J, Brok J, Imberger G, Gluud C. User manual for trial sequential analysis (TSA). Copenhagen Trial Unit, Centre for Clinical Intervention research, Copenhagen, Denmark. 2011: 1–115 available from www.ctu.dk/tsa.
Wetterslev J, Thorlund K, Brok J, Gluud C. Estimating required information size by quantifying diversity in a randomeffects metaanalysis. BMC Med Res Methodol. 2009;9:86.
Thorlund K, Engstrøm J, Wetterslev J, Brok J, Imberger G, Gluud C. Software for trial sequential analysis (TSA) ver. 0.9.5.5 Beta. Copenhagen Trial Unit, Centre for Clinical Intervention Research, Copenhagen, Denmark, freeware available at www.ctu.dk/tsa.
Young C, Horton R. Putting clinical trials into context. Lancet. 2005;366:107–8.
Clarke M, Horton R. Bringing it all together: LancetCochrane collaborate on systematic reviews. Lancet. 2001;357:1728.
Clarke M, Hopewell S, Chalmers I. Clinical trials should begin and end with systematic reviews of relevant evidence: 12 years and waiting. Lancet. 2010;376:20–1.
Nielsen N, Friberg H, Gluud C, Wetterslev J. Hypothermia after cardiac arrest should be further evaluated—a systematic review of randomised trials with metaanalysis and trial sequential analysis. Int J Cardiol. 2011;151:333–41.
Nielsen N, Wetterslev J, Cronberg T, Erlinge D, Gasche Y, Hassager C, Horn J, Hovdenes J, Kjaergaard J, Kuiper M, Pellis T, Stammet P, Wanscher M, Wise MP, Åneman A, AlSubaie N, Boesgaard S, BroJeppesen J, Brunetti I, Bugge JF, Hingston CD, Juffermans NP, Koopmans M, Køber L, Langørgen J, Lilja G, Møller JE, Rundgren M, Rylander C, Smid O, Werer C, Winkel P, Friberg H, TTM Trial Investigators. Targeted temperature management at 33°C versus 36°C after cardiac arrest. N Engl J Med. 2013;369:2197–206.
Nielsen N, Wetterslev J, alSubaie N, Andersson B, BroJeppesen J, Bishop G, Brunetti I, Cranshaw J, Cronberg T, Edqvist K, Erlinge D, Gasche Y, Glover G, Hassager C, Horn J, Hovdenes J, Johnsson J, Kjaergaard J, Kuiper M, Langørgen J, Macken L, Martinell L, Martner P, Pellis T, Pelosi P, Petersen P, Persson S, Rundgren M, Saxena M, Svensson R, Stammet P, Thorén A, Undén J, Walden A, Wallskog J, Wanscher M, Wise MP, Wyon N, Aneman A, Friberg H. Target temperature management after outofhospital cardiac arrest – a randomised, parallelgroup, assessorblinded clinical trial – rationale and design. Am Heart J. 2012;163:541–8.
Lan KKG, DeMets DL. Discrete sequential boundaries for clinical trials. Biometrika. 1983;70:659–63.
Peberdy MA, Callaway CW, Neumar RW, Geocadin RG, Zimmerman JL, Donnino M, Gabrielli A, Silvers SM, Zaritsky AL, Merchant R, Vanden Hoek TL, Kronick SL, American Heart Association. Part 9: postcardiac arrest care: American Heart Association Guidelines for Cardiopulmonary Resuscitation and Emergency Cardiovascular Care. Circulation. 2010;122 suppl 3:S768–86.
Armitage P, McPherson CK, Rowe BC. Repeated significance tests on accumulating data. J Royal Stat Soc Series A (General). 1969;132:235–44.
Pocock SJ. Group sequential methods in the design and analysis of clinical trials. Biometrika. 1977;64:191–9.
Berkey CS, Mosteller F, Lau J, Antman EM. Uncertainty of the time of first significance in random effects cumulative metaanalysis. Control Clin Trials. 1996;17:357–71.
Imberger G, Vejlby AD, Hansen SB, Møller AM, Wetterslev J. Statistical multiplicity in systematic reviews of anaesthesia interventions: a quantification and comparison between Cochrane and nonCochrane reviews. PLoS One. 2011;6:e28422.
Wald A. Contributions to the theory of statistical estimation and testing hypotheses. Ann Math Stat. 1939;10:299–326.
Wald A. Sequential tests of statistical hypotheses. Ann Math Stat. 1945;16:117–86.
Wald A, Wolfowitz J. Bayes solutions of sequential decision problems. Proc Natl Acad Sci U S A. 1949;35:99–102.
Winkel P, Zhang NF. Statistical development of quality in medicine. Chichester, West Sussex: Wiley; 2007. p. 1–224.
Armitage P. The evolution of ways of deciding when clinical trials should stop recruiting. James Lind Library Bulletin 2013. www.jameslindlibrary.org.
Dunn OJ. Multiple comparisons among means. J Am Stat Assoc. 1961;56:52–64.
Peto R, Pike MC, Armitage P, Breslow NE, Cox DR, Howard SV, Mantel N, McPherson K, Peto J, Smith PG. Design and analysis of randomised clinical trials requiring prolonged observation of each patient. I. Introduction and design. Br J Cancer. 1976;34:585–612.
O’Brien PC, Fleming TR. A multiple testing procedure for clinical trials. Biometrics. 1979;35:549–56.
ICH Harmonised Tripartite Guideline. Statistical principles for clinical trials. International Conference on Harmonisation E9 Expert Working Group. Stat Med. 1999;18:1905–42.
Kim K, DeMets DL. Confidence intervals following group sequential tests in clinical trials. Biometrics. 1987;43:857–64.
DeMets DL. Group sequential procedures: calendar versus information time. Stat Med. 1989;8:1191–8.
Jennison C, Turnbull BW. Group sequential methods with applications to clinical trials. Boca Raton: Chapman & Hall/CRC Press; 2000.
Grant AM, Altman DG, Babiker AB, Campbell MK, Clemens FJ, Darbyshire JH, Elbourne DR, McLeer SK, Parmar MK, Pocock SJ, Spiegelhalter DJ, Sydes MR, Walker AE, Wallace SA, DAMOCLES Study Group. Issues in data monitoring and interim analysis of trials. Health Technol Assess. 2005;9:1–238. iii–iv.
Chow S, Shao J, Wang H. Sample size calculations in clinical research. Taylor & Francis/CRC: Boca Raton; 2003.
Reboussin DM, DeMets DL, Kim KM, Lan KK. Computations for group sequential boundaries using the LanDeMets spending function method. Control Clin Trials. 2000;21:190–207.
DerSimonian R, Laird N. Metaanalysis in clinical trials. Control Clin Trials. 1986;7:177–88.
Deeks JJ, Higgins JPT. Statistical algorithms in Review Manager ver. 5.3. On behalf of the Statistical Methods Group of The Cochrane Collaboration. 2010.
Higgins JP, Thompson SG. Quantifying heterogeneity in a metaanalysis. Stat Med. 2002;21:1539–58.
Kulinskaya E, Wood J. Trial sequential methods for metaanalysis. Res Synth Methods. 2014;5:212–20.
Thorlund K, Devereaux PJ, Wetterslev J, Guyatt G, Ioannidis JP, Thabane L, Gluud LL, AlsNielsen B, Gluud C. Can trial sequential monitoring boundaries reduce spurious inferences from metaanalyses? Int J Epidemiol. 2009;38:276–86.
Imberger G, Thorlund K, Gluud C, Wetterslev J. False positive findings in cumulative metaanalysis with and without application of trial sequential analysis: an empirical review. BMJ Open. 2016;6(8):e011890.
Imberger G, Gluud C, Boylan J, Wetterslev J. Systematic reviews of anesthesiologic interventions reported as statistically significant: problems with power, precision, and type 1 error protection. Anesth Analg. 2015;121:1611–22.
Mascha EJ. Alpha, beta, meta: guidelines for assessing power and type I error in metaanalyses. Anesth Analg. 2015;121:1430–3.
Turner RM, Davey J, Clarke MJ, Thompson SG, Higgins JP. Predicting the extent of heterogeneity in metaanalysis, using empirical data from the Cochrane Database of Systematic Reviews. Int J Epidemiol. 2012;41:818–27.
Thorlund K, Imberger G, Johnston BC, Walsh M, Awad T, Thabane L, Gluud C, Devereaux PJ, Wetterslev J. Evolution of heterogeneity (I^{2}) estimates and their 95% confidence intervals in large metaanalyses. PLoS One. 2012;7:e39471.
Brok J, Thorlund K, Gluud C, Wetterslev J. Trial sequential analysis reveals insufficient information size and potentially false positive results in many metaanalyses. J Clin Epidemiol. 2008;61:763–9.
Higgins JPT, Green S. red. Cochrane Handbook for Systematic Reviews of Interventions Version 5.1.0. The Cochrane Collaboration, 2011. www.cochranehandbook.org.
Keus F, Wetterslev J, Gluud C, van Laarhoven CJ. Evidence at a glance: error matrix approach for overviewing available evidence. BMC Med Res Methodol. 2010;10:90.
Garattini S, Jakobsen JC, Wetterslev J, Berthele’ V, Banzi R, Rath A, Neugebauer E, Laville M, Maisson Y, Hivert Y, Eickermann M, Aydin B, Ngwabyt S, Martinho C, Giradi C, Szmigielski C, DemotesMaynard J, Gluud C. Evidencebased clinical practice: overview of threats to the validity of evidence. Eur J Intern Med. 2016;32:13–21.
Kjaergard LL, Villumsen J, Gluud C. Reported methodological quality and discrepancies between large and small randomised trials in metaanalyses. Ann Intern Med. 2001;135:982–9. err 2008;149:219.
Savović J, Jones HE, Altman DG, Harris RJ, Jüni P, Pildal J, AlsNielsen B, Balk EM, Gluud C, Gluud LL, Ioannidis JP, Schulz KF, Beynon R, Welton NJ, Wood L, Moher D, Deeks JJ, Sterne JA. Influence of reported study design characteristics on intervention effect estimates from randomised, controlled trials. Ann Intern Med. 2012;157:429–38.
Lundh A, Sismondo S, Lexchin J, Busuioc OA, Bero L. Industry sponsorship and research outcome. Cochrane Database Syst Rev. 2012;12:MR000033.
Chan AW, Hróbjartsson A, Haahr MT, Gøtzsche PC, Altman DG. Empirical evidence for selective reporting of outcomes in randomised trials: comparison of protocols to published articles. JAMA. 2004;291:2457–65.
Andrews JC, Schünemann HJ, Oxman AD, Pottie K, Meerpohl JJ, Coello PA, Rind D, Montori VM, Brito JP, Norris S, Elbarbary M, Post P, Nasser M, Shukla V, Jaeschke R, Brozek J, Djulbegovic B, Guyatt G. GRADE guidelines: 15. Going from evidence to recommendationdeterminants of a recommendation’s direction and strength. J Clin Epidemiol. 2013;66:726–35.
The Fermi paradox. http://en.wikipedia.org/wiki/Fermi_paradox. Accessed 27 Feb 2017.
Roberts I, Ker K, Edwards P, Beecher D, Manno D, Sydenham E. The knowledge system underpinning healthcare is not fit for purpose and must change. BMJ. 2015;350:h2463. doi:10.1136/bmj.h2463.
Bolland MJ, Grey A, Gamble GD, Reid IR. The effect of vitamin D supplementation on skeletal, vascular, or cancer outcomes: a trial sequential metaanalysis. Lancet Diabetes Endocrinol. 2014;2(4):307–20. doi:10.1016/S22138587(13)702122.
Tovey DI, Bero L, Farquhar C, Lasserson T, MacLehose H, Macdonald G, et al. A response to Ian Roberts and his colleagues. Rapid response. BMJ. 2015;350:h2463.
Wetterslev J, Engstrøm J, Gluud C, Thorlund K. Trial sequential analysis: methods and software for cumulative metaanalyses. Cochrane Methods Cochrane Database Syst Rev. 2012;2 suppl 1:29–31.
Higgins JPT. Comment on “Trial sequential analysis: methods and software for cumulative metaanalyses”. Cochrane Methods Cochrane Database Syst Rev. 2012;2 suppl 1:32–3.
Wetterslev J, Engstrøm J, Gluud C, Thorlund K. Response to “Comment by Higgins”. Cochrane Methods Cochrane Database Syst Rev. 2012;2 suppl 1:33–5.
Higgins JP, Whitehead A, Simmonds M. Sequential methods for randomeffects metaanalysis. Stat Med. 2011;30:903–21.
Fleisher LA, Beckman JA, Brown KA, Calkins H, Chaikof EL, Fleischmann KE, Freeman WK, Froehlich JB, Kasper EK, Kersten JR, Riegel B, Robb JF, Smith Jr SC, Jacobs AK, Adams CD, Anderson JL, Antman EM, Buller CE, Creager MA, Ettinger SM, Faxon DP, Fuster V, Halperin JL, Hiratzka LF, Hunt SA, Lytle BW, Nishimura R, Ornato JP, Page RL, Riegel B, Tarkington LG, Yancy CW. ACC/AHA 2007 guidelines on perioperative cardiovascular evaluation and care for noncardiac surgery: executive summary: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines (Writing Committee to Revise the 2002 Guidelines on Perioperative Cardiovascular Evaluation for Noncardiac Surgery). J Am Coll Cardiol. 2007;50:1707–32.
Popper KR. Logik der Forschung. Vienna: Springer; 1959.
Bangalore S, Wetterslev J, Pranesh S, Sawhney S, Gluud C, Messerli FH. Perioperative beta blockers in patients having noncardiac surgery: a metaanalysis. Lancet. 2008;372:1962–76.
Jakobsen JC, Wetterslev J, Winkel P, Lange T, Gluud C. The threshold for statistical and clinical significance in systematic reviews with metaanalytic methods. Med Res Methodol. 2014;14:120.
Sterne JA. Teaching hypothesis tests – time for significant change? Stat Med 2002;21: 985–94, 995–9, 1001.
Jakobsen JC, Gluud C, Winkel P, Lange T, Wetterslev J. The thresholds for statistical and clinical significance – a fivestep procedure for evaluation of intervention effects in randomised clinical trials. BMC Med Res Methodol. 2014;14:34.
Roloff V, Higgins JP, Sutton AJ. Planning future studies based on the conditional power of a metaanalysis. Stat Med. 2013;32:11–24.
IntHout J, Ioannidis JP, Borm GF. Obtaining evidence by a single wellpowered trial or several modestly powered trials. Stat Methods Med Res. 2016;25(2):538–52.
Valentine JC, Pigott TD, Rothstein HR. How many studies do you need? A primer on statistical power for metaanalysis. J Educ Behav Stat. 2010;35(2):215–47.
Borenstein M, Hedges LV, Higgins JPT, Rothstein HR. Introduction to metaanalysis. Chichester: John Wiley & Sons Ltd.; 2009.
Higgins JP, Spiegelhalter DJ. Being sceptical about metaanalyses: a Bayesian perspective on magnesium trials in myocardial infarction. Int J Epidemiol. 2002;31:96–104.
Spiegelhalter DJ, Abrams KR, Myles JP. Bayesian approaches to clinical trials and healthcare evaluation. Statistics in practice. Chichester: John Wiley & Sons Ltd; 2004.
Higgins JP, Thompson SG, Spiegelhalter DJ. A reevaluation of randomeffects metaanalysis. J R Stat Soc Ser A Stat Soc. 2009;172:137–59.
Jennison C, Turnbull BW. Efficient group sequential designs when there are several effect sizes under consideration. Stat Med. 2006;25:917–32.
Pereira TV, Horwitz RI, Ioannidis JP. Empirical evaluation of very large treatment effects of medical interventions. JAMA. 2012;308:1676–84.
Lindley DV. A statistical paradox. Biometrika. 1957;44:187–92.
Ioannidis JPA. Why most published research findings are false. PLoS Med. 2005;2:e124.
Fisher R. Statistical methods and scientific induction. J R Stat Soc Ser B. 1955;17:69–78.
Johnson EV. Revised standards for statistical evidence. PNAS. 2013, 110:48:19313–19317. Accessed Dec 2016. www.pnas.org/cgi/doi/10.1073/pnas.1313476110.
Acknowledgements
We thank Ema Erkocevic Petersen, M.Sci. Biomedical Engineering and Informatics, Data Manager at Copenhagen Trial Unit for the excellent art work on figures in this article. We thank Dimitrinka Nikolova, M.A. in English Philology, Coordinating Editor of The Cochrane HepatoBiliary Group, for thorough reading and correcting the language of the article text.
Funding
The work on this article has exclusively been funded by the Copenhagen Trial Unit with no other funding involved.
Availability of data and materials
TSA software and Manual are available at http://www.ctu.dk/tsa/for free.
Data for the motivating example are available from (16–18).
Authors’ contribution
JW, JCJ, and CG conceived the idea for this manuscript. JW wrote the first draft of the manuscript. JCJ and CG critically revised the manuscript. All authors read and approved the final version of the manuscript.
Competing interest
None of the authors have financial interests related to this article. JCJ does not have any other known competing interests. JW and CG are members of the task force to develop theory and software for Trial Sequential Analysis at the Copenhagen Trial Unit.
Consent for publication
Not applicable since no individual patient data are presented.
Ethics approval and consent to participate
Not applicable.
Author information
Authors and Affiliations
Corresponding author
Appendix
Appendix
Sample size in a single randomised trial
Dichotomous outcome
When the intervention and the control group are equalsized groups, the sample size (SS) is calculated as [38]:
In this formula, P _{ E } and P _{ C } are the frequencies of the outcome in the experimental group and the control group, ϴ is the intervention effect that the trial wants to address with μ = P _{ C } – P _{ E }, α and β are the maximally allowed type I and type II error risks, and ν is the variance of the outcome difference between the two groups.
The test statistic Z _{ i } at the ith interim analysis for the dichotomous outcome measure:
with P _{ Ei } − P _{ Ci } being the estimate of the intervention effect at the ith interim analysis, and Var(P _{ Ei } − P _{ Ci }) is the variance of the estimate [36].
Continuous outcome
If the anticipated difference between the means in the control group and the experimental group is assumed to be X _{1}–X _{2} with the standard deviation SD, α and β being the type I and type II error risks, the SS is calculated as [38]:
The cumulative Z _{ i }value at the ith interim analysis is calculated as:
X _{ 1i } – X _{ 2i } being the estimate of the intervention effect at the ith interim analysis, and SD _{ (X1i –X2i) } its standard deviation [36].
The required information size in a metaanalysis
The required information size (RIS) in a metaanalysis of a dichotomous outcome can be calculated as (see the definition of μ and ν above) [11]:
Therefore, RIS emerges as the sample size estimated for a single trial with corresponding power to detect or reject an anticipated intervention effect μ, multiplied with a factor adjusting for the final expected or the present diversity (D ^{2}) in the metaanalysis. P_{C} can be estimated as the pooled unweighted proportion of outcomes in the control groups in the included trials [11, 39–41].
Alternatively, but with equal result, the RIS for a randomeffects model can be calculated as:
Where ν _{ random } is the variance in the randomeffects model. This is the model variance based calculation of RIS.
When the required information size has been estimated, the metaanalysis can be analysed in relation to the trial sequential monitoring boundaries, constructed as per LanDeMets’ critical monitoring boundaries, analogous to the sequential analysis of a single randomised trial [17, 42]. Similar to the Z _{ i }value for the ith cumulative metaanalysis, we use here the ratio between the logarithm of the MantelHaenszels weighted relative risik (RR _{ iMH }) and its standard error SE[ln(RR _{ iMH })] in a randomeffects model, e.g., as suggested by DerSimonian and Laird [44, 45]:
The required information size and the number of required extra trials in a metaanalysis
As per Elena Kulinskaya and John Wood [43], using the same mathematical notation, the sample size of future K ^{'} equally sized trials with equally sized (n _{ i }) intervention and control groups is:
where σ ^{2} is the variance in these future trials and τ ^{2} the between trial variance in the first K trials, and c is a constant when the results from previous trials in the randomeffects metaanalysis are known being:
V ^{− 1}_{ R(K)} is the reciprocal of the variance in the randomeffects model of the pooled estimate from the first K trials which is equal to SD ^{− 2}_{ R(K)} in a metaanalysis of a continuous outcome. In the examples, given by Kulinskaya and Wood, σ ^{2} is determined by the simple unweighted average of the within trial variances in the first K trials [43]. However, as it may be more appropriate to use an estimate of the variance in the future trials which is the weighted average of the pooled squared standard deviations from the fixedeffect model, SD ^{2}_{ F(K)} , of the first K trials for continuous outcomes, we propose:
where N is the total number of participants in the metaanalysis of the first K trials and SE ^{2}_{ F(K)} is the squared pooled standard error in the fixedeffect model. If N ^{’} is the sum of sample sizes from the equally sized K ^{'} future trials, the formula (1) for the sample size in future trials can be rewritten to:
from which it is evident that K’ has to be an integer greater than c ⋅ τ ^{2} (the minimum required number of extra trials) and that N’ converges to 4 · σ ^{2} ⋅ c when K’ approaches the infinite as:
meaning that N’ can never be less than 4 · σ ^{2} · c, despite the number of future trials being carried out. Because the randomeffects model variance of the first K trials is always greater than the fixedeffect model variance of the K first trials, we have:
and as: c = ([Z _{1 − α/2} + Z _{1 − β}]/θ)^{2} − SD ^{− 2}_{ R(K)} , we get
If RIS is defined as the required information size for a metaanalysis to have as many participants as an equally powered trial to address a hypothesis of the same intervention effect. with the randomeffects model variance experienced [11] so far, we will have:
and thereby the equation (2) can be rewritten as:
which shows that RIS is always greater than the minimum required extra participants in the final metaanalysis. If V ^{− 1}_{ R(K)} = SD ^{− 2}_{ R(K)} is small (close to zero when abundant information has already been collected), the statement that RIS > 4 · σ ^{2} · c is merely the trivial that the required information size will be greater than the sample size in one future trial with an anticipated variance of σ ^{2}.
Two scenarios now cover all possible situations:

A)
If: RIS − AIS > 4 ⋅ SD ^{2}_{ F(K)} ⋅ c, then we will be able to calculate the corresponding required number (K’ _{ RISAIS }) of future trials to cover the information gap indicated by the RISAIS participants according to formula (1), having:
$$ K{'}_{RIS AIS} = \frac{c\cdot {\tau}^2}{1\frac{4\cdot {\sigma}^2\cdot c}{\left( RIS AIS\right)}} $$ 
B)
If RIS − AIS ≤ 4 ⋅ SD ^{2}_{ F(K)} ⋅ c, then the situation is that the extra required information size is less than what is required in one extra trial. The question of whether there is an intervention effect greater or equal to θ may be answered if RIS has been achieved in the minimal number of required trials.
Moreover, if for example:
with RIS − AIS = 3 ⋅ 4 ⋅ σ ^{2} ⋅ c in the formula for K ’ _{ RIS − AIS }, then this leads to:
and we have:
meaning that the extra number of required trials, K ’ _{ RIS − AIS }, for the randomeffects metaanalysis to be adequately powered using RISAIS, would be 50% higher than the minimal required number of trials. In contrast, the minimal required number of trials, c ⋅ τ ^{2}, demands a very large total information size.
The law of reverse implication in mathematical logic
If A and B are statements that can be determined to be true or false, it follows from mathematical logic that the statement: (from A follows B), is equivalent to the statement: (from nonB follows nonA). That is: the statement (»from the null hypothesis being true follows that the result of the trial is unlikely«), is equivalent with the statement: (»from the trial result being likely follows that the null hypothesis is false«). If the probability that the trial result has emerged from a true null hypothesis is regarded as »unlikely«, when the Pvalue is lower than a certain level of significance, the null hypothesis can be rejected.
The statistical concepts used in this article
The null hypothesis (H _{0}) is the hypothesis that the mean or the occurence of a certain outcome is precisely the same in the groups compared.
The alternative hypothesis (H _{A}) is the assumed magnitude of the difference between the means or the occurence of a certain outcome in the groups compared.
The Pvalue is the probability (P) that a specific dataset (D), or something even more extreme, would appear if the null hypothesis is true, that is: P = P(DH _{0}) which shold be read as: the Pvalue is equal to the probability of getting the dataset, given the null hypothesis is true.
The Pvalue is intuitively difficult to understand and is often wrongly interpreted. In Bayesian statistic, one can obtain results much easier to understand, and the clinical significance of a result from a trial or a metaanalysis may, therefore, be easier to achieve. One can present, e.g., the »reverse« probability, being how large is the probability that the null hypothesis is true given a specific dataset (D) from a trial or a metaanalysis: P(H _{0}D). Furthermore, how large is the probability that an alternative hypothesis is true given a specific dataset from a trial or a metaanalysis can be expressed as the ‘alternative’ Pvalue: P(H _{A}D).
Bayes factor (BF) is the ratio between the Pvalue and the alternative Pvalue: BF = P(DH _{ 0 })/P(DH _{ A }).
Due to Bayes theorem, we have:
A cumulative metaanalysis is a metaanalysis where the result data from the last conducted trial are added to the result data from all the previous trials and tested on the accumulated number of participants.
The Zvalue is the general name for the test statistics, e.g., the tvalue in a student’s ttest or a χ ^{2}test statistic in a χ ^{2}test or the MantelHaenszels test statistic. When the theoretical distribution of the test statistic is known, the Zvalue can be converted into a Pvalue.
A cumulative Z _{ i }value is the test statistic calculated after addition of the data from the ith trial.
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.
About this article
Cite this article
Wetterslev, J., Jakobsen, J.C. & Gluud, C. Trial Sequential Analysis in systematic reviews with metaanalysis. BMC Med Res Methodol 17, 39 (2017). https://doi.org/10.1186/s1287401703157
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287401703157