Experimental design and optimal information
Choose informative predictor settings for a stated model, while preserving randomization, replication, and checks against misspecification.
On this page
The question for today
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.
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.
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.
Work it through
Spend four runs on an intercept and a slope
For the line model 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).
Both predictor sums are zero, so the intercept and slope columns are perpendicular. For A, ; for B it is . Both designs have rank 2.
The determinants are 16 and 8. The slope variances are respectively and , while both intercept variances equal . A has twice B's determinant and half its slope variance.
For any four allowed settings, . 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.
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.
The living notebook
Explore the model
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.
- Area scale · det A
- 1.4
- Singular values · σmax, σmin
- 1.597, 0.876
- Condition number · κ₂(A)
- 1.823
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.
Pause & reason
Check your understanding
Your turn
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 , with determinant 6 and slope variance . For D it is , with determinant 8 and inverse ; its slope variance is . 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.
Keep the idea
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.