Math 2.0

Statistics Definition

Optimal experimental design

Builds on Design of experimentsGeneral linear modelDeterminantEigenvalueCovariance matrix

Definition

Let 𝒳\mathcal{X} be a compact set, the design region, and let f ⁣:𝒳→ℝpf \colon \mathcal{X} \to \R^p be continuous. Consider the linear model y=f(x)Tβ+εy = f(x)^T \beta + \varepsilon with unknown β∈ℝp\beta \in \R^p and uncorrelated errors of mean 00 and common variance σ2\sigma^2. An approximate design is a probability measure ξ\xi on 𝒳\mathcal{X} with finite support, and its information matrix is M(ξ)=∫𝒳f(x)f(x)T dξ(x)=∑x∈supp⁡ξξ({x})f(x)f(x)T.M(\xi) = \int_{\mathcal{X}} f(x) f(x)^T \, d\xi(x) = \sum_{x \in \operatorname{supp} \xi} \xi(\set{x}) f(x) f(x)^T. A design ξ∗\xi^* is:

  1. D-optimal if it maximises det⁡M(ξ)\det M(\xi) over all approximate designs;
  2. A-optimal if M(ξ∗)M(\xi^*) is invertible and it minimises tr⁡M(ξ)−1\tr M(\xi)^{-1} over the designs with M(ξ)M(\xi) invertible;
  3. E-optimal if it maximises the smallest eigenvalue of M(ξ)M(\xi);
  4. G-optimal if it minimises max⁡x∈𝒳f(x)TM(ξ)−1f(x)\max_{x \in \mathcal{X}} f(x)^T M(\xi)^{-1} f(x) over the designs with M(ξ)M(\xi) invertible.

Remarks

For an exact design with nn runs at x1,…,xnx_1, \dots, x_n and design matrix XX with rows f(xi)Tf(x_i)^T, one has M=1nXTXM = \frac{1}{n} X^T X. If MM is invertible, the least squares estimator has covariance matrix σ2nM−1\frac{\sigma^2}{n} M^{-1}. So a D-optimal design minimises the volume of the confidence ellipsoid for β\beta under normal errors.

Example: 𝒳=[−1,1]\mathcal{X} = [-1, 1] and f(x)=(1,x)Tf(x) = (1, x)^T, a straight line. With moments m1=∫x dξm_1 = \int x \, d\xi and m2=∫x2 dξm_2 = \int x^2 \, d\xi, one gets det⁡M(ξ)=m2−m12\det M(\xi) = m_2 - m_1^2, the variance of xx under ξ\xi. On [−1,1][-1, 1] this variance is at most 11, with equality exactly for the design with mass 1/21/2 at −1-1 and 1/21/2 at 11, which is therefore D-optimal.

Non-example: the design with all mass at x=0x = 0 has M=(1000)M = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}, which is singular, so det⁡M=0\det M = 0 and the slope cannot be estimated; it fails clause 1.

Reading path

To understand this page from the ground up, read these 6 pages in order. The foundations come first.

Show the path as a graph
Each page uses the pages above it. The colour of the bar on a box shows how strong the evidence is, and the box also gives it as text. The list has the same pages.
  • no evidence level
  1. Covariance matrix concept
  2. Design of experiments concept
  3. Determinant concept
  4. Eigenvalue concept
  5. General linear model concept
  6. Optimal experimental design definition

Further reading

  1. Web https://en.wikipedia.org/wiki/Optimal_experimental_design
About this page
TagMKAMVQJXBP
URIm2://MKAMVQJXBP
Kinddefinition
WritingW0 stub
MSC62K05
DefinesOptimal experimental design
UsesDesign of experiments, General linear model, Determinant, Eigenvalue, Covariance matrix

Markdown · JSON · Edit on GitHub

What this page needs for W1

  • an inventory row: the row of the inventory that the node accounts for (0 of 1)
  • a source with a locator (0 of 1)

The writing levels say how complete the text of a page is. They do not measure its truth.

What links here

No page links to this page yet.