Mitchell Aaron Schepps, Jérémy Seurat, France Mentré, Weng Kee Wong
{"title":"Design optimization of longitudinal studies using metaheuristics: Application to lithium pharmacokinetics.","authors":"Mitchell Aaron Schepps, Jérémy Seurat, France Mentré, Weng Kee Wong","doi":"10.1177/09622802251350262","DOIUrl":"https://doi.org/10.1177/09622802251350262","url":null,"abstract":"<p><p>Lithium is recommended as a first line treatment for patients with bipolar disorder. However, only certain patients show a good response to the drug, and the variability and tolerability of lithium response are poorly understood. Greater precision in the early identification of individuals who are likely to respond well to lithium is a significant unmet clinical need. We create optimal designs to better understand the pharmacokinetic exposition of lithium for patients with and without a genetic covariate. From a Fisher information matrix based method, we find different optimal designs for estimating various parameters in a complicated pharmacokinetics/pharmacodynamics nonlinear mixed effects model with multiple physician specified constraints. Our approach uses flexible state-of-the-art metaheuristics to find various types of efficient designs, including multiple-objective optimal designs that can balance the competitiveness of the objectives and deliver higher efficiencies for more important objectives. Results from this article will be used as part of a broader study to implement efficient designs to better understand the exposition of sustained-release lithium in patients with bipolar disorder.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"9622802251350262"},"PeriodicalIF":1.6,"publicationDate":"2025-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144326884","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Isabelle Wilson, Steven Julious, Christina Yap, Susan Todd, Munyaradzi Dimairo
{"title":"Response adaptive randomisation in clinical trials: Current practice, gaps and future directions.","authors":"Isabelle Wilson, Steven Julious, Christina Yap, Susan Todd, Munyaradzi Dimairo","doi":"10.1177/09622802251348183","DOIUrl":"https://doi.org/10.1177/09622802251348183","url":null,"abstract":"<p><p><b>Introduction:</b> Adaptive designs (ADs) offer clinical trials flexibility to modify design aspects based on accumulating interim data. Response adaptive randomisation (RAR) adjusts treatment allocation according to interim results, favouring promising treatments. Despite scientific appeal, RAR adoption lags behind other ADs. Understanding methods and applications could provide insights and resources and reveal future research needs. This study examines RAR application, trial results and achieved benefits, reporting gaps, statistical tools and concerns, while highlighting examples of effective practices. <b>Methods:</b> RAR trials with comparative efficacy, effectiveness or safety objectives, classified at least phase I/II, were identified via statistical literature, trial registries, statistical resources and researcher-knowledge. Search spanned until October 2023, including results until February 2024. Analysis was descriptive and narrative. <b>Results:</b> From 652 articles/trials screened, 65 planned RAR trials (11 platform trials) were identified, beginning in 1985 and gradually increasing through to 2023. Most trials were in oncology (25%) and drug-treatments (80%), with 63% led by US teams. Predominantly Phase II (62%) and multi-arm (63%), 85% used Bayesian methods, testing superiority hypotheses (86%). Binary outcomes appeared in 55%, with a median observation of 56 days. Bayesian RAR algorithms were applied in 83%. However, 71% of all trials lacked clear details on statistical implementation. Subgroup-level RAR was seen in 23% of trials. Allocation was restricted in 51%, and 88% was included a burn-in period. Most trials (85%) planned RAR alongside other adaptations. Of trials with results, 92% used RAR, but over 50% inadequately reported allocation changes. A mean 22% reduction in sample size was seen, with none over-allocating to ineffective arms. <b>Conclusion:</b> RAR has shown benefits in conditions like sepsis, COVID-19 and cancer, enhancing effective treatment allocation and saving resources. However, complexity, costs and simulation need limit wider adoption. This review highlights RAR's benefits and suggests enhancing statistical tools to encourage wider adoption in clinical research.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"9622802251348183"},"PeriodicalIF":1.6,"publicationDate":"2025-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144317876","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lara Maleyeff, Fan Li, Sebastien Haneuse, Rui Wang
{"title":"Permutation tests for detecting treatment effect heterogeneity in cluster randomized trials.","authors":"Lara Maleyeff, Fan Li, Sebastien Haneuse, Rui Wang","doi":"10.1177/09622802251348999","DOIUrl":"https://doi.org/10.1177/09622802251348999","url":null,"abstract":"<p><p>Cluster randomized trials are widely used in healthcare research for the evaluation of intervention strategies. Beyond estimating the average treatment effect, it is often of interest to assess whether the treatment effect varies across subgroups. While conventional methods based on tests of interaction terms between treatment and covariates can be used to detect treatment effect heterogeneity in cluster randomized trials, they typically rely on parametric assumptions that may not hold in practice. Adapting existing permutation tests from individually randomized trials, however, requires conceptual clarification and modification due to the multiple possible interpretations of treatment effect heterogeneity in the cluster randomized trial context. In this work, we develop variations of permutation tests and clarify key causal definitions in order to assess treatment effect heterogeneity in cluster randomized trials. Our procedure enables investigators to simultaneously test for effect modification across a large number of covariates, while maintaining nominal type I error rates and reasonable power in simulation studies. In the Pain Program for Active Coping and Training (PPACT) study, the proposed methods are able to detect treatment effect heterogeneity that was not identified by conventional methods assessing treatment-covariate interactions.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"9622802251348999"},"PeriodicalIF":1.6,"publicationDate":"2025-06-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144317875","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Partial areas under the curve of the cumulative distribution function as a new composite estimand for randomized clinical trials.","authors":"Masataka Taguri, Kenichi Hayashi","doi":"10.1177/09622802251314195","DOIUrl":"10.1177/09622802251314195","url":null,"abstract":"<p><p>Clinical trials often face the challenge of post-randomization events, such as the initiation of rescue therapy or the premature discontinuation of randomized treatment. Such events, called \"intercurrent events\" (ICEs) in ICH E9(R1), may influence the estimation and interpretation of treatment effects. According to ICH E9(R1), there are five strategies for handling ICEs. This study focuses on the composite strategy, which incorporates ICEs in the outcome of interest and defines the treatment effects using composite endpoints that combine the measured continuous variables and ICEs. An advantage of this strategy is that it avoids the occurrence of missing data because they are defined as part of the outcome of interest. In this study, we propose a new composite estimand: the difference in the partial areas under the curves (pAUCs) of the cumulative distribution function. While the pAUC is closely related to the trimmed mean approach proposed by Permutt and Li, it offers the advantage of allowing pre-specification of the cutoff value for a \"good\" response based on clinical considerations. This ensures that the pAUC can be calculated irrespective of the proportion of ICEs. We describe the causal interpretation of our method and its relationship with two other strategies (treatment policy and hypothetical strategies) using a potential outcome framework. We present simulation results in which our method performs reasonably well compared to several existing approaches in terms of type I error, power, and the proportion of undefined test statistics.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1097-1113"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144062239","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analysis of recurrent events in cluster randomised trials: The PLEASANT trial case study.","authors":"Kelly Grant, Steven A Julious","doi":"10.1177/09622802251316972","DOIUrl":"10.1177/09622802251316972","url":null,"abstract":"<p><p>Recurrent events for many clinical conditions, such as asthma, can indicate poor health outcomes. Recurrent events data are often analysed using statistical methods such as Cox regression or negative binomial regression, suffering event or time information loss. This article re-analyses the preventing and lessening exacerbations of asthma in school-age children associated with a new term (PLEASANT) trial data as a case study, investigating the utility, extending recurrent events survival analysis methods to cluster randomised trials. A conditional frailty model is used, with the frailty term at the general practitioner practice level, accounting for clustering. A rare events bias adjustment is applied if few participants had recurrent events and truncation of small event risk sets is explored, to improve model accuracy. Global and event-specific estimates are presented, alongside a mean cumulative function plot to aid interpretation. The conditional frailty model global results are similar to PLEASANT results, but with greater precision (include time, recurrent events, within-participant dependence, and rare events adjustment). Event-specific results suggest an increasing risk reduction in medical appointments for the intervention group, in September-December 2013, as medical contacts increase over time. The conditional frailty model is recommended when recurrent events are a study outcome for clinical trials, including cluster randomised trials, to help explain changes in event risk over time, assisting clinical interpretation.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1079-1096"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12209553/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144080727","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Permutation-based global rank test with adaptive weights for multiple primary endpoints.","authors":"Satoshi Yoshida, Yusuke Yamaguchi, Kazushi Maruo, Masahiko Gosho","doi":"10.1177/09622802251334886","DOIUrl":"10.1177/09622802251334886","url":null,"abstract":"<p><p>Multiple efficacy endpoints are investigated in clinical trials, and selecting the appropriate primary endpoints is key to the study's success. The global test is an analysis approach that can handle multiple endpoints without multiplicity adjustment. This test, which aggregates the statistics from multiple primary endpoints into a single statistic using weights for the statistical comparison, has been gaining increasing attention. A key consideration in the global test is determination of the weights. In this study, we propose a novel global rank test in which the weights for each endpoint are estimated based on the current study data to maximize the test statistic, and the permutation test is applied to control the type I error rate. Simulation studies conducted to compare the proposed test with other global tests show that the proposed test can control the type I error rate at the nominal level, regardless of the number of primary endpoints and correlations between endpoints. Additionally, the proposed test offers higher statistical powers when the efficacy is considerably different between endpoints or when endpoints are moderately correlated, such as when the correlation coefficient is greater than or equal to 0.5.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1254-1266"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144080760","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Rank-based estimators of global treatment effects for cluster randomized trials with multiple endpoints on different scales.","authors":"Emma Davies Smith, Vipul Jairath, Guangyong Zou","doi":"10.1177/09622802251338387","DOIUrl":"10.1177/09622802251338387","url":null,"abstract":"<p><p>Cluster randomized trials commonly employ multiple endpoints. When a single summary of treatment effects across endpoints is of primary interest, global methods represent a common analysis strategy. However, specification of the required joint distribution is non-trivial, particularly when the endpoints have different scales. We develop rank-based interval estimators for a global treatment effect referred to here as the \"global win probability, or the mean of multiple Wilcoxon Mann-Whitney probabilities, and interpreted as the probability that a treatment individual responds better than a control individual on average. Using endpoint-specific ranks among the combined sample and within each arm, each individual-level observation is converted to a \"win fraction\" which quantifies the proportion of wins experienced over every observation in the comparison arm. An individual's multiple observations are then replaced with a single \"global win fraction\" by averaging win fractions across endpoints. A linear mixed model is applied directly to the global win fractions to obtain point, variance, and interval estimates adjusted for clustering. Simulation demonstrates our approach performs well concerning confidence interval coverage and type I error, and methods are easily implemented using standard software. A case study using public data is provided with corresponding R and SAS code.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1267-1289"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144080761","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jonathan W Bartlett, Camila Olarte Parra, Emily Granger, Ruth H Keogh, Erik W van Zwet, Rhian M Daniel
{"title":"G-formula with multiple imputation for causal inference with incomplete data.","authors":"Jonathan W Bartlett, Camila Olarte Parra, Emily Granger, Ruth H Keogh, Erik W van Zwet, Rhian M Daniel","doi":"10.1177/09622802251316971","DOIUrl":"10.1177/09622802251316971","url":null,"abstract":"<p><p>G-formula is a popular approach for estimating the effects of time-varying treatments or exposures from longitudinal data. G-formula is typically implemented using Monte-Carlo simulation, with non-parametric bootstrapping used for inference. In longitudinal data settings missing data are a common issue, which are often handled using multiple imputation, but it is unclear how G-formula and multiple imputation should be combined. We show how G-formula can be implemented using Bayesian multiple imputation methods for synthetic data, and that by doing so, we can impute missing data and simulate the counterfactuals of interest within a single coherent approach. We describe how this can be achieved using standard multiple imputation software and explore its performance using a simulation study and an application from cystic fibrosis.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1130-1143"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12209542/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143753917","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Duncan T Wilson, Andrew Hall, Julia M Brown, Rebecca Ea Walwyn
{"title":"Optimising error rates in programmes of pilot and definitive trials using Bayesian statistical decision theory.","authors":"Duncan T Wilson, Andrew Hall, Julia M Brown, Rebecca Ea Walwyn","doi":"10.1177/09622802251322987","DOIUrl":"10.1177/09622802251322987","url":null,"abstract":"<p><p>Pilot trials are often conducted in advance of definitive trials to assess their feasibility and to inform their design. Although pilot trials typically collect primary endpoint data, preliminary tests of effectiveness have been discouraged given their typically low power. Power could be increased at the cost of a higher type I error rate, but there is little methodological guidance on how to determine the optimal balance between these operating characteristics. We consider a Bayesian decision-theoretic approach to this problem, introducing a utility function and defining an optimal pilot and definitive trial programme as that which maximises expected utility. We base utility on changes in average primary outcome, the cost of sampling, treatment costs, and the decision-maker's attitude to risk. We apply this approach to re-design OK-Diabetes, a pilot trial of a complex intervention with a continuous primary outcome with known standard deviation. We then examine how optimal programme characteristics vary with the parameters of the utility function. We find that the conventional approach of not testing for effectiveness in pilot trials can be considerably sub-optimal.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1162-1177"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12209550/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143754028","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Betrand Fesuh Nono, Georges Nguefack-Tsague, Martin Kegnenlezom, Eugène-Patrice N Nguéma
{"title":"An iterative matrix uncertainty selector for high-dimensional generalized linear models with measurement errors.","authors":"Betrand Fesuh Nono, Georges Nguefack-Tsague, Martin Kegnenlezom, Eugène-Patrice N Nguéma","doi":"10.1177/09622802251316963","DOIUrl":"10.1177/09622802251316963","url":null,"abstract":"<p><p>Measurement error is a prevalent issue in high-dimensional generalized linear regression that existing regularization techniques may inadequately address. Most require estimating error distributions, which can be computationally prohibitive or unrealistic. We introduce an error distribution-free approach for variable selection called the Iterative Matrix Uncertainty Selector (IMUS). IMUS employs the matrix uncertainty selector framework for linear models, which is known for its selection consistency properties. It features an efficient iterative algorithm easily implemented for any generalized linear model within the exponential family. Empirically, we demonstrate that IMUS performs well in simulations and on three microarray gene expression datasets, achieving effective covariate selection with smoother convergence and clearer elbow criteria compared to other error distribution free methods. Notably, simulation studies in logistic and Poisson regression showed that IMUS exhibited smoother convergence and clearer elbow criteria, performing comparably to the Generalized Matrix Uncertainty Selector (GMUS) and Generalized Matrix Uncertainty Lasso (GMUL) in covariate selection. In many scenarios, IMUS had smaller estimation errors than GMUL and GMUS, measured by both the 1- and 2-norms. In applications to three microarray datasets with noisy measurements, GMUS faced convergence issues, while GMUL converged but lacked well-defined elbows for two datasets. In contrast, IMUS converged with well-defined elbows for all datasets, providing a potentially effective solution for high dimensional regression problems involving measurement errors.</p>","PeriodicalId":22038,"journal":{"name":"Statistical Methods in Medical Research","volume":" ","pages":"1114-1129"},"PeriodicalIF":1.6,"publicationDate":"2025-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143664493","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}