UpSkill Sprint Consulting logo UpSkill Sprint Consulting
Browse lessons

Advanced statistics · Minitab workflow

Finding the Best Predictive Regression Model in Minitab

Build a strong candidate term pool, let Minitab search it without violating hierarchy, validate the search on unseen data, compare interpretable nonlinear alternatives, and select a model an engineer can defend.

Advanced 30–40 minutes Interactive Continuous response

Learning objectives

What you will be able to do

By the end of this lesson, you will be able to configure Minitab to search a controlled set of linear, interaction, and polynomial terms—and distinguish a model that merely fits the current data from one that is likely to predict new production data.

01

Construct

Build a candidate model with main effects, squares, and defensible interactions without creating an uncontrolled term explosion.

02

Configure

Use centering, hierarchical forward selection, and cross-validation correctly in Minitab.

03

Compare

Challenge polynomial regression with MARS and automated machine-learning models while preserving interpretability.

04

Diagnose

Review residual patterns, VIF, leverage, Cook’s distance, time order, and extrapolation risk.

05

Select

Choose the simplest model whose validation performance is practically equivalent to the strongest competitor.

06

Defend

Document the model search, validation design, limitations, and engineering rationale for technical review.

The real problem

You are not searching for the equation with the largest R²

A model can fit historical data extremely well and still fail on the next heat, coil, batch, customer, or production campaign. The objective is to find a model that balances predictive accuracy, stability, interpretability, and process credibility.

Y = β₀ + β₁X₁ + β₂X₂ + β₃X₁² + β₄X₁X₂ + ε

Track 1 — Interpretable polynomial regression

Search a controlled pool of main effects, square terms, and two-way interactions. This produces a conventional equation, coefficient table, VIF values, residual plots, and familiar engineering interpretations.

Track 2 — Interpretable nonlinear regression

Use MARS when the process contains thresholds, slope changes, or operating regions that a single quadratic curve cannot represent well.

Critical correction: in Minitab’s forward selection with validation, the next term is selected primarily by the stepwise significance rule; validation is then used to identify the best model size or step. It is not a brute-force search that directly tests every possible equation and adds whichever term gives the largest validation improvement.
Do not use the ordinary training R² as the selection criterion. It nearly always increases as terms are added—even when the added terms reduce generalization.

Step 1

Build a large enough—but still defensible—candidate pool

Start with process physics and measurement logic. Give Minitab plausible terms to consider; do not ask it to manufacture meaning from every mathematical combination.

Candidate-term builder

See how quickly model complexity grows as you add squares, interactions, and individual cubic terms.

Try it

14 terms

candidate terms, excluding the intercept

Recommended first-pass candidate model

Main effects + individual squares + two-way interactions

This is exactly what Minitab’s Terms through order: 2 option creates for selected continuous predictors: powers through order 2 and interactions through order 2. If categorical predictors are selected, Minitab also creates the applicable interactions with those predictors.

Why not start with “Terms through order 3”? It adds more than individual cubes. It can also generate mixed terms such as X₁²X₂, X₁X₂², and X₁X₂X₃. The search space, multicollinearity, and interpretation burden can grow sharply.

Use individual cubic terms only when justified

Add Xᵢ³ selectively when process knowledge or residual curvature suggests asymmetric behavior that a quadratic cannot represent. For more complicated nonlinear structure, MARS is usually a better challenger than a high-order global polynomial.

Hierarchy challenge

Select every lower-order term required to support the chosen higher-order term.

Practice
Choose the necessary lower-order terms, then check your answer.

Step 2

Design validation before model selection

Validation is part of the model-building design—not a final decorative statistic. Decide what constitutes an independent future prediction before opening the regression dialog.

Use random K-fold cross-validation only when rows are reasonably independent and exchangeable. For smaller datasets, 5-fold or 10-fold cross-validation is a practical default. For larger datasets, a separate test set may be faster and easier to audit.

Prevent group leakage

If multiple rows come from the same heat, coil, pipe, batch, customer, production campaign, or test specimen, keep that entire group in one fold. Otherwise, the validation set contains near-duplicates of information already seen during training.

Random rows vs grouped folds

Toggle the split method and observe whether the same production group appears in both training and validation folds.

Visualize
Important implementation detail: prepare a fold column such as CV_Fold containing fold labels 1–5 or 1–10. Assign every row from the same heat or coil to the same fold label. Do not use a raw identifier with hundreds of unique values unless those values are intentionally the folds you want Minitab to use.

Recommended worksheet structure

ColumnExampleRole
ResponseYield StrengthContinuous outcome to predict
Continuous predictorsTemperature, CE, GaugeNumeric inputs
Categorical predictorsGrade, Mill, ShiftFactor levels coded by Minitab
Group IDHeat_IDDefines the unit that must stay together
CV fold1, 2, 3, 4, 5Used by Minitab’s “Assign rows of each fold by ID column” option
Time orderProduction timestampUsed to detect drift, cycles, and serial dependence

