Abstract:Quality control in industrial processes is increasingly making use of prior scientific knowledge, often encoded in physical models that require numerical approximation. Statistical prediction, and subsequent optimization, is key to ensuring the process output meets a specification target. However, the numerical expense of approximating the models poses computational challenges to the identification of combinations of the process factors where there is confidence in the quality of the response. Recent work in Bayesian computation and statistical approximation (emulation) of expensive computational models is exploited to develop a novel strategy for optimizing the posterior probability of a process meeting specification. The ensuing methodology is motivated by, and demonstrated on, a chemical synthesis process to manufacture a pharmaceutical product, within which an initial set of substances evolve according to chemical reactions, under certain process conditions, into a series of new substances. One of these substances is a target pharmaceutical product and two are unwanted by-products. The aim is to determine the combinations of process conditions and amounts of initial substances that maximize the probability of obtaining sufficient target pharmaceutical product whilst ensuring unwanted by-products do not exceed a given level. The relationship between the factors and amounts of substances of interest is theoretically described by the solution to a system of ordinary differential equations incorporating temperature dependence. Using data from a small experiment, it is shown how the methodology can approximate the multivariate posterior predictive distribution of the pharmaceutical target and by-products, and therefore identify suitable operating values. 1 1 Materials to replicate the analysis can be found at www.github.com/amo105/chemicalkinetics.

Bayesian hierarchical modeling of drug stability data

A Modern Approach to Stability Studies via Bayesian Linear Mixed Models Incorporating Auxiliary Effects

Bayesian estimation for longitudinal data in a joint model with HPCs

Phase I Trial Design for Drug Combinations with Bayesian Model Averaging.

Bayesian Nonparametric Risk Assessment in Developmental Toxicity Studies with Ordinal Responses

Bayesian aggregation of average data: An application in drug development

Stability Modeling Methodologies to Enable Earlier Patient Access

Optimal Bayesian Hierarchical Model to Accelerate the Development of Tissue-Agnostic Drugs and Basket Trials.

Estimating Shelf Life Through Tolerance Intervals Extended to Nonlinear Response Trends

Development, Comparison, and Qualification of Predictive Pharmaceutical Stability Models Using Independent Statistical Tools

Bayesian joint modeling of chemical structure and dose response curves

An Inventory Model for Deteriorating Drugs with Stochastic Lead Time

A simple Bayesian model to estimate proportions and ratios from count data with a hierarchical error structure with an application to droplet digital PCR experiments

Hierarchical Stochastic Model in Bayesian Inference: Theoretical Implications and Efficient Approximation

Flexible Bayesian Modeling for Longitudinal Binary and Ordinal Responses

Statistical Modeling of Within-Laboratory Precision Using a Hierarchical Bayesian Approach

Bayesian prediction for physical models with application to the optimization of the synthesis of pharmaceutical products using chemical kinetics

A method for Bayesian regression modelling of composition data

Identifying treatment effect heterogeneity in dose-finding trials using Bayesian hierarchical models

Bayesian Analysis of Population Health Data