Pharmaceutical Statistics最新文献

Drug-drug interaction (DDI) trials are an important part of drug development as they provide evidence on the benefits and risks when two or more drugs are taken concomitantly. Sample size calculation is typically recommended to be based on the existence of clinically justified no-effect boundaries but these are challenging to define in practice, while the default no-effect boundaries of 0.8-1.25 are known to be overly conservative requiring a large sample size. In addition, no-effect boundaries are of little use when there is prior pharmacological evidence that a mild or moderate interaction between two drugs may be present, in which case effect boundaries would be more useful. We introduce precision-based sample size calculation that accounts for both the stochastic nature of the pharmacokinetic parameters and the anticipated width of (no-)effect boundaries, should these exist. The methodology is straightforward, requires considerably less sample size and has favorable operating characteristics. A case study on statins is presented to illustrate the ideas.

The Bayesian logistic regression method (BLRM) is a widely adopted and flexible design for finding the maximum tolerated dose in oncology phase I studies. However, the BLRM design has been criticized in the literature for being overly conservative due to the use of the overdose control rule. Recently, a discussion paper titled "Improving the performance of Bayesian logistic regression model with overall control in oncology dose-finding studies" in Statistics in Medicine has proposed an overall control rule to address the "excessive conservativeness" of the standard BLRM design. In this short communication, we discuss the relative conservativeness of the standard BLRM design and also suggest a dose-switching rule to further enhance its performance.

When the distributions of treatment effect modifiers differ between a randomized trial and an external target population, the sample average treatment effect in the trial may be substantially different from the target population average treatment, and accurate estimation of the latter requires adjusting for the differential distribution of effect modifiers. Despite the increasingly rich literature on transportability, little attention has been devoted to methods for transporting trial results to estimate counterfactual survival functions in target populations, when the primary outcome is time to event and subject to right censoring. In this article, we study inverse probability weighting and doubly robust estimators to estimate counterfactual survival functions and the target average survival treatment effect in the target population, and provide their respective approximate variance estimators. We focus on a common scenario where the target population information is observed only through a complex survey, and elucidate how the survey weights can be incorporated into each estimator we considered. Simulation studies are conducted to examine the finite-sample performances of the proposed estimators in terms of bias, efficiency and coverage, under both correct and incorrect model specifications. Finally, we apply the proposed method to assess transportability of the results in the Action to Control Cardiovascular Risk in Diabetes-Blood Pressure (ACCORD-BP) trial to all adults with Diabetes in the United States.

Motivated by the need to model dose-response or dose-toxicity curves in clinical trials, we develop a new horseshoe-based prior for Bayesian isotonic regression modeling a binary outcome against an ordered categorical predictor, where the probability of the outcome is assumed to be monotonically non-decreasing with the predictor. The set of differences between outcome probabilities in consecutive categories of the predictor is equipped with a multivariate prior having support over simplex. The Dirichlet distribution, which can be derived from a normalized sum of independent gamma-distributed random variables, is a natural choice of prior, but using mathematical and simulation-based arguments, we show that the resulting posterior is prone to underflow and other numerical instabilities, even under simple data configurations. We propose an alternative prior based on horseshoe-type shrinkage that is numerically more stable. We show that this horseshoe-based prior is not subject to the numerical instability seen in the Dirichlet/gamma-based prior and that the horseshoe-based posterior can estimate the underlying true curve more efficiently than the Dirichlet-based one. We demonstrate the use of this prior in a model predicting the occurrence of radiation-induced lung toxicity in lung cancer patients as a function of dose delivered to normal lung tissue. Our methodology is implemented in the R package isotonicBayes and therefore suitable for use in the design of dose-finding studies or other dose-response modeling contexts.

As an alternative to the Frequentist p-value, the Bayes factor (or ratio of marginal likelihoods) has been regarded as one of the primary tools for Bayesian hypothesis testing. In recent years, several researchers have begun to re-analyze results from prominent medical journals, as well as from trials for FDA-approved drugs, to show that Bayes factors often give divergent conclusions from those of p-values. In this paper, we investigate the claim that Bayes factors are straightforward to interpret as directly quantifying the relative strength of evidence. In particular, we show that for nested hypotheses with consistent priors, the Bayes factor for the null over the alternative hypothesis is the posterior mean of the likelihood ratio. By re-analyzing 39 results previously published in the New England Journal of Medicine, we demonstrate how the posterior distribution of the likelihood ratio can be computed and visualized, providing useful information beyond the posterior mean alone.

