Beyond p-values: Rethinking Statistical Inference and Strategic Portfolio Stewardship in Clinical Development

For decades, drug development has operated under a implicit compromise. We built clinical trials around fixed-sample sizes, rigid hypothesis testing, and binary p-value thresholds (p < 0.05). While this frequentist foundation provided a clean mathematical framework for an earlier computational era, it often forced us to evaluate dynamic biological questions through static, blunt instruments.

In modern oncology and precision medicine, where therapy is increasingly targeted to small, molecularly defined sub-populations, the limits of legacy frameworks are no longer just academic concerns—they directly affect patient outcomes and portfolio survival. Transitioning from conventional trial execution to adaptive, Bayesian architectures is not simply a technical upgrade in software; it represents a fundamental paradigm shift in how we learn from clinical data and allocate R&D capital.

Rethinking Early Escalation: The Shift from MTD to Optimal Biological Dose

In early-phase oncology, the traditional 3+3 rule-based escalation algorithm was designed around a toxicological assumption: that higher doses inevitably yield higher efficacy, and that our primary goal is finding the Maximum Tolerated Dose (MTD).

With targeted therapies, antibody-drug conjugates (ADCs), and immuno-oncology combinations, that assumption frequently falls apart. Therapeutic efficacy often reaches a biological plateau long before unacceptable toxicity occurs. Pushing cohorts up to the MTD exposes patients to unnecessary toxicity without generating additional clinical benefit.

FDA’s Project Optimus formalizes this reality, catalyzing a transition from traditional toxicological MTD estimation to characterizing the Optimal Biological Dose (OBD). Modern model-assisted and model-based designs—such as Bayesian Optimal Interval (BOIN), modified toxicity probability interval (mTPI, mTPI2), Continual Reassessment Method (CRM), and Bayesian Logistic Regression Modeling (BLRM)—allow us to evaluate joint dose-toxicity utility contours dynamically. By incorporating the modern estimand framework (e.g. designing continuous exposure-response modeling and randomized expansion cohorts with composite endpoints of interest early in development), we establish true biological saturation before committing assets to confirmatory Phase 3 trials.

Resolving the Small Subgroup Dilemma via Dynamic Borrowing

Precision medicine presents a frequentist dilemma when evaluating multiple cohorts. Consider an oncology basket trial testing a targeted agent across six distinct disease indications with small cohort sizes (N = 6 per arm).

Under an unpooled frequentist analysis, each arm is evaluated in isolation. Due to high sampling variance at N = 6, statistical power to detect an Objective Response Rate (ORR > 20%) drops to roughly 42%, leaving active treatment signals vulnerable to false negatives and premature pipeline abandonment. Conversely, lumping all patients into a complete pooling model distorts true biological signals by diluting a highly active cohort with non-responding arms (details upon request).

Bayesian Hierarchical Modeling (BHM) resolves this tension through dynamic partial borrowing. By estimating the between-cohort heterogeneity variance, BHM dynamically adjusts information sharing:

  • When cohort profiles demonstrate clinical consistency (between-cohort variance -> 0), the model borrows precision across active clusters, expanding the Effective Sample Size (ESS) for responding arms (e.g., expanding N = 6 to ESS = 18, a +200% sample size gain) and boosting statistical power from 42% to 88% (details upon request).

  • If a cohort diverges as a non-responder (e.g., 0/6 responses), the heterogeneity parameter automatically severs borrowing, isolating the non-responding arm and triggering an early futility stop.

This ability to share information when biologically justified—and disconnect when cohorts conflict—ensures we neither miss targeted efficacy signals nor subsidize non-performing indications.

Real-World Evidence, Synthetic Control Arms, and AI Governance

In rare indications, biomarker-stratified sub-populations, and advanced cell therapy settings, enrolling concurrent randomized control arms is increasingly unfeasible or ethically challenging. Incorporating Real-World Evidence (RWE) through Synthetic Control Arms (SCA) offers a vital path forward, but only if selection biases and historical control drift are rigorously controlled.

By combining Propensity Score Matching (PSM) or Inverse Probability of Treatment Weighting (IPTW) with Meta-Analytic Predictive (MAP) mixture priors, we can translate matched EHR or registry data into informative control priors for active studies. The key lies in robustification: assigning a heavy-tailed non-informative mixture component to the MAP prior ensures that if unexpected drift occurs between historical controls and concurrent trial data, the model automatically discounts the external data to preserve Type I error control.

