Peter Müller,
Ph.D. Many Bayesian design problems
require the solution of analytically intractable integration and optimization
problems. We will discuss three related approaches. The first is optimal design
by curve fitting of Robert E. Kass,
Ph.D. One of the most important
techniques in learning about brain function has involved examining neuronal
activity in laboratory animals under varying experimental conditions. Neural information is represented and communicated
through series of action potentials, or spike trains, and a fundamental
question in neuroscience is precisely how this is accomplished, that is, what
physiological significance should be attached to a particular neuron firing
pattern. My colleagues and I have framed scientific questions in terms of point
process intensity functions, and have used Bayesian methods to fit the point
process models to neuronal data. I will describe the neurophysiological setting, and then use it as background
to discuss a general approach to curve fitting with free-knot splines and reversible-jump MCMC, which may be applied in
the point process setting. With this analytical
foundation in place I will outline the progress we’ve made and the substantive
problems we hope to address in the next few years. Andy P. Grieve, Ph.D. The context of our case study will be a dose-response study
in the treatment of acute stroke. Such studies are an extremely important part
of the drug development process as knowledge of the relationship between response
and dose is an essential requirement for making informed decisions about
dosage. One potential difficulty in the
use of a limited number of doses in a parallel group design to investigate
dose response is the danger that the steep part of the dose-response curve may
fall between two doses and little is learned. The use of a large number of
doses in order to circumvent this problem is potentially wasteful in its use of
patients because a large number of patients will either be receiving doses
which are little different from placebo or, at the other extreme, receiving doses which
have a greater potential for causing adverse side effects. Ideally, the vast
majority of patients should receive doses in the steepest part of the dose-response
function. The design decided upon was a sequential Bayesian adaptive design in
which knowledge of the dose-response curve was updated on an on-going basis in
order to inform two decisions. First, the dose that should be
allocated to the next patient; second whether the study should continue or should
be stopped. There have been two major hindrances to the use of Bayesian
methods in pharmaceutical R&D, one of which is practical, the other more
philosophical. The practical constraint
has been the lack of availability of methods and software for their
implementation. The philosophical constraint has been a perceived antipathy
from regulators to the use of priors. We will discuss aspects of these issues in the practical implementation
of the proposed design. In particular we
will consider:
Gary L. Rosner,
Sc.D. I will highlight various Bayesian methods
we use to design and analyze population-based studies of anticancer therapy.
In particular, we use Bayesian nonparametric models and combine
information across studies and patients within studies by hierarchical
modeling. The specific
application concerns predicting the optimal dose for a leukemia patient
undergoing high doses of chemotherapy, followed by bone marrow transplantation. The optimization
seeks to find the dose that minimizes the expected loss, where the loss
function associates costs when the area under the concentration-time curve
(AUC), a measure of systemic exposure, is either below or above the limits of a
target range. By first giving a patient
a sub-therapeutic test dose of the anticancer drug, we estimate the
patient-specific pharmacokinetics and use this information to predict the
patient’s AUC as a function of dose. The
study design uses data from the following sources. First, we have historical data on leukemia
patients who underwent the same high-dose chemotherapy. These data consist of pharmacokinetics and
clinical outcomes. A subsequent study
collected pharmacokinetic information on patients who received a fixed low dose
and a non-individualized high-dose of the drug.
In the third study, patients receive a low, test dose of the drug, the
same low dose as in the second study. We
fit a pharmacokinetic model to the concentrations of the drug measured in the
patient after administration of the test dose to infer the patient-specific
parameters in the model of drug disposition.
We determine the optimal dose by averaging a loss function with respect
to the predictive distribution for the patient’s AUC as a function of
dose. The patient will then receive the
optimal dose, which is the dose that is associated with the smallest expected
loss for the patient. Our design
incorporates historical information from the other two studies, along with the
current patient’s data, borrowing strength to improve the precision of the
prediction. Lurdes Y. T. Inoue, Ph.D. A sequential Bayesian phase II/III design is proposed for
comparative clinical trials. The design
is based on both survival time and discrete early events that may be related to
survival through a parametric mixture model.
Phase II involves a small number of centers. Patients are randomized between treatments
throughout, and sequential decisions are based on predictive probabilities of
concluding superiority of the experimental treatment. Whether to stop early, continue, or shift into phase III is assessed
repeatedly in phase II. Phase III begins
when additional institutions are incorporated into the ongoing phase II trial. Using
simulation studies in the context of a non-small-cell lung cancer trial, we will
show that the proposed method maintains overall size and power while usually
requiring substantially smaller sample size and shorter trial duration when
compared to conventional group-sequential phase III designs. Peter F. Thall, Ph.D. Outcome-adaptive
decision-making during an ongoing experiment uses the data that become
available at successive times as a basis for deciding what to do next. This is especially useful in clinical trials,
where the decisions may be what dose or treatment to give the next patient,
whether to drop a treatment arm, or whether to continue or terminate the
trial. The Bayesian paradigm provides a
natural basis for this process. The
posterior is updated repeatedly as new data become available, and decisions are
based on posterior or predictive probabilities.
This talk will consist of several examples of oncology trials using
adaptive Bayesian methods, including (1) a trial of allogeneic
donor lymphocytes for treating relapsed acute leukemia patients in which
adaptive randomization is used to optimize the lymphocyte infusion time of each
patient, (2) a new method for dose-finding in phase I trials where the doses of
two different agents used in combination are varied, and (3) a trial to
determine whether Gleevec has substantive
anti-disease activity in sarcoma that uses a hierarchical model to account for
multiple disease subtypes. Dalene Stangl, Ph.D. Conflicting conclusions between
large clinical trials and meta-analyses fuel a debate about the usefulness of meta-analysis. Much of the controversy is overstated,
stemming from the disregard of naturally occurring variation between
studies. This course will discuss the
controversy, review some basic concepts, and introduce recent methodological
developments in meta-analysis. New
developments will include methods for incorporating 1) inconsistency between designs
and outcomes across studies, 2) study-level covariates (including measures of
study quality), and 3) adjustments for publication bias. Methods will be taught
through the presentation of examples from clinical and community trials,
epidemiology, and health policy. Giovanni Parmigiani,
Ph.D. Over the last decade Bayesian hierarchical
models have been increasingly used in numerous areas,
including clinical trials and epidemiological studies. This is now a well
established methodology for handling study-to-study heterogeneity, small sample
sizes, heterogeneous study designs, publication bias, and other
complexities. The necessary computations
for fitting Bayesian hierarchical models in a wide range of situations can be
carried out conveniently using standard software packages such as BUGS. As results from Bayesian hierarchical models are
increasingly used to support clinical and policy decision making, the issue
arises of whether they provide a sound way for comparing treatments. In this
case study we will consider a meta-analysis of 2x2 tables, each arising from a
study comparing adverse event counts for a treatment arm and a control arm. Our
analysis will highlight strengths as well as important potential limitations of
Bayesian hierarchical approaches, and will emphasize approaches that ensure
robustness of conclusions to modeling choices such as parameterization,
distributional assumptions and prior hyperparameters. Bradley P. Carlin, Ph.D. Survival models have a long
history in the biomedical and biostatistical literature, and are enormously
popular in the analysis of time-to-event data.
Very often these data will be grouped into strata, such as clinical
sites, geographic regions, and so on.
Such data will often be available over multiple time periods, and for multiple
diseases. In this talk we will consider
spatial survival models from two general points of view: {\em
geostatistical approaches}, where we use the exact
geographic locations (e.g., latitude and longitude) of the strata, and {\em lattice approaches}, where we use only the positions of
the strata relative to each other (e.g., which counties neighbor which
others). We will compare these
approaches in the context of a data set on infant mortality in We will then consider
hierarchical spatial process models for multivariate survival data sets which
are spatio-temporally arranged. Such models must
account for correlations between survival rates in neighboring spatial regions,
adjacent time periods, and similar diseases (say, different forms of
cancer). We will investigate Cox semiparametric survival modeling approaches, adding spatial
and temporal effects in a hierarchical structure. Here, a multivariate lattice model for the
region-specific frailties is most convenient. Exemplification will be provided
using time-to-event data for various cancers from the National Cancer
Institute’s Surveillance, Epidemiology, and End Results (SEER) database. Donald A. Berry, Ph.D.
A Bayesian statistical approach to decision making from a medical perspective can be helpful
in pharmaceutical company decision making. I will show how the Bayesian approach can be used to
develop efficient designs for clinical trials. I will apply a decision-analytic
approach to the question of proceeding to the next phase of drug development
and address the optimal design of future trials. I will present case studies of
trial designs that have the goal of maximizing information while minimizing
cost, where cost is measured in terms of time, money and patient resources. |