Next Article in Journal
“I Have Eight Different Files at Eight Different Places”: Perspectives of Youths and Their Family Caregivers on Transitioning from Pediatric to Adult Rehabilitation and Community Services
Previous Article in Journal
Predictability of Radiologically Measured Psoas Muscle Area for Intraoperative Hypotension in Older Adult Patients Undergoing Femur Fracture Surgery
Previous Article in Special Issue
Multidisciplinary Guidelines for the Rational Use of Topical Non-Steroidal Anti-Inflammatory Drugs for Musculoskeletal Pain (2022)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methods for the Inclusion of Real-World Evidence in a Rare Events Meta-Analysis of Randomized Controlled Trials

1
Chinese Evidence-Based Medicine Center and MAGIC China Center, West China Hospital, Sichuan University, Chengdu 610041, China
2
NMPA Key Laboratory for Real World Data Research and Evaluation in Hainan, Chengdu 610041, China
3
Sichuan Center of Technology Innovation for Real World Data, Chengdu 610041, China
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2023, 12(4), 1690; https://doi.org/10.3390/jcm12041690
Submission received: 29 November 2022 / Revised: 10 February 2023 / Accepted: 17 February 2023 / Published: 20 February 2023
(This article belongs to the Special Issue Advances in Clinical Trials and Evidence Synthesis (CTES))

Abstract

:
Background: Many rare events meta-analyses of randomized controlled trials (RCTs) have lower statistical power, and real-world evidence (RWE) is becoming widely recognized as a valuable source of evidence. The purpose of this study is to investigate methods for including RWE in a rare events meta-analysis of RCTs and the impact on the level of uncertainty around the estimates. Methods: Four methods for the inclusion of RWE in evidence synthesis were investigated by applying them to two previously published rare events meta-analyses: the naïve data synthesis (NDS), the design-adjusted synthesis (DAS), the use of RWE as prior information (RPI), and the three-level hierarchical models (THMs). We gauged the effect of the inclusion of RWE by varying the degree of confidence placed in RWE. Results: This study showed that the inclusion of RWE in a rare events meta-analysis of RCTs could increase the precision of the estimates, but this depended on the method of inclusion and the level of confidence placed in RWE. NDS cannot consider the bias of RWE, and its results may be misleading. DAS resulted in stable estimates for the two examples, regardless of whether we placed high- or low-level confidence in RWE. The results of the RPI approach were sensitive to the confidence level placed in RWE. The THM was effective in allowing for accommodating differences between study types, while it had a conservative result compared with other methods. Conclusion: The inclusion of RWE in a rare events meta-analysis of RCTs could increase the level of certainty of the estimates and enhance the decision-making process. DAS might be appropriate for inclusion of RWE in a rare event meta-analysis of RCTs, but further evaluation in different scenarios of empirical or simulation studies is still warranted.

1. Introduction

Well-conducted randomized controlled trials (RCTs) have been considered the gold standard, and as such are often used in meta-analyses to evaluate the effects of healthcare interventions [1]. For rare events meta-analyses of RCTs, the outcomes of the included trials may be very sparse, with no events even being observed in some trials [2,3], resulting in lower statistical power [4]. Recently, there has been a growing interest in the use of real-world evidence (RWE) in clinical drug assessments and health care evaluations [5,6,7]. RWE is clinical evidence from multiple sources outside the typical clinical research setting, including electronic health records and billing databases [5]. Regulatory bodies such as the US Food and Drug Administration (FDA) have been increasingly using RWE to monitor adverse post-market safety events and make regulatory decisions [5,8]. An important potential advantage of RWE is that it can represent real-life clinical practice across a broader and more diverse distribution of patients than RCTs [9]. The results from RWE studies, such as those obtained from epidemiology databases, health surveys, and electronic medical records, typically include larger sample sizes and have longer follow-up periods, increasing the probability of finding rare events [10,11,12].
RCTs and RWE studies are potentially valuable evidence sources for assessing the effect of rare events, and integrating RWE studies and RCTs to assess the effects of rare events may help increase the overall level of certainty of evidence [13,14]. However, the inclusion of RWE in a rare events meta-analysis of RCTs is not a straightforward issue, as the estimates obtained from RWE studies may be subject to bias due to the potential selection and information biases of real-world data [15]. Therefore, RCTs and RWE studies should not be directly integrated without considering the bias of the RWE studies. Generalized evidence synthesis statistical approaches for integrating RCTs and RWE studies can provide not only an integrated result but also a quantitative analysis of the influence of RWE on the integrated evidence under different bias assumptions for the RWE studies, which is more attractive for decision-making [16,17,18,19]. Relevant methods have been published in the statistical literature but have not been systematically verified or widely used in empirical analyses [16].
Therefore, the aim of this study is to investigate whether these statistical methods, including naïve data synthesis (NDS), design-adjusted synthesis (DAS), the use of RWE as prior information (RPI), and three-level hierarchical models (THMs), could be used for integrating RWE in a rare events meta-analysis of RCTs by considering the potential inherent biases of RWE studies and their impact on the level of uncertainty around the estimates.

2. Materials and Methods

2.1. Study Design Overview

In this study, we used four approaches to combine RCTs and RWE studies [18,20]: (1) NDS, where data from all studies, regardless of the study design, were directly combined; (2) DAS, where RCTs and RWE studies were synthesized and where the information from the RWE studies was adjusted to reflect the confidence in the study findings; (3) the RPI, which included RWE via the prior for the treatment effect and combined this with a likelihood based only on the data from RCTs; (4) the THM, which was used to simultaneously model the between-study heterogeneity of treatment effects within each study design (RCT or RWE studies) and across study designs. The methodologies were applied to two recently published meta-analyses concerning the risk of diabetic ketoacidosis (DKA) among patients using sodium/glucose cotransporter 2 (SGLT-2) inhibitors compared with active comparators [21] and the effectiveness of 23-valent pneumococcal polysaccharide vaccine (PPV23) vaccination against invasive pneumococcal disease (IPD) in elderly patients [22]. Here, we describe how these methods can be performed and how the bias of the RWE studies can be considered. Finally, we discuss the pros and cons of each method.

