Pharmaceutical Statistics

Immunoassays play an important role in drug development of products targeting the immune system. Consistent quality of the results from an immunoassay is essential to make unbiased and accurate claims about the drug product during preclinical and clinical development stages. Assay qualification and validation shed light on the performance of the assay. It is the first evaluation and the verification, respectively, of the assay's performance. This tutorial explains and illustrates the calculation methodology for important assay qualification parameters including precision, relative accuracy, linearity, the lower limit of quantification (LLOQ), the upper limit of quantification (ULOQ), the assay range and dilutability...

It is well known that medication adherence is critical to patient outcomes and can decrease patient mortality. The Pharmacy Quality Alliance (PQA) has recognized and identified medication adherence as an important indicator of medication-use quality. Hence, there is a need to use the right methods to assess medication adherence. The PQA has endorsed the proportion of days covered (PDC) as the primary method of measuring adherence. Although easy to calculate, the PDC has however several drawbacks as a method of measuring adherence...

We present the motivation, experience, and learnings from a data challenge conducted at a large pharmaceutical corporation on the topic of subgroup identification. The data challenge aimed at exploring approaches to subgroup identification for future clinical trials. To mimic a realistic setting, participants had access to 4 Phase III clinical trials to derive a subgroup and predict its treatment effect on a future study not accessible to challenge participants. A total of 30 teams registered for the challenge with around 100 participants, primarily from Biostatistics organization...

In recently conducted phase III trials in a rare disease area, patients received monthly treatment at a high dose of the drug, which targets to lower a specific biomarker level, closely associated with the efficacy endpoint, to around 10% across patients. Although this high dose demonstrated strong efficacy, treatments were withheld due to the reports of serious adverse events. Dosing in these studies were later resumed at a reduced dosage which targets to lower the biomarker level to 15%-35% across patients...

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...

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...

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...

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...

No abstract text is available yet for this article.

The pharmaceutical industry is plagued with long, costly development and high risk. Therefore, a company's effective management and optimisation of a portfolio of projects is critical for success. Project metrics such as the probability of success enable modelling of a company's pipeline accounting for the high uncertainty inherent within the industry. Making portfolio decisions inherently involves managing risk, and statisticians are ideally positioned to champion not only the derivation of metrics for individual projects, but also advocate decision-making at a broader portfolio level...

Since the publication of ICH E9 (R1), "Addendum to statistical principles for clinical trials: on choosing appropriate estimands and defining sensitivity analyses in clinical trials," there has been a lot of debate about the hypothetical strategy for handling intercurrent events. Arguments against the hypothetical strategy are twofold: (1) the clinical question has limited clinical/regulatory interest; (2) the estimation may need strong statistical assumptions. In this article, we provide an example of a hypothetical strategy handling use of rescue medications in the acute pain setting...

We propose a novel frailty model with change points applying random effects to a Cox proportional hazard model to adjust the heterogeneity between clusters. In the specially focused eight Empowered Action Group (EAG) states in India, there are problems with different survival curves for children up to the age of five in different states. Therefore, when analyzing the survival times for the eight EAG states, we need to adjust for the effects among states (clusters). Because the frailty model includes random effects, the parameters are estimated using the expectation-maximization (EM) algorithm...

Count outcomes are collected in clinical trials for new drug development in several therapeutic areas and the event rate is commonly used as a single primary endpoint. Count outcomes that are greater than the mean value are termed overdispersion; thus, count outcomes are assumed to have a negative binomial distribution. However, in clinical trials for treating asthma and chronic obstructive pulmonary disease (COPD), a regulatory agency has suggested that a continuous endpoint related to lung function must be evaluated as a primary endpoint in addition to the event rate...

The delayed treatment effect is a common feature of immunotherapy, characterized by a gradual onset of action ranging from no effect to full effect. In this study, we propose a generalized delayed treatment effect function to depict the delayed effective process precisely and flexibly. To reduce potential power loss caused by the delayed treatment effect in a group sequential trial, we employ the maximin efficiency robust test, which enhances power robustness across a range of possible delays. We present novel approaches based on the Markov chain method for determining group sequential boundaries, calculating the power function, and estimating the maximum sample size through iterative regressions between the square root of the maximum sample size and the normal quantile of power...

The sum of the longest diameter (SLD) of the target lesions is a longitudinal biomarker used to assess tumor response in cancer clinical trials, which can inform about early treatment effect. This biomarker is semicontinuous, often characterized by an excess of zeros and right skewness. Conditional two-part joint models were introduced to account for the excess of zeros in the longitudinal biomarker distribution and link it to a time-to-event outcome. A limitation of the conditional two-part model is that it only provides an effect of covariates, such as treatment, on the conditional mean of positive biomarker values, and not an overall effect on the biomarker, which is often of clinical relevance...

We compare the performance of nonparametric estimators for the mean number of recurrent events and provide a systematic overview for different recurrent event settings. The mean number of recurrent events is an easily interpreted marginal feature often used for treatment comparisons in clinical trials. Incomplete observations, dependencies between successive events, terminating events acting as competing risk, or gaps between at risk periods complicate the estimation. We use survival multistate models to represent different complex recurrent event situations, profiting from recent advances in nonparametric estimation for non-Markov multistate models, and explain several estimators by using multistate intensity processes, including the common Nelson-Aalen-type estimators with and without competing mortality...

With the advent of cancer immunotherapy, some special features including delayed treatment effect, cure rate, diminishing treatment effect and crossing survival are often observed in survival analysis. They violate the proportional hazard model assumption and pose a unique challenge for the conventional trial design and analysis strategies. Many methods like cure rate model have been developed based on mixture model to incorporate some of these features. In this work, we extend the mixture model to deal with multiple non-proportional patterns and develop its geometric average hazard ratio (gAHR) to quantify the treatment effect...

Cross-over designs are commonly used in randomized clinical trials to estimate efficacy of a new treatment. They have received a lot of attention, particularly in connection with regulatory requirements for new drugs. The main advantage of using cross-over designs over conventional parallel designs is increased precision, thanks to within-subject comparisons. In the statistical literature, more recent developments are discussed in the analysis of cross-over trials, in particular regarding repeated measures...

Multiregional clinical trials (MRCTs) have become increasingly common during the development of new drugs to obtain simultaneous drug approvals worldwide. When planning MRCTs, a major statistical challenge is determination of the regional sample size. In general, the regional sample size must be determined as the sample size such that the regional consistency probability, defined as the probability of meeting the regional consistency criterion, is greater than a prespecified value. The Japanese Ministry of Health, Labour and Welfare proposed two criteria for regional consistency...

Matching reduces confounding bias in comparing the outcomes of nonrandomized patient populations by removing systematic differences between them. Under very basic assumptions, propensity score (PS) matching can be shown to eliminate bias entirely in estimating the average treatment effect on the treated. In practice, misspecification of the PS model leads to deviations from theory and matching quality is ultimately judged by the observed post-matching balance in baseline covariates. Since covariate balance is the ultimate arbiter of successful matching, we argue for an approach to matching in which the success criterion is explicitly specified and describe an evolutionary algorithm to directly optimize an arbitrary metric of covariate balance...