Article ID: | iaor20164717 |
Volume: | 63 |
Issue: | 6 |
Start Page Number: | 1528 |
End Page Number: | 1546 |
Publication Date: | Dec 2015 |
Journal: | Operations Research |
Authors: | Bjarnadttir Margrt V, Zenios Stefanos A, Goh Joel, Bayati Mohsen |
Keywords: | quality & reliability, medicine |
Postmarketing drug surveillance is the process of monitoring the adverse events of pharmaceutical or medical devices after they are approved by the appropriate regulatory authorities. Historically, such surveillance was based on voluntary reports by medical practitioners, but with the widespread adoption of electronic medical records and comprehensive patient databases, surveillance systems that utilize such data are of considerable interest. Unfortunately, existing methods for analyzing the data in such systems ignore the open‐ended exploratory nature of such systems that requires the assessment of multiple possible adverse events. In this article, we propose a method, SEQMEDS, that assesses the effect of a single drug on multiple adverse events by analyzing data that accumulate sequentially and explicitly captures interdependencies among the multiple events. The method continuously monitors a vector‐valued test‐statistic derived from the cumulative number of adverse events. It flags a potential adverse event once the test‐statistic crosses a stopping boundary. We employ asymptotic analysis that assumes a large number of observations in a given window of time to show how to compute the stopping boundary by solving a convex optimization problem that achieves a desired Type I error and minimizes the expected time to detection under a pre‐specified alternative hypothesis. We apply our method to a model in which the interdependency among the multiple adverse events is captured by a Cox proportional hazards model with time‐dependent covariates and demonstrate that it provides an approximation of a fully sequential test for the maximum hazard ratio of the drug over multiple adverse events. A numerical study verifies that our method delivers Type I /II errors that are below pre‐specified levels and is robust to distributional assumptions and parameter values.