1
|
Hartley BF, Lunn D, Mander AP. Efficient Study Design and Analysis of Longitudinal Dose-Response Data Using Fractional Polynomials. Pharm Stat 2024. [PMID: 39073285 DOI: 10.1002/pst.2425] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Revised: 06/12/2024] [Accepted: 06/24/2024] [Indexed: 07/30/2024]
Abstract
Correctly characterising the dose-response relationship and taking the correct dose forward for further study is a critical part of the drug development process. We use optimal design theory to compare different designs and show that using longitudinal data from all available timepoints in a continuous-time dose-response model can substantially increase the efficiency of estimation of the dose-response compared to a single timepoint model. We give theoretical results to calculate the efficiency gains for a large class of these models. For example, a linearly growing Emax dose-response in a population with a between/within-patient variance ratio ranging from 0.1 to 1 measured at six visits can be estimated with between 1.43 and 2.22 times relative efficiency gain, or equivalently, with 30% to a 55% reduced sample size, compared to a single model of the final timepoint. Fractional polynomials are a flexible way to incorporate data from repeated measurements, increasing precision without imposing strong constraints. Longitudinal dose-response models using two fractional polynomial terms are robust to mis-specification of the true longitudinal process while maintaining, often large, efficiency gains. These models have applications for characterising the dose-response at interim or final analyses.
Collapse
Affiliation(s)
| | - Dave Lunn
- Department of Biostatistics, GSK Research and Development, Brentford, UK
| | - Adrian P Mander
- Department of Biostatistics, GSK Research and Development, Brentford, UK
| |
Collapse
|
2
|
Yan Z, Yang M. Statistical considerations in model-based dose finding for binary responses under model uncertainty. Stat Med 2024; 43:2472-2485. [PMID: 38605556 DOI: 10.1002/sim.10082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2023] [Revised: 02/21/2024] [Accepted: 04/01/2024] [Indexed: 04/13/2024]
Abstract
The statistical methodology for model-based dose finding under model uncertainty has attracted increasing attention in recent years. While the underlying principles are simple and easy to understand, developing and implementing an efficient approach for binary responses can be a formidable task in practice. Motivated by the statistical challenges encountered in a phase II dose finding study, we explore several key design and analysis issues related to the hybrid testing-modeling approaches for binary responses. The issues include candidate model selection and specifications, optimal design and efficient sample size allocations, and, notably, the methods for dose-response testing and estimation. Specifically, we consider a class of generalized linear models suited for the candidate set and establish D-optimal designs for these models. Additionally, we propose using permutation-based tests for dose-response testing to avoid asymptotic normality assumptions typically required for contrast-based tests. We perform trial simulations to enhance our understanding of these issues.
Collapse
Affiliation(s)
- Zhiwu Yan
- Biostatistics Department, 89bio, Inc., San Francisco, California, USA
| | - Min Yang
- Department of Mathematics, Statistics, and Computer Science, University of Illinois at Chicago, Chicago, Illinois
| |
Collapse
|
3
|
Baksh RA, Strydom A, Carter B, Carriere I, Ritchie K. Toward the right treatment at the right time: Modeling the trajectory of cognitive decline to identify the earliest age of change in people with Alzheimer's disease. ALZHEIMER'S & DEMENTIA (AMSTERDAM, NETHERLANDS) 2024; 16:e12563. [PMID: 38463041 PMCID: PMC10921067 DOI: 10.1002/dad2.12563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Revised: 01/26/2024] [Accepted: 02/03/2024] [Indexed: 03/12/2024]
Abstract
Introduction Age is the greatest risk factor for Alzheimer's disease (AD). A limitation of randomized control trials in AD is a lack of specificity in the age ranges of participants who are enrolled in studies of disease-modifying therapies. We aimed to apply Emax (i.e., maximum effect) modeling as a novel approach to identity ideal treatment windows. Methods Emax curves were fitted to longitudinal cognitive data of 101 participants with AD and 1392 healthy controls. We included the Mini-Mental State Examination (MMSE) and tests of verbal fluency and executive functioning. Results In people with AD, the earliest decline in the MMSE could be detected in the 67-71 age band while verbal fluency declined from the 41-45 age band. In healthy controls, changes in cognition showed a later trajectory of decline. Discussion Emax modeling could be used to design more efficient trials which has implications for randomized control trials targeting the earlier stages of AD.
Collapse
Affiliation(s)
- R. Asaad Baksh
- Department of Forensic and Neurodevelopmental Sciences, Institute of Psychiatry, Psychology, and NeuroscienceKing's College LondonDenmark HillLondonUK
- The LonDownS ConsortiumDenmark HillLondonUK
| | - André Strydom
- Department of Forensic and Neurodevelopmental Sciences, Institute of Psychiatry, Psychology, and NeuroscienceKing's College LondonDenmark HillLondonUK
- The LonDownS ConsortiumDenmark HillLondonUK
- South London and Maudsley NHS Foundation TrustMichael Rutter CentreLondonUK
| | - Ben Carter
- Department of Biostatistics and Health InformaticsInstitute of Psychiatry, Psychology and NeuroscienceKing's College LondonLondonUK
| | - Isabelle Carriere
- INSERM, Institut de Neurosciences de Montpellier INMMontpellierFrance
| | - Karen Ritchie
- INSERM, Institut de Neurosciences de Montpellier INMMontpellierFrance
| |
Collapse
|
4
|
Liu F, Zhao Q, Rodgers AJ, Mehrotra DV. Calculation of Phase 2 dose-finding study sample size for reliable Phase 3 dose selection. Pharm Stat 2023; 22:1076-1088. [PMID: 37550963 DOI: 10.1002/pst.2330] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Revised: 07/09/2023] [Accepted: 07/14/2023] [Indexed: 08/09/2023]
Abstract
Sample sizes of Phase 2 dose-finding studies, usually determined based on a power requirement to detect a significant dose-response relationship, will generally not provide adequate precision for Phase 3 target dose selection. We propose to calculate the sample size of a dose-finding study based on the probability of successfully identifying the target dose within an acceptable range (e.g., 80%-120% of the target) using the multiple comparison and modeling procedure (MCP-Mod). With the proposed approach, different design options for the Phase 2 dose-finding study can also be compared. Due to inherent uncertainty around an assumed true dose-response relationship, sensitivity analyses to assess the robustness of the sample size calculations to deviations from modeling assumptions are recommended. Planning for a hypothetical Phase 2 dose-finding study is used to illustrate the main points. Codes for the proposed approach is available at https://github.com/happysundae/posMCPMod.
Collapse
Affiliation(s)
- Fang Liu
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| | - Qing Zhao
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| | - Anthony J Rodgers
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| | - Devan V Mehrotra
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| |
Collapse
|
5
|
Vishwakarma GK, Bhattacharjee A, Tank F, Pashchenko AF. Subgroup identification of targeted therapy effects on biomarker for time to event data. Cancer Biomark 2023; 38:413-424. [PMID: 37980650 DOI: 10.3233/cbm-230181] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
BACKGROUND The initiation biomarker-driven trials have revolutionized oncology drug development by challenging the traditional phased approach and introducing basket studies. Notable successes in non-small cell lung cancer (NSCLC) with ALK, ALK/ROS1, and EGFR inhibitors have prompted the need to expand this approach to other cancer sites. OBJECTIVES This study explores the use of dose response modeling and time-to-event algorithms on the biomarker molecular targeted agent (MTA). By simulating subgroup identification in MTA-related time-to-event data, the study aims to develop statistical methodology supporting biomarker-driven trials in oncology. METHODS A total of n patients are selected assigned for different doses. A dataset is prepared to mimic the situation on Subgroup Identification of MTA for time to event data analysis. The response is measured through MTA. The MTA value is also measured through ROC. The Markov Chain Monte Carlo (MCMC) techniques are prepared to perform the proposed algorithm. The analysis is carried out with a simulation study. The subset selection is performed through the Threshold Limit Value (TLV) by the Bayesian approach. RESULTS The MTA is observed with range 12-16. It is expected that there is a marginal level shift of the MTA from pre to post-treatment. The Cox time-varying model can be adopted further as causal-effect relation to establishing the MTA on prolonging the survival duration. The proposed work in the statistical methodology to support the biomarker-driven trial for oncology research. CONCLUSION This study extends the application of biomarker-driven trials beyond NSCLC, opening possibilities for implementation in other cancer sites. By demonstrating the feasibility and efficacy of utilizing MTA as a biomarker, the research lays the foundation for refining and validating biomarker use in clinical trials. These advancements aim to enhance the precision and effectiveness of cancer treatments, ultimately benefiting patients.
Collapse
Affiliation(s)
| | | | | | - Alexander F Pashchenko
- Laboratory of Intellectual Control Systems and Simulation, V. A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences, Moscow, Russia
| |
Collapse
|
6
|
A parallel sampling framework for model averaging: Application to dose response studies. Contemp Clin Trials 2022; 123:106957. [PMID: 36228983 DOI: 10.1016/j.cct.2022.106957] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 09/23/2022] [Accepted: 10/01/2022] [Indexed: 01/27/2023]
|
7
|
Kirby S, Chuang-Stein C. The Acute Stroke Therapy by Inhibition of Neutrophils study - Key features and impact. Pharm Stat 2022; 21:720-728. [PMID: 35819119 DOI: 10.1002/pst.2218] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Revised: 01/25/2022] [Accepted: 01/25/2022] [Indexed: 11/05/2022]
Abstract
The Acute Stroke Therapy by Inhibition of Neutrophils (ASTIN) study, initiated in November of the year 2000, is now widely recognized as having been a landmark study in the history of clinical trials. We look at why this is the case by considering its key features and impact. These key features are: the use of Bayesian design and analysis; the use of the normal dynamic linear model; the response adaptive nature of the study; the use of real-time dosing decisions; and the use of an integrated model to predict 90-day response on the Scandinavian Stroke Scale. Our overall conclusion is that the ASTIN study's main impact came from showing the clinical trial community the feasibility of the novel design and analysis used when most of these key features were rarely used in industry trials, let alone used together in one trial in a disease area with a tremendous unmet medical need.
Collapse
|
8
|
Duda JC, Kappenberg F, Rahnenführer J. Model selection characteristics when using MCP-Mod for dose-response gene expression data. Biom J 2022; 64:883-897. [PMID: 35187701 DOI: 10.1002/bimj.202000250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 11/11/2021] [Accepted: 01/12/2022] [Indexed: 11/10/2022]
Abstract
We extend the scope of application for MCP-Mod (Multiple Comparison Procedure and Modeling) to in vitro gene expression data and assess its characteristics regarding model selection for concentration gene expression curves. Precisely, we apply MCP-Mod on single genes of a high-dimensional gene expression data set, where human embryonic stem cells were exposed to eight concentration levels of the compound valproic acid (VPA). As candidate models we consider the sigmoid E max $E_{\max }$ (four-parameter log-logistic), linear, quadratic, E max $E_{\max }$ , exponential, and beta model. Through simulations we investigate the impact of omitting one or more models from the candidate model set to uncover possibly superfluous models and to evaluate the precision and recall rates of selected models. Each model is selected according to Akaike information criterion (AIC) for a considerable number of genes. For less noisy cases the popular sigmoid E max $E_{\max }$ model is frequently selected. For more noisy data, often simpler models like the linear model are selected, but mostly without relevant performance advantage compared to the second best model. Also, the commonly used standard E max $E_{\max }$ model has an unexpected low performance.
Collapse
Affiliation(s)
- Julia C Duda
- Department of Statistics, TU Dortmund University, Dortmund, Germany
| | | | | |
Collapse
|
9
|
Pedder H, Dias S, Boucher M, Bennetts M, Mawdsley D, Welton NJ. Methods to assess evidence consistency in dose-response model based network meta-analysis. Stat Med 2021; 41:625-644. [PMID: 34866221 DOI: 10.1002/sim.9270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2021] [Revised: 08/06/2021] [Accepted: 11/08/2021] [Indexed: 11/10/2022]
Abstract
Network meta-analysis (NMA) simultaneously estimates multiple relative treatment effects based on evidence that forms a network of treatment comparisons. Heterogeneity in treatment definitions, such as dose, can lead to a violation of the consistency assumption that underpins NMA. Model-based NMA (MBNMA) methods have been proposed that allow functional dose-response relationships to be estimated within an NMA, which avoids lumping different doses together and thereby reduces the likelihood of inconsistency. Dose-response MBNMA relies on appropriate specification of the dose-response relationship as well as consistency of relative effects. In this article we describe methods to check for inconsistency in dose-response MBNMA models. Global and local (node-splitting) tests for inconsistency are described that account for studies with ≥3 arms that are typical in dose-finding trials. We show that consistency needs to be assessed with respect to the choice of dose-response function. We illustrate the methods using a network comparing biologics for moderate-to-severe psoriasis. By comparing results from an Emax and an exponential dose-response function we show that failure to correctly characterise the dose-response can introduce apparent inconsistency. The number of comparisons for which node-splitting is possible is also shown to be dependent on the complexity of the selected dose-response function. We highlight that the nature of dose-finding studies, which typically compare multiple doses of the same agent, provide limited scope to assess inconsistency, but these study designs help guard against inconsistency in the first place. We demonstrate the importance of assessing consistency to obtain robust relative effects to inform drug-development and policy decisions.
Collapse
Affiliation(s)
- Hugo Pedder
- Population Health Sciences, Canynge Hall, University of Bristol, Bristol, UK
| | - Sofia Dias
- Centre for Reviews and Dissemination, University of York, York, UK
| | | | - Meg Bennetts
- Pharmacometrics, Pfizer Ltd., Sandwich, Kent, UK
| | - David Mawdsley
- Population Health Sciences, Canynge Hall, University of Bristol, Bristol, UK
| | - Nicky J Welton
- Population Health Sciences, Canynge Hall, University of Bristol, Bristol, UK
| |
Collapse
|
10
|
Nalley EM, Tuttle LJ, Barkman AL, Conklin EE, Wulstein DM, Richmond RH, Donahue MJ. Water quality thresholds for coastal contaminant impacts on corals: A systematic review and meta-analysis. THE SCIENCE OF THE TOTAL ENVIRONMENT 2021; 794:148632. [PMID: 34323749 DOI: 10.1016/j.scitotenv.2021.148632] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2021] [Revised: 06/15/2021] [Accepted: 06/19/2021] [Indexed: 06/13/2023]
Abstract
Reduced water quality degrades coral reefs, resulting in compromised ecosystem function and services to coastal communities. Increasing management capacity on reefs requires prioritization of the development of data-based water-quality thresholds and tipping points. To meet this urgent need of marine resource managers, we conducted a systematic review and meta-analysis that quantified the effects on scleractinian corals of chemical pollutants from land-based and atmospheric sources. We compiled a global dataset addressing the effects of these pollutants on coral growth, mortality, reproduction, physiology, and behavior. The resulting quantitative review of 55 articles includes information about industrial sources, modes of action, experimentally tested concentrations, and previously identified tolerance thresholds of corals to 13 metals, 18 pesticides, 5 polycyclic aromatic hydrocarbons (PAHs), a polychlorinated biphenyl (PCB), and a pharmaceutical. For data-rich contaminants, we make more robust threshold estimates by adapting models for Bayesian hierarchical meta-analysis that were originally developed for biopharmaceutical application. These models use information from multiple studies to characterize the dose-response relationships (i.e., Emax curves) between a pollutant's concentration and various measures of coral health. Metals used in antifouling paints, especially copper, have received a great deal of attention to-date, thus enabling us to estimate the cumulative impact of copper across coral's early life-history. The effects of other land-based pollutants on corals are comparatively understudied, which precludes more quantitative analysis. We discuss opportunities to improve future research so that it can be better integrated into quantitative assessments of the effects of more pollutant types on sublethal coral stress-responses. We also recommend that managers use this information to establish more conservative water quality thresholds that account for the synergistic effects of multiple pollutants on coral reefs. Ultimately, active remediation of local stressors will improve the resistance, resilience, and recovery of individual reefs and reef ecosystems facing the global threat of climate change.
Collapse
Affiliation(s)
- Eileen M Nalley
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, 46-007 Lilipuna Road, Kāne'ohe, HI 96744, USA.
| | - Lillian J Tuttle
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, 46-007 Lilipuna Road, Kāne'ohe, HI 96744, USA; NOAA Pacific Islands Regional Office, Honolulu, HI 96860, USA
| | - Alexandria L Barkman
- Kewalo Marine Laboratory, Pacific Biosciences Research Center, University of Hawai'i at Mānoa, 41 Ahui Street, Honolulu, HI 96813, USA
| | - Emily E Conklin
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, 46-007 Lilipuna Road, Kāne'ohe, HI 96744, USA
| | - Devynn M Wulstein
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, 46-007 Lilipuna Road, Kāne'ohe, HI 96744, USA
| | - Robert H Richmond
- Kewalo Marine Laboratory, Pacific Biosciences Research Center, University of Hawai'i at Mānoa, 41 Ahui Street, Honolulu, HI 96813, USA
| | - Megan J Donahue
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, 46-007 Lilipuna Road, Kāne'ohe, HI 96744, USA
| |
Collapse
|
11
|
VanBuren JM, Casper TC, Nishijima DK, Kuppermann N, Lewis RJ, Dean JM, McGlothlin A. The design of a Bayesian adaptive clinical trial of tranexamic acid in severely injured children. Trials 2021; 22:769. [PMID: 34736498 PMCID: PMC8567588 DOI: 10.1186/s13063-021-05737-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Accepted: 10/20/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Trauma is the leading cause of death and disability in children in the USA. Tranexamic acid (TXA) reduces the blood transfusion requirements in adults and children during surgery. Several studies have evaluated TXA in adults with hemorrhagic trauma, but no randomized controlled trials have occurred in children with trauma. We propose a Bayesian adaptive clinical trial to investigate TXA in children with brain and/or torso hemorrhagic trauma. METHODS/DESIGN We designed a double-blind, Bayesian adaptive clinical trial that will enroll up to 2000 patients. We extend the traditional Emax dose-response model to incorporate a hierarchical structure so multiple doses of TXA can be evaluated in different injury populations (isolated head injury, isolated torso injury, or both head and torso injury). Up to 3 doses of TXA (15 mg/kg, 30 mg/kg, and 45 mg/kg bolus doses) will be compared to placebo. Equal allocation between placebo, 15 mg/kg, and 30 mg/kg will be used for an initial period within each injury group. Depending on the dose-response curve, the 45 mg/kg arm may open in an injury group if there is a trend towards increasing efficacy based on the observed relationship using the data from the lower doses. Response-adaptive randomization allows each injury group to differ in allocation proportions of TXA so an optimal dose can be identified for each injury group. Frequent interim stopping periods are included to evaluate efficacy and futility. The statistical design is evaluated through extensive simulations to determine the operating characteristics in several plausible scenarios. This trial achieves adequate power in each injury group. DISCUSSION This trial design evaluating TXA in pediatric hemorrhagic trauma allows for three separate injury populations to be analyzed and compared within a single study framework. Individual conclusions regarding optimal dosing of TXA can be made within each injury group. Identifying the optimal dose of TXA, if any, for various injury types in childhood may reduce death and disability.
Collapse
Affiliation(s)
- John M. VanBuren
- Department of Pediatrics, University of Utah School of Medicine, 295 Chipeta Way, Salt Lake City, UT 84108 USA
| | - T. Charles Casper
- Department of Pediatrics, University of Utah School of Medicine, 295 Chipeta Way, Salt Lake City, UT 84108 USA
| | - Daniel K. Nishijima
- Department of Emergency Medicine, University of California, Davis School of Medicine, Sacramento, CA 95817 USA
| | - Nathan Kuppermann
- Department of Emergency Medicine, University of California, Davis School of Medicine, Sacramento, CA 95817 USA
- Department of Pediatrics, University of California, Davis School of Medicine, Sacramento, CA 95817 USA
| | - Roger J. Lewis
- Department of Emergency Medicine, Harbor-UCLA Medical Center, Torrance, CA 90509 USA
- Berry Consultants, LLC, Austin, TX 78746 USA
| | - J. Michael Dean
- Department of Pediatrics, University of Utah School of Medicine, 295 Chipeta Way, Salt Lake City, UT 84108 USA
| | | | - For the TIC-TOC Collaborators of the Pediatric Emergency Care Applied Research Network (PECARN)
- Department of Pediatrics, University of Utah School of Medicine, 295 Chipeta Way, Salt Lake City, UT 84108 USA
- Department of Emergency Medicine, University of California, Davis School of Medicine, Sacramento, CA 95817 USA
- Department of Pediatrics, University of California, Davis School of Medicine, Sacramento, CA 95817 USA
- Department of Emergency Medicine, Harbor-UCLA Medical Center, Torrance, CA 90509 USA
- Berry Consultants, LLC, Austin, TX 78746 USA
| |
Collapse
|
12
|
Tang F, Carlson S, Wick J, Gajewski BJ. Bayesian EMAX model with a mixture of normal distributions for dose-response in clinical trials. Contemp Clin Trials 2021; 110:106571. [PMID: 34555517 DOI: 10.1016/j.cct.2021.106571] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2021] [Revised: 08/24/2021] [Accepted: 09/16/2021] [Indexed: 11/30/2022]
Abstract
When a dose-response relationship is monotonic, the EMAX model has been shown to provide a good empirical fit for designing and analyzing dose-response data across a wide range of pharmaceutical studies. However, the EMAX model has never been applied to a finite mixture distribution. Motivated by a proposal investigating DHA dose effect on preterm birth (PTB, <37 weeks gestation) rate, we developed a Bayesian EMAX mixture model incorporating the three normal components finite mixture model into the EMAX framework. The proposed Bayesian EMAX mixture model analyzes gestational age as a continuous variable, which allows for statistically efficient estimates of PTB rate using various cut point with the same parsimonious model. For example, we can estimate the rate of early PTB (ePTB, <34 weeks gestation), PTB (<37 weeks gestation), and late-term birth (>41 weeks gestation) using the same model. We compared our proposed EMAX mixture model with an EMAX logistic model and an independent doses logistic model for a dichotomized endpoint using extensive simulations. Across the scenarios under consideration, the EMAX mixture model achieved higher power than the EMAX logistic model and the independent doses logistic model in detecting the effect of DHA supplementation on the PTB rate. The EMAX mixture model also resulted in smaller mean squared errors (MSE) in PTB rate estimates.
Collapse
Affiliation(s)
- Fengming Tang
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Kansas City, KS 66160, United States of America; Saint Luke's Health System, Kansas City, MO 64111, United States of America.
| | - Susan Carlson
- Department of Dietetics and Nutrition, University of Kansas Medical Center, Kansas City, KS 66160, United States of America
| | - Jo Wick
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Kansas City, KS 66160, United States of America
| | - Byron J Gajewski
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Kansas City, KS 66160, United States of America.
| |
Collapse
|
13
|
Thomas M, Bornkamp B, Ickstadt K. Identifying treatment effect heterogeneity in dose-finding trials using Bayesian hierarchical models. Pharm Stat 2021; 21:17-37. [PMID: 34258861 DOI: 10.1002/pst.2150] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Revised: 04/09/2021] [Accepted: 06/14/2021] [Indexed: 11/12/2022]
Abstract
An important task in drug development is to identify patients, which respond better or worse to an experimental treatment. Identifying predictive covariates, which influence the treatment effect and can be used to define subgroups of patients, is a key aspect of this task. Analyses of treatment effect heterogeneity are however known to be challenging, since the number of possible covariates or subgroups is often large, while samples sizes in earlier phases of drug development are often small. In addition, distinguishing predictive covariates from prognostic covariates, which influence the response independent of the given treatment, can often be difficult. While many approaches for these types of problems have been proposed, most of them focus on the two-arm clinical trial setting, where patients are given either the treatment or a control. In this article we consider parallel groups dose-finding trials, in which patients are administered different doses of the same treatment. To investigate treatment effect heterogeneity in this setting we propose a Bayesian hierarchical dose-response model with covariate effects on dose-response parameters. We make use of shrinkage priors to prevent overfitting, which can easily occur, when the number of considered covariates is large and sample sizes are small. We compare several such priors in simulations and also investigate dependent modeling of prognostic and predictive effects to better distinguish these two types of effects. We illustrate the use of our proposed approach using a Phase II dose-finding trial and show how it can be used to identify predictive covariates and subgroups of patients with increased treatment effects.
Collapse
Affiliation(s)
- Marius Thomas
- Clinical Development and Analytics, Novartis Pharma AG, Basel, Switzerland
| | - Björn Bornkamp
- Clinical Development and Analytics, Novartis Pharma AG, Basel, Switzerland
| | - Katja Ickstadt
- Faculty of Statistics, TU Dortmund University, Dortmund, Germany
| |
Collapse
|
14
|
Aouni J, Bacro JN, Toulemonde G, Colin P, Darchy L. Utility-Based Dose Selection for Phase II Dose-Finding Studies. Ther Innov Regul Sci 2021; 55:818-840. [PMID: 33851358 DOI: 10.1007/s43441-021-00273-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Accepted: 02/26/2021] [Indexed: 10/21/2022]
Abstract
BACKGROUND AND OBJECTIVES Dose selection is a key feature of clinical development. Poor dose selection has been recognized as a major driver of development failure in late phase. It usually involves both efficacy and safety criteria. The objective of this paper is to develop and implement a novel fully Bayesian statistical framework to optimize the dose selection process by maximizing the expected utility in phase III. METHODS The success probability is characterized by means of a utility function with two components, one for efficacy and one for safety. Each component refers to a dose-response model. Moreover, a sequential design (with futility and efficacy rules at the interim analysis) is compared to a fixed design in order to allow one to hasten the decision to perform the late phase study. Operating characteristics of this approach are extensively assessed by simulations under a wide range of dose-response scenarios. RESULTS AND CONCLUSIONS Simulation results illustrate the difficulty of simultaneously estimating two complex dose-response models with enough accuracy to properly rank doses using an utility function combining the two. The probability of making the good decision increases with the sample size. For some scenarios, the sequential design has good properties: with a quite large probability of study termination at interim analysis, it enables to reduce the sample size while maintaining the properties of the fixed design.
Collapse
Affiliation(s)
- Jihane Aouni
- Sanofi, Research and Development, 91385, Chilly-Mazarin, France. .,IMAG, Univ Montpellier, CNRS, Montpellier, France.
| | | | - Gwladys Toulemonde
- IMAG, Univ Montpellier, CNRS, Montpellier, France.,Lemon, INRIA, Montpellier Cedex 5, France
| | - Pierre Colin
- Sanofi, Research and Development, 91385, Chilly-Mazarin, France
| | - Loic Darchy
- Sanofi, Research and Development, 91385, Chilly-Mazarin, France
| |
Collapse
|
15
|
Ursino M, Röver C, Zohar S, Friede T. Random-effects meta-analysis of Phase I dose-finding studies using stochastic process priors. Ann Appl Stat 2021. [DOI: 10.1214/20-aoas1390] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Affiliation(s)
- Moreno Ursino
- Centre de Recherche des Cordeliers, INSERM, Sorbonne Université, USPC, Université de Paris
| | - Christian Röver
- Department of Medical Statistics, University Medical Center Göttingen
| | - Sarah Zohar
- Centre de Recherche des Cordeliers, INSERM, Sorbonne Université, USPC, Université de Paris
| | - Tim Friede
- Department of Medical Statistics, University Medical Center Göttingen
| |
Collapse
|
16
|
Pedder H, Dias S, Bennetts M, Boucher M, Welton NJ. Joining the Dots: Linking Disconnected Networks of Evidence Using Dose-Response Model-Based Network Meta-Analysis. Med Decis Making 2021; 41:194-208. [PMID: 33448252 PMCID: PMC7879230 DOI: 10.1177/0272989x20983315] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Accepted: 11/30/2020] [Indexed: 12/20/2022]
Abstract
BACKGROUND Network meta-analysis (NMA) synthesizes direct and indirect evidence on multiple treatments to estimate their relative effectiveness. However, comparisons between disconnected treatments are not possible without making strong assumptions. When studies including multiple doses of the same drug are available, model-based NMA (MBNMA) presents a novel solution to this problem by modeling a parametric dose-response relationship within an NMA framework. In this article, we illustrate several scenarios in which dose-response MBNMA can connect and strengthen evidence networks. METHODS We created illustrative data sets by removing studies or treatments from an NMA of triptans for migraine relief. We fitted MBNMA models with different dose-response relationships. For connected networks, we compared MBNMA estimates with NMA estimates. For disconnected networks, we compared MBNMA estimates with NMA estimates from an "augmented" network connected by adding studies or treatments back into the data set. RESULTS In connected networks, relative effect estimates from MBNMA were more precise than those from NMA models (ratio of posterior SDs NMA v. MBNMA: median = 1.13; range = 1.04-1.68). In disconnected networks, MBNMA provided estimates for all treatments where NMA could not and were consistent with NMA estimates from augmented networks for 15 of 18 data sets. In the remaining 3 of 18 data sets, a more complex dose-response relationship was required than could be fitted with the available evidence. CONCLUSIONS Where information on multiple doses is available, MBNMA can connect disconnected networks and increase precision while making less strong assumptions than alternative approaches. MBNMA relies on correct specification of the dose-response relationship, which requires sufficient data at different doses to allow reliable estimation. We recommend that systematic reviews for NMA search for and include evidence (including phase II trials) on multiple doses of agents where available.
Collapse
Affiliation(s)
- Hugo Pedder
- Department of Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
| | - Sofia Dias
- Centre for Reviews and Dissemination, University of York, York, North Yorkshire, UK
| | | | | | - Nicky J. Welton
- Department of Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
| |
Collapse
|
17
|
Günhan BK, Meyvisch P, Friede T. Shrinkage Estimation for Dose–Response Modeling in Phase II Trials With Multiple Schedules. Stat Biopharm Res 2020. [DOI: 10.1080/19466315.2020.1850519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Burak Kürsad Günhan
- Department of Medical Statistics, University Medical Center Göttingen, Göttingen, Germany
| | | | - Tim Friede
- Department of Medical Statistics, University Medical Center Göttingen, Göttingen, Germany
| |
Collapse
|
18
|
Bergougnan L, Andersen G, Plum-Mörschel L, Evaristi MF, Poirier B, Tardat A, Ermer M, Herbrand T, Arrubla J, Coester HV, Sansone R, Heiss C, Vitse O, Hurbin F, Boiron R, Benain X, Radzik D, Janiak P, Muslin AJ, Hovsepian L, Kirkesseli S, Deutsch P, Parkar AA. Endothelial-protective effects of a G-protein-biased sphingosine-1 phosphate receptor-1 agonist, SAR247799, in type-2 diabetes rats and a randomized placebo-controlled patient trial. Br J Clin Pharmacol 2020; 87:2303-2320. [PMID: 33125753 PMCID: PMC8247405 DOI: 10.1111/bcp.14632] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Revised: 10/20/2020] [Accepted: 10/24/2020] [Indexed: 12/12/2022] Open
Abstract
Aims SAR247799 is a G‐protein‐biased sphingosine‐1 phosphate receptor‐1 (S1P1) agonist designed to activate endothelial S1P1 and provide endothelial‐protective properties, while limiting S1P1 desensitization and consequent lymphocyte‐count reduction associated with higher doses. The aim was to show whether S1P1 activation can promote endothelial effects in patients and, if so, select SAR247799 doses for further clinical investigation. Methods Type‐2 diabetes patients, enriched for endothelial dysfunction (flow‐mediated dilation, FMD <7%; n = 54), were randomized, in 2 sequential cohorts, to 28‐day once‐daily treatment with SAR247799 (1 or 5 mg in ascending cohorts), placebo or 50 mg sildenafil (positive control) in a 5:2:2 ratio per cohort. Endothelial function was assessed by brachial artery FMD. Renal function, biomarkers and lymphocytes were measured following 5‐week SAR247799 treatment (3 doses) to Zucker diabetic fatty rats and the data used to select the doses for human testing. Results The maximum FMD change from baseline vs placebo for all treatments was reached on day 35; mean differences vs placebo were 0.60% (95% confidence interval [CI] −0.34 to 1.53%; P = .203) for 1 mg SAR247799, 1.07% (95% CI 0.13 to 2.01%; P = .026) for 5 mg SAR247799 and 0.88% (95% CI −0.15 to 1.91%; P = .093) for 50 mg sildenafil. Both doses of SAR247799 were well tolerated, did not affect blood pressure, and were associated with minimal‐to‐no lymphocyte reduction and small‐to‐moderate heart rate decrease. Conclusion These data provide the first human evidence suggesting endothelial‐protective properties of S1P1 activation, with SAR247799 being as effective as the clinical benchmark, sildenafil. Further clinical testing of SAR247799, at sub‐lymphocyte‐reducing doses (≤5 mg), is warranted in vascular diseases associated with endothelial dysfunction.
Collapse
Affiliation(s)
- Luc Bergougnan
- Sanofi R&D, 1 Avenue Pierre Brossolette, Chilly Mazarin, France
| | | | | | | | - Bruno Poirier
- Sanofi R&D, 1 Avenue Pierre Brossolette, Chilly Mazarin, France
| | - Agnes Tardat
- Sanofi R&D, 371 Rue du Professeur Blayac, Montpellier, France
| | | | | | | | | | - Roberto Sansone
- Division of Cardiology, Pulmonary diseases and Vascular medicine, University Hospital Düsseldorf, Düsseldorf, Germany
| | - Christian Heiss
- Department of Clinical and Experimental Medicine, University of Surrey, Stag Hill, Guildford, UK
| | - Olivier Vitse
- Sanofi R&D, 371 Rue du Professeur Blayac, Montpellier, France
| | - Fabrice Hurbin
- Sanofi R&D, 371 Rue du Professeur Blayac, Montpellier, France
| | - Rania Boiron
- Sanofi R&D, 1 Avenue Pierre Brossolette, Chilly Mazarin, France
| | - Xavier Benain
- Sanofi R&D, 371 Rue du Professeur Blayac, Montpellier, France
| | - David Radzik
- Sanofi R&D, 1 Avenue Pierre Brossolette, Chilly Mazarin, France
| | - Philip Janiak
- Sanofi R&D, 1 Avenue Pierre Brossolette, Chilly Mazarin, France
| | | | | | | | | | | |
Collapse
|
19
|
Wadsworth I, Hampson LV, Bornkamp B, Jaki T. Exposure-response modelling approaches for determining optimal dosing rules in children. Stat Methods Med Res 2020; 29:2583-2602. [PMID: 32050840 PMCID: PMC7528535 DOI: 10.1177/0962280220903751] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Within paediatric populations, there may be distinct age groups characterised by different exposure-response relationships. Several regulatory guidance documents have suggested general age groupings. However, it is not clear whether these categorisations will be suitable for all new medicines and in all disease areas. We consider two model-based approaches to quantify how exposure-response model parameters vary over a continuum of ages: Bayesian penalised B-splines and model-based recursive partitioning. We propose an approach for deriving an optimal dosing rule given an estimate of how exposure-response model parameters vary with age. Methods are initially developed for a linear exposure-response model. We perform a simulation study to systematically evaluate how well the various approaches estimate linear exposure-response model parameters and the accuracy of recommended dosing rules. Simulation scenarios are motivated by an application to epilepsy drug development. Results suggest that both bootstrapped model-based recursive partitioning and Bayesian penalised B-splines can estimate underlying changes in linear exposure-response model parameters as well as (and in many scenarios, better than) a comparator linear model adjusting for a categorical age covariate with levels following International Conference on Harmonisation E11 groupings. Furthermore, the Bayesian penalised B-splines approach consistently estimates the intercept and slope more accurately than the bootstrapped model-based recursive partitioning. Finally, approaches are extended to estimate Emax exposure-response models and are illustrated with an example motivated by an in vitro study of cyclosporine.
Collapse
Affiliation(s)
- Ian Wadsworth
- Department of Mathematics & Statistics, Fylde College, Lancaster University, Lancaster, UK
- Phastar, Macclesfield, UK
| | - Lisa V Hampson
- Advanced Methodology & Data Science, Novartis Pharma AG, Basel, Switzerland
| | - Björn Bornkamp
- Advanced Methodology & Data Science, Novartis Pharma AG, Basel, Switzerland
| | - Thomas Jaki
- Department of Mathematics & Statistics, Fylde College, Lancaster University, Lancaster, UK
| |
Collapse
|
20
|
Huang X, Gajewski BJ. Comparison of hierarchical EMAX and NDLM models in dose-response for early phase clinical trials. BMC Med Res Methodol 2020; 20:194. [PMID: 32690004 PMCID: PMC7370408 DOI: 10.1186/s12874-020-01071-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Accepted: 07/01/2020] [Indexed: 11/29/2022] Open
Abstract
Background Phase II clinical trials primarily aim to find the optimal dose and investigate the relationship between dose and efficacy relative to standard of care (control). Therefore, before moving forward to a phase III confirmatory trial, the most effective dose is needed to be identified. Methods The primary endpoint of a phase II trial is typically a binary endpoint of success or failure. The EMAX model, ubiquitous in pharmacology research, was fit for many compounds and described the data well, except for a single compound, which had nonmonotone dose–response (Thomas et al., Stat Biopharmaceutical Res. 6:302-317 2014). To mitigate the risk of nonmonotone dose response one of the alternative options is a Bayesian hierarchical EMAX model (Gajewski et al., Stat Med. 38:3123-3138 2019). The hierarchical EMAX adapts to its environment. Results When the dose-response curve is monotonic it enjoys the efficiency of EMAX. When the dose-response curve is non-monotonic the additional random effect hyperprior makes the hierarchical EMAX model more adjustable and flexible. However, the normal dynamic linear model (NDLM) is a useful model to explore dose-response relationships in that the efficacy at the current dose depends on the efficacy of the previous dose(s). Previous research has compared the EMAX to the hierarchical EMAX (Gajewski et al., Stat Med. 38:3123-3138 2019) and the EMAX to the NDLM (Liu et al., BMC Med Res Method 17:149 2017), however, the hierarchical EMAX has not been directly compared to the NDLM. Conclusions The focus of this paper is to compare these models and discuss the relative merit for each of their uses for an ongoing early phase dose selection study.
Collapse
Affiliation(s)
- Xiaqing Huang
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Mail Stop 1026, 3901 Rainbow Blvd., Kansas City, KS, 66160, USA
| | - Byron J Gajewski
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Mail Stop 1026, 3901 Rainbow Blvd., Kansas City, KS, 66160, USA.
| |
Collapse
|
21
|
Möllenhoff K, Bretz F, Dette H. Equivalence of regression curves sharing common parameters. Biometrics 2019; 76:518-529. [DOI: 10.1111/biom.13149] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Accepted: 09/04/2019] [Indexed: 11/28/2022]
Affiliation(s)
| | - Frank Bretz
- Novartis Pharma AGBasel Switzerland
- Center for Medical Statistics, Informatics and Intelligent SystemsMedical University of Vienna Vienna Austria
| | - Holger Dette
- Department of MathematicsRuhr‐Universität Bochum Bochum Germany
| |
Collapse
|
22
|
Gajewski BJ, Meinzer C, Berry SM, Rockswold GL, Barsan WG, Korley FK, Martin RH. Bayesian hierarchical EMAX model for dose-response in early phase efficacy clinical trials. Stat Med 2019; 38:3123-3138. [PMID: 31070807 DOI: 10.1002/sim.8167] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2018] [Revised: 03/14/2019] [Accepted: 03/14/2019] [Indexed: 11/07/2022]
Abstract
A primary goal of a phase II dose-ranging trial is to identify a correct dose before moving forward to a phase III confirmatory trial. A correct dose is one that is actually better than control. A popular model in phase II is an independent model that puts no structure on the dose-response relationship. Unfortunately, the independent model does not efficiently use information from related doses. One very successful alternate model improves power using a pre-specified dose-response structure. Past research indicates that EMAX models are broadly successful and therefore attractive for designing dose-response trials. However, there may be instances of slight risk of nonmonotone trends that need to be addressed when planning a clinical trial design. We propose to add hierarchical parameters to the EMAX model. The added layer allows information about the treatment effect in one dose to be "borrowed" when estimating the treatment effect in another dose. This is referred to as the hierarchical EMAX model. Our paper compares three different models (independent, EMAX, and hierarchical EMAX) and two different design strategies. The first design considered is Bayesian with a fixed trial design, and it has a fixed schedule for randomization. The second design is Bayesian but adaptive, and it uses response adaptive randomization. In this article, a randomized trial of patients with severe traumatic brain injury is provided as a motivating example.
Collapse
Affiliation(s)
- Byron J Gajewski
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Kansas City, Kansas
| | - Caitlyn Meinzer
- Department of Public Health Sciences, Medical University of South Carolina, Charleston, South Carolina
| | - Scott M Berry
- Department of Biostatistics & Data Science, University of Kansas Medical Center, Kansas City, Kansas
- Berry Consultants, LLC, Austin, Texas
| | | | - William G Barsan
- Department of Emergency Medicine, University of Michigan, Ann Arbor, Michigan
| | - Frederick K Korley
- Department of Emergency Medicine, University of Michigan, Ann Arbor, Michigan
| | - Renee' H Martin
- Department of Public Health Sciences, Medical University of South Carolina, Charleston, South Carolina
| |
Collapse
|
23
|
Buttgereit F, Strand V, Lee EB, Simon-Campos A, McCabe D, Genet A, Tammara B, Rojo R, Hey-Hadavi J. Fosdagrocorat (PF-04171327) versus prednisone or placebo in rheumatoid arthritis: a randomised, double-blind, multicentre, phase IIb study. RMD Open 2019; 5:e000889. [PMID: 31168411 PMCID: PMC6525626 DOI: 10.1136/rmdopen-2018-000889] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 02/27/2019] [Accepted: 03/21/2019] [Indexed: 02/06/2023] Open
Abstract
Objectives Glucocorticoids have anti-inflammatory, transrepression-mediated effects, although adverse events (AEs; transactivation-mediated effects) limit long-term use in patients with rheumatoid arthritis (RA). We evaluated the efficacy and safety of fosdagrocorat (PF-04171327), a dissociated agonist of the glucocorticoid receptor, versus prednisone or placebo. Methods In this 12-week, phase II, randomised controlled trial, 323 patients with moderate to severe RA were randomised 1:1:1:1:1:1:1 to fosdagrocorat (1 mg, 5 mg, 10 mg or 15 mg), prednisone (5 mg or 10 mg) or placebo, once daily. The primary endpoints (week 8) were American College of Rheumatology 20% improvement criteria (ACR20) responses, and percentage changes from baseline in biomarkers of bone formation (procollagen type 1 N-terminal peptide [P1NP]) and resorption (urinary N-telopeptide to urinary creatinine ratio [uNTx:uCr]). Safety was assessed. Results ACR20 responses with fosdagrocorat 10 mg and 15 mg were superior to placebo, and fosdagrocorat 15 mg was non-inferior to prednisone 10 mg (week 8 model-predicted ACR20 responses: 47%, 61%, 69% and 73% vs 51%, 71% and 37% with fosdagrocorat 1 mg, 5 mg, 10 mg and 15 mg vs prednisone 5 mg, 10 mg and placebo, respectively). Percentage changes from baseline in P1NP with fosdagrocorat 1 mg, 5 mg and 10 mg met non-inferiority criteria to prednisone 5 mg. Corresponding changes in uNTx:uCr varied considerably. All fosdagrocorat doses reduced glycosylated haemoglobin levels. AEs were similar between groups; 63 (19.5%) patients reported treatment-related AEs; 9 (2.8%) patients reported serious AEs. No patients had adrenal insufficiency, treatment-related significant infections or laboratory abnormalities. No deaths were reported. Conclusion In patients with RA, fosdagrocorat 10 mg and 15 mg demonstrated efficacy similar to prednisone 10 mg and safety similar to prednisone 5 mg. Trial registration number NCT01393639
Collapse
Affiliation(s)
- Frank Buttgereit
- Rheumatology and Clinical Immunology, Charité University Medicine Berlin (CCM), Berlin, Germany
| | - Vibeke Strand
- Stanford University School of Medicine, Palo Alto, California, USA
| | - Eun Bong Lee
- Internal Medicine, Seoul National University College of Medicine, Seoul, Republic of Korea
| | | | | | | | | | | | | |
Collapse
|
24
|
Smith CL, Jin Y, Raddad E, McNearney TA, Ni X, Monteith D, Brown R, Deeg MA, Schnitzer T. Applications of Bayesian statistical methodology to clinical trial design: A case study of a phase 2 trial with an interim futility assessment in patients with knee osteoarthritis. Pharm Stat 2018; 18:39-53. [DOI: 10.1002/pst.1906] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2017] [Revised: 06/11/2018] [Accepted: 08/22/2018] [Indexed: 12/16/2022]
Affiliation(s)
| | - Yan Jin
- Eli Lilly and Company; Indianapolis IN USA
| | | | | | - Xiao Ni
- Eli Lilly and Company; Indianapolis IN USA
- Novartis Institutes for Biomedical Research; Cambridge MA USA
| | - David Monteith
- Eli Lilly and Company; Indianapolis IN USA
- Xenon Pharmaceuticals Inc.; Burnaby BC Canada
| | | | - Mark A. Deeg
- Eli Lilly and Company; Indianapolis IN USA
- Regulus Therapeutics Inc; San Diego CA USA
| | - Thomas Schnitzer
- Feinberg School of Medicine; Northwestern University; Chicago IL USA
| |
Collapse
|
25
|
Novick S, Ho S, Best N. Data-Driven Prior Distributions for A Bayesian Phase-2 COPD Dose-Finding Clinical Trial. Stat Biopharm Res 2018. [DOI: 10.1080/19466315.2018.1462728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Steven Novick
- Department of Advanced Biostatistics and Data Analytics, GlaxoSmithKline, Uxbridge, Middlesex, UK
- Department of Statistical Sciences, MedImmune, Gaithersburg, MD
| | - Shuyen Ho
- Department of Advanced Biostatistics and Data Analytics, GlaxoSmithKline, Uxbridge, Middlesex, UK
- Department of Statistical Sciences & Innovation, UCB BioSciences, Inc, Raleigh, NC
| | - Nicky Best
- Department of Advanced Biostatistics and Data Analytics, GlaxoSmithKline, Uxbridge, Middlesex, UK
| |
Collapse
|
26
|
Buatois S, Ueckert S, Frey N, Retout S, Mentré F. Comparison of Model Averaging and Model Selection in Dose Finding Trials Analyzed by Nonlinear Mixed Effect Models. AAPS JOURNAL 2018; 20:56. [DOI: 10.1208/s12248-018-0205-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2017] [Accepted: 02/16/2018] [Indexed: 11/30/2022]
|
27
|
Thomas M, Bornkamp B, Seibold H. Subgroup identification in dose-finding trials via model-based recursive partitioning. Stat Med 2018; 37:1608-1624. [DOI: 10.1002/sim.7594] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Revised: 10/24/2017] [Accepted: 11/29/2017] [Indexed: 12/15/2022]
Affiliation(s)
| | | | - Heidi Seibold
- Universität Zürich; Hirschengraben 84 Zürich CH-8001 Switzerland
| |
Collapse
|
28
|
Arnold LM, McCarberg BH, Clair AG, Whalen E, Thomas N, Jorga A, Pauer L, Vissing R, Park PW. Dose–response of pregabalin for diabetic peripheral neuropathy, postherpetic neuralgia, and fibromyalgia. Postgrad Med 2017; 129:921-933. [DOI: 10.1080/00325481.2017.1384691] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Lesley M. Arnold
- Department of Psychiatry and Behavioral Neuroscience, Women's Health Research Program, University of Cincinnati College of Medicine, Cincinnati, OH, USA
| | - Bill H. McCarberg
- Department of Family Medicine, University of California at San Diego School of Medicine, San Diego, CA, USA
| | | | - Ed Whalen
- Statistics, Pfizer, New York, NY, USA
| | | | | | - Lynne Pauer
- Global Product Development - Clinical Sciences & Operations, Pfizer, Groton, CT, USA
| | - Richard Vissing
- Neuroscience and Pain Division, Pfizer Inc, Louisville, KY, USA
| | - Peter W. Park
- North America Medical Affairs, Pfizer Inc, New York, NY, USA
| |
Collapse
|
29
|
Liu F, Walters SJ, Julious SA. Design considerations and analysis planning of a phase 2a proof of concept study in rheumatoid arthritis in the presence of possible non-monotonicity. BMC Med Res Methodol 2017; 17:149. [PMID: 28969588 PMCID: PMC5625783 DOI: 10.1186/s12874-017-0416-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Accepted: 08/31/2017] [Indexed: 05/10/2023] Open
Abstract
Background It is important to quantify the dose response for a drug in phase 2a clinical trials so the optimal doses can then be selected for subsequent late phase trials. In a phase 2a clinical trial of new lead drug being developed for the treatment of rheumatoid arthritis (RA), a U-shaped dose response curve was observed. In the light of this result further research was undertaken to design an efficient phase 2a proof of concept (PoC) trial for a follow-on compound using the lessons learnt from the lead compound. Methods The planned analysis for the Phase 2a trial for GSK123456 was a Bayesian Emax model which assumes the dose-response relationship follows a monotonic sigmoid “S” shaped curve. This model was found to be suboptimal to model the U-shaped dose response observed in the data from this trial and alternatives approaches were needed to be considered for the next compound for which a Normal dynamic linear model (NDLM) is proposed. This paper compares the statistical properties of the Bayesian Emax model and NDLM model and both models are evaluated using simulation in the context of adaptive Phase 2a PoC design under a variety of assumed dose response curves: linear, Emax model, U-shaped model, and flat response. Results It is shown that the NDLM method is flexible and can handle a wide variety of dose-responses, including monotonic and non-monotonic relationships. In comparison to the NDLM model the Emax model excelled with higher probability of selecting ED90 and smaller average sample size, when the true dose response followed Emax like curve. In addition, the type I error, probability of incorrectly concluding a drug may work when it does not, is inflated with the Bayesian NDLM model in all scenarios which would represent a development risk to pharmaceutical company. The bias, which is the difference between the estimated effect from the Emax and NDLM models and the simulated value, is comparable if the true dose response follows a placebo like curve, an Emax like curve, or log linear shape curve under fixed dose allocation, no adaptive allocation, half adaptive and adaptive scenarios. The bias though is significantly increased for the Emax model if the true dose response follows a U-shaped curve. Conclusions In most cases the Bayesian Emax model works effectively and efficiently, with low bias and good probability of success in case of monotonic dose response. However, if there is a belief that the dose response could be non-monotonic then the NDLM is the superior model to assess the dose response. Electronic supplementary material The online version of this article (10.1186/s12874-017-0416-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Feng Liu
- GlaxoSmithKline, Inc, 1250 South Collegeville Road, PO Box 5089, Collegeville, PA, 19426-0989, USA. .,Medical Statistics Group, University of Sheffield, Sheffield, UK.
| | | | - Steven A Julious
- Medical Statistics Group, University of Sheffield, Sheffield, UK
| |
Collapse
|
30
|
Feller C, Schorning K, Dette H, Bermann G, Bornkamp B. Optimal designs for dose response curves with common parameters. Ann Stat 2017. [DOI: 10.1214/16-aos1520] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
31
|
Understanding MCP-MOD dose finding as a method based on linear regression. Stat Med 2017; 36:4401-4413. [DOI: 10.1002/sim.7424] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2016] [Revised: 06/20/2017] [Accepted: 06/26/2017] [Indexed: 11/07/2022]
|
32
|
Musuamba FT, Manolis E, Holford N, Cheung S, Friberg LE, Ogungbenro K, Posch M, Yates J, Berry S, Thomas N, Corriol-Rohou S, Bornkamp B, Bretz F, Hooker AC, Van der Graaf PH, Standing JF, Hay J, Cole S, Gigante V, Karlsson K, Dumortier T, Benda N, Serone F, Das S, Brochot A, Ehmann F, Hemmings R, Rusten IS. Advanced Methods for Dose and Regimen Finding During Drug Development: Summary of the EMA/EFPIA Workshop on Dose Finding (London 4-5 December 2014). CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2017; 6:418-429. [PMID: 28722322 PMCID: PMC5529745 DOI: 10.1002/psp4.12196] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/21/2016] [Revised: 03/27/2017] [Accepted: 03/27/2017] [Indexed: 02/05/2023]
Abstract
Inadequate dose selection for confirmatory trials is currently still one of the most challenging issues in drug development, as illustrated by high rates of late‐stage attritions in clinical development and postmarketing commitments required by regulatory institutions. In an effort to shift the current paradigm in dose and regimen selection and highlight the availability and usefulness of well‐established and regulatory‐acceptable methods, the European Medicines Agency (EMA) in collaboration with the European Federation of Pharmaceutical Industries Association (EFPIA) hosted a multistakeholder workshop on dose finding (London 4–5 December 2014). Some methodologies that could constitute a toolkit for drug developers and regulators were presented. These methods are described in the present report: they include five advanced methods for data analysis (empirical regression models, pharmacometrics models, quantitative systems pharmacology models, MCP‐Mod, and model averaging) and three methods for study design optimization (Fisher information matrix (FIM)‐based methods, clinical trial simulations, and adaptive studies). Pairwise comparisons were also discussed during the workshop; however, mostly for historical reasons. This paper discusses the added value and limitations of these methods as well as challenges for their implementation. Some applications in different therapeutic areas are also summarized, in line with the discussions at the workshop. There was agreement at the workshop on the fact that selection of dose for phase III is an estimation problem and should not be addressed via hypothesis testing. Dose selection for phase III trials should be informed by well‐designed dose‐finding studies; however, the specific choice of method(s) will depend on several aspects and it is not possible to recommend a generalized decision tree. There are many valuable methods available, the methods are not mutually exclusive, and they should be used in conjunction to ensure a scientifically rigorous understanding of the dosing rationale.
Collapse
Affiliation(s)
- F T Musuamba
- EMA Modelling and Simulation Working Group, London, UK.,Federal Agency for Medicines and Health Products, Brussels, Belgium.,UMR850 INSERM, Université de Limoges, Limoges, France
| | - E Manolis
- EMA Modelling and Simulation Working Group, London, UK.,European Medicines Agency, London, UK
| | - N Holford
- Department of Pharmacology & Clinical Pharmacology, University of Auckland, Auckland, New Zealand
| | | | | | | | - M Posch
- Center for Medical Statistics, Informatics and Intelligent Systems, Medical University of Vienna, Vienna, Austria
| | | | - S Berry
- Berry consultants, Austin, Texas, USA
| | | | | | | | - F Bretz
- Center for Medical Statistics, Informatics and Intelligent Systems, Medical University of Vienna, Vienna, Austria.,Novartis, London, UK
| | | | - P H Van der Graaf
- Leiden Academic Centre for Drug Research, Leiden, The Netherlands.,Certara QSP, Canterbury, UK
| | - J F Standing
- EMA Modelling and Simulation Working Group, London, UK.,University College London, London, UK
| | - J Hay
- EMA Modelling and Simulation Working Group, London, UK.,Medicines and Healthcare Products Regulatory Agency, London, UK
| | - S Cole
- EMA Modelling and Simulation Working Group, London, UK.,Medicines and Healthcare Products Regulatory Agency, London, UK
| | - V Gigante
- EMA Modelling and Simulation Working Group, London, UK.,Agenzia Italiana del Farmaco, Roma, Italy
| | - K Karlsson
- EMA Modelling and Simulation Working Group, London, UK.,Medical Products Agency, Uppsala, Sweden
| | | | - N Benda
- EMA Modelling and Simulation Working Group, London, UK.,Bundesinstitut für Arzneimittel und Medizinprodukte, Bonn, Germany
| | - F Serone
- EMA Modelling and Simulation Working Group, London, UK.,Agenzia Italiana del Farmaco, Roma, Italy
| | - S Das
- AstraZeneca UK Limited, London, UK
| | | | - F Ehmann
- European Medicines Agency, London, UK
| | - R Hemmings
- Medicines and Healthcare Products Regulatory Agency, London, UK
| | - I Skottheim Rusten
- EMA Modelling and Simulation Working Group, London, UK.,Norvegian Medicines Agency, Oslo, Norway
| |
Collapse
|
33
|
Thomas N, Roy D. Analysis of Clinical Dose–Response in Small-Molecule Drug Development: 2009–2014. Stat Biopharm Res 2017. [DOI: 10.1080/19466315.2016.1256229] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
| | - Dooti Roy
- Boehringer-Ingelheim, Ridgefield, CT
| |
Collapse
|
34
|
Rekowski J, Köllmann C, Bornkamp B, Ickstadt K, Scherag A. Phase II dose-response trials: A simulation study to compare analysis method performance under design considerations. J Biopharm Stat 2017; 27:885-901. [PMID: 28362145 DOI: 10.1080/10543406.2017.1293078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
Phase II trials are intended to provide information about the dose-response relationship and to support the choice of doses for a pivotal phase III trial. Recently, new analysis methods have been proposed to address these objectives, and guidance is needed to select the most appropriate analysis method in specific situations. We set up a simulation study to evaluate multiple performance measures of one traditional and three more recent dose-finding approaches under four design options and illustrate the investigated analysis methods with an example from clinical practice. Our results reveal no general recommendation for a particular analysis method across all design options and performance measures. However, we also demonstrate that the new analysis methods are worth the effort compared to the traditional ANOVA-based approach.
Collapse
Affiliation(s)
- Jan Rekowski
- a Institute for Medical Informatics, Biometry and Epidemiology , University of Duisburg-Essen , Germany
| | | | - Björn Bornkamp
- c Statistical Methodology , Novartis Pharma AG , Basel , Switzerland
| | | | - André Scherag
- d Clinical Epidemiology, Center for Sepsis Control and Care , University Hospital Jena , Germany
| |
Collapse
|
35
|
Hampson LV, Fisch R, Van LM, Jaki T. Asymmetric inner wedge group sequential tests with applications to verifying whether effective drug concentrations are similar in adults and children. Stat Med 2017; 36:426-441. [PMID: 27859519 DOI: 10.1002/sim.7154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2015] [Revised: 09/19/2016] [Accepted: 10/03/2016] [Indexed: 11/08/2022]
Abstract
Extrapolating from information available on one patient group to support conclusions about another is common in clinical research. For example, the findings of clinical trials, often conducted in highly selective patient cohorts, are routinely extrapolated to wider populations by policy makers. Meanwhile, the results of adult trials may be used to support conclusions about the effects of a medicine in children. For example, if the effective concentration of a drug can be assumed to be similar in adults and children, an appropriate paediatric dosing rule may be found by 'bridging', that is, by matching the adult effective concentration. However, this strategy may result in children receiving an ineffective or hazardous dose if, in fact, effective concentrations differ between adults and children. When there is uncertainty about the equality of effective concentrations, some pharmacokinetic-pharmacodynamic data may be needed in children to verify that differences are small. In this paper, we derive optimal group sequential tests that can be used to verify this assumption efficiently. Asymmetric inner wedge tests are constructed that permit early stopping to accept or reject an assumption of similar effective drug concentrations in adults and children. Asymmetry arises because the consequences of under- and over-dosing may differ. We show how confidence intervals can be obtained on termination of these tests and illustrate the small sample operating characteristics of designs using simulation. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Lisa V Hampson
- Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics and Statistics, Lancaster University, Lancaster, LA1 4YF, U.K
| | - Roland Fisch
- Biostatistical Science and Pharmacometrics, Novartis Pharma AG, Basel, CH-4002, Switzerland
| | - Linh M Van
- Biostatistical Science and Pharmacometrics, Novartis Pharmaceutical, Cambridge, 02139, MA, U.S.A
| | - Thomas Jaki
- Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics and Statistics, Lancaster University, Lancaster, LA1 4YF, U.K
| |
Collapse
|
36
|
Wu J, Banerjee A, Jin B, Menon SM, Martin SW, Heatherington AC. Clinical dose-response for a broad set of biological products: A model-based meta-analysis. Stat Methods Med Res 2017; 27:2694-2721. [PMID: 28067121 DOI: 10.1177/0962280216684528] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Characterizing clinical dose-response is a critical step in drug development. Uncertainty in the dose-response model when planning a dose-ranging study can often undermine efficiency in both the design and analysis of the trial. Results of a previous meta-analysis on a portfolio of small molecule compounds from a large pharmaceutical company demonstrated a consistent dose-response relationship that was well described by the maximal effect model. Biologics are different from small molecules due to their large molecular sizes and their potential to induce immunogenicity. A model-based meta-analysis was conducted on the clinical efficacy of 71 distinct biologics evaluated in 91 placebo-controlled dose-response studies published between 1995 and 2014. The maximal effect model, arising from receptor occupancy theory, described the clinical dose-response data for the majority of the biologics (81.7%, n = 58). Five biologics (7%) with data showing non-monotonic trend assuming the maximal effect model were identified and discussed. A Bayesian model-based hierarchical approach using different joint specifications of prior densities for the maximal effect model parameters was used to meta-analyze the whole set of biologics excluding these five biologics ( n = 66). Posterior predictive distributions of the maximal effect model parameters were reported and they could be used to aid the design of future dose-ranging studies. Compared to the meta-analysis of small molecules, the combination of fewer doses, narrower dosing ranges, and small sample sizes further limited the information available to estimate clinical dose-response among biologics.
Collapse
Affiliation(s)
- Joseph Wu
- 1 Biometrics and Data Management, Global Product Development, Groton, CT, USA
| | - Anindita Banerjee
- 2 Early Clinical Development, Worldwide Research & Development, Cambridge, MA, USA
| | - Bo Jin
- 2 Early Clinical Development, Worldwide Research & Development, Cambridge, MA, USA
| | - Sandeep M Menon
- 3 Statistical Research Consulting Center, Global Product Development, Cambridge, MA, USA
| | - Steven W Martin
- 4 Pharmacometrics, Global Product Development, Cambridge, MA, USA
| | - Anne C Heatherington
- 2 Early Clinical Development, Worldwide Research & Development, Cambridge, MA, USA
| |
Collapse
|
37
|
Affiliation(s)
- Bo Jin
- Early Oncology Development and Clinical Research, Pfizer, Inc., Cambridge, MA, USA
| | - Kerry B. Barker
- Early Oncology Development and Clinical Research, Pfizer, Inc., Cambridge, MA, USA
| |
Collapse
|
38
|
Crippa A, Orsini N. Dose-response meta-analysis of differences in means. BMC Med Res Methodol 2016; 16:91. [PMID: 27485429 PMCID: PMC4971698 DOI: 10.1186/s12874-016-0189-0] [Citation(s) in RCA: 80] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2015] [Accepted: 07/13/2016] [Indexed: 11/10/2022] Open
Abstract
Background Meta-analytical methods are frequently used to combine dose-response findings expressed in terms of relative risks. However, no methodology has been established when results are summarized in terms of differences in means of quantitative outcomes. Methods We proposed a two-stage approach. A flexible dose-response model is estimated within each study (first stage) taking into account the covariance of the data points (mean differences, standardized mean differences). Parameters describing the study-specific curves are then combined using a multivariate random-effects model (second stage) to address heterogeneity across studies. Results The method is fairly general and can accommodate a variety of parametric functions. Compared to traditional non-linear models (e.g. Emax, logistic), spline models do not assume any pre-specified dose-response curve. Spline models allow inclusion of studies with a small number of dose levels, and almost any shape, even non monotonic ones, can be estimated using only two parameters. We illustrated the method using dose-response data arising from five clinical trials on an antipsychotic drug, aripiprazole, and improvement in symptoms in shizoaffective patients. Using the Positive and Negative Syndrome Scale (PANSS), pooled results indicated a non-linear association with the maximum change in mean PANSS score equal to 10.40 (95 % confidence interval 7.48, 13.30) observed for 19.32 mg/day of aripiprazole. No substantial change in PANSS score was observed above this value. An estimated dose of 10.43 mg/day was found to produce 80 % of the maximum predicted response. Conclusion The described approach should be adopted to combine correlated differences in means of quantitative outcomes arising from multiple studies. Sensitivity analysis can be a useful tool to assess the robustness of the overall dose-response curve to different modelling strategies. A user-friendly R package has been developed to facilitate applications by practitioners.
Collapse
Affiliation(s)
- Alessio Crippa
- Department of Public Health Sciences, Karolinska Institutet, Stockholm, Sweden.
| | - Nicola Orsini
- Department of Public Health Sciences, Karolinska Institutet, Stockholm, Sweden
| |
Collapse
|
39
|
Thomas N, Harel O, Little RJA. Analyzing clinical trial outcomes based on incomplete daily diary reports. Stat Med 2016; 35:2894-906. [PMID: 26888661 DOI: 10.1002/sim.6890] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2015] [Revised: 12/19/2015] [Accepted: 01/04/2016] [Indexed: 11/07/2022]
Abstract
A case study is presented assessing the impact of missing data on the analysis of daily diary data from a study evaluating the effect of a drug for the treatment of insomnia. The primary analysis averaged daily diary values for each patient into a weekly variable. Following the commonly used approach, missing daily values within a week were ignored provided there was a minimum number of diary reports (i.e., at least 4). A longitudinal model was then fit with treatment, time, and patient-specific effects. A treatment effect at a pre-specified landmark time was obtained from the model. Weekly values following dropout were regarded as missing, but intermittent daily missing values were obscured. Graphical summaries and tables are presented to characterize the complex missing data patterns. We use multiple imputation for daily diary data to create completed data sets so that exactly 7 daily diary values contribute to each weekly patient average. Standard analysis methods are then applied for landmark analysis of the completed data sets, and the resulting estimates are combined using the standard multiple imputation approach. The observed data are subject to digit heaping and patterned responses (e.g., identical values for several consecutive days), which makes accurate modeling of the response data difficult. Sensitivity analyses under different modeling assumptions for the data were performed, along with pattern mixture models assessing the sensitivity to the missing at random assumption. The emphasis is on graphical displays and computational methods that can be implemented with general-purpose software. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
| | - Ofer Harel
- University of Connecticut, Mansfield, CT, U.S.A
| | | |
Collapse
|
40
|
Schorning K, Bornkamp B, Bretz F, Dette H. Model selection versus model averaging in dose finding studies. Stat Med 2016; 35:4021-40. [PMID: 27226147 DOI: 10.1002/sim.6991] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2015] [Revised: 04/13/2016] [Accepted: 04/17/2016] [Indexed: 11/08/2022]
Abstract
A key objective of Phase II dose finding studies in clinical drug development is to adequately characterize the dose response relationship of a new drug. An important decision is then on the choice of a suitable dose response function to support dose selection for the subsequent Phase III studies. In this paper, we compare different approaches for model selection and model averaging using mathematical properties as well as simulations. We review and illustrate asymptotic properties of model selection criteria and investigate their behavior when changing the sample size but keeping the effect size constant. In a simulation study, we investigate how the various approaches perform in realistically chosen settings. Finally, the different methods are illustrated with a recently conducted Phase II dose finding study in patients with chronic obstructive pulmonary disease. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Kirsten Schorning
- Fakultät für Mathematik, Ruhr-Universität Bochum, Bochum, 44780, Germany
| | - Björn Bornkamp
- Novartis Pharma AG, Lichtstrasse 35, 4002, Basel, Switzerland
| | - Frank Bretz
- Novartis Pharma AG, Lichtstrasse 35, 4002, Basel, Switzerland
| | - Holger Dette
- Fakultät für Mathematik, Ruhr-Universität Bochum, Bochum, 44780, Germany
| |
Collapse
|