Abstract:Clinical Trials, Ahead of Print. 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.

Extending Bayesian analysis of circular data to comparison of multiple groups

Bayesian Analysis of Circular Data Using Wrapped Distributions

Bringing the analysis of animal orientation data full circle: model-based approaches with maximum likelihood

Bayesian estimation for longitudinal data in a joint model with HPCs

Bayesian model selection on linear mixed-effects models for comparisons between multiple treatments and a control

Recommendations for Bayesian hierarchical model specifications for case-control studies in mental health

A comparison of Bayesian sampling algorithms for high-dimensional particle physics and cosmology applications

Bayesian analysis of multivariate mixed longitudinal ordinal and continuous data

Hierarchical Bayesian Clustering Design of Multiple Biomarker Subgroups (HCOMBS)

A Bayesian Phase II Proof-of-concept Design for Clinical Trials with Longitudinal Endpoints.

Parallel Markov Chain Monte Carlo for Bayesian Hierarchical Models with Big Data, in Two Stages

Bayesian joint analysis of heterogeneous- and skewed-longitudinal data and a binary outcome, with application to AIDS clinical studies

Bayesian Solutions for Assessing Differential Effects in Biomarker Positive and Negative Subgroups

Finding our Way in the Dark: Approximate MCMC for Approximate Bayesian Methods

A comparison of computational algorithms for the Bayesian analysis of clinical trials

Investigating the causal effects of multiple treatments using longitudinal data: a simulation study

Statistical Inference in Circular Structural Model and Fitting Circles to Noisy Data

Nonparametric multimodal regression for circular data

Parallel Markov Chain Monte Carlo for Non-Gaussian Posterior Distributions

Bayesian Sequentially Monitored Multi-arm Experiments with Multiple Comparison Adjustments

Bayesian Nonparametric Cost-Effectiveness Analyses: Causal Estimation and Adaptive Subgroup Discovery