2.2. Illustrative Example Dataset

The first example was a recent meta-analysis conducted by Falkenhorst et al. [22]. In this meta-analysis, the authors used the evidence from RCTs and RWE studies to investigate the effectiveness of PPV23 vaccination against IPD and pneumococcal pneumonia in adults aged ≥ 60 years. Here, we focus on IPD. This study included four RCTs, five cohort studies, three case-control studies, and five case–case studies. We excluded one RCT (for reporting no events in both treatment groups) and five case–case studies (for IPD not being reported) in our study. The outcomes of these studies are displayed in Table 1. The authors assessed the risk of bias of individual studies using the Cochrane Risk of Bias tool [23] for RCTs and the Newcastle-Ottawa Scale for observational studies [24]. A total of 16 IPD events were reported across three RCTs, which randomized a total of 41,992 patients, and the incidence rate was 0.04%. We observed 165 IPD events among 575,454 patients in five cohort studies.
The second example showed how to integrate RCTs and RWE studies when assessing the effects of rare adverse events. A typical example of rare adverse events is the risk of DKA caused by SGLT-2 inhibitors. The data used in this study were obtained from the study conducted by Alkabbani et al. [21]. This study included twelve placebo-controlled RCTs, seven active comparator RCTs, and seven active comparator RWE studies. All were retrospective propensity score-matched cohort studies. Our primary concern was whether SGLT-2 inhibitors increase the risk of DKA compared with active comparators. We excluded one RWE study because its control was not an active comparator [36]. In this study, the authors assessed the risk of bias in each study using the checklist proposed by Downs et al. [37] for RCTs and RWE studies. Table 2 shows the basic characteristics of the included RCTs and RWE studies. A total of 8 DKA events were reported across all RCTs, which randomized 8100 patients, and the incidence rate was 0.1%. We observed 2016 DKA events among 1,072,992 patients across all RWE studies.

2.3. Methods for Incorporating RWE in a Rare Events Meta-Analysis of RCTs

In the following, we describe how these methods can be performed and how the bias of the RWE studies can be considered. We used j = 1, 2, … k to denote study j. The odds ratio (OR) metric was used as the effect measure because rare events are often binary outcomes. These methods are also suitable for other measures, such as the relative effect (RR) and risk difference (RD).

2.3.1. The Naïve Data Synthesis

The NDS is identical to the conventional random effects meta-analysis model, namely the normal–normal hierarchical model. The model was then written as follows:
d j ~ N ( θ j , s j 2 )
θ j ~ N ( θ , τ 2 )
where d j and s j 2 are the study-specific empirical effect and the within-study variances for study j, respectively. Both the treatment effect ( d i ) and variance ( s i 2 ) are on the log OR scale. Here, θ j denotes the true underlying log-OR for the jth study, θ is the summary log-OR, and τ 2 is the between-study variance. We assigned a weakly informative prior for θ , i.e., θ ~ N ( 0 , 2.82 2 ) , following Günhan et al. [49]. For the heterogeneity parameter τ , we used a half-normal prior with a scale of 0.5 ( τ ~ H N ( 0.5 ) ) , following Friede et al. [50].

2.3.2. The Design-Adjusted Synthesis

The DAS approach was originally used for integrating RCTs and RWE studies in network meta-analysis. We simplified this model to fit a pairwise meta-analysis. This model gives less weight to RWE studies compared to RCTs. This down-weighting is achieved by inflating the variance of the mean effect of RWE studies. The general conventional meta-analysis model was extended by including a down weighting factor ω j in the variance of the estimated treatment effect of study j:
d j ~ N ( θ j , s j 2 ω j )
θ j ~ N ( θ , τ 2 )
where ω j is the variance inflation factor for study j, and the other parameters are the same as the DAS; ω j was equal to 1 when the study type was an RCT, since the estimates obtained from RCTs are often considered unbiased. For the RWE studies, the ω j values ranged from zero to one, and the variance of the treatment effect was inflated so that the weight in the meta-analysis was decreased. We evaluated varying levels of trust in the RWE studies, with ω j ~beta (4, 1), ω j ~beta (1.5, 1), and ω j ~beta (0.25, 1) standing for high, medium, and lower levels of confidence in RWE studies. The prior specifications for θ and τ were identical to those of the NDS.

2.3.3. The Use of Real-World Evidence as Prior Information

This approach included RWE via the prior for the treatment effect and combined this with a likelihood based only on the data from RCTs. We first synthesized the RWE studies to generate θ ^ R W E with a corresponding variance V ^ R W E . Then, we centered the informative prior for the population mean ( θ ) on θ ^ R W E but use an inflated variance of the mean effect of RWE; that is, θ ~ N ( θ ^ R W E , V ^ R W E ω ) , where ω is the variance inflation factor. The prior specification for τ is identical to the NDS. We also evaluated varying levels of trust in the RWE studies, with ω j ~beta (4, 1), ω j ~beta (1.5, 1), and ω j ~beta (0.25, 1).

2.3.4. The Three-Level Hierarchical Model

For the THM, the evidence from the RCTs and RWE studies were modeled separately at the within-study and within-design levels, and then the estimates for RCTs and RWE studies were combined in an overall measure of the treatment effect using the random-effects model, assuming i = 1, 2, where 1 represents the RCTs and 2 represents the RWE studies. The model is presented as follows:
d j , i ~ N ( θ j , i , s j , i 2 )
θ j , i ~ N ( θ i , τ i 2 )
θ i ~ N ( θ , τ 2 )
where d j , i and s j 2 are the study-specific empirical effect and the within-study variances for study j and design i, respectively; θ j , i denotes the true underlying effect for the jth study and ith design; θ i and τ i 2 are the summary effect and the variance for the ith design, respectively; and the other parameters are the same as the NDS. The RWE studies could be down-weighted by inflating the variance; that is, θ 2 ~ N ( θ , τ 2 ω ) . The prior specifications for θ and τ were identical to those of the NDS. For the study-design variance τ i 2 , we also assigned a half-normal prior with a scale of 0.5. The prior specifications for ω were also set in three different scenarios, namely ω j ~beta (4, 1), ω j ~beta (1.5, 1), and ω j ~beta (0.25, 1).

