Fact-checked by Grok 2 weeks ago

Central composite design

Central composite design (CCD) is a statistical experimental design used in (RSM) to efficiently fit second-order () models to experimental data, enabling the estimation of and interactions in response surfaces without the need for a complete three-level . It combines a two-level or with additional axial (or star) points positioned at a specified distance from the center and multiple center points replicated to assess pure error and improve model stability. Introduced by and Kenneth B. Wilson in their seminal 1951 paper, CCD originated as a foundational in RSM for attaining optimal conditions in experimental processes, building on earlier steepest ascent techniques to systematically explore and optimize response variables. The design always includes twice as many star points as the number of factors (k), with the axial distance parameter α typically chosen to ensure rotatability—providing uniform prediction variance across the design space—or for unbiased estimates. CCDs are categorized into three primary types based on the positioning of axial points relative to constraints: the circumscribed () design, which extends beyond the for broader exploration (requiring five levels per ); the inscribed (CCI) design, a scaled version of fitted within predefined limits; and the face-centered () design, where axial points lie on the faces of the (using only three levels per ). This flexibility allows CCDs to accommodate sequential experimentation, starting from screening designs and augmenting them as needed. In practice, CCDs are valued for their efficiency in estimating first- and second-order model coefficients with fewer runs than full factorials, supporting orthogonal blocking to control for nuisance factors and facilitating the identification of optimal process settings. They are extensively applied in fields such as for reaction optimization, for process improvement (e.g., adjusting and in injection molding), pharmaceuticals for development, and other areas involving multivariable optimization where curvature in responses is anticipated.

Background

Definition and Purpose

The central composite design (CCD) is a statistical experimental design used in (RSM) to fit second-order models that capture curvature in the response variable. It augments a two-level full or with additional center points and axial (also known as star) points, allowing for the estimation of quadratic terms and interactions among factors. This structure enables efficient modeling of nonlinear relationships between input variables and the output response. The primary purpose of the is to facilitate the exploration and optimization of response surfaces in experimental settings where linear approximations are insufficient, such as in and product development. By incorporating points that extend beyond the , the design detects and quantifies effects that might otherwise go unnoticed, supporting sequential experimentation to refine models iteratively. This approach is particularly valuable in RSM for identifying optimal operating conditions with minimal resource expenditure. Compared to full three-level designs, which require 3^k runs for k factors and can become prohibitively large, the offers key advantages in efficiency, typically needing only 2^k factorial points plus 2k axial points and several center points to estimate pure experimental error and improve model robustness. For example, for three factors, a requires 20 runs compared to 27 for a full three-level , while retaining the ability to assess interactions and nonlinearity, makes the a practical choice for applied research aiming to balance information gain with cost.

Historical Context

The foundations of central composite design (CCD) trace back to the development of factorial designs in the early 20th century, pioneered by Ronald A. Fisher and Frank Yates at the Rothamsted Experimental Station in the 1920s and 1930s. These designs enabled the efficient study of multiple factors and their interactions in agricultural experiments, laying the groundwork for more advanced sequential optimization techniques in experimental design. Central composite design was invented by and Kenneth B. Wilson in as a key component of (RSM), aimed at modeling and optimizing processes through second-order approximations. In their seminal paper, they introduced the CCD alongside the method of steepest ascent to explore response surfaces systematically, building on factorial designs to include axial and center points for curvature estimation. This innovation marked a shift toward efficient designs for industrial optimization, particularly in chemical processes. Subsequent refinements focused on enhancing the statistical properties of CCDs, notably by Box and J. Stuart Hunter in 1957, who developed criteria for rotatability to ensure uniform prediction variance across the design space. Their work on multifactor experimental designs for response surfaces standardized the axial distance parameter (alpha) to achieve this property, making CCDs more practical for second-order modeling. The adoption of has profoundly influenced modern (DOE) practices, integrating into software tools such as and Design-Expert for automated generation and analysis of rotatable designs. In industries like , has become a standard for process optimization and , as evidenced by its widespread use in RSM applications since the mid-20th century.

Design Construction

Core Components

