Bayesian experimental design


Bayesian experimental design

Bayesian experimental design provides a general probability-theoretical framework from which other theories on experimental design can be derived. It is based on Bayesian inference to interpret the observations/data acquired during the experiment. This allows accounting for both any prior knowledge on the parameters to be determined as well as uncertainties in observations.

The theory of Bayesian experimental design is to a certain extent based on the theory for making optimal decisions under uncertainty. The aim when designing an experiment is to maximize the expected utility of the experiment outcome. The utility is most commonly defined in terms of a measure of the accuracy of the information provided by the experiment (e.g. the Shannon information or the negative variance), but may also involve factors such as the financial cost of performing the experiment. What will be the optimal experiment design depends on the particular utility criterion chosen.

Contents

Relations to more specialized optimal design theory

Linear theory

If the model is linear, the prior probability density function (PDF) is homogeneous and observational errors are normally distributed, the theory simplifies to the classical optimal experimental design theory.

Approximate normality

In numerous publications on Bayesian experimental design, it is (often implicitly) assumed that all posterior PDFs will be approximately normal. This allows for the expected utility to be calculated using linear theory, averaging over the space of model parameters, an approach reviewed in Chaloner & Verdinelli (1995). Caution must however be taken when applying this method, since approximate normality of all possible posteriors is difficult to verify, even in cases of normal observational errors and uniform prior PDF.

Mathematical formulation

Notation
\theta\, parameters to be determined
y\, observation or data
\xi\, design
p(y|\theta,\xi)\, PDF for making observation y, given parameter values θ and design ξ
p(\theta)\, prior PDF
p(y|\xi)\, marginal PDF in observation space
p(\theta | y, \xi)\,    posterior PDF
U(\xi)\,    utility of the design ξ
U(y, \xi)\,    utility of the experiment outcome after observation y with design ξ

Given a vector θ of parameters to determine, a prior PDF p(θ) over those parameters and a PDF p(y | θ,ξ) for making observation y, given parameter values θ and an experiment design ξ, the posterior PDF can be calculated using Bayes' theorem

p(\theta | y, \xi) = \frac{p(y | \theta, \xi) p(\theta)}{p(y | \xi)}  \, ,

where p(y | ξ) is the marginal probability density in observation space

p(y|\xi) = \int{p(\theta)p(y|\theta,\xi)d\theta}\, .

The expected utility of an experiment with design ξ can then be defined

U(\xi)=\int{p(y|\xi)U(y,\xi)dy}\, ,

where U(y,ξ) is some real-valued functional of the posterior PDF p(θ | y,ξ) after making observation y using an experiment design ξ.

Gain in Shannon information as utility

If the utility is defined as the prior-posterior gain in Shannon information

 U(y, \xi) = \int{\log(p(\theta | y, \xi))p(\theta | y, \xi)d\theta} - \int{\log(p(\theta))p(\theta)d\theta} \, .

Lindley (1956) noted that the expected utility will then be coordinate-independent and can be written in two forms


\begin{alignat}{2}
 U(\xi) & = \int{\int{\log(p(\theta | y,\xi))p(\theta, y | \xi)d\theta}dy} - \int{\log(p(\theta))p(\theta)d\theta} \\
      & = \int{\int{\log(p(y | \theta,\xi))p(\theta, y | \xi)dy}d\theta} - \int{\log(p(y| \xi))p(y| \xi)dy} ,
\end{alignat}
\,

of which the latter can be evaluated without the need for evaluating individual posterior PDFs p(θ | y,ξ) for all possible observations y. Worth noting is that the first term on the second equation line will not depend on the design ξ, as long as the observational uncertainty doesn't. On the other hand, the integral of p(θ)log p(θ) in the first form is constant for all ξ, so if the goal is to choose the design with the highest utility, the term need not be computed at all. Several authors have considered numerical techniques for evaluating and optimizing this criterion, e.g. van den Berg, Curtis & Trampert (2003) and Ryan (2003).

See also

References

  • Ryan, K. J. (2003), "Estimating Expected Information Gains for Experimental Designs With Application to the Random Fatigue-Limit Model", Journal of Computational and Graphical Statistics 12 (3): 585–603, doi:10.1198/1061860032012 

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Design of experiments — In general usage, design of experiments (DOE) or experimental design is the design of any information gathering exercises where variation is present, whether under the full control of the experimenter or not. However, in statistics, these terms… …   Wikipedia

  • Design matrix — In statistics, a design matrix is a matrix of explanatory variables, often denoted by X, that is used in certain statistical models, e.g., the general linear model.[1][2] It can contain indicator variables (ones and zeros) that indicate group… …   Wikipedia

  • Optimal design — This article is about the topic in the design of experiments. For the topic in optimal control theory, see shape optimization. Gustav Elfving developed the optimal design of experiments, and so minimized surveyors need for theodolite measurements …   Wikipedia

  • Randomized block design — In the statistical theory of the design of experiments, blocking is the arranging of experimental units in groups (blocks) that are similar to one another. Typically, a blocking factor is a source of variability that is not of primary interest to …   Wikipedia

  • Completely randomized design — In the design of experiments, completely randomized designs are for studying the effects of one primary factor without the need to take other nuisance variables into account. This article describes completely randomized designs that have one… …   Wikipedia

  • Combinatorial design — theory is the part of combinatorial mathematics that deals with the existence and construction of systems of finite sets whose intersections have specified numerical properties. For instance, a balanced incomplete block design (usually called for …   Wikipedia

  • Regression discontinuity design — In statistics, econometrics, epidemiology and related disciplines, a regression discontinuity design (RDD) is a design that elicits the causal effects of interventions by exploiting a given exogenous threshold determining assignment to treatment …   Wikipedia

  • Clinical study design — is the formulation of trials and experiments in medical and epidemiological research, sometimes known as clinical trials. Many of the considerations here are shared under the more general topic of design of experiments but there can be others, in …   Wikipedia

  • Oscar Kempthorne — Born January 31, 1911(1911 01 31) St. Tudy, Cornwall Died November 15, 2000(2000 11 15) (aged …   Wikipedia

  • List of statistics topics — Please add any Wikipedia articles related to statistics that are not already on this list.The Related changes link in the margin of this page (below search) leads to a list of the most recent changes to the articles listed below. To see the most… …   Wikipedia