2.4. Implementation and Model Fit

For RCTs, the effect of the single-zero event trial was estimated by the method proposed by Greenland et al. [51]. All models were performed using Hamiltonian Monte Carlo algorithms, as implemented in the RStan package (version 2.21.2) in the R statistical environment (version 4.0.3, R Foundation for Statistical Computing, Vienna, Austria). As a comparison, we estimated the results of RCTs using a Bayesian random effects meta-analysis. We fitted four chains for each model, each with 5000 iterations, the first half of which were considered warmups and discarded. Convergence was judged to have occurred when R ^ (the potential scale reduction factor) was no greater than 1.1 for all parameters [52].

3. Results

3.1. Effectiveness of PPV23 Vaccination against IPD in Elderly Patients

The Bayesian meta-analysis including only RCTs had a posterior median of 0.29 and a 95% credibility interval (CrI) [0.06, 1.45], which indicated that based on the RCT evidence, the PPV23 is not effective against IPD. The posterior median and 95% CrI based on the NDS were 0.42 and (0.29, 0.58), which indicated the effectiveness of PPV23 against IPD in elderly patients. Table 3 shows the posterior median and 95% CrI of the effectiveness of PPV23 vaccination against IPD in elderly patients using the DAS, RPI, and THM. The models, except for the THM, consistently showed the effectiveness of PPV23 vaccination against IPD in elderly patients. By comparing the length of the intervals between the NDS and the other three models, we found that the length of the interval of the NDS was much smaller than that of the other three models. We also observed that the length of the interval in the estimates obtained from the four models was much smaller than the estimates from only RCTs.

3.2. Risk of DKA among Users Receiving SGLT-2 Inhibitors Versus Active Comparators

We did not observe an increased risk of DKA when the results were pooled from the RCTs (OR = 1.02, 95% CrI: 0.37–2.86), while we observed an increased risk of DKA when directly combining the RCTs and RWE studies using NDS (OR = 1.58, 95% CrI: 1.22–2.04). The estimates of each model under different variance inflation factors, which represent different levels of confidence in RWE studies, are presented in Table 4. The results of the different models were inconsistent. The DAS showed that SGLT2 inhibitors were associated with an increased risk of DKA under different confidence levels placed in RWE studies. For the RPI, we observed an increased risk when we placed a high level of confidence in RWE studies (OR = 1.56, 95% CrI: 1.12–2.12); however, we did not observe an increased risk when we placed a medium level of confidence (OR = 1.52, 95% CrI: 0.98–2.19) or a low level (OR = 1.41, 95% CrI: 0.56–2.45) of confidence in RWE studies. We also did not observe an increased risk of DKA by using the THM, even though we placed a high level of confidence in RWE studies (OR = 1.42, 95% CrI: 0.50–3.75). As with the first example, the interval lengths of the four model estimates were smaller than the outcomes estimated based on RCTs alone.

4. Discussion

In this study, we discussed approaches for incorporating RWE in a rare events meta-analysis of RCTs, with particular interest in the effectiveness of PPV23 against IPD in elderly patients and the risk of DKA among patients using SGLT-2 inhibitors compared with active comparators. We did not observe the effectiveness of PPV23 and the risk of DKA among patients using SGLT-2 inhibitors when the data were only from RCTs. In the effectiveness example of PPV 23, all methods, except for the THM, found that the inclusion of RWE studies in rare events meta-analysis of RCTs could complement the evidence of the findings of the RCT-only analysis. However, for the risk of DKA among patients using SGLT-2 inhibitors, the results of the different models were inconsistent. For instance, for the RPI, an obvious finding was that we did not observe an increased risk when low and medium confidence was placed in the RWE studies. The results showed that the inclusion of RWE studies during the evidence synthesis process could increase the certainty of the estimates when the rare events meta-analysis of RCTs could provide enough evidence, but this depended on the method of inclusion and the confidence level placed in the RWE studies.
We do not recommend using NDS as the main method of analysis to combine rare events data. A recent scoping review found that the NDS approach was the most frequently used in empirical analyses [16]. However, NDS ignores differences in study designs and cannot consider the potential bias of RWE studies [20]. Furthermore, compared with RCTs, the results of RWE studies often show a large effect because of some uncontrolled confounding bias factors [53], and their interval estimates are much smaller because the events and the sample size are usually much larger [21]. Therefore, with the inclusion of RWE studies in a rare events meta-analysis of RCTs using NDS, not only can the bias of RWE studies not be adjusted but it would also give a larger weight than that of RCTs. Our two illustrative examples were also confirmed. This indicated that addressing the bias of RWE studies leads to more informed decision-making when integrating RCTs and RWE to assess the effects of rare events.
For the DAS, the data from the RWE studies are down-weighted based on the confidence levels of their credibility. An obvious feature of this model is that it needs to separate assessments of bias for each RWE study, which could be done by eliciting expert opinions regarding the bias parameters [54]. Another useful method is to gauge the impact on the model estimates by varying the amount of confidence placed in RWE studies [20]. For the RPI, there are two key differences with the DAS: one is in the estimation of heterogeneity that is performed separately for RCTs and RWE studies, and the other is that it only needs to set the magnitude of possible bias for the total RWE. Therefore, this model may be more attractive in empirical analyses, because one can not only set the prior distribution of heterogeneity parameters according to the characteristics of RCT and RWE studies (such as patient follow-up and statistical analysis methods), but can also set the magnitude of bias based on the risk of bias evaluation results for the whole RWE rather than for each RWE study.
Although the advantage of the THM is that it can obtain the estimates of each study design to increase its precision, the results of the total combined effect are conservative. In these two examples, compared to the results of the other models, the most notable difference was that the interval length was much larger when using the THM, even though we placed a high level of confidence in the RWE studies. This was because the THM explicitly considers the heterogeneity between study designs, allowing for additional variability across studies [20]. The main goal of the inclusion of RWE studies in a rare events meta-analysis of RCTs is to increase the power in testing whether the true effect exists. Rare events are often associated with safety results such as serious adverse events, indicating that the THM as the main method of analysis may not be suitable and may increase the chances of patient exposure to unnecessary danger [55]. The estimates can be improved by using expert opinions and meta-epidemiology to set a WIP for the combined effects [56]. For instance, constructing a WIP works via the consideration of the prior expected range of the treatment effect [57] and the heterogeneity values [58].
There were some limitations of this study that need to be recognized. First, this study only utilized two illustrative examples, the results may differ in other outcomes, and broader analyses in further empirical analyses are necessary. Second, this study only investigated methods for the inclusion of RWE in rare event meta-analyses of RCTs, and the performance of these methods was not explored. There is no clear rule for choosing the method to combine RCTs and RWE studies; when different methods lead to contradicting conclusions, researchers can choose those methods that lead to the desired outcome. Thus, a further evaluation of these methods in different scenarios, including the use of comprehensive simulation studies, is warranted. Third, meta-regression was not considered in this study. Although meta-regression can explain some between-study heterogeneity, it may be limited by the information on available covariates of each study or the number of studies in a meta-analysis. Fourth, although we used the OR as the effect measure, these methods are generalizable to other measures of association frequently used in meta-analyses, such as the RR, RD, and weighted mean differences. The performance of these methods in other measures needs to be further explored.