The central composite design (CCD) is constructed around a core structure consisting of a factorial portion and center points, which together provide the foundational data for estimating linear and interaction effects in . The factorial portion comprises $2^k points, where k is the number of factors, positioned at the levels \pm 1 in a coded scale for each factor; this embedded two-level full or allows for the estimation of main effects and two-factor interactions through orthogonal contrasts. Center points, replicated multiple times at the (0, 0, \dots, 0) in the coded units, serve to estimate the pure experimental error and detect any in the response surface by comparing variability against the factorial points. Typically, 5 to 6 center point replicates are included to ensure sufficient for error estimation and model stability, balancing efficiency with reliability in second-order modeling. These core components support the fitting of a second-order model of the form y = \beta_0 + \sum_{i=1}^k \beta_i x_i + \sum_{i=1}^k \beta_{ii} x_i^2 + \sum_{i<j} \beta_{ij} x_i x_j + \epsilon, where the factorial points primarily inform the linear terms \beta_i and interaction terms \beta_{ij}, while center points contribute to the intercept \beta_0 and help assess quadratic curvature via \beta_{ii}. The total number of runs in this core structure is N = 2^k + n_c, with n_c denoting the number of center points; this is often augmented with axial points to fully resolve the quadratic terms.

Axial Points and Alpha Selection

In central composite designs (CCD), axial points, also known as star points, are additional experimental runs positioned along the axes of the factor space to facilitate the estimation of quadratic effects in the response surface model. These points are located at coordinates (±α, 0, ..., 0) for each factor, with permutations across all k factors, resulting in 2k axial points. By placing these points beyond the factorial portion of the design, they enable the fitting of second-order terms without relying solely on the linear factorial points, which alone cannot capture curvature effectively. This structure was introduced as part of the foundational CCD framework by to explore response surfaces systematically. The parameter α represents the distance from the design center to each axial point, scaled in coded units where the factorial points lie at ±1. To prevent overlap or confounding between linear and quadratic estimates, α is typically selected greater than 1, ensuring the axial points extend outside the hypercube defined by the . When α = 1, the points lie on the faces of the cube, which may introduce some confounding but is suitable for constrained experimental regions. The choice of α critically influences the geometric properties and statistical efficiency of the , as detailed in subsequent developments by . Several methods exist for selecting α, each tailored to specific design objectives such as feasibility, variance properties, or blocking. In face-centered CCD (CCF), α = 1 confines all points within the [-1, 1] bounds, requiring only three levels per factor and simplifying implementation in cuboidal regions, though it sacrifices some precision for higher-order terms. For rotatable CCD (CCC), α is chosen to achieve uniform prediction variance at points equidistant from the center, given by the equation \alpha = (2^k)^{1/4} = (n_f)^{1/4}, where k is the number of factors and n_f = 2^k is the number of factorial points in a full factorial design; this value ensures the design sphere has constant variance, ideal for spherical experimental regions. For orthogonal CCDs, particularly in blocked experiments, α is selected to ensure uncorrelated estimates, often approximately \sqrt{2^k} for standard two-block designs with balanced center points, but generally computed based on the specific blocking structure and number of center points to orthogonalize blocks to model terms. These selections stem from optimality criteria outlined by Box and Hunter, balancing estimability and efficiency. The value of α directly impacts the resolution and aliasing characteristics of the CCD. Larger α values, as in rotatable or orthogonal designs, improve the resolution for quadratic terms by better isolating them from higher-order interactions, particularly when augmenting a resolution V fractional factorial base, which minimizes aliasing of main effects and two-factor interactions with quadratic terms. Conversely, smaller α (e.g., 1 in face-centered designs) may lead to increased aliasing with cubic or higher effects, reducing the design's ability to distinguish subtle curvature but allowing feasible experimentation within bounded regions. Optimal α selection thus depends on the anticipated response behavior and experimental constraints, as analyzed in standard DOE references.

Full Design Matrix

