Back to List

From Factor Effects to Prediction and Exploration: Diagnosing a DOE Model

Connects full-factorial responses to effects and ANOVA, hierarchical regression models, residual diagnostics, natural-unit equations, and Profiler-based condition exploration.

Advanced
|
32min
|
Verified (2026-08-14)
model hierarchypoolingcoded unitnatural unitPrediction ProfilerdesirabilityresidualJMP
Progress0/28 (0%)

Filling a design table with response values is not the end. Only after calculating effects, expressing a model, and checking whether residuals support that candidate model can you explore conditions.

This unit asks one question.


How do we turn design-table responses into a diagnosable model and candidate next conditions?
model hierarchyIf there is an interaction, the related main effects are kept together.
coded unitDesign coordinates changed to −1 and +1
natural unitReal units such as °C, minutes, and concentration
ProfilerExpression for exploring model predictions according to input changes

Effects and regression coefficients express the same information on different scales

For a 2² design, the coded model is ŷ=b₀+bAA+bBB+bABAB. An effect is high mean minus low mean, while a ±1 coded coefficient is effect/2. In a balanced design, b₀ is the grand mean.

U19 · Figure 01
Effects become model terms, and the model makes condition-specific predictions.
A80B50AB60
When reducing terms, don't just chase the p-value, but maintain a hierarchy of interactions and relevant main effects.

ANOVA partitions SS explained by the model and residual SS. With replication, residuals can be divided further into pure error and lack of fit. Before scanning each effect p-value, inspect the complete candidate model, degrees of freedom, and residual information.

Model hierarchy preserves an interpretable coordinate system

When including AB interaction, retain A and B main effects; when including A², retain A. This is the hierarchy principle. Removing lower-order terms because their p-values are large changes the model's meaning when the coding origin changes and breaks interaction interpretation.

In a small DOE, fitting every term leaves almost no residual degrees of freedom. Effect plots, a prior sparsity assumption, and scientific knowledge can help, but repeatedly deleting terms through automatic stepwise p-value rules creates selection bias. A pooled term is not declared exactly zero; it is a modelling judgment to combine it into error estimation.

An attractive effect plot does not exempt residual diagnostics

Time drift, missing curvature, unequal variance, and influential points can distort effect size. Inspect actual by predicted, residual by predicted, residual by run order, and raw data together.

Coded and natural-unit equations have different roles

Coded −1/+1 units make effect comparisons and orthogonal structure clear. A natural-unit equation in real °C, pH, or rpm is needed for operational interpretation and communication.

For low X_L and high X_H, coded x is (X−center)/half-range. Transforming to natural units changes coefficient numbers and units, but must give the same predictions within the same range. Preserve the transformation, units, and valid range together.

The Profiler moves a cross-section of the fitted model

Prediction Profiler shows a prediction while moving one factor and holding the others at their current values. With interaction, the slope for A changes with B's setting, so panels are not independent. A contour shows a two-factor combination at once.

Desirability turns “maximize,” “target,” or “range” into 0–1 functions and can combine responses. Maximum desirability is optimal only for the supplied goals; it does not automatically define biological importance or manufacturing robustness. Changing function shapes or weights changes the candidate.

In-Silico Lab: turn effects into an equation and move coordinates

  1. Compare B-dependence of the A profiler at AB=0 and AB=6.
  2. Move A and B from −1 to +1 and read the current prediction.
  3. Explain why values outside the design range are unavailable.
  4. List range information needed to convert coded coefficients to natural units.
In-Silico Lab · U19

Connect effects with model formulas and conditional exploration

Move A·B at ±1 coded coordinates to read the predictions of the interaction model and see items to check before natural unit interpretation.

If this is your first time: What should I press?
  1. 1. Read the question firstIn the Lab title, check the one thing you will compare this time.
  2. 2. Change just one conditionInitially, change only one of the inputs: n, effect, or spread.
  3. 3. new composite specimen pressureNew synthetic data is created. The same conditions may vary depending on the sample.
  4. 4. Pictures and calculation results CompareWrite in one sentence what moves and what stays the same before and after the change.

If it gets stuckresetGo back to see the default results and change just one condition. This Lab is not a correct answer tester but a pattern observation tool.

Synthetic observations of the same settings

A coded, predicted cross section with B fixed at the current value.

calculation result

intercept70.0
A/B/AB coefficients8 / 5 / 6
current forecast70.0
design scope−1 ≤ A,B ≤ +1

The maximum desirability point is the maximum number of goals entered by the user. It does not automatically guarantee scientific significance and confirmatory testing.

educational synthetic modelbjs-factorial-sequence-v1. Actual research judgments must separately reflect experimental units, missingness, distribution, multiplicity, pre-planning, and domain criteria.

Read JMP analysis in judgment order, not output order

Effect / ANOVA

See which terms explain the response variation, but maintain the hierarchy.

Residual Diagnostics

Look for model failures behind good-looking effects charts.

Profiler / Contour

It explores natural unit conditions and goals, but does not replace confirmatory experiments.

Recommended reasoning order:

  1. verify actual runs, missingness, and deviations;
  2. inspect effect directions and interaction plots;
  3. define candidate model and hierarchy;
  4. review ANOVA and effect uncertainty;
  5. diagnose residuals, influence, and run order;
  6. reconcile coded and natural equations;
  7. explore candidates with Profiler and contour;
  8. record confirmation experiments and prediction error.

An optimum is a candidate for the next experiment

The model is estimated from finite design points and error. At a predicted maximum, perform an independent confirmation run and compare observed values with the prediction and CI or PI. A boundary optimum may suggest a wider range is useful, but never extend automatically beyond safety or physical limits.

Example result statement

We hierarchically fitted the pre-specified candidate model A+B+AB. Its coded equation was 70+8A+5B+6AB, with no time trend in residuals by run order. A=+1 and B=+1 maximized prediction in the Profiler, but this was a design boundary, not a declared final optimum. We planned an independent confirmation batch to check prediction error.

Takeaways

  • A factorial effect and a ±1 coded coefficient differ by a factor of two.
  • Retain related lower-order terms with interactions and squared terms.
  • Pooling and automatic term deletion require modelling judgment.
  • Residuals and run order must not trail behind the effect table.
  • Attach units and a valid range to natural-unit equations.
  • A Profiler optimum is a candidate requiring confirmation.
The effect table is not the end. Conditional exploration becomes the basis for a hierarchical model, residual diagnosis, natural unit analysis, and confirmation experiments.

The next unit is U20. It examines which effects become aliased and indistinguishable when runs from a full design are reduced, then stops automatically.

Official supplementary resources

The model and Lab in this article are educational synthetic material.

💬 Questions & Comments

0 comments

You can post without signing in. Guest comments cannot be edited or deleted by their author.

0/2000

Loading...