5. Conclusions

In summary, including RWE in a rare events meta-analysis has the potential to corroborate findings from RCTs, increase precision, and enhance the decision-making process, but this depends on the method of inclusion and the assumption for the magnitude of bias risk of RWE. The statistical performance of these approaches requires further evaluation in different scenarios of empirical or simulation studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm12041690/s1, supplementary file: R code for the methods included this study.

Author Contributions

X.S., L.L. and M.Y. conceived and designed the study. X.S., M.Y. and L.L. acquired the funding. X.S. and M.Y. drafted the manuscript. M.Y., Y.W. and F.M. conducted the data analysis. X.S., M.Y., L.L., F.M., Y.W. and K.Z. critically revised the article. X.S. is the guarantor. All authors have read and agreed to the published version of the manuscript.

Funding

The National Natural Science Foundation of China (Grant No. 72204173, 82274368, and 71904134), National Science Fund for Distinguished Young Scholars (Grant No. 82225049), Sichuan Provincial Central Government Guides Local Science and Technology Development Special Project (Grant No. 2022ZYD0127), Fundamental Research Funds for the Central Public Welfare Research Institutes (Grant No. 2020YJSZX-3).

Data Availability Statement

The computing code used in the study is provided in the Supplementary Materials.

Conflicts of Interest

The authors declare no competing interest.