Step 3

Exact Minitab workflow: polynomial regression with validated selection

The sequence below is designed for a continuous response and a mix of continuous and categorical predictors.

1Define the model purpose and data boundaryPrediction target, operating range, grouping unit, and holdout logic

State what the model will predict, for which population, and within what process range. Decide whether the validation target is a new row, a new heat, a new coil, a future campaign, or a later time period.

  • Exclude post-outcome variables and data created after the prediction decision.
  • Flag impossible values, missingness, duplicate records, and data-quality rules before fitting.
  • Reserve a true future or external test set when the decision is high consequence.
2Open Fit Regression ModelAssign the response and predictor types correctly
  • Enter the continuous outcome in Response.
  • Enter numeric process inputs in Continuous predictors.
  • Enter grade, mill, product family, shift, operator, or other factors in Categorical predictors.
  • Verify the categorical reference level under Coding so coefficients compare against a meaningful baseline.
3Create the candidate term poolMain effects, squares, and two-way interactions
  1. Select the continuous predictors you want to expand.
  2. Choose Terms through order.
  3. Enter 2.
  4. Click Add.
  5. Review the term list and remove combinations that are physically impossible, redundant, or unsupported.

Terms through order 2 adds individual squares for continuous predictors and all interactions through order 2 among the selected predictors.

4Center continuous predictorsReduce nonessential collinearity from powers and interactions

Centering changes each predictor from X to X − X̄. This reduces the artificial correlation between X, X², and interaction terms and makes lower-order coefficients meaningful near the center of the observed operating region.

Xcentered = X − X̄

Do not confuse centering with dividing by the standard deviation. Scaling can be useful for coefficient comparison, but it changes the unit interpretation.

5Configure forward selection with validationLet validation choose the model size
  • Method: Forward selection with validation
  • Validation: K-fold cross-validation for smaller data; test-set validation for larger data
  • Folds: 5 or 10 are practical starting points; Minitab’s regression default is 10
  • Group-aware folds: Assign rows by a prepared fold-ID column
  • Display: Include model-selection details
  • Graph: Display R² versus model-selection step
How it works: within each fold, Minitab performs forward selection. It then evaluates the sequence across folds and identifies the step with the best K-fold stepwise R². Finally, it runs forward selection on the full dataset and stops at that selected step.
6Enforce hierarchyKeep higher-order terms statistically coherent
  • Hierarchical model: Require a hierarchical model at each step
  • Require hierarchy for: All terms
  • Terms entering: Extra terms can enter to maintain hierarchy

For X₁X₂, retain X₁ and X₂. For X₁³, retain X₁ and X₁². For X₁²X₂, retain the applicable lower-order components required by the hierarchy rule.

7Request diagnostics and stored statisticsDo not stop at the selected equation
  • Display the four-in-one residual plots.
  • Display VIF values and unusual observations.
  • Store fits, residuals, standardized residuals, leverage, Cook’s distance, and DFITS when deeper investigation is needed.
  • Plot residuals against production order or timestamp when the process is time-dependent.
8Interpret the validation curve—not only the peakPrefer a stable plateau over a fragile maximum

If the highest validation R² occurs at step 12 but steps 7–10 are nearly identical, prefer the simpler model unless the additional terms provide a material, repeatable improvement.

  • Compare training R² with validation R².
  • Look for a widening training–validation gap.
  • Check whether selected terms are stable under alternate fold assignments.
  • Confirm coefficient directions and interactions make engineering sense.
Stepwise selection is exploratory. It cannot incorporate every practical consideration, can choose one of several correlated predictors arbitrarily, and can produce different “best” models under different criteria. Treat selected-term p-values cautiously and validate with new data whenever possible.

Step 4

Challenge the polynomial with MARS regression

MARS—Multivariate Adaptive Regression Splines—fits connected piecewise functions. It is especially useful when the process has thresholds, breakpoints, or different slopes in different operating regions.

Ŷ = β₀ + β₁·max(0, X − c) + β₂·max(0, c − X)

A “knot” at c allows the slope to change. This can be more interpretable than forcing one cubic curve across the full process range.

526.3 MPa

fictional predicted response

Minitab path and starting settings

Validation

Use K-fold cross-validation or a separate test set. Use a fold-ID column when groups must stay together.

Interactions

Start with an additive model or allow interactions up to order 2. Higher orders rapidly reduce interpretability.

Selection criterion

Compare maximum validation R² with minimum validation MAD. MAD is less influenced by the worst-fitting observations.

Interpretation trade-off: MARS can produce an equation-like representation and variable-importance output. However, partial dependence plots are not available in Minitab when the MARS analysis includes interactions.

Step 5

Use Discover Best Model as a challenger—not an automatic final decision

With Minitab’s Predictive Analytics Module, Discover Best Model for a continuous response compares five model families: multiple regression, MARS, CART, TreeNet, and Random Forests.

Model-selection decision lab

Select a model, then apply the “simplest model within practical equivalence” rule.

Decide
ModelValidation R²Validation RMSEMADInterpretabilityComplexity
Polynomial regression84.6%11.28.2High8 terms
MARS88.9%9.46.8Good11 basis functions
CART82.1%12.18.9Medium14 terminal nodes
TreeNet89.3%9.16.6Limited300 trees
Random Forests88.5%9.66.9Limited500 trees
Select a model row or apply the rule to see the recommended decision.
Recommended principle: choose the simplest interpretable model whose validation performance is practically equivalent to the strongest model. “Practically equivalent” should be defined before seeing the final ranking using engineering error tolerance—not improvised afterward.

Step 6

The selected model is only a candidate until it passes diagnostics

Validation performance answers “how accurately did this workflow predict held-out observations?” Diagnostics answer “is the model structurally credible, stable, and safe to use?”

Residuals versus fits

Look for random scatter around zero. Curvature suggests missing terms; fanning suggests nonconstant variance.

Residuals versus order

Trends, cycles, or shifts indicate dependence, process drift, omitted time effects, or a changing measurement system.

Normal probability plot

Strong curvature or distant points can make p-values and confidence intervals unreliable, even when point predictions remain useful.

Leverage and Cook’s distance

Investigate observations with unusual predictor combinations or disproportionate influence. Never delete a valid point only because it weakens the model.

Required final review

Prediction

Validation R², RMSE, MAD, and the gap between training and validation performance.

Complexity

Number of retained terms or basis functions; stability across alternate folds.

Collinearity

VIF values, coefficient instability, and redundant predictors.

Residual adequacy

Linearity, constant variance, normality when inference matters, and independence.

Influence

Leverage, Cook’s distance, DFITS, data quality, and sensitivity analyses.

Engineering validity

Logical signs, plausible thresholds, operating-range coverage, and no unsafe extrapolation.

Step 7

A defensible model-selection rule

Do not allow a single metric to decide. Use sequential gates so that a model cannot “win” on accuracy while failing interpretability, stability, or process logic.

Gate 1: Valid designNo leakage, correct grouping, correct prediction horizon, and no post-outcome predictors.
Gate 2: Competitive errorValidation error meets the practical tolerance and is close to the best challenger.
Gate 3: Adequate diagnosticsNo unresolved structural residual pattern, severe instability, or dominant data defect.
Gate 4: Engineering utilityInterpretation, measurement burden, deployment feasibility, and safe operating range are acceptable.
Recommended final decision statement: “MARS was selected because its grouped five-fold validation RMSE was within 3% of the best model, while providing substantially greater interpretability and stable predictor effects. The model passed residual, influence, and operating-range checks. TreeNet was retained as a benchmark but not selected for deployment.”

When to reject the automated winner

  • The improvement is smaller than normal measurement or process variation.
  • The model relies on unavailable, delayed, expensive, or unreliable predictors.
  • The validation split leaked information across the same production units.
  • The model extrapolates beyond the observed operating region.
  • The predictor effects contradict known physics without a plausible explanation.
  • The model is unstable when folds, time windows, or influential points change.

Knowledge check

Can you run the workflow defensibly?

8-question applied quiz

Choose the strongest answer. Feedback is shown after submission.

1. What does “Terms through order 2” create?
2. How does forward selection with validation work?
3. Why center continuous predictors?
4. Several rows come from each heat. What is the strongest validation design?
5. If X₁³ is retained, which lower-order terms should ordinarily remain?
6. Which model is especially useful for interpretable thresholds?
7. TreeNet has 0.4 percentage points more validation R² than MARS. What should you do?
8. Residuals show a strong trend over production order. What does this imply?

Field checklist

Your model-selection standard work

Use this checklist each time you search for a predictive regression model in Minitab.

Define the independent prediction unit.New row, heat, coil, customer, campaign, or future time period.
Audit data before modeling.Missing values, impossible values, duplicate records, leakage, and measurement integrity.
Start with a controlled second-order pool.Main effects, individual squares, and defensible two-way interactions.
Center continuous predictors.Use Coding → Subtract the mean when powers and interactions are present.
Use forward selection with validation.Require hierarchy and display selection details plus the validation R² curve.
Protect grouped observations.Prepare a fold-ID column that keeps each group entirely within one fold.
Challenge the polynomial.Compare MARS and, where available, Discover Best Model.
Review diagnostics.Residuals, VIF, leverage, Cook’s distance, time order, and extrapolation.
Select by practical equivalence.Prefer the simpler interpretable model when the predictive difference is immaterial.
Document limitations and monitor.State the operating range, expected error, retraining trigger, and model owner.

Verified sources

Official Minitab documentation

  1. Specify model terms for Fit Regression Model
  2. Specify coding and standardization
  3. Perform stepwise regression
  4. Specify the validation method
  5. Using stepwise and best subsets regression
  6. Interpret residual plots
  7. Overview of MARS Regression
  8. MARS validation methods
  9. Discover Best Model for a continuous response