For novel immuno-oncology therapies, the primary purpose of a dose-finding trial is to identify an optimal dose (OD), defined as the tolerable dose having adequate efficacy and immune response under the unpredictable dose-outcome (toxicity, efficacy, and immune response) relationships. In addition, the multiple low or moderate-grade toxicities rather than dose-limiting toxicities (DLTs) and multiple levels of efficacy should be evaluated differently in dose-finding to determine true OD for developing novel immuno-oncology therapies. We proposed a generalized Bayesian optimal interval design for immunotherapy, simultaneously considering efficacy and toxicity grades and immune response outcomes. The proposed design, named gBOIN-ETI design, is model-assisted and easy to implement to develop immunotherapy efficiently. The operating characteristics of the gBOIN-ETI are compared with other dose-finding trial designs in oncology by simulation across various realistic settings. Our simulations show that the gBOIN-ETI design could outperform the other available approaches in terms of both the percentage of correct OD selection and the average number of patients allocated to the OD across various realistic trial settings.

Subgroup analysis may be used to investigate treatment effect heterogeneity among subsets of the study population defined by baseline characteristics. Several methodologies have been proposed in recent years and with these, statistical issues such as multiplicity, complexity, and selection bias have been widely discussed. Some methods adjust for one or more of these issues; however, few of them discuss or consider the stability of the subgroup assignments. We propose exploring the stability of subgroups as a sensitivity analysis step for stratified medicine to assess the robustness of the identified subgroups besides identifying possible factors that may drive this instability. After applying Bayesian credible subgroups, a nonparametric bootstrap can be used to assess stability at subgroup-level and patient-level. Our findings illustrate that when the treatment effect is small or not so evident, patients are more likely to switch to different subgroups (jumpers) across bootstrap resamples. In contrast, when the treatment effect is large or extremely convincing, patients generally remain in the same subgroup. While the proposed subgroup stability method is illustrated through Bayesian credible subgroups method on time-to-event data, this general approach can be used with other subgroup identification methods and endpoints.

Single-arm trials (SATs), while not preferred, remain in use throughout the drug development cycle. They may be accepted by regulators in particular contexts (e.g., in oncology or rare diseases) when the potential effects of new treatments are very large and placebo treatment is unethical. However, in the postregulatory space, SATs are common, and perhaps even more poorly suited to address the questions of interest. In this manuscript, we review regulatory and HTA positions on SATs; challenges posed by SATs to address research questions beyond regulators, evolving statistical methods to provide context for SATs, case studies where SATs could and could not address questions of interest, and communication strategies to influence decision making and optimize study design to address evidence needs.

Recombinant adeno-associated virus (AAV) has become a popular platform for many gene therapy applications. The strength of AAV-based products is a critical quality attribute that affects the efficacy of the drug and is measured as the concentration of vector genomes, or physical titer. Because the dosing of patients is based on the titer measurement, it is critical for manufacturers to ensure that the measured titer of the drug product is close to the actual concentration of the batch. Historically, dosing calculations have been performed using the measured titer, which is reported on the drug product label. However, due to recent regulatory guidance, sponsors are now expected to label the drug product with nominal or "target" titer. This new expectation for gene therapy products can pose a challenge in the presence of process and analytical variability. In particular, the manufacturer must decide if a dilution of the drug substance is warranted at the drug product stage to bring the strength in line with the nominal value. In this paper, we present two straightforward statistical methods to aid the manufacturer in the dilution decision. These approaches use the understanding of process and analytical variability to compute probabilities of achieving the desired drug product titer. We also provide an approach for determining an optimal assay replication strategy for achieving the desired probability of meeting drug product release specifications.

Delayed outcome is common in phase I oncology clinical trials. It causes logistic difficulty, wastes resources, and prolongs the trial duration. This article investigates this issue and proposes the time-to-event 3 + 3 (T3 + 3) design, which utilizes the actual follow-up time for at-risk patients with pending toxicity outcomes. The T3 + 3 design allows continuous accrual without unnecessary trial suspension and is costless and implementable with pretabulated dose decision rules. Besides, the T3 + 3 design uses the isotonic regression to estimate the toxicity rates across dose levels and therefore can accommodate for any targeted toxicity rate for maximum tolerated dose (MTD). It dramatically facilitates the trial preparation and conduct without intensive computation and statistical consultation. The extension to other algorithm-based phase I dose-finding designs (e.g., i3 + 3 design) is also studied. Comprehensive computer simulation studies are conducted to investigate the performance of the T3 + 3 design under various dose-toxicity scenarios. The results confirm that the T3 + 3 design substantially shortens the trial duration compared with the conventional 3 + 3 design and yields much higher accuracy in MTD identification than the rolling six design. In summary, the T3 + 3 design addresses the delayed outcome issue while keeping the desirable features of the 3 + 3 design, such as simplicity, transparency, and costless implementation. It has great potential to accelerate early-phase drug development.