Math 2.0

Statistics Design of experiments Theorem

Kiefer–Wolfowitz equivalence theorem

Builds on Optimal experimental designLinear regressionFisher informationDeterminantProbability measure

Statement

Let 𝒳\mathcal{X} be a compact metric space, the design region. Let f ⁣:𝒳→ℝpf \colon \mathcal{X} \to \R^p be continuous, and consider the linear regression model E[y(x)]=θTf(x)E[y(x)] = \theta^T f(x) with θ∈ℝp\theta \in \R^p and uncorrelated errors of equal variance. A design is a probability measure ξ\xi on the Borel sets of 𝒳\mathcal{X}. Its information matrix is M(ξ)=∫𝒳f(x)f(x)T dξ(x)M(\xi) = \int_{\mathcal{X}} f(x) f(x)^T \, d\xi(x). When M(ξ)M(\xi) is nonsingular, its standardized variance function is d(x,ξ)=f(x)TM(ξ)−1f(x)d(x, \xi) = f(x)^T M(\xi)^{-1} f(x). Assume that some design has a nonsingular information matrix. For a design ξ∗\xi^* with M(ξ∗)M(\xi^*) nonsingular, the following three conditions are equivalent:

  1. ξ∗\xi^* is D-optimal: it maximizes det⁡M(ξ)\det M(\xi) over all designs ξ\xi.
  2. ξ∗\xi^* is G-optimal: it minimizes max⁡x∈𝒳d(x,ξ)\max_{x \in \mathcal{X}} d(x, \xi) over all designs ξ\xi with M(ξ)M(\xi) nonsingular.
  3. max⁡x∈𝒳d(x,ξ∗)=p\max_{x \in \mathcal{X}} d(x, \xi^*) = p.

Moreover, all D-optimal designs have the same information matrix, and d(x,ξ∗)=pd(x, \xi^*) = p at every support point xx of a D-optimal design ξ∗\xi^*.

Remarks

For every design with nonsingular information matrix, ∫d(x,ξ) dξ(x)=tr⁡(M(ξ)−1M(ξ))=p\int d(x, \xi) \, d\xi(x) = \tr(M(\xi)^{-1} M(\xi)) = p. So max⁡xd(x,ξ)≥p\max_x d(x, \xi) \ge p always holds, and this explains the number pp in condition 3. The function ξ↦log⁡det⁡M(ξ)\xi \mapsto \log \det M(\xi) is concave on the convex set of designs, and its directional derivative at ξ\xi towards the point mass δx\delta_x equals d(x,ξ)−pd(x, \xi) - p. So condition 3 says that no direction increases log⁡det⁡M\log \det M, which for a concave function means a maximum. Example: for the line f(x)=(1,x)Tf(x) = (1, x)^T on 𝒳=[−1,1]\mathcal{X} = [-1, 1], the design ξ∗=12δ−1+12δ1\xi^* = \frac12 \delta_{-1} + \frac12 \delta_{1} has M(ξ∗)=I2M(\xi^*) = I_2. Then d(x,ξ∗)=1+x2≤2=pd(x, \xi^*) = 1 + x^2 \le 2 = p, so ξ∗\xi^* is D-optimal and G-optimal. The theorem was proved by Jack Kiefer and Jacob Wolfowitz in 1960.

Going deeper

For a student

The equivalence theorem turns a hard optimization over probability measures into a checkable condition: compute d(x,ξ)d(x, \xi) and see whether it stays below pp. It is the basis of the Fedorov–Wynn type algorithms, which add weight at the point where d(x,ξ)d(x, \xi) is largest. Similar equivalence theorems hold for other concave criteria of optimal experimental design, with the variance function replaced by the corresponding directional derivative.

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. Determinant concept
  2. Fisher information concept
  3. Linear regression concept
  4. Optimal experimental design concept
  5. Probability measure concept
  6. Kiefer–Wolfowitz equivalence theorem statement

Further reading

  1. Web https://en.wikipedia.org/wiki/Optimal_experimental_design
About this page
Tag71XBFB49HE
URIm2://71XBFB49HE
Kindstatement, theorem
WritingW0 stub
MSC62K05
StatusThis copy of the site has no build data. Run m2 status 71XBFB49HE.
AboutOptimal experimental design
UsesOptimal experimental design, Linear regression, Fisher information, Determinant, Probability measure

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.