References

  1. Zabor, E.C.; Kaizer, A.M.; Hobbs, B.P. Randomized Controlled Trials. Chest 2020, 158, S79–S87. [Google Scholar] [CrossRef] [PubMed]
  2. Xu, C.; Li, L.; Lin, L.; Chu, H.; Thabane, L.; Zou, K.; Sun, X. Exclusion of studies with no events in both arms in meta-analysis impacted the conclusions. J. Clin. Epidemiol. 2020, 123, 91–99. [Google Scholar] [CrossRef]
  3. Hodkinson, A.; Kontopantelis, E. Applications of simple and accessible methods for meta-analysis involving rare events: A simulation study. Stat. Methods Med. Res. 2021, 30, 1589–1608. [Google Scholar] [CrossRef]
  4. Jia, P.; Lin, L.; Kwong, J.; Xu, C. Many meta-analyses of rare events in the Cochrane Database of Systematic Reviews were underpowered. J. Clin. Epidemiol. 2020, 131, 113–122. [Google Scholar] [CrossRef] [PubMed]
  5. Sherman, R.E.; Anderson, S.A.; Dal Pan, G.J.; Gray, G.W.; Gross, T.; Hunter, N.L.; LaVange, L.; Marinac-Dabic, D.; Marks, P.W.; Robb, M.A.; et al. Real-World Evidence—What Is It and What Can It Tell Us? N. Engl. J. Med. 2016, 375, 2293–2297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sun, X.; Tan, J.; Tang, L.; Guo, J.J.; Li, X. Real world evidence: Experience and lessons from China. BMJ 2018, 360, j5262. [Google Scholar] [CrossRef] [Green Version]
  7. Brown, J.P.; Wing, K.; Evans, S.J.; Bhaskaran, K.; Smeeth, L.; Douglas, I.J. Use of real-world evidence in postmarketing medicines regulation in the European Union: A systematic assessment of European Medicines Agency referrals 2013–2017. BMJ Open 2019, 9, e028133. [Google Scholar] [CrossRef]
  8. Bolislis, W.R.; Fay, M.; Kühler, T.C. Use of Real-world Data for New Drug Applications and Line Extensions. Clin. Ther. 2020, 42, 926–938. [Google Scholar] [CrossRef]
  9. Wu, J.; Wang, C.; Toh, S.; Pisa, F.E.; Bauer, L. Use of real-world evidence in regulatory decisions for rare diseases in the United States-Current status and future directions. Pharm. Drug Saf. 2020, 29, 1213–1218. [Google Scholar] [CrossRef]
  10. Douros, A.; Lix, L.M.; Fralick, M.; Dell’Aniello, S.; Shah, B.R.; Ronksley, P.E.; Tremblay, É.; Hu, N.; Alessi-Severini, S.; Fisher, A.; et al. Sodium-Glucose Cotransporter-2 Inhibitors and the Risk for Diabetic Ketoacidosis: A Multicenter Cohort Study. Ann. Intern. Med. 2020, 173, 417–425. [Google Scholar] [CrossRef]
  11. Li, L.; Li, S.; Deng, K.; Liu, J.; Vandvik, P.O.; Zhao, P.; Zhang, L.; Shen, J.; Bala, M.M.; Sohani, Z.N.; et al. Dipeptidyl peptidase-4 inhibitors and risk of heart failure in type 2 diabetes: Systematic review and meta-analysis of randomised and observational studies. BMJ 2016, 352, i610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hutton, B.; Joseph, L.; Fergusson, D.; Mazer, C.D.; Shapiro, S.; Tinmouth, A. Risks of harms using antifibrinolytics in cardiac surgery: Systematic review and network meta-analysis of randomised and observational studies. BMJ 2012, 345, e5798. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cuello-Garcia, C.A.; Santesso, N.; Morgan, R.L.; Verbeek, J.; Thayer, K.; Ansari, M.T.; Meerpohl, J.; Schwingshackl, L.; Katikireddi, S.V.; Brozek, J.L.; et al. GRADE Guidance 24. Optimizing the integration of randomized and non-randomized studies of interventions in evidence syntheses and health guidelines. J. Clin. Epidemiol. 2022, 142, 200–208. [Google Scholar] [CrossRef] [PubMed]
  14. Jenkins, D.A.; Hussein, H.; Martina, R.; Dequen-O’Byrne, P.; Abrams, K.R.; Bujkiewicz, S. Methods for the inclusion of real-world evidence in network meta-analysis. BMC Med. Res. Methodol. 2021, 21, 207. [Google Scholar] [CrossRef] [PubMed]
  15. Fang, Y.; He, W.; Wang, H.; Wu, M. Key considerations in the design of real-world studies. Contemp. Clin. Trials. 2020, 96, 106091. [Google Scholar] [CrossRef]
  16. Zhang, K.; Arora, P.; Sati, N.; Béliveau, A.; Troke, N.; Veroniki, A.A.; Rodrigues, M.; Rios, P.; Zarin, W.; Tricco, A.C. Characteristics and methods of incorporating randomized and nonrandomized evidence in network meta-analyses: A scoping review. J. Clin. Epidemiol. 2019, 113, 1–10. [Google Scholar] [CrossRef]
  17. Sarri, G.; Patorno, E.; Yuan, H.; Guo, J.J.; Bennett, D.; Wen, X.; Zullo, A.R.; Largent, J.; Panaccio, M.; Gokhale, M.; et al. Framework for the synthesis of non-randomised studies and randomised controlled trials: A guidance on conducting a systematic review and meta-analysis for healthcare decision making. BMJ Evid. Based Med. 2020, 27, 109–119. [Google Scholar] [CrossRef]
  18. Verde, P.E.; Ohmann, C. Combining randomized and non-randomized evidence in clinical research: A review of methods and applications. Res. Synth Methods 2015, 6, 45–62. [Google Scholar] [CrossRef]
  19. Verde, P.E. A bias-corrected meta-analysis model for combining, studies of different types and quality. Biom. J. 2021, 63, 406–422. [Google Scholar] [CrossRef]
  20. Efthimiou, O.; Mavridis, D.; Debray, T.P.; Samara, M.; Belger, M.; Siontis, G.C.; Leucht, S.; Salanti, G. GetReal Work Package 4. Combining randomized and non-randomized evidence in network meta-analysis. Stat. Med. 2017, 36, 1210–1226. [Google Scholar] [CrossRef]
  21. Alkabbani, W.; Pelletier, R.; Gamble, J.M. Sodium/Glucose Cotransporter 2 Inhibitors and the Risk of Diabetic Ketoacidosis: An Example of Complementary Evidence for Rare Adverse Events. Am. J. Epidemiol. 2021, 190, 1572–1581. [Google Scholar] [CrossRef]
  22. Falkenhorst, G.; Remschmidt, C.; Harder, T.; Hummers-Pradier, E.; Wichmann, O.; Bogdan, C. Effectiveness of the 23-Valent Pneumococcal Polysaccharide Vaccine (PPV23) against Pneumococcal Disease in the Elderly: Systematic Review and Meta-Analysis. PLoS ONE 2017, 12, e0169368. [Google Scholar] [CrossRef] [PubMed]
  23. Higgins, J.P.; Altman, D.G.; Gøtzsche, P.C.; Jüni, P.; Moher, D.; Oxman, A.D.; Savovic, J.; Schulz, K.F.; Weeks, L.; Sterne, J.A. The Cochrane Collaboration’s tool for assessing risk of bias in randomised trials. BMJ 2011, 343, d5928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wells, G.A.; Shea, B.; O’Connell, D.; Peterson, J.; Welch, V.; Losos, M.; Tugwell, P.; Ga, S.W. The Newcastle-Ottawa Scale (NOS) for Assessing the Quality of Nonrandomised Studies in Meta-Analyses. Available online: https://www.ohri.ca/programs/clinical_epidemiology/oxford.asp (accessed on 10 February 2023).
  25. Honkanen, P.O.; Keistinen, T.; Miettinen, L.; Herva, E.; Sankilampi, U.; Läärä, E.; Leinonen, M.; Kivelä, S.L.; Mäkelä, P.H. Incremental effectiveness of pneumococcal vaccine on simultaneously administered influenza vaccine in preventing pneumonia and pneumococcal pneumonia among persons aged 65 years or older. Vaccine 1999, 17, 2493–2500. [Google Scholar] [CrossRef] [PubMed]
  26. Maruyama, T.; Taguchi, O.; Niederman, M.S.; Morser, J.; Kobayashi, H.; Kobayashi, T.; D’Alessandro-Gabazza, C.; Nakayama, S.; Nishikubo, K.; Noguchi, T.; et al. Efficacy of 23-valent pneumococcal vaccine in preventing pneumonia and improving survival in nursing home residents: Double blind, randomised and placebo controlled trial. BMJ 2010, 340, c1004. [Google Scholar] [CrossRef] [Green Version]
  27. Ortqvist, A.; Hedlund, J.; Burman, L.A.; Elbel, E.; Höfer, M.; Leinonen, M.; Lindblad, I.; Sundelöf, B.; Kalin, M. Randomised trial of 23-valent pneumococcal capsular polysaccharide vaccine in prevention of pneumonia in middle-aged and elderly people. Swedish Pneumococcal Vaccination Study Group. Lancet 1998, 351, 399–403. [Google Scholar] [CrossRef]
  28. Hechter, R.C.; Chao, C.; Jacobsen, S.J.; Slezak, J.M.; Quinn, V.P.; Van Den Eeden, S.K.; Tseng, H.F. Clinical effectiveness of pneumococcal polysaccharide vaccine in men: California Men’s Health Study. Vaccine 2012, 30, 5625–5630. [Google Scholar] [CrossRef]
  29. Jackson, L.A.; Neuzil, K.M.; Yu, O.; Benson, P.; Barlow, W.E.; Adams, A.L.; Hanson, C.A.; Mahoney, L.D.; Shay, D.K.; Thompson, W.W. Effectiveness of pneumococcal polysaccharide vaccine in older adults. N. Engl. J. Med. 2003, 348, 1747–1755. [Google Scholar] [CrossRef] [Green Version]
  30. Ochoa-Gondar, O.; Vila-Corcoles, A.; Rodriguez-Blanco, T.; Gomez-Bertomeu, F.; Figuerola-Massana, E.; Raga-Luria, X.; Hospital-Guardiola, I. Effectiveness of the 23-valent pneumococcal polysaccharide vaccine against community-acquired pneumonia in the general population aged ≥ 60 years: 3 years of follow-up in the CAPAMIS study. Clin. Infect Dis. 2014, 58, 909–917. [Google Scholar] [CrossRef]
  31. Tsai, Y.H.; Hsieh, M.J.; Chang, C.J.; Wen, Y.W.; Hu, H.C.; Chao, Y.N.; Huang, Y.C.; Yang, C.T.; Huang, C.C. The 23-valent pneumococcal polysaccharide vaccine is effective in elderly adults over 75 years old--Taiwan’s PPV vaccination program. Vaccine 2015, 33, 2897–2902. [Google Scholar] [CrossRef]
  32. Vila-Córcoles, A.; Ochoa-Gondar, O.; Hospital, I.; Ansa, X.; Vilanova, A.; Rodríguez, T.; Llor, C.; EVAN Study Group. Protective effects of the 23-valent pneumococcal polysaccharide vaccine in the elderly population: The EVAN-65 study. Clin. Infect Dis. 2006, 43, 860–868. [Google Scholar] [CrossRef] [PubMed]
  33. Dominguez, A.; Salleras, L.; Fedson, D.S.; Izquierdo, C.; Ruiz, L.; Ciruela, P.; Fenoll, A.; Casal, J. Effectiveness of pneumococcal vaccination for elderly people in Catalonia, Spain: A case-control study. Clin. Infect Dis. 2005, 40, 1250–1257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Leventer-Roberts, M.; Feldman, B.S.; Brufman, I.; Cohen-Stavi, C.J.; Hoshen, M.; Balicer, R.D. Effectiveness of 23-valent pneumococcal polysaccharide vaccine against invasive disease and hospital-treated pneumonia among people aged ≥ 65 years: A retrospective case-control study. Clin. Infect Dis. 2015, 60, 1472–1480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Vila-Corcoles, A.; Salsench, E.; Rodriguez-Blanco, T.; Ochoa-Gondar, O.; de Diego, C.; Valdivieso, A.; Hospital, I.; Gomez-Bertomeu, F.; Raga, X. Clinical effectiveness of 23-valent pneumococcal polysaccharide vaccine against pneumonia in middle-aged and older adults: A matched case-control study. Vaccine 2009, 27, 1504–1510. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, Y.; Desai, M.; Ryan, P.B.; DeFalco, F.J.; Schuemie, M.J.; Stang, P.E.; Berlin, J.A.; Yuan, Z. Incidence of diabetic ketoacidosis among patients with type 2 diabetes mellitus treated with SGLT2 inhibitors and other antihyperglycemic agents. Diabetes Res. Clin. Pract. 2017, 128, 83–90. [Google Scholar] [CrossRef] [Green Version]
  37. Downs, S.H.; Black, N. The feasibility of creating a checklist for the assessment of the methodological quality both of randomised and non-randomised studies of health care interventions. J. Epidemiol. Community Health 1998, 52, 377–384. [Google Scholar] [CrossRef] [Green Version]
  38. Lavalle-González, F.J.; Januszewicz, A.; Davidson, J.; Tong, C.; Qiu, R.; Canovatchel, W.; Meininger, G. Efficacy and safety of canagliflozin compared with placebo and sitagliptin in patients with type 2 diabetes on background metformin monotherapy: A randomised trial. Diabetologia 2013, 56, 2582–2592. [Google Scholar] [CrossRef] [Green Version]
  39. Roden, M.; Merker, L.; Christiansen, A.V.; Roux, F.; Salsali, A.; Kim, G.; Stella, P.; Woerle, H.J.; Broedl, U.C.; EMPA-REG EXTEND™ MONO investigators. Safety, tolerability and effects on cardiometabolic risk factors of empagliflozin monotherapy in drug-naïve patients with type 2 diabetes: A double-blind extension of a Phase III randomized controlled trial. Cardiovasc. Diabetol. 2015, 14, 154. [Google Scholar] [CrossRef] [Green Version]
  40. Haering, H.U.; Merker, L.; Christiansen, A.V.; Roux, F.; Salsali, A.; Kim, G.; Meinicke, T.; Woerle, H.J.; Broedl, U.C.; EMPA-REG EXTEND™ METSU investigators. Empagliflozin as add-on to metformin plus sulphonylurea in patients with type 2 diabetes. Diabetes Res. Clin. Pract. 2015, 110, 82–90. [Google Scholar] [CrossRef]
  41. Frías, J.P.; Guja, C.; Hardy, E.; Ahmed, A.; Dong, F.; Öhman, P.; Jabbour, S.A. Exenatide once weekly plus dapagliflozin once daily versus exenatide or dapagliflozin alone in patients with type 2 diabetes inadequately controlled with metformin monotherapy (DURATION-8): A 28 week, multicentre, double-blind, phase 3, randomised controlled trial. Lancet Diabetes Endocrinol. 2016, 4, 1004–1016. [Google Scholar]
  42. Hollander, P.; Liu, J.; Hill, J.; Johnson, J.; Jiang, Z.W.; Golm, G.; Huyck, S.; Terra, S.G.; Mancuso, J.P.; Engel, S.S. Ertugliflozin Compared with Glimepiride in Patients with Type 2 Diabetes Mellitus Inadequately Controlled on Metformin: The VERTIS SU Randomized Study. Diabetes Ther. 2018, 9, 193–207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Pratley, R.E.; Eldor, R.; Raji, A.; Golm, G.; Huyck, S.B.; Qiu, Y.; Sunga, S.; Johnson, J.; Terra, S.G.; Mancuso, J.P. Ertugliflozin plus sitagliptin versus either individual agent over 52 weeks in patients with type 2 diabetes mellitus inadequately controlled with metformin: The VERTIS FACTORIAL randomized trial. Diabetes Obes. Metab. 2018, 20, 1111–1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gallo, S.; Charbonnel, B.; Goldman, A.; Shi, H.; Huyck, S.; Darekar, A.; Lauring, B.; Terra, S.G. Long-term efficacy and safety of ertugliflozin in patients with type 2 diabetes mellitus inadequately controlled with metformin monotherapy: 104-week VERTIS MET trial. Diabetes Obes. Metab. 2019, 21, 1027–1036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Fralick, M.; Schneeweiss, S.; Patorno, E. Risk of Diabetic Ketoacidosis after Initiation of an SGLT2 Inhibitor. N. Engl. J. Med. 2017, 376, 2300–2302. [Google Scholar] [CrossRef] [PubMed]
  46. Kim, Y.G.; Jeon, J.Y.; Han, S.J.; Kim, D.J.; Lee, K.W.; Kim, H.J. Sodium-glucose co-transporter-2 inhibitors and the risk of ketoacidosis in patients with type 2 diabetes mellitus: A nationwide population-based cohort study. Diabetes Obes. Metab. 2018, 20, 1852–1858. [Google Scholar] [CrossRef]
  47. Ueda, P.; Svanström, H.; Melbye, M.; Eliasson, B.; Svensson, A.M.; Franzén, S.; Gudbjörnsdottir, S.; Hveem, K.; Jonasson, C.; Pasternak, B. Sodium glucose cotransporter 2 inhibitors and risk of serious adverse events: Nationwide register based cohort study. BMJ 2018, 363, k4365. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, L.; Voss, E.A.; Weaver, J.; Hester, L.; Yuan, Z.; DeFalco, F.; Schuemie, M.J.; Ryan, P.B.; Sun, D.; Freedman, A.; et al. Diabetic ketoacidosis in patients with type 2 diabetes treated with sodium glucose co-transporter 2 inhibitors versus other antihyperglycemic agents: An observational study of four US administrative claims databases. Pharm. Drug Saf. 2019, 28, 1620–1628. [Google Scholar]
  49. Günhan, B.; Röver, C.; Friede, T. Random-effects meta-analysis of few studies involving rare events. Res. Synth. Methods 2020, 11, 74–90. [Google Scholar] [CrossRef] [Green Version]
  50. Friede, T.; Röver, C.; Wandel, S.; Neuenschwander, B. Meta-analysis of few small studies in orphan diseases. Res. Synth. Methods 2017, 8, 79–91. [Google Scholar] [CrossRef]
  51. Greenland, S.; Mansournia, M.A.; Altman, D.G. Sparse data bias: A problem hiding in plain sight. BMJ 2016, 352, i1981. [Google Scholar] [CrossRef] [Green Version]
  52. Gelman, A.; Hill, J. Data Analysis Using Regression and Multilevel/Hierarchical Models; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
  53. Hong, Y.D.; Jansen, J.P.; Guerino, J.; Berger, M.L.; Crown, W.; Goettsch, W.G.; Mullins, C.D.; Willke, R.J.; Orsini, L.S. Comparative effectiveness and safety of pharmaceuticals assessed in observational studies compared with randomized controlled trials. BMC Med. 2021, 19, 307. [Google Scholar] [CrossRef] [PubMed]
  54. Schnell-Inderst, P.; Iglesias, C.P.; Arvandi, M.; Ciani, O.; Matteucci Gothe, R.; Peters, J.; Blom, A.W.; Taylor, R.S.; Siebert, U. A bias-adjusted evidence synthesis of RCT and observational data: The case of total hip replacement. Health Econ. 2017, 26 (Suppl. 1), 46–69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Li, L.; Deng, K.; Busse, J.W.; Zhou, X.; Xu, C.; Liu, Z.; Ren, Y.; Zou, K.; Sun, X. A systematic survey showed important limitations in the methods for assessing drug safety among systematic reviews. J. Clin. Epidemiol. 2020, 123, 80–90. [Google Scholar] [CrossRef]
  56. Anglemyer, A.; Horvath, H.T.; Bero, L. Healthcare outcomes assessed with observational study designs compared with those assessed in randomized trials. Cochrane Database Syst. Rev. 2014, 2014, Mr000034. [Google Scholar] [CrossRef] [PubMed]
  57. Greenland, S. Bayesian perspectives for epidemiological research: I. Foundations and basic methods. Int. J. Epidemiol. 2006, 35, 765–775. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Turner, R.M.; Davey, J.; Clarke, M.J.; Thompson, S.G.; Higgins, J.P. Predicting the extent of heterogeneity in meta-analysis, using empirical data from the Cochrane Database of Systematic Reviews. Int. J. Epidemiol. 2012, 41, 818–827. [Google Scholar] [CrossRef] [Green Version]
