 Research article
 Open Access
 Published:
Considerations for pooling realworld data as a comparator cohort to a single arm trial: a simulation study on assessment of heterogeneity
BMC Medical Research Methodology volume 23, Article number: 193 (2023)
Abstract
Background
Novel precision medicine therapeutics target increasingly granular, genomicallydefined populations. Rare subgroups make it challenging to study within a clinical trial or single realworld data (RWD) source; therefore, pooling from disparate sources of RWD may be required for feasibility. Heterogeneity assessment for pooled data is particularly complex when contrasting a pooled realworld comparator cohort (rwCC) with a singlearm clinical trial (SAT), because the individual comparisons are not independent as all compare a rwCC to the same SAT. Our objective was to develop a methodological framework for pooling RWD focused on the rwCC use case, and simulate novel approaches of heterogeneity assessment, especially for small datasets.
Methods
We present a framework with the following steps: prespecification, assessment of dataset eligibility, and outcome analyses (including assessment of outcome heterogeneity). We then simulated heterogeneity assessments for a binary response outcome in a SAT compared to two rwCCs, using standard methods for metaanalysis, and an Adjusted Cochran’s Q test, and directly comparing the individual participant data (IPD) from the rwCCs.
Results
We found identical power to detect a true difference for the adjusted Cochran’s Q test and the IPD method, with both approaches superior to a standard Cochran’s Q test. When assessing the impact of heterogeneity in the null scenario of no difference between the SAT and rwCCs, a lack of statistical power led to Type 1 error inflation. Similarly, in the alternative scenario of a true difference between SAT and rwCCs, we found substantial Type 2 error, with underpowered heterogeneity testing leading to underestimation of the treatment effect.
Conclusions
We developed a methodological framework for pooling RWD sources in the context of designing a rwCC for a SAT. When testing for heterogeneity during this process, the adjusted Cochran’s Q test matches the statistical power of IPD heterogeneity testing. Limitations of quantitative heterogeneity testing in protecting against Type 1 or Type 2 error indicate these tests are best used descriptively, and after careful selection of datasets based on clinical/data considerations. We hope these findings will facilitate the rigorous pooling of RWD to unlock insights to benefit oncology patients.
Background
In the era of precision medicine, generating insights for rare genomicallydefined subpopulations often requires pooling of data from multiple sources [1]. Within oncology, an increasingly common use case is the pooling of realworld data (RWD) from different sources to serve as a realworld comparator cohort (rwCC) for a singlearm clinical trial (SAT) in rare subpopulations for research and to support regulatory decisionmaking [2,3,4].
General considerations for combining data from multiple clinical research studies have been described, with pooling individual participant data (IPD) [5, 6] as the preferred method; studylevel metaanalysis approaches are also frequently used [7,8,9]. Testing for heterogeneity in outcomes between studies is a standard step in the pooling process [7, 10], and is critical when creating a rwCC from disparate datasets, as a fundamental assumption of a pooled analyses is that comparison of a SAT to two different rwCC should yield statistically identical results (i.e., no heterogeneity).
However, there remain challenges for establishing best practices for pooling heterogeneous RWD into a rwCC. First, using RWD as a data source for clinical research requires consideration of nuances that may not be as critical with clinical trial data (e.g., differences in quality and organization of the datasets) [11, 12]. Yet, there is no overarching framework, specifically for pooling RWD, that addresses both questions of which studies should be merged (‘what to pool?’), as well as how to integrate them (‘how to pool?’) [13, 14]. Second, standard methods for conducting a heterogeneity assessment are insufficient, as they require unmet assumptions for the rwCC use case; notably, the effect estimates pooled in the metaanalysis are not independent because the SAT is included in all comparisons. Third, the small sample sizes often observed within genomicallydefined subpopulations [2] can negatively affect all methods for assessing heterogeneity, with implications for the role of this testing in data analyses, and the interpretation of its results.
To address a gap in published guidance for addressing the above issues, the following provides a framework for the methodologic considerations regarding pooling RWD and specifically for the rwCC use case from the perspective of the sponsor of a novel treatment. Supporting simulations then explore various methods of assessing heterogeneity in outcomes between rwCCs, including an approach that solves issues of nonindependence. Conclusions are drawn regarding the strengths and limitations of various heterogeneity tests in different scenarios, and how their results should be interpreted in the context of a rwCC analysis.
A framework for pooling datasets into a realworld comparator cohort
In the following sections, the steps of the proposed framework describing the pooling of datasets into a rwCC are described, with component steps summarized in Table 1. This framework was developed based on the authors’ review of the literature and experience with data pooling projects.
Prespecification of research question and pooling processes
First, the clinical research question should be prespecified [15], followed by preparation of the associated statistical analysis plan. The data requirements for these analyses can then inform the prespecification of the processes for both assessing the eligibility of the candidate datasets for inclusion in the systematic review [16], and whether pooling is appropriate.
Assess dataset eligibility
Qualitatively assess metadata for all datasets & variables for relevance, reliability, and harmonizability
Before undertaking any assessments related to pooling, each individual dataset that has been identified for detailed assessment by meeting basic eligibility criteria [16] should first be evaluated for its suitability as a potential rwCC given the prespecified requirements. This could include assessment of their fitnessforuse (i.e., relevance and reliability) [17], and suitability as a rwCC to a specific SAT [18, 19]. Interpretation of the assessment results may require complex clinical judgment, however, published guidance is sparse [20].
Once candidate datasets are identified that individually could serve as rwCCs, the feasibility of pooling them must be determined. For this, the prespecified key variables must be harmonizable across datasets [17, 21]. This includes coding them to have comparable values, which typically involves using the lowest common denominator coding scheme, e.g., collapsing packyears of smoking into ever/never. Such reduction of the granularity of a confounder could lead to residual confounding, in which case investigators may need to weigh the tradeoff between bias and increased sample size [21]. In other cases, the values may be comparable but with different definitions and/or ascertainment, in which case pooling could introduce bias. For example, the outcome variable realworld response, for which patients without a response assessment are sometimes counted as nonresponders and other times excluded.
Quantitatively assess nonoutcome characteristics and sample size after cohort selection
After applying the relevant inclusion/exclusion (I/E) criteria to each dataset, descriptive analyses can assess the sample size of each dataset, as well as the distributions and missingness of key variables in order to assess their comparability to each other and comparison with a priori expectations (e.g., published literature). If desired and allowable given the prespecification, the tradeoff can be assessed between the strictness of the I/E criteria and the sample size. Prior to conducting this analysis, any participants included in multiple datasets should be identified and deduplicated prior to pooling [22], to avoid doublecounting. An a priori hierarchy should determine which data source to use for these participants (e.g., selecting the dataset that optimizes criteria such as data quality and quantity/followup time, either for a particular participant or overall).
Finally, after these assessments are complete and the resulting information is synthesized, the eligibility for the final data analysis of each dataset should be determined according to the prespecified process.
Outcomes analysis
Primary analyses
If it is determined based on the preceding steps that pooling for a particular endpoint is appropriate, then analysis can proceed, treating the realworld datasets as a combined rwCC. If the previous steps provide evidence the individual datasets are not poolable, the sponsor will need to investigate the reasons for any differences [23], as well as analyze and interpret the pooled results appropriately.
Heterogeneity assessment
Assessment of heterogeneity is a standard feature of metaanalyses [7]. Extended to the rwCC use case, the fundamental assumption is that comparison of a SAT to two different rwCC should yield statistically identical results (i.e., no heterogeneity). There has not been indepth consideration of optimal methodologies for testing for heterogeneity in the rwCC use case; this topic is explored further in the remainder of this manuscript. Statistical considerations (described further in the Results) suggest not making the pooling decision solely based on heterogeneity testing, but rather assessing heterogeneity descriptively after the pooling decision has been made based on the considerations for assessing dataset eligibility described above.
Sensitivity analyses
A primary benefit of pooling is that it facilitates certain sensitivity analyses that would be infeasible due to sample size limitations if carried out in each dataset separately. These analyses may include assessment of whether intervention effects vary by participant characteristics [24] (e.g., comparing the effect size between different groups, or restricting to a subset of participants after applying additional criteria that increase similarity with the SAT). In principle, heterogeneity should be described separately for each such sensitivity analysis, but if sample size is reduced in each individual dataset after application of these more stringent criteria, then the power of tests of heterogeneity would be decreased.
Methods
The fundamental assumption of a rwCC analysis is that, after appropriate adjustment for confounders, the SATrwCC comparison yields–up to sampling variability–the same estimate that would have resulted from a randomized controlled trial. Given this assumption, comparison of a SAT to two different rwCC should yield identical results up to sampling variability. However, if heterogeneity is detected, this underlying assumption is likely violated, and pooled results should be interpreted with caution. If there are many rwCCs, a metaanalysis using randomeffects could be used to summarize the control arms before comparing to the SAT, or to summarize the results of the SATrwCC comparisons (although the validity of generating a single summary measure in the presence of heterogeneity is under debate [25]). Here we assume there are only a few rwCCs (e.g., as low as 2 or 3), not enough to make such a metaanalysis feasible [26].
Methods for testing heterogeneity
We consider two ways in which heterogeneity of outcomes across SAT to rwCC comparisons can be assessed.
First, heterogeneity could be assessed by directly comparing the IPD from the rwCCs (the “IPD” method). Before comparing the rwCCs, it is desirable to make the distribution of prognostic covariates similar between each of these rwCCs and the SAT. This could be accomplished using weighting. Each rwCC would be weighted to the SAT using, for example, average treatment effect on the treated weights. Then the weighted rwCC could be compared to each other using appropriate statistical methods (e.g., using logistic regression if outcomes are binary).
Second, heterogeneity among the effect estimates separately comparing the SAT to each rwCC could be assessed; we call this the aggregate method since it requires aggregate statistics from each SATrwCC comparison. This is a traditional nonIPD metaanalysis approach in which the SAT would be compared to each rwCC separately, e.g., using a matching or weighting estimator. Then, the resulting statistics can be assessed for heterogeneity, typically using Cochran’s Q test.
A methodology that accounts for dependence: the Adjusted Qtest.
Applying the aggregate method is challenging because the effect estimates being pooled in the metaanalysis will not be independent, since the SAT is included in all comparisons. Therefore, standard methods for assessing heterogeneity, like Cochran’s Q test for heterogeneity, which rely on independence, will be invalid. However, a modified version of Cochran’s Q test can account for this dependence. We call this adjusted test, which is similar to standard methods in the literature for metaanalysis of dependent results [27, 28], the “Adjusted Qtest”. It requires three steps to calculate:

1.
In the first step, the covariance matrix Σ of the effect estimates (e.g., risk difference, odds ratio, hazard ratio) from the various SATrwCC comparisons is calculated. The effect estimates from the various SATrwCC comparisons will be correlated, since each comparison is to the same SAT. One way to calculate this covariance matrix (with dimensions kbyk, where k is the number of SATrwCC comparisons) is to use the bootstrap. A common randomization seed could be used to generate bootstrap samples of the SAT. Then, for each bootstrap sample of the SAT, effect estimates and associated standard errors for each SATrwCC comparison could be calculated, using a corresponding bootstrap sample of each rwCC. Then these could be used to calculate the covariance matrix of the effect estimates for the rwCCSAT comparisons.

2.
In the second step, we multiply the effect estimates from the SATrwCC comparisons by the inverse Cholesky transformation Σ^{−½}. This is a socalled “whitening transformation” which renders the effect estimates uncorrelated, since if the covariance of the effect estimates is Σ, after multiplying those effect estimates by Σ^{−½}, their covariance matrix will now be the identity matrix [29].

3.
In the third step, Cochran’s Q can be calculated using the resulting statistics, whose assumed covariance matrix is now the identity matrix. Q is then referred to the chisquare distribution with I1 degrees of freedom, where I is the number of rwCC datasets. This is in accordance with normal practice, although it should be noted that, when, as here, the covariance of the statistics is estimated, using the chisquare distribution is only an approximation to the true distribution of Q and may inflate Type I error. Moreover, the distribution of Q depends on the effect measure i.e., hazard ratios or risk differences [30]. We note that bootstrap methods can be used instead of the chisquare distribution for a potentially more accurate test (see Joshi et al.) [31]. Further work could combine the bootstrap methods of Joshi et al. with the use of the bootstrap above in Step 1.
Detecting heterogeneity
We wish to evaluate the impact of the effect sizes on Cochran’s Q and Adjusted Q when assessing heterogeneity of two SATrwCC comparisons using a common SAT. We simulate data from two rwCCs and one SAT, using a binary response outcome. We assume the true response rate (RR) is 50% in the SAT and 50% in one of the rwCC (i.e., no difference). However, due to bias from unmeasured confounding, the observed RR in the other rwCC may be different from 50% (i.e., heterogeneity). We calculate the probability that Cochran’s Q test and Adjusted Q will reject the null hypothesis that the effect sizes versus the SAT are the same (i.e., probability of detecting heterogeneity). We compare to the IPD method by also showing the probability that a logistic regression model rejects the null hypothesis that the rwCC arms have the same RR. There are 100 subjects in the SAT and in each rwCC, and we run 1000 replicates for each scenario. For the Adjusted Q test, we generate 1000 bootstrap samples for each replicate.
Impact of heterogeneity testing
Although potentially inadvisable, one could adopt a conditional approach, where if the above methods did not reveal heterogeneity, then one would proceed with a pooled analysis. This is similar to the testthenpool approach described by Viele et al. [32]. To explore the operating characteristics of this approach, we consider two hypothetical scenarios of interest. In both scenarios, we assume that a SATrwCC comparison is designed using the assumption of an 80% RR in the SAT and 50% RR in the rwCC. To attain 90% power to detect a difference of this size at an alpha level of 5%, roughly 50 participants are required in the SAT and in the rwCC. The sponsor aims to pool two rwCC, each with 25 participants, to achieve the desired power.
Scenario 1: impact of heterogeneity testing in null scenario
The first scenario of interest is a null scenario, where the true RR in the SAT and in one of the rwCC is 50%. However, due to some bias (e.g., missing data on patient responses), the RR in the second rwCC may be lower than 50%. To protect against the Type I error inflation that might result from pooling in this scenario, the sponsor pools only if the pvalue from an IPD heterogeneity assessment of the two rwCCs is greater than 0.10 (0.10 as opposed to 0.05 is used to increase the power to detect heterogeneity).
Scenario 2: impact of heterogeneity testing when treatment improves outcomes
In the second scenario, we consider the RR in the SAT is 80%. The RR in the first rwCC is 50% but due to some bias, the RR in the second rwCC could be higher than 50% (e.g., an unmeasured confounder), potentially leading to Type 2 error.
Results
Methods for testing heterogeneity
Detecting heterogeneity
Figure 1 shows the probability that Cochran’s Q test and Adjusted Q will reject the null hypothesis that the effect sizes versus the SAT are the same (i.e., probability of detecting heterogeneity). When rwCCs with identical true RRs are compared to the same SAT, the Type I error of the Q test is lower than the nominal 5% alpha level since there is no variability in the SAT arm across the comparisons. However, Type I error is close to the nominal level with the Adjusted Q test. In addition, the power functions of the Adjusted Q test and the IPD method are nearly identical, and substantially higher than for the traditional Cochran’s Q. The Adjusted Q test therefore properly accounts for the dependence in test statistics across the SATrwCC comparisons, unlike the traditional Cochran’s Q.
Impact of heterogeneity testing
Scenario 1: impact of heterogeneity testing in null scenario
In Table 2, we show the expected operating characteristics with a range of different RRs in the second rwCC. In the null scenario, where because of some bias, the RR in one of the two rwCC may be lower than the true value of 50%. As expected, when the response difference between the rwCCs increases, the Type 1 error, observed response difference, and probability of detecting heterogeneity also increase. However, given the low sample sizes in the rwCC, power is relatively low for detecting heterogeneity even at the relaxed pvalue threshold of 0.10. A 30% RR difference is detected in only 64% of trials; the Type I error rate is 30% in the remaining 36% of trials where no heterogeneity is detected, leading to an overall Type I error rate of nearly 11% (30% * 36%). Thus, although heterogeneity testing successfully identifies some highbias scenarios, in many reasonable scenarios it will be underpowered, leading to Type 1 error inflation (albeit less than if we always pool).
In these simulations, the RR in the SAT and in the first rwCC is 50%, and the RR in the second rwCC varies. The pooling probability is the probability that the heterogeneity test fails to reject (if it rejects, we do not pool). The Type I error rate is calculated under two different decision rules: we pool only when the heterogeneity test fails to reject (‘Conditional on pooling’), or we always pool (‘If always pool’). The average estimated response difference is also calculated (last column) and compared to the expected response difference when always pooling. Monte Carlo standard errors are in parentheses, for 100,000 replicates.
Note that in each row of Table 2, the Type I error and response difference conditional on pooling are approximately equal to those of if we always pooled. This counterintuitive finding is because the average response across the two rwCC depends only weakly on whether the observed RRs are similar to each other or not, conditional on the true RRs in the two rwCC. For example, if the true RRs are 50% in rwCC #1 and 40% in rwCC #2, a typical case when similar RRs would be observed might be a rate of 45% in both datasets, while a typical case when different RRs would be observed might be a RR of 50% in one dataset and 40% in the other. In either case, the average rwCC RR is 45%, leading to Type I error inflation and bias in estimation of the response rate difference.
Scenario 2: impact of heterogeneity testing when treatment improves outcomes
In Scenario 2, we examine a situation where there should be a true RR difference between the SAT and both rwCCs, but the RR in one of the rwCCs is biased upward leading to Type 2 error. As shown in Table 3, the low power to detect heterogeneity means that the datasets will be inappropriately pooled often, leading to potentially severe underestimation of the treatment effect. For example, with a 20% RR difference between the rwCC, datasets will be pooled 67% of the time, and the average estimate of the treatment effect will be 20% when pooling, a third lower than the true RR difference of 30%.
In these simulations, the RR in the SAT is 80% and in the first rwCC is 50%, and the RR in the second rwCC varies. Monte Carlo standard errors are in parentheses, for 1,000,000 replicates.
Discussion
This article makes several contributions to the existing literature on data pooling [6, 21, 23, 33]. First, we present a framework for data pooling focused on RWD, and specifically rwCC analyses. The three steps of this framework are: prespecification, dataset selection, and outcomes analyses (including analysis of heterogeneity).
Second, we evaluate how best to detect heterogeneity in the rwCC use case and demonstrate that standard approaches to heterogeneity assessment are biased. To address this gap, we describe an Adjusted Cochran’s Q test, similar to existing methods in the literature for metaanalysis of dependent results, which accounts for correlation in the treatment estimates given the common SAT, and matches the statistical power of IPD heterogeneity testing. In scenarios in which all IPD can be combined, heterogeneity can be assessed either by directly comparing the rwCCs, or using the Adjusted Q. The former is the preferable approach; however, a use case in which only the Adjusted Q is feasible is when IPD may be available but not fully shareable across study sites. This can occur when IPD from the SAT can be combined with IPD from each rwCC, but IPD from different rwCCs cannot be combined (i.e., if one or more of the rwCCs is derived from a registry) [34]. If the sponsor wished to combine these results, this could be accomplished using a fixed effect metaanalysis approach, although adjustments (e.g., similar to those proposed for the Adjusted Q test) should be made to standard errors to account for nonindependence of the results being metaanalyzed (see the Appendix for further discussion).
Third, our simulation results highlight the limitations of quantitative heterogeneity testing in protecting against Type 1 or Type 2 error especially in the presence of small sample sizes. While heterogeneity testing does identify scenarios with higher bias, its lack of statistical power when there are only a few rwCC (especially with the low sample sizes that would lead a sponsor to consider pooling) make it an insufficient mitigation on its own. This is in line with previous research illustrating the dangers of twostep procedures where fixed effects or random effects metaanalysis is selected based on a preliminary test of heterogeneity [35, 36]. We note that if the goal of pooling RWD were to achieve greater generalizability rather than a rwCC use case, heterogeneity would be expected (and potentially even desirable), and thus, testing would not be needed. The latter use case is out of scope for the current manuscript.
Together, these results underscore the importance of careful selection of datasets (‘what to pool’) based on clinical/data considerations, to minimize heterogeneity prior to any statistical testing. This remains an area in which further research is needed, given the multitude of considerations and difficulty testing the underlying assumptions [14]. Subsequent heterogeneity test results should be used descriptively to contextualize study findings. If statistical heterogeneity is detected or if qualitatively large but nonsignificant differences are seen, the investigators should present datasetspecific results, and conduct sensitivity analyses to investigate potential sources of heterogeneity [23].
Conclusions
In conclusion, we have provided a methodological framework for pooling separate RWD sources, particularly applicable to the use case of designing a rwCC for a single arm trial in a rare subpopulation. Additionally, we provide simulations demonstrating the performance of a heterogeneity assessment methodology and describe considerations for this essential component of pooling. As researchers and regulators seek to benefit patients with novel insights from RWD pooled from multiple available sources, we hope the methods presented herein can enhance the rigor of such analyses.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Abbreviations
 RWD:

Realworld data
 rwCC:

Realworld comparator cohort
 SAT:

Singlearm clinical trial
 IPD:

Individual participant data
 RR:

Response rate
References
Vega DM, Nishimura KK, Zariffa N, Thompson JC, Hoering A, Cilento V, et al. Changes in Circulating Tumor DNA Reflect Clinical Benefit Across Multiple Studies of Patients With Non–SmallCell Lung Cancer Treated With Immune Checkpoint Inhibitors. JCO Precis Oncol. 2022;6:e2100372.
Minchom A, Viteri S, Bazhenova L, Gadgeel SM, Ou SI, Trigo J, et al. Amivantamab compared with realworld therapies in patients with advanced nonsmall cell lung cancer harboring EGFR exon 20 insertion mutations who progressed after platinumbased chemotherapy. Lung Cancer (Amsterdam, Netherlands). 2022;168:74–82.
Patel J, Ling YL, Tang J, Nakasato A, Kang B, Balu S. PCN239 Development of the NOBLE (NOvartis Braf+ meLanoma patients obsErvational) Dataset: a Novel Approach to Creating a Longitudinal Melanoma Data Source by Harmonizing Multiple Electronic Health Record Databases. Value in health. 2021;24:S64–5.
Jagannath S, Lin Y, Goldschmidt H, Reece D, Nooka A, Senin A, et al. KarMMaRW: comparison of idecabtagene vicleucel with realworld outcomes in relapsed and refractory multiple myeloma. Blood Cancer J. 2021;11(6):116.
Blettner M, Sauerbrei W, Schlehofer B, Scheuchenpflug T, Friedenreich C. Traditional reviews, metaanalyses and pooled analyses in epidemiology. Int J Epidemiol. 1999;28(1):1–9.
Tierney JF, Vale C, Riley R, Smith CT, Stewart L, Clarke M, et al. Individual Participant Data (IPD) Metaanalyses of Randomised Controlled Trials: Guidance on Their Use. PLoS Med. 2015;12(7): e1001855.
Deeks JJ, Higgins JPT, Altman DG, on behalf of the Cochrane Statistical, Methods Group. Analysing data and undertaking metaanalyses. Cochrane Handbook for Systematic Reviews of Interventions. 2019. p. 241–84.
Integrated Summary of Effectiveness; Guidance for Industry; Availability. The Federal Register / FIND. 2015;80(195):60916.
European Agency for the Evaluation of Medicinal Products. Points to consider on application with 1. Metaanalysis; 2. One pivotal study. London: European Agency for the Evaluation of Medicinal Products; 2001.
Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in metaanalyses. BMJ. 2003;327(7414):557–60.
Hersh WR, Weiner MG, Embi PJ, Logan JR, Payne PR, Bernstam EV, et al. Caveats for the use of operational electronic health record data in comparative effectiveness research. Med Care. 2013;51(8 Suppl 3):30.
Bastarache L, Brown JS, Cimino JJ, Dorr DA, Embi PJ, Payne PRO, et al. Developing realworld evidence from realworld data: Transforming raw data into analytical datasets. Learn Health Syst. 2021;6(1):e10293.
Bangdiwala SI, Bhargava A, O’Connor DP, Robinson TN, Michie S, Murray DM, et al. Statistical methodologies to pool across multiple intervention studies. Transl Behav Med. 2016;6(2):228–35.
Demets DL. Methods for combining randomized clinical trials: strengths and limitations. Stat Med. 1987;6(3):341–50.
Kahan BC, Forbes G, Cro S. How to design a prespecified statistical analysis approach to limit phacking in clinical trials: the PreSPEC framework. BMC Med. 2020;18(1):253.
Moher D, Liberati A, Tetzlaff J, Altman DG. Preferred Reporting Items for Systematic Reviews and MetaAnalyses: The PRISMA Statement. J Clin Epidemiol. 2009;62(10):1006–12.
Mahendraratnam N, Silcox C, Mercon K, Kroetsch A, Romine M, Harrison N, et al. Determining RealWorld Data’s Fitness for Use and the Role of Reliability. 2019.
Seeger JD, Davis KJ, Iannacone MR, Zhou W, Dreyer N, Winterstein AG, et al. Methods for external control groups for single arm trials or longterm uncontrolled extensions to randomized clinical trials. Pharmacoepidemiol Drug Saf. 2020;29(11):1382–92.
Tan WK, Segal BD, Curtis MD, Baxi SS, Capra WB, GarrettMayer E, et al. Augmenting control arms with RealWorld Data for cancer trials: Hybrid control arm methods and considerations. arXiv preprint arXiv:2108.07335 2021.
Chen J, Ho M, Lee K, Song Y, Fang Y, Goldstein BA, et al. The Current Landscape in Biostatistics of RealWorld Data and Evidence: Clinical Study Design and Analysis. Stat Biopharm Res. 2023;15(10):29–42.
Lesko CR, Jacobson LP, Althoff KN, Abraham AG, Gange SJ, Moore RD, et al. Collaborative, pooled and harmonized study designs for epidemiologic research: challenges and opportunities. Int J Epidemiol. 2018;47(2):654–68.
Datavant. (Datavant) Overview of Datavant's DeIdentification and Linking Technology for Structured Data. https://datavant.com/wpcontent/uploads/dlm_uploads/2018/09/WhitePaper_DeIdentifyingandLinkingStructuredData.pdf. Accessed 22 Aug 2023.
Thompson SG. Why sources of heterogeneity in metaanalysis should be investigated. BMJ. 1994;309(6965):1351–5.
Fisher DJ, Copas AJ, Tierney JF, Parmar MK. A critical review of methods for the assessment of patientlevel interactions in individual participant data metaanalysis of randomized trials, and guidance for practitioners. J Clin Epidemiol. 2011;64(9):949–67.
Petitti DB. Metaanalysis, decision analysis, and costeffectiveness analysis. 2. Aufl. ed. New York: Oxford University Press; 2000.
Friede T, Röver C, Wandel S, Neuenschwander B. Metaanalysis of few small studies in orphan diseases. Res Synth Methods. 2017;8(1):79–91.
Sera F, Armstrong B, Blangiardo M, Gasparrini A. An extended mixedeffects framework for metaanalysis. Stat Med. 2019;38(29):5429–44.
Gleser LJ, Olkin I. Stochastically dependent effect sizes. New York: Russell Sage Foundation; 1994. p. 339–55.
Graybill FA. Theory and application of the linear model. 1. printing ed. North Scituate, Mass: Duxbury; 1976.
Kulinskaya E, Dollinger MB, Bjørkestøl K. On the moments of Cochran’s Q statistic under the null hypothesis, with application to the metaanalysis of risk difference. Res Synth Methods. 2011;2(4):254–70.
Joshi M, Pustejovsky JE, Beretvas SN. Cluster wild bootstrapping to handle dependent effect sizes in metaanalysis with a small number of studies. Res Synth Methods. 2022;13(4):457–77.
Viele K, Berry S, Neuenschwander B, Amzal B, Chen F, Enas N, et al. Use of historical control data for assessing treatment effects in clinical trials. Pharm Stat. 2014;13(1):41–54.
Higgins JPT, Thomas J, Chandler J, et al editors. Cochrane Handbook for Systematic Reviews of Interventions. 6.2nd ed. Newark: John Wiley & Sons, Ltd; 2021.
Costa LJ, Lin Y, Cornell RF, Martin T, Chhabra S, Usmani SZ, et al. Comparison of Ciltacel, an AntiBCMA CART Cell Therapy, Versus Conventional Treatment in Patients With Relapsed/Refractory Multiple Myeloma. Clin Lymphoma Myeloma Leuk. 2022;22(5):326–35.
Shuster JJ, Jones LS, Salmon DA. Fixed vs random effects metaanalysis in rare event studies: the rosiglitazone link with myocardial infarction and cardiac death. Stat Med. 2007;26(24):4375–85.
Borenstein M, Hedges L, Higgins J, Rothstein H. Introduction to MetaAnalysis. 1st ed. Chichster: Wiley; 2009.
Acknowledgements
The authors would like to thank Hannah Gilham and Darren Johnson of Flatiron Health, Inc. for publication management support, and Shrujal Baxi and Aracelis Torres of Flatiron Health, Inc. for their scientific input and support of this research.
Funding
This study was collaboratively sponsored by Flatiron Health, Inc., which is an independent member of the Roche group, and Janssen.
Author information
Authors and Affiliations
Contributions
Preparation: DB, TR, JP, MS, OH. Analysis: DB, OH. Writing and interpretation: DB, TR, JP, MS, OH.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not Applicable.
Consent for publication
Not Applicable.
Competing interests
At the time of this study, TR, MS, and OH report employment at Flatiron Health, Inc., which is an independent member of the Roche Group, and stock ownership in Roche. DB and JP are employees of Johnson & Johnson and hold Johnson & Johnson stock or stock options. OH reports current employment at Regeneron Pharmaceuticals and stock ownership in Roche and Regeneron Pharmaceuticals (performed study while previously employed at Flatiron Health). MS reports current employment at Calico Life Sciences and stock ownership in Roche and Calico Life Sciences (performed study while previously employed at Flatiron Health). Flatiron Health, Inc. and Janssen were involved in the study design and administration, as well as collection, management, analysis, and interpretation of the data. The parent company of Flatiron Health, Inc. (Roche group) was not involved in the design or any aspect of the conduct of the study, nor were other Flatiron Health or Janssen employees aside from the authors. Daniel Backenroth and Olivier Humblet had full access to the data in the study and are accountable for the integrity of this report. All the authors collaborated in the preparation, review, and approval of the manuscript and agreed to the decision to submit this manuscript for publication.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Backenroth, D., Royce, T., Pinheiro, J. et al. Considerations for pooling realworld data as a comparator cohort to a single arm trial: a simulation study on assessment of heterogeneity. BMC Med Res Methodol 23, 193 (2023). https://doi.org/10.1186/s12874023020027
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12874023020027
Keywords
 Pooling
 Metaanalysis
 Realworld data
 Oncology
 Small population
 Heterogeneity
 Singlearm trial
 Realworld comparator cohort