Modeling Monotone Nonlinear Disease Progression and Checking the Correctness of the Associated Software

Donald B. Rubin,Samantha Rachel Cook

2004-01-01

Abstract:Management of chronic diseases often involves monitoring one or more outcomes over time. Modeling long term disease progression in untreated patients can be useful for drug development as well as patient monitoring. The outcomes of interest in chronic diseases often have expected monotone or approximately monotone progression over time, possibly preceded by a period during which the outcome measurement is relatively constant. Many diseases have some treatment available for current patients, if only to relieve some symptoms or to slow progression, making disease modeling in untreated current patients difficult. Then historical data can often be useful for developing models for disease progression without treatment. The first part of this thesis presents a model for disease progression in untreated patients suffering from Fabry disease, an X-linked recessive disorder. We develop a Bayesian hierarchical random effects regression model for serum creatinine, a critically important blood measurement for Fabry patients, using a historical data base compiled from medical records and patient registries. We then use this model and results from the historical patients' data to aid in multiply imputing missing placebo data from a clinical trial for a new drug specifically developed to treat Fabry disease; data are missing when placebo patients drop out or switch to the new drug, and thus imputation using only the current untreated patients would rely solely on extrapolation. The second and third parts of this thesis expand upon two topics encountered in developing the disease progression models considered in Part One. Part Two outlines a general strategy for testing the correctness of software for fitting Bayesian models, capitalizing on properties of Bayesian posterior distributions. Software testing becomes a critically important enterprise with models fit using Markov Chain Monte Carlo methods. Part Three presents a method for constructing vague prior distributions for complex Bayesian models, based on the asymptotic normality of the likelihood function. Throughout, we illustrate concepts and techniques using the disease progression model developed in Part One.