A comparison of computational algorithms for the Bayesian analysis of clinical trials.
Background: Clinical trials are increasingly using Bayesian methods for their design and analysis. Inference in Bayesian trials typically uses simulation-based approaches such as Markov Chain Monte Carlo methods. Markov Chain Monte Carlo has high computational cost and can be complex to implement. T...
| Publicado en: | Clinical Trials Vol. 21; no. 6; pp. 689 - 701 |
|---|---|
| Autores principales: | , , , , , , , , , , , , , |
| Formato: | equations & formulas research tables/charts Journal Article |
| Publicado: |
Sage Publications, Ltd.
Dec2024
|
| Acceso en línea: | Ver este registro en EBSCOhost |
| fields | @attributes: recordID: 1 pdfLink: plink: https://search.ebscohost.com/login.aspx?direct=true&db=ccm&AN=180624273&site=ehost-live header: @attributes: shortDbName: ccm uiTerm: 180624273 longDbName: CINAHL Complete uiTag: AN controlInfo: bkinfo: dissinfo: jinfo: jid: 17407745 AB6 jtl: Clinical Trials issn: 17407745 maglogo: N pubinfo: dt: Dec2024 vid: 21 iid: 6 pid: 33180 pub: Sage Publications, Ltd. place: <Blank> artinfo: ui: 180624273 177265500 180624273 180624273 10.1177/17407745241247334 180624273 ppf: 689 ppct: 12 formats: tig: atl: A comparison of computational algorithms for the Bayesian analysis of clinical trials. aug: au: Chen, Ziming Berger, Jeffrey S Castellucci, Lana A Farkouh, Michael Goligher, Ewan C Hade, Erinn M Hunt, Beverley J Kornblith, Lucy Z Lawler, Patrick R Leifer, Eric S Lorenzi, Elizabeth Neal, Matthew D Zarychanski, Ryan Heath, Anna affil: The Hospital for Sick Children, Toronto, ON, Canada sug: subj: Algorithms Probability Anticoagulants Therapeutic Use Heparin Therapeutic Use COVID-19 Clinical Trials Human Multimethod Studies Survival Analysis Length of Stay Descriptive Statistics Models, Statistical Comparative Studies Systems Analysis Funding Source ab: Background: Clinical trials are increasingly using Bayesian methods for their design and analysis. Inference in Bayesian trials typically uses simulation-based approaches such as Markov Chain Monte Carlo methods. Markov Chain Monte Carlo has high computational cost and can be complex to implement. The Integrated Nested Laplace Approximations algorithm provides approximate Bayesian inference without the need for computationally complex simulations, making it more efficient than Markov Chain Monte Carlo. The practical properties of Integrated Nested Laplace Approximations compared to Markov Chain Monte Carlo have not been considered for clinical trials. Using data from a published clinical trial, we aim to investigate whether Integrated Nested Laplace Approximations is a feasible and accurate alternative to Markov Chain Monte Carlo and provide practical guidance for trialists interested in Bayesian trial design. Methods: Data from an international Bayesian multi-platform adaptive trial that compared therapeutic-dose anticoagulation with heparin to usual care in non-critically ill patients hospitalized for COVID-19 were used to fit Bayesian hierarchical generalized mixed models. Integrated Nested Laplace Approximations was compared to two Markov Chain Monte Carlo algorithms, implemented in the software JAGS and stan, using packages available in the statistical software R. Seven outcomes were analysed: organ-support free days (an ordinal outcome), five binary outcomes related to survival and length of hospital stay, and a time-to-event outcome. The posterior distributions for the treatment and sex effects and the variances for the hierarchical effects of age, site and time period were obtained. We summarized these posteriors by calculating the mean, standard deviations and the 95% equitailed credible intervals and presenting the results graphically. The computation time for each algorithm was recorded. Results: The average overlap of the 95% credible interval for the treatment and sex effects estimated using Integrated Nested Laplace Approximations was 96% and 97.6% compared with stan, respectively. The graphical posterior densities for these effects overlapped for all three algorithms. The posterior mean for the variance of the hierarchical effects of age, site and time estimated using Integrated Nested Laplace Approximations are within the 95% credible interval estimated using Markov Chain Monte Carlo but the average overlap of the credible interval is lower, 77%, 85.6% and 91.3%, respectively, for Integrated Nested Laplace Approximations compared to stan. Integrated Nested Laplace Approximations and stan were easily implemented in clear, well-established packages in R, while JAGS required the direct specification of the model. Integrated Nested Laplace Approximations was between 85 and 269 times faster than stan and 26 and 1852 times faster than JAGS. Conclusion: Integrated Nested Laplace Approximations could reduce the computational complexity of Bayesian analysis in clinical trials as it is easy to implement in R, substantially faster than Markov Chain Monte Carlo methods implemented in JAGS and stan, and provides near identical approximations to the posterior distributions for the treatment effect. Integrated Nested Laplace Approximations was less accurate when estimating the posterior distribution for the variance of hierarchical effects, particularly for the proportional odds model, and future work should determine if the Integrated Nested Laplace Approximations algorithm can be adjusted to improve this estimation. pubtype: Academic Journal doctype: equations & formulas research tables/charts Journal Article ougenre: Article language: English refInfo: holdings: @attributes: islocal: N |
|---|