The full design matrix for a central composite design (CCD) is constructed by integrating the factorial points, axial points, and center points into a single matrix \mathbf{X}, where all coordinates are in coded units scaled to the interval [-1, 1] for factorial and center points, and to \pm \alpha along individual axes for axial points. This assembly allows for efficient estimation of second-order response surface models by providing points that capture linear, interaction, and quadratic effects. The total number of runs is $2^k + 2k + n_c, where k is the number of factors, $2^k are the factorial points (or a fraction thereof), $2k are the axial points, and n_c (typically 5–6 or more) are the replicated center points used to estimate pure error and assess lack of fit. Coded units standardize the design space, transforming actual (uncoded) variable levels via the formula
x_i = \frac{Z_i - \bar{Z}_i}{0.5 (Z_{i,\max} - Z_{i,\min})},
where Z_i is the actual value of factor i, \bar{Z}_i is its midpoint (average of low and high levels), and the denominator represents half the range. This scaling ensures factorial points lie at \pm 1 and facilitates model interpretation, with uncoded values recoverable by reversing the transformation. For instance, if temperature ranges from 50°C to 100°C, the midpoint is 75°C and half-range is 25°C, so a coded value of 1 corresponds to 100°C. A concrete example for k=2 factors (e.g., temperature and pressure) uses a full $2^2 factorial with \alpha = \sqrt{2} \approx 1.414 for rotatability, yielding 4 factorial points, 4 axial points, and typically 1–6 center points (here shown with 1 for brevity, but replications are recommended). The matrix is:
Runx_1x_2
1-1-1
21-1
3-11
411
500
6-1.4140
71.4140
80-1.414
901.414
Additional center points can be appended as repeated rows of (0, 0) to improve precision. For k=3 factors (e.g., temperature, pressure, and catalyst concentration), a full factorial CCD with \alpha \approx 1.682 (for rotatability) and 6 center points totals 20 runs: 8 factorial, 6 axial, and 6 centers. The matrix, randomized in practice but shown here systematically, is:
Runx_1x_2x_3Type
1-1-1-1Factorial
21-1-1Factorial
3-11-1Factorial
411-1Factorial
5-1-11Factorial
61-11Factorial
7-111Factorial
8111Factorial
91.68200Axial
10-1.68200Axial
1101.6820Axial
120-1.6820Axial
13001.682Axial
1400-1.682Axial
15–20000Center (×6)
Note that while \alpha = 2 may be selected in some non-rotatable or orthogonal variants for simplicity (e.g., when prioritizing uniform precision over rotatability), the value here aligns with standard rotatable properties. Replicated designs emphasize multiple center points to quantify experimental error, while blocked designs separate factorial and axial points into distinct blocks to control for time drifts or nuisance factors; for example, in a two-block setup, factorial points occupy one block and axial points the other, with centers distributed across blocks for balance. This approach maintains estimability of model parameters without confounding.

Statistical Properties

Rotatability and Uniform Precision

Rotatability in experimental design refers to the property where the prediction variance of the fitted response surface model remains constant for all points equidistant from the design center. This symmetry ensures that rotating the coordinate system around the center does not alter the variance, making it particularly desirable for response surface methodology (RSM) as it provides isotropic prediction accuracy across the design space. Introduced by in their foundational work on second-order designs, rotatability facilitates reliable estimation of quadratic effects without directional bias in variance. In central composite designs (CCD), rotatability is achieved by selecting the axial distance parameter \alpha such that \alpha = (2^k)^{1/4}, where k is the number of factors; this value, detailed further in the context of axial point placement, positions the star points to create spherical variance contours. For instance, with k=2, \alpha \approx 1.414, and for k=3, \alpha \approx 1.682, ensuring the design's prediction variance depends solely on the radial distance from the origin in the coded factor space. This configuration applies to rotatable variants like the circumscribed (CCC) and inscribed (CCI) CCDs, while face-centered (CCF) designs lack this property due to constrained \alpha = 1. Uniform precision builds on rotatability by emphasizing minimized and consistent variance throughout the coded design region, often enhanced by including multiple center points to stabilize estimates near the origin. In rotatable CCDs, this results in a design where prediction reliability is balanced between the factorial hypercube and the axial extensions, avoiding inflated variances at the periphery that could occur in non-rotatable setups. The prediction variance for a point x in the design space is given by \text{Var}(\hat{y}) = \sigma^2 \left[ x' (X'X)^{-1} x + 1 \right], where \sigma^2 is the error variance, X is the design matrix, and the term x' (X'X)^{-1} x captures the scaled distance effect. For rotatable designs, this simplifies such that \text{Var}(\hat{y}) is a function only of the Euclidean distance r = \|x\| from the center, yielding spherical iso-variance contours rather than ellipsoidal or irregular shapes. Variance plots illustrate the superiority of rotatable CCDs: in a two-factor rotatable design, contours form perfect spheres centered at the origin, maintaining low and uniform variance (e.g., scaling linearly with r^2) up to the star points, whereas non-rotatable versions like CCF exhibit directional stretching and higher peripheral variance, potentially leading to biased optimizations. These visualizations, as shown in comparative figures for CCC versus CCF, highlight how rotatability enhances the design's efficiency for RSM applications.

Orthogonality Conditions

In central composite designs (CCDs), orthogonality refers to the property where the columns of the design matrix X are orthogonal, resulting in X'X being a diagonal matrix. This condition ensures that the least-squares estimates of the regression coefficients \beta in the second-order response surface model are unbiased and uncorrelated with each other. For a CCD to satisfy orthogonality, particularly in blocked experiments where the factorial points form one block and the axial points another, the axial distance parameter \alpha must be selected as \alpha = \sqrt{2^k}, where k is the number of factors and the factorial portion consists of a full $2^k design. This choice prevents confounding between the linear, quadratic, and interaction terms by ensuring that block effects are orthogonal to the model parameters. In the resulting moment matrix X'X, the off-diagonal elements are zero under this condition, while the diagonal elements for the quadratic coefficients \beta_{ii} take the form $1 + \alpha^2 in normalized units (assuming no center points for simplicity and scaling the factorial contributions to unity). The primary benefits of orthogonality in CCDs include simplified computation of least-squares estimates, as the inverse of the diagonal X'X matrix is straightforward, and the absence of bias in directions such as steepest ascent, where linear term estimates remain unaffected by quadratic or interaction effects. However, selecting \alpha = \sqrt{2^k} for orthogonality often conflicts with achieving rotatability, which requires a different value of \alpha (typically (2^k)^{1/4}) to ensure uniform prediction variance; thus, experimenters must prioritize based on whether parameter independence or prediction uniformity is more critical.

Variance Properties

The variance properties of central composite designs (CCDs) are evaluated through metrics that assess prediction accuracy and estimation efficiency for second-order models. A key measure of efficiency is the average prediction variance at the design points, which for many designs approximates \sigma^2 \frac{p}{N}, where p is the number of model parameters and N the number of runs; lower values indicate more efficient designs for parameter estimation and prediction. This metric highlights how CCDs balance run economy with variance minimization, particularly when compared to less efficient alternatives. In terms of overall efficiency, CCDs require significantly fewer experimental runs than a full three-level factorial design; for instance, a three-factor CCD typically uses 20 runs versus 27 for the full three-level factorial, and for five factors, ≈48 versus 243. This reduction stems from augmenting a two-level factorial with axial and center points, yielding comparable model precision with reduced resource demands. The power of CCDs to detect quadratic effects relies on the axial distance \alpha and the number of center points n_c, with the minimum detectable curvature determined by the standard error of quadratic coefficients β_{ii}, which is \sigma / (\alpha \sqrt{2}); larger \alpha values enhance detection sensitivity for curvature. Larger n_c further improves power by stabilizing error estimates, allowing reliable identification of nonlinear responses at significance levels like \alpha = 0.05. CCDs exhibit moderate sensitivity to outliers, particularly at center points, where a single aberrant observation can inflate pure error estimates and bias lack-of-fit tests; however, multiple replicated center points mitigate this by enabling robust pure error estimation via the mean square among replicates, independent of the model fit. This feature provides an internal check for experimental variability, enhancing design reliability in noisy environments. The D-efficiency metric, which measures the relative volume of the confidence ellipsoid for parameters compared to an optimal design, is given by \left( \frac{\det(X'X)}{\det(X_{\text{opt}}'X_{\text{opt}})} \right)^{1/p}, where p is the number of model parameters; higher values signify better design optimality, with rotatable CCDs often achieving 80-95% D-efficiency relative to theoretical optima for k = 3-5 factors. Variations in \alpha, such as arithmetic means, can boost D-efficiency by 5-10% over standard cuboidal settings in constrained regions.

Analysis and Implementation

Model Fitting Process

The model fitting process for a (CCD) involves estimating the parameters of a second-order polynomial response surface model using ordinary least squares regression, applied to the design matrix X that incorporates the factorial, axial, and center points, including columns for linear terms (x_i), quadratic terms (x_i^2), and cross-product terms (x_i x_j). The model is typically expressed as y = \beta_0 + \sum \beta_i x_i + \sum \beta_{ii} x_i^2 + \sum \beta_{ij} x_i x_j + \epsilon, where the parameter vector \beta is estimated via \hat{\beta} = (X^T X)^{-1} X^T y, with y as the vector of observed responses. This approach allows for the quantification of linear, quadratic, and interaction effects while accounting for experimental error. The process begins with fitting a first-order linear model using data from the factorial portion of the CCD (or fractional factorial) and the center points. A statistical test for curvature is then performed by comparing the average response at the center points to the average at the factorial points, often via a t-test or analysis of variance (ANOVA) on the difference; if significant (e.g., p-value < 0.05), this indicates nonlinearity and justifies proceeding to the quadratic model. Replicated center points are essential here, as they provide an estimate of pure experimental error for the test. If curvature is detected, axial (star) points are incorporated, and the full second-order model is refit using least squares on the complete dataset. Model adequacy is assessed through diagnostics, including the lack-of-fit test, which compares the model's residual sum of squares to the pure error sum of squares derived from replicated center points via an F-test; a non-significant p-value (e.g., > 0.05) supports the model's fit. Additionally, the measures the proportion of total variability explained by the model, while the adjusted R^2 accounts for the number of predictors to avoid overestimation in smaller designs. In cases of non-normality or heteroscedasticity in residuals (detected via plots or tests like Shapiro-Wilk), a Box-Cox power transformation may be applied to the response variable to stabilize variance and achieve approximate normality, selecting the optimal \lambda parameter (e.g., via maximum likelihood) before refitting the model. Finally, intervals for the coefficients \beta are computed using the estimated variance-covariance (X^T X)^{-1} \hat{\sigma}^2, where \hat{\sigma}^2 is the , and for predicted responses on the surface, joint regions or prediction intervals are derived to quantify uncertainty.

Sequential Experimentation

Central composite designs (CCDs) are well-suited for sequential experimentation in , enabling the iterative construction of second-order models by building upon initial screening experiments to efficiently explore and optimize response surfaces. This approach, pioneered by Box and Wilson, begins with a such as $2^{k-p} to identify significant factors among the k potential variables, typically requiring a III or higher to estimate main effects and low-order interactions with minimal runs. Following factor screening, a path of steepest ascent is followed using the fitted model to shift the experimental region toward the anticipated optimum, with step sizes scaled to the factor units (e.g., increasing by 5 and concentration by 0.25 g/L per unit change in the coded model). Upon detecting nonlinearity—often through replicated center points showing a response plateau or decline—this transitions to modeling by augmenting the into a CCD. Augmentation leverages the existing points as the base block of the , adding $2kaxial points at a distance\alphafrom the center (commonly\alpha = \sqrt{2^k}for rotatability in full factorials) and 3 to 6 additional center points to estimate pure error and pure [quadratic](/page/Quadratic) effects. This modular assembly allows flexibility, as a2^{k-p}$ fractional can similarly form the factorial portion if does not confound terms. Decision rules for augmentation rely on analysis of variance (ANOVA) from the initial , where a significant lack-of-fit test (e.g., via replicated centers) or normal probability plots indicating prompts the addition; for instance, if observed center responses deviate substantially from linear predictions (e.g., $688 vs. &#36;670), terms are warranted. The adaptive nature of this strategy reduces total experimental runs—often from 30+ for a standalone to 20-25 by reusing screening data—while systematically approaching the optimum, as illustrated in process optimization workflows starting with 8-16 screening trials, 4-8 ascent steps, and 8-12 augmentations. Model fitting then follows standard least-squares to estimate the second-order . A key limitation is the assumption of negligible higher-order interactions or effects beyond second degree, potentially requiring redesign if the fitted surface exhibits poor adequacy or if the process involves categorical factors that cannot be continuously ascended.

Software Tools

Several packages facilitate the generation and analysis of central composite designs (CCDs) in . The rsm package (version 2.10.6 as of March 2025) provides functions for creating CCDs, including and support for or fractional blocking structures, while also enabling plots for response surface visualization. Similarly, the DoE.base package, often used in conjunction with DoE.wrapper, supports CCD construction through orthogonal arrays and designs, allowing for flexible specification of factors and replications. Commercial software offers robust platforms tailored for (DOE), with specialized CCD capabilities. JMP's DOE platform (version 19 as of September 2025) includes templates for CCDs up to eight factors, supporting options like rotatable, orthogonal, and uniform precision variants, along with blocking and replication controls. Minitab's Response Surface Design tool (version 22.3.0 as of 2025) generates CCDs for 2–10 factors, incorporating blocks, center points, and to model effects. Design-Expert (version 25.0.4 as of October 2025) emphasizes rotatable CCDs, providing menus for alpha value selection (with rotatable as default for up to five factors) and tools for optimizing models. In , the pyDOE library constructs matrices, supporting central, axial, and factorial points with customizable alpha values for . For model fitting post-design generation, the statsmodels package handles quadratic of data, including estimation of coefficients and prediction intervals. Key features across these tools include interactive menus for alpha selection to achieve rotatability or , options for random blocking to mitigate experimental bias, and export functionalities to generate design matrices in tabular formats for further analysis or simulation. A recommended when using software for CCDs is to verify after design generation, typically by examining the correlation matrix of the model terms or using built-in diagnostics to ensure unbiased coefficient estimates.

Applications

Response Surface Optimization

Once a response surface model has been fitted using data from a central composite design (CCD), optimization begins with techniques to identify potential optima. plots project the response surface onto a two-dimensional plane, displaying lines of constant response values to reveal regions of high or low performance, while three-dimensional surface plots provide a more intuitive depiction of the curvature and interactions between factors. These visualizations, generated from the fitted polynomial equation, help practitioners locate stationary points or desirable regions within the experimental space. For initial optimization, the method of steepest ascent is applied when a model indicates a favorable direction, moving along the path of maximum response increase proportional to the regression coefficients until the response plateaus or declines. For models derived from , more advanced techniques such as grid search—evaluating the response at discrete points across the factor space—or gradient-based methods iteratively adjust factors toward the defined by setting partial derivatives to zero. These approaches efficiently navigate the curved surface to approximate the optimum. In cases involving multiple responses, the desirability function approach, proposed by and Suich, transforms each individual response into a desirability score d_i (ranging from 0 to 1) based on predefined goals (e.g., maximization, minimization, or target), then combines them into an overall desirability D. The composite function is given by D = \left( \sum_{i=1}^k w_i d_i^p \right)^{1/p}, where k is the number of responses, w_i are weights summing to 1 reflecting response priorities, p determines the function's shape (often p=1 for linear or p=2 for quadratic), and optimization maximizes D over the factor settings. This method balances trade-offs without requiring a single composite objective. Canonical analysis further aids interpretation by rotating the to principal axes using the eigenvectors of the form's , yielding a simplified model in variables that reveals the surface's nature—such as a maximum (all negative eigenvalues), minimum (all positive), or (mixed signs)—and the distances to these points along each . This highlights the and , guiding decisions on whether the optimum lies within the experimental region. Finally, predicted optima are validated through confirmation runs, additional experiments at the suggested factor levels, to assess agreement between observed and modeled responses, accounting for experimental error and confirming practical achievability. Multiple replicates in these runs provide estimates of process variability at the optimum.

Industrial Case Studies

In , one of the earliest applications of central composite design (CCD) was demonstrated in the seminal work by Box and Wilson. In their 1951 paper, they illustrated the by optimizing the of a laboratory-scale (A + B + C → D + other products), using sequential designs including a CCD to explore factors such as reactant proportions, concentration, and reaction time. This approach led to substantial improvements in product by systematically exploring the of the response, highlighting CCD's efficiency in industrial process optimization over exhaustive full designs. In the , has been widely used for formulating oral solid to balance critical quality attributes like tablet hardness and drug dissolution rates. A representative case involved the development of fast-disintegrating tablets using as a , where a face-centered optimized the concentrations of , sodium starch glycolate (superdisintegrant), and spray-dried . The design assessed effects on tablet hardness and disintegration time, resulting in an optimized with balanced strength and rapid drug release profiles. This multi-response optimization underscores 's role in ensuring consistent drug release in tablet . In , particularly injection molding, CCD facilitates tuning of process parameters like melt and injection pressure to enhance product and operational efficiency. For example, a face-centered CCD has been applied to optimize three key parameters—melt , packing pressure, and cooling time—for producing parts, using 20 runs to minimize warpage and reduce cycle time, thereby enhancing production efficiency and meeting dimensional tolerances. This demonstrates CCD's practical value in reducing defects and boosting throughput in polymer processing. In , CCD supports sensory and nutritional optimization of baked goods through multi-response desirability functions. For cookie formulation, a central composite design was used to optimize composite flour cookies from germinated , , and flours, adjusting ingredient proportions to enhance , (spread ratio, crispness), and overall sensory acceptability. This illustrates CCD's utility in balancing , , and benefits in product development. Key lessons from these industrial applications emphasize 's efficiency in run reduction compared to traditional multi-level factorials; for instance, a three-factor (k=3) rotatable requires about 20 runs versus 27 for a full 3^k design, enabling faster experimentation while providing robust second-order estimates for optimization. This sequential approach allows practitioners to confirm first-order models before augmenting to , minimizing resource use in resource-constrained industrial settings.