Table 1. Characteristics of included studies in the meta-analysis of the efficacy and effectiveness of pneumococcal polysaccharide vaccine (PPV23) vaccination against invasive pneumococcal disease (IPD) in elderly patients.
Table 1. Characteristics of included studies in the meta-analysis of the efficacy and effectiveness of pneumococcal polysaccharide vaccine (PPV23) vaccination against invasive pneumococcal disease (IPD) in elderly patients.
First Author (Year)DesignCases/ParticipantsLog (OR)SERisk of Bias
Honkanen (1999) [25]RCT7/38,037−0.970.84Unclear
Maruyama (2010) [26]RCT3/2289−1.941.51Low
Ortqvist (1998) [27]RCT6/1666−1.521.1Low
Hechter (2012) [28]Cohort study9/31,282−1.040.89High
Jackson (2003) [29]Cohort study61/47,365−0.580.26Low
Ochoa-Gondar (2014) [30]Cohort study16/27,204−0.960.75Low
Tsai (2015) [31]Cohort study57/458,362−1.420.33High
Vila-Corcoles (2006) [32]Cohort study22/11,241−0.510.51Low
Dominguez (2005) [33]Case-control study149/596−1.180.26Low
Leventer-Roberts (2015) [34]Case-control study212/1060−0.540.17Low
Vila-Corcoles (2009) [35] Case-control study94/282−1.080.34Low
Table 2. Characteristics of studies included in the meta-analysis of the risk of diabetic ketoacidosis (DKA) among patients using sodium/glucose cotransporter 2 (SGLT-2) inhibitors compared with active comparators.
Table 2. Characteristics of studies included in the meta-analysis of the risk of diabetic ketoacidosis (DKA) among patients using sodium/glucose cotransporter 2 (SGLT-2) inhibitors compared with active comparators.
First Author (Year)DesignCases/ParticipantsLog (OR)SEStudy Quality
Lavalle-González (2013) [38]RCT1/1284−0.791.23Good
Roden (2015) [39]RCT1/6800.491.24Good
Haering (2015) [40]RCT2/2702−0.231.01Good
Frías (2016) [41]RCT1/463−0.701.23Good
Hollander (2018) [42]RCT1/13610.501.24Fair
Pratley (2018) [43]RCT1/12320.311.27Good
Gallo (2019) [44]RCT1/4140.701.23Good
Fralick (2017) [45]Cohort study81/76,0900.790.24Good
Wang (2017) [36]Cohort study55/60,9320.650.38Fair
Kim (2018) [46]Cohort study63/112,650−0.050.25Good
Ueda (2018) [47]Cohort study30/34,4260.760.38Good
Wang-CCAE (2019) [48]Cohort study668/220,5040.340.10Good
Wang-MDCD (2019) [48]Cohort study155/20,5320.170.20Good
Wang-MDCR (2019) [48]Cohort study80/27,7640.980.34Good
Wang-Optum (2019) [48]Cohort study379/115,7220.250.14Good
Douros (2020) [10]Cohort study505/404,3721.050.18Good
Table 3. Estimated effects and 95% credibility intervals of the effectiveness of pneumococcal polysaccharide vaccine (PPV23) vaccination in preventing invasive pneumococcal disease (IPD) in elderly patients.
Table 3. Estimated effects and 95% credibility intervals of the effectiveness of pneumococcal polysaccharide vaccine (PPV23) vaccination in preventing invasive pneumococcal disease (IPD) in elderly patients.
Variance Inflation Factor (w)DASRPITHM
w~beta (0.25, 1)0.40 (0.22–0.67)0.38 (0.12–0.91)0.37 (0.09–1.32)
w~beta (1.5, 1)0.42 (0.27–0.62)0.40 (0.23–0.64)0.39 (0.12–1.20)
w~beta (4, 1)0.42 (0.28–0.59)0.41 (0.26–0.61)0.39 (0.14–1.14)
Table 4. Estimated effects and 95% credibility intervals of the risk of diabetic ketoacidosis (DKA) among users receiving sodium/glucose cotransporter 2 (SGLT-2) inhibitors versus active comparators.
Table 4. Estimated effects and 95% credibility intervals of the risk of diabetic ketoacidosis (DKA) among users receiving sodium/glucose cotransporter 2 (SGLT-2) inhibitors versus active comparators.
Variance Inflation Factor (w)DASRPITHM
w~beta (0.25, 1)1.46 (1.06–2.07)1.41 (0.56–2.45)1.32 (0.36–3.81)
w~beta (1.5, 1)1.53 (1.19–2.00)1.52 (0.98–2.19)1.40 (0.47–3.73)
w~beta (4, 1)1.56 (1.21–2.03)1.56 (1.12–2.12)1.42 (0.50–3.75)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yao, M.; Wang, Y.; Mei, F.; Zou, K.; Li, L.; Sun, X. Methods for the Inclusion of Real-World Evidence in a Rare Events Meta-Analysis of Randomized Controlled Trials. J. Clin. Med. 2023, 12, 1690. https://doi.org/10.3390/jcm12041690

AMA Style

Yao M, Wang Y, Mei F, Zou K, Li L, Sun X. Methods for the Inclusion of Real-World Evidence in a Rare Events Meta-Analysis of Randomized Controlled Trials. Journal of Clinical Medicine. 2023; 12(4):1690. https://doi.org/10.3390/jcm12041690

Chicago/Turabian Style

Yao, Minghong, Yuning Wang, Fan Mei, Kang Zou, Ling Li, and Xin Sun. 2023. "Methods for the Inclusion of Real-World Evidence in a Rare Events Meta-Analysis of Randomized Controlled Trials" Journal of Clinical Medicine 12, no. 4: 1690. https://doi.org/10.3390/jcm12041690

APA Style

Yao, M., Wang, Y., Mei, F., Zou, K., Li, L., & Sun, X. (2023). Methods for the Inclusion of Real-World Evidence in a Rare Events Meta-Analysis of Randomized Controlled Trials. Journal of Clinical Medicine, 12(4), 1690. https://doi.org/10.3390/jcm12041690

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop