core28 min studyLesson 260 in Mathematics

Experimental design and optimal information

Choose informative predictor settings for a stated model, while preserving randomization, replication, and checks against misspecification.

On this page

What you’ll understand

  • Relate a design's information matrix to coefficient uncertainty.
  • Compare determinant, trace, and target-contrast criteria.
  • Explain why a mathematically optimal design can still omit essential scientific checks.

Regression analysis often begins after predictors have already been observed. Experimental design asks an earlier question: where should measurements be taken so that the intended comparisons are informative? A useful answer specifies the mean model, error structure, feasible settings, budget, and target. Optimizing a matrix criterion can then improve precision under that specification. It cannot decide which scientific assumptions are credible, guarantee a causal interpretation, or discover curvature at settings that the chosen design never visits. Those responsibilities belong in the design process too.

Information depends on the model and covariance

For a fixed full-rank linear Gaussian model with known common error variance, Fisher information for the coefficient vector is the Gram matrix divided by that variance. Even without Gaussianity, the same inverse Gram matrix controls the least-squares covariance under spherical errors. Thus selecting predictor rows changes which coefficient directions are measured precisely. Repeated observations at a setting add information under independent errors; highly correlated repeats can add much less, requiring a covariance-aware information matrix instead.

A design criterion compares feasible matrices subject to an explicit run budget and factor restrictions. Rank must be checked first: a singular information matrix leaves some parameter directions unidentified. Scaling matters when criteria compare particular coordinate uncertainties. Changing the model basis consistently leaves determinant comparisons unchanged up to a common factor, but a trace criterion can change under rescaling because it sums coordinate-specific variances. State the parameterization and target rather than calling one design intrinsically informative without context.

I(β)=1σ2XTX,Cov⁡(β^)=σ2(XTX)−1,Var⁡(cTβ^)=σ2cT(XTX)−1c\mathcal I(\beta)=\frac1{\sigma^2}X^{\mathsf T}X,\qquad \operatorname{Cov}(\widehat\beta)=\sigma^2(X^{\mathsf T}X)^{-1},\qquad \operatorname{Var}(c^{\mathsf T}\widehat\beta)=\sigma^2c^{\mathsf T}(X^{\mathsf T}X)^{-1}c
The Fisher-information identity assumes a Gaussian linear model with known σ². The covariance identities require the correct full-rank mean model and spherical covariance, but not Gaussianity.↔ Scroll wide equations horizontally.

Different objectives compress information differently

D-optimality maximizes the determinant of the information matrix, equivalently minimizing the determinant of coefficient covariance for fixed variance. In a Gaussian confidence ellipsoid, volume is proportional to the square root of that covariance determinant, with confidence level and residual scale held fixed. A-optimality minimizes the trace of inverse information, which averages coordinate variances in the chosen scaling. A c-optimal design instead minimizes uncertainty in one prespecified contrast c. These criteria need not select the same measurements.

A large determinant is a product-of-directions criterion, not a promise that every direction is individually as precise as possible. Subject to a bounded predictor region and a simple line model, endpoint measurements can be very efficient for slope estimation. If a quadratic term is scientifically plausible, endpoints alone cannot distinguish its effect from the intercept when squared predictor values are constant. A design optimized for a misspecified model can therefore be precise about an inadequate representation. Reserve opportunities to assess that representation.

D: maximize det⁡(XTX),A: minimize tr⁡ ⁣((XTX)−1),c: minimize cT(XTX)−1c\text{D: maximize }\det(X^{\mathsf T}X),\qquad \text{A: minimize }\operatorname{tr}\!\left((X^{\mathsf T}X)^{-1}\right),\qquad \text{c: minimize }c^{\mathsf T}(X^{\mathsf T}X)^{-1}c
These fixed-variance criteria compare full-rank designs for a common stated model, feasible region, and run budget. A and c criteria depend on the chosen parameter scaling or scientific target.↔ Scroll wide equations horizontally.

Randomization and model checking serve different purposes

Randomizing run order helps prevent treatment settings from being systematically aligned with uncontrolled time trends. Blocking can account for known nuisance groups, while replication provides information about variability. Center points or additional factor levels can expose curvature that endpoint-only designs conceal. These choices should be planned alongside precision criteria, not treated as cosmetic adjustments after an optimizer has consumed every run. Randomization supports a design-based argument but does not automatically make measurement errors Gaussian or remove every form of interference.

Computer searches over a finite candidate set often return a good design rather than a proved global optimum. Record the candidate set, restrictions, model terms, optimization criterion, and any optimality certificate or approximation claim. For nonlinear models or GLMs, information may depend on unknown coefficients, making a locally optimal design conditional on a working parameter value. Sensitivity analysis, sequential updating, or robust design criteria can address that dependence, but the original local optimum should not be advertised as universally optimal.

Spend four runs on an intercept and a slope

For the line model E(Y)=β0+β1xE(Y)=\beta_0+\beta_1x with independent common-variance errors and x restricted to [−1,1], compare four-run designs A: (−1,−1,1,1) and B: (−1,0,0,1).

  1. Both predictor sums are zero, so the intercept and slope columns are perpendicular. For A, XATXA=diag⁡(4,4)X_A^{\mathsf T}X_A=\operatorname{diag}(4,4); for B it is diag⁡(4,2)\operatorname{diag}(4,2). Both designs have rank 2.

  2. The determinants are 16 and 8. The slope variances are respectively σ2/4\sigma^2/4 and σ2/2\sigma^2/2, while both intercept variances equal σ2/4\sigma^2/4. A has twice B's determinant and half its slope variance.

  3. For any four allowed settings, det⁡(XTX)=4∑xi2−(∑xi)2≤16\det(X^{\mathsf T}X)=4\sum x_i^2-(\sum x_i)^2\le16. A attains 16 by putting all observations at endpoints with balanced signs. However, its squared predictor column is constant, so it cannot separately identify a quadratic term.

The result

A is D-optimal for this fixed four-run line model on [−1,1]. B sacrifices line-model precision to include center settings that can help investigate curvature when replication and a richer model are considered.

Explore the model

Interactive lab

Interactive lab · matrix

Compare information volume with directional imbalance

Pause and Reset; set a=1 and b=0. Predict determinant 1 and equal singular values. For a separate run Reset and set a=1 and b=2. The determinant stays 1 while singular values become unequal. If this matrix is viewed only as a hypothetical design factor, its Gram determinant also stays 1, yet directional precision changes. The instrument enforces no run budget or candidate set and does not optimize an experimental design.

Compare information volume with directional imbalanceThe selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA determinant summarizes a product of information scales; scientific targets and model checks still determine which design is useful.

Check your understanding

Why can an endpoint-only design optimized for a straight line be inadequate when curvature is plausible?

Try it yourself

For a three-run intercept-and-slope model, compare designs C: (−1,0,1) and D: (−1,−1,1). Compute each Gram determinant and slope variance under independent errors with common variance σ².

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Use the two-by-two Gram matrix with entries n, sum x, and sum x². The slope variance is the lower-right entry of its inverse times σ².

Compare with the solution

For C the Gram matrix is diag⁡(3,2)\operatorname{diag}(3,2), with determinant 6 and slope variance σ2/2\sigma^2/2. For D it is (3−1−13)\begin{pmatrix}3&-1\\-1&3\end{pmatrix}, with determinant 8 and inverse 18(3113)\frac18\begin{pmatrix}3&1\\1&3\end{pmatrix}; its slope variance is 3σ2/83\sigma^2/8. D gives greater determinant and smaller slope variance for the line model, while C includes a center setting useful for a different model-checking objective.

Take it with you

  • Information criteria compare designs under a stated model, covariance, and budget.
  • Determinant, average-variance, and contrast-specific objectives can disagree.
  • Randomization, replication, and curvature checks supply protections that a matrix optimum alone cannot provide.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.