As AI-driven algorithms and synthetic data generation accelerate statistical workflows, regulatory agencies are actively establishing dedicated guidance frameworks to oversee their application in clinical submissions. To maintain audit-readiness under 21 CFR Part 11 and align with FDA Good Machine Learning Practice (GMLP) principles, sponsors are establishing complete data provenance, validating end-to-end CDISC ADaM data lineage, and implementing human-in-the-loop verification for AI-assisted statistical workflows prior to database lock.

Biostatistics as Strategic Capital Allocation

Ultimately, modern biostatistics is more than a retroactive reporting function designed to compute p-values at study lock. When structured effectively through master protocols (Basket, Umbrella, Platform) and continuous interim probability monitoring, Biometrics becomes a primary driver of portfolio ROI and capital allocation.

Using Predictive Probability of Success (PPoS) frameworks, executive governance boards can evaluate the conditional probability that a program will meet its target product profile before committing capital to prospective Phase 3 studies. High-probability assets receive immediate acceleration, while failing candidates are halted early through objective futility boundaries. The capital saved from early futility terminations—often tens of millions of dollars per program—can be re-allocated directly into accelerating substantially de-risked assets.

By uniting mathematical rigor, operational discipline, and forward-looking regulatory alignment, we fulfill our ultimate obligation to both clinical science and the patients waiting for novel therapies: making decisions faster, clearer, and with uncompromised integrity.

Simulation Based Clinical Trial Design

Modern clinical trial designs often rely on commercial software platforms. They are easy to use, can present results in a straightforward manner, and bring nice visualizations. The major drawbacks for many may be threefold: 1) the underlying methological assumptions can be like a black box to users, which may lead to the second drawback, 2) lack of flexibility and interpretability from the strategic perspective, and 3) is prone to error without a proper validation process, particularly in high-pressure, fast-paced environments. As the design paradigm steadily shifts to the Bayesian framework, simulation can be a powerful tool to assist advanced clinical trial design and enhance communication and collaboration between stakeholders.

Simulation lies on the distribution assumption of data. Data that follow distributions, such as normal, binomial or Poisson distributions, can be generated by a single function of a typical statistical programming language. When it comes to simulating survival data, a frequently used endpoint of many of today’s clinical trials, it may take a set of comprehensive steps, involving generating two variables, a time and a censoring variable. For the survival data with a mixture distribution or a piecewise hazard function, the steps are beyond simply generating two variables, each following a basic distribution. There are a few choices in R that can be found to assist generating such types of survival data. I have recently run a few design scenarios with the simsurv package, made a side-by-side comparison with the gsDesign package, the gsDesign2 package as well as the EAST® Software by Cytel, and found they have a high degree of alignment. The simsurv R package creators also published a tutorial paper in the Journal of Statistical Software that showed cases of simulating survival data from standard parametric distributions (e.g. exponential, Weibull) to complex situations, such as a user-defined or time-dependent hazard function. These situations may facilitate survival data generated under extreme distributions that may fit in design scenarios for pseudo-progression, delayed treatment effects, treatment crossovers, rare diseases or such.

Once there is a way to generate the desired survival data (time-to-event endpoints), a trial can be designed under either a frequentist or a Bayesian framework. In addition, the trial outcome can be anticipated with the simulated data prior to its initiation, and also be updated as the trial is at an interim look.

Simulated survival data with a piecewise hazard assumption for the control treatment, and non-proportional hazard ratios for the two experimental treatments against the control

Simulated survival data with a piecewise hazard assumption for the control treatment, and non-proportional hazard ratios for the two experimental treatments against the control.


A Real-world, Comparative Study of FDA- approved Diagnostic Assays PD-L1 IHC 28-8 and 22C3 in Lung Cancer and Other Malignancies

Between different types of diagnostic tests, different versions of a diagnostic assay, statistical agreement assessments are often needed, in order to bridge the assay performance or the clinical utility from one test to another, or to upgrade an outdated diagnostic test. The statistical methods involved in agreement assessments include Positive Percent Agreement (PPA), Negative Percent Agreement (NPA), Overall Percent Agreement (OPA) for qualitative assays, and Bland-Altman plots, Passing Bablok or Deming regression for quantitative assays, etc.

A paper that addressed a high degree of agreement between 2 FDA-approved PD-L1 diagnostic assays, 28-8 and 22C3 showed some examples of statistical agreement assessments.