DOE Pro XL includes 2 Level Full and Fractional Factorial designs, 3 Level Full Factorial, Taguchi, Plackett-Burman, Central Composite, Box - Behnken, and Custom Designs. If you need help, use the Computer Aided design selection feature to pick the best design for your problem. Stat-ease design-expert 11. 4 free download. Disintegration time measured in box-behnken design runs. Design-expert® download design-expert is a software for design. Applications of plackett–burman and central composite design for. Colorado matters: npr. Software downloads. What's new and exciting in v10 of design-expert® software.
Digging Into DOE
Selecting the right central composite design
for response surface methodology applications
_________________
by Richard Verseput
Research and development are the heart and soul of improvement efforts in manufacturing, and it's fast becoming standard practice to employ design of experiments methods in industrial R&D. In the early stages of their work, experimenters typically use screening experiment designs that normally consist of trials run at the extreme lower- and upper-bound level setting combinations of the variable study ranges. They provide information on the direct additive effects of the study variables and on pairwise (two-variable) interaction effects. Screening designs enable experimenters to select the best materials and equipment from available alternatives and to focus on the correct variables and ranges for further study.
In the later stages of the experimental work, the goal shifts from screening to product and process optimization. The statistical experiment designs most widely used in optimization experiments are termed 'response surface designs.' In addition to trials at the extreme level settings of the variables, response surface designs contain trials in which one or more of the variables is set at the midpoint of the study range (other levels in the interior of the range may also be represented). Thus, these designs provide information on direct effects, pairwise interaction effects and curvilinear variable effects. Response surface methodology, one approach to product and process optimization work, derives its name from the use of these widely used optimization experiment designs.
Most practitioners of RSM now generate their experiment designs and analyze their data using a statistical software program running on a personal computer. Many of these software programs can generate many classes of RSM designs and, in some cases, offer several varieties of each class. However, the central composite design is the most popular of the many classes of RSM designs due to the following three properties:
A CCD can be run sequentially. It can be naturally partitioned into two subsets of points; the first subset estimates linear and two-factor interaction effects while the second subset estimates curvature effects. The second subset need not be run when analysis of the data from the first subset points indicates the absence of significant curvature effects.
CCDs are very efficient, providing much information on experiment variable effects and overall experimental error in a minimum number of required runs.
CCDs are very flexible. The availability of several varieties of CCDs enables their use under different experimental regions of interest and operability.
Three main varieties of CCD are available in most statistical software programs: face-centered, rotatable and inscribed. Many experimenters are uncertain about the best variety to use in a given study. To make the right selection, the experimenter must first understand the differences between these varieties in terms of the experimental region of interest and region of operability, according to the following definitions.
Region of interest --a geometric region defined by lower and upper limits on study-variable level setting combinations that are of interest to the experimenter
Region of operability --a geometric region defined by lower and upper limits on study-variable level setting combinations that can be operationally achieved with acceptable safety and that will output a testable product
The face-centered CCD
Figure 1: Two-Variable Face-Centered CCD |
Note that the radius, designated a, determines the geometry of the design region. An a of 1.0 defines a square design geometry (a cube for three variables, a hypercube for four or more variables, etc.). As the value of a increases, the axial points extend beyond the faces of the square and the design region becomes more spherical.
To successfully carry out any designed experiment, the region of operability must encompass the region of interest. In this example, this means that the process must be able to be operated at all time and temperature level setting combinations in the region of interest described by the variable ranges. Therefore, the first step in selecting the right variety of CCD is to compare the region of interest to the region of operability. If the process can't be operated at one or more of the extremes of the region (in this case the corners of the square), then a face-centered CCD is inappropriate. This leaves the experimenter with two options: reduce the variable ranges and generate a new face-centered CCD, or generate an inscribed CCD. The process for making this decision will be discussed in the section on inscribed designs.
Figure 2: Two-Variable Rotatable CCD |
Figure 2 presents a rotatable CCD for time and temperature using the previously defined lower and upper variable bounds. As the figure shows, the rotatable CCD uses an a value of 1.4 to describe a circular design geometry (a sphere for three variables, a hypersphere for four or more variables, etc.). You may recognize this design as the one used in Example 14-2 of Design and Analysis of Experiments by Montgomery.1
For a given number of variables, the a value required to achieve rotatability is computed as a = 4 ¥&F where F is the number of points in the 2k factorial design (k is the number of study variables). Rotatability refers to the uniformity of prediction error. In rotatable designs, all points at the same radial distance (r) from the center point have the same magnitude of prediction error. Note that in most rotatable designs the factorial points and axial points lie on different concentric geometric spheres and thus have different magnitudes of prediction error. For example, in the 23 rotatable CCD, the factorial points are located on a sphere of radius 1.73, while the axial points are located on a sphere of radius 1.68.
Setting the a to ¥k changes the design from a rotatable CCD to a spherical CCD. In spherical CCDs, all design points occur on the same geometric sphere. These designs are not exactly rotatable, but they are near-rotatable. As stated in Myers and Montgomery,2 'If the desired region of the design is spherical, the CCD that is most effective from a variance point of view is to use a = ¥k and three to five center runs.' Generally a spherical design region is desired when the experimenter strongly suspects that the response will express a simple maxima or minima within the study variable ranges. As previously stated, it's required when the process can't be operated at one or more of the extremes of the design region.
Central Composite Design software, free download
The value of making rotatability (or near-rotatability) a design goal becomes clear when you consider that the experimenter doesn't know the location of the optimum point within the region of interest before the experiment is conducted, so it's desirable that all points a given distance (r) from the center point in any direction have the same magnitude of prediction error.
Relative to the a value of 1.0 in the face-centered design, the rotatable design a value of 1.4 has extended the design region beyond the defined variable bounds. Thus, predicted responses at or near the axial points, which would have been extrapolations in a face-centered design, are now within the design region. This is an important consequence, because the magnitude of prediction error increases geometrically with distance outside the design region.
However, extending the axial points beyond the defined variable bounds when using rotatable CCDs requires operating the process at five level settings of each variable (as will be shown later, the same is true of inscribed CCDs). The face-centered CCD, in contrast, requires operating the process at only three level settings of each variable. This makes the face-centered CCD a simpler design to carry out.
As previously stated, a face-centered CCD can be used when the region of operability encompasses the region of interest as defined by the variable bounds. To decide between a face-centered CCD and a rotatable (or spherical) CCD in this case, the experimenter must decide whether the rotatable design benefits of uniform prediction error and extension of the design region adequately offset the added complexity of operating the process at two additional level settings of each variable.
The experimenter should not take this added complexity lightly. Even when reconfiguring the process is not expensive and time-consuming --and it's normally both --the increased number of reconfigurations provides a much greater opportunity for sources of experimental error associated with setup and operation to express themselves. It's my experience that unexpectedly large experimental error is the most common cause of failure in studies employing designed experiments, and it makes sense to select options that reduce the chance of error expression. For the following two reasons, the benefits of rotatable designs do not offset the added complexity and associated risk in most cases, making the face-centered CCD more than adequate for most experiments.
Reason 1: Insufficient reduction in prediction error
|
For the Montgomery data, the 95-percent PI for the axial points is ±0.802 (recall that the axial points are design points in the rotatable CCD). For a face-centered CCD, predictions at the axial points would be extrapolations. However, the 95-percent PI for the axial points is ±0.991 for the face-centered CCD. In practical terms, a reduction from ±1.0 to ±0.8 in the 95-percent P.I. for a predicted yield of, say, 80 percent is of no real consequence. The same lack of practical consequence is true for reduction from ±14.9 to ±12.1 given a 5-percent RSD.
For the k = 3 designs, the reductions in the 95-percent PI seem more consequential. However, one must remember that the rotation a value increases from 1.4 in the k = 2 CCD to 1.68 in the k = 3 CCD. This means that the experimental ranges have been extended by 1.68 times the original ranges defined by the experimenter, for example from 10 minutes to 16.8 minutes for time.
Reason 2: Insufficient difference in effects estimation
Table 2: Regression model Coefficient Data--Exanmple 14-2 From Montgomery
Table 2 presents the estimated regression model coefficients and corresponding ±95-percent confidence intervals from Example 14-2 of Design and Analysis of Experiments by Montgomery for the yield response. The last column in the table presents the estimated coefficients using a 13-point face-centered CCD for the same variables. The actual yield data from Example 14-2 were used for the center points and factorial design points that the two designs have in common. Predicted yield data were used for the face-centered axial design points. The predictions were obtained using the regression model from Example 14-2 with a correctly propagated experimental error added to each predicted value. As the Table 2 data show, the face-centered CCD coefficients are all within the rotatable CCD coefficient ±95-percent confidence intervals. In fact, a comparison of the estimated coefficients for the two design types clearly shows that they are either identical or nearly so in all cases.
The inscribed CCD
|
The inscribed option is a convenient way to generate a rotatable CCD that enables the experimenter to study the full ranges of the experiment variables while excluding nonallowable operating conditions at one or more of the extremes of the design region. The excluded portion of the original face-centered CCD region is shown in gray in Figure 3, as are the excluded face-centered CCD points.
However, for the k = 3 inscribed CCD, the 95-percent PI for the eight excluded factorial points is ±1.074, as opposed to ±0.788 for the face-centered CCD that contains these points. This means that the prediction error will be 36-percent higher (1.074/0.788) for the extreme points that are excluded by the inscribing but are within the region of operability of the face-centered CCD. When only one extreme point is nonallowable, eliminating that run from a face-centered CCC only increases the prediction error for four design points --by 7.4 percent for three axial points and by 10.8 percent for one extreme point. It's better in this case to generate a face-centered CCD and simply not run the nonallowable point rather than suffer a 36-percent increase in prediction error for the seven extreme points excluded by inscribing but within the region of operability.
General guidelines for CCD selection
A rotatable or near-rotatable CCD should be used when nonallowable operating conditions exist at two or more of the extremes of the design region. The inscribed design may be the best rotatable design option in this circumstance, because experimenters commonly think of setting individual variable bounds in terms of practical or operational limits. The inscribed design lets the experimenter define the variable bounds in these terms while excluding nonallowable operating conditions at the extremes of the design region.
A face-centered CCD, on the other hand, should be used whenever the region of operability encompasses the full region of interest as described by the variable ranges or when nonallowable operating conditions exist at only one of the extremes of the design region.
In rotatable and near-rotatable CCDs, each experiment variable is represented at five levels. The face-centered CCD requires only three levels of each experiment variable, making it the simplest variety of CCD to carry out as well as the least prone to corruption due to sources of experimental error associated with setup and operation.
Compared to face-centered CCDs, rotatable or near-rotatable CCDs offer reduced prediction error for, and improved estimation of, quadratic (curvature) effects. However, given a reasonable magnitude of overall experimental error, say a 5-percent RSD, these benefits do not outweigh the added complexity of requiring each variable to be run at five levels.
References
Central Composite Design Example
1. Douglas C. Montgomery. Design and Analysis of Experiments , 4th Edition, John Wiley and Sons. New York: 1997.
2. Raymond H. Myers and Douglas C. Montgomery. Response Surface Methodology, John Wiley and Sons. New York: 1995.
About the author
Richard Verseput is president of S-Matrix Corp. ( www.s-matrix-corp.com). He has published many articles on applied DOE, and has created and taught applied DOE courses in the aerospace, automotive, biotechnology, chemicals and pharmaceutical industries. E-mail him at rverseput@qualitydigest.com .
Here's a representation of a classic central composite design for 2 factors.
- Four corners of the square represent the factorial (+/- 1) design points
- Four star points represent the axial (+/- alpha) design points
- Replicated center point
Numeric Factors: The number of numeric factors involved in the experiment
Categoric Factors: The number of categoric factors involved in the experiment
Note
For the Montgomery data, the 95-percent PI for the axial points is ±0.802 (recall that the axial points are design points in the rotatable CCD). For a face-centered CCD, predictions at the axial points would be extrapolations. However, the 95-percent PI for the axial points is ±0.991 for the face-centered CCD. In practical terms, a reduction from ±1.0 to ±0.8 in the 95-percent P.I. for a predicted yield of, say, 80 percent is of no real consequence. The same lack of practical consequence is true for reduction from ±14.9 to ±12.1 given a 5-percent RSD.
For the k = 3 designs, the reductions in the 95-percent PI seem more consequential. However, one must remember that the rotation a value increases from 1.4 in the k = 2 CCD to 1.68 in the k = 3 CCD. This means that the experimental ranges have been extended by 1.68 times the original ranges defined by the experimenter, for example from 10 minutes to 16.8 minutes for time.
Reason 2: Insufficient difference in effects estimation
Table 2: Regression model Coefficient Data--Exanmple 14-2 From Montgomery
Table 2 presents the estimated regression model coefficients and corresponding ±95-percent confidence intervals from Example 14-2 of Design and Analysis of Experiments by Montgomery for the yield response. The last column in the table presents the estimated coefficients using a 13-point face-centered CCD for the same variables. The actual yield data from Example 14-2 were used for the center points and factorial design points that the two designs have in common. Predicted yield data were used for the face-centered axial design points. The predictions were obtained using the regression model from Example 14-2 with a correctly propagated experimental error added to each predicted value. As the Table 2 data show, the face-centered CCD coefficients are all within the rotatable CCD coefficient ±95-percent confidence intervals. In fact, a comparison of the estimated coefficients for the two design types clearly shows that they are either identical or nearly so in all cases.
The inscribed CCD
|
The inscribed option is a convenient way to generate a rotatable CCD that enables the experimenter to study the full ranges of the experiment variables while excluding nonallowable operating conditions at one or more of the extremes of the design region. The excluded portion of the original face-centered CCD region is shown in gray in Figure 3, as are the excluded face-centered CCD points.
However, for the k = 3 inscribed CCD, the 95-percent PI for the eight excluded factorial points is ±1.074, as opposed to ±0.788 for the face-centered CCD that contains these points. This means that the prediction error will be 36-percent higher (1.074/0.788) for the extreme points that are excluded by the inscribing but are within the region of operability of the face-centered CCD. When only one extreme point is nonallowable, eliminating that run from a face-centered CCC only increases the prediction error for four design points --by 7.4 percent for three axial points and by 10.8 percent for one extreme point. It's better in this case to generate a face-centered CCD and simply not run the nonallowable point rather than suffer a 36-percent increase in prediction error for the seven extreme points excluded by inscribing but within the region of operability.
General guidelines for CCD selection
A rotatable or near-rotatable CCD should be used when nonallowable operating conditions exist at two or more of the extremes of the design region. The inscribed design may be the best rotatable design option in this circumstance, because experimenters commonly think of setting individual variable bounds in terms of practical or operational limits. The inscribed design lets the experimenter define the variable bounds in these terms while excluding nonallowable operating conditions at the extremes of the design region.
A face-centered CCD, on the other hand, should be used whenever the region of operability encompasses the full region of interest as described by the variable ranges or when nonallowable operating conditions exist at only one of the extremes of the design region.
In rotatable and near-rotatable CCDs, each experiment variable is represented at five levels. The face-centered CCD requires only three levels of each experiment variable, making it the simplest variety of CCD to carry out as well as the least prone to corruption due to sources of experimental error associated with setup and operation.
Compared to face-centered CCDs, rotatable or near-rotatable CCDs offer reduced prediction error for, and improved estimation of, quadratic (curvature) effects. However, given a reasonable magnitude of overall experimental error, say a 5-percent RSD, these benefits do not outweigh the added complexity of requiring each variable to be run at five levels.
References
Central Composite Design Example
1. Douglas C. Montgomery. Design and Analysis of Experiments , 4th Edition, John Wiley and Sons. New York: 1997.
2. Raymond H. Myers and Douglas C. Montgomery. Response Surface Methodology, John Wiley and Sons. New York: 1995.
About the author
Richard Verseput is president of S-Matrix Corp. ( www.s-matrix-corp.com). He has published many articles on applied DOE, and has created and taught applied DOE courses in the aerospace, automotive, biotechnology, chemicals and pharmaceutical industries. E-mail him at rverseput@qualitydigest.com .
Here's a representation of a classic central composite design for 2 factors.
- Four corners of the square represent the factorial (+/- 1) design points
- Four star points represent the axial (+/- alpha) design points
- Replicated center point
Numeric Factors: The number of numeric factors involved in the experiment
Categoric Factors: The number of categoric factors involved in the experiment
Note
For CCD's the number of runs required for the numeric factors will be multiplied by the number of categoric combinations. If two categoric factors having three levels each are added to this design the number of runs will be multiplied by nine. Consider optimal designs when you have both numeric and categoric factors.
Enter factor ranges in terms of …
+/-1 to define the limits for the area of interest where the optimum is believed to exist. Axial points will typically be outside this limit. Enter the limits in the Low and High columns.
alphas to ensure that even the extreme axial runs are within the area of operability. The area of interest must be within the area of operability. Enter the limits in the -alpha and +alpha columns.
Type: There will be different options available in this pull down based on how many factors are to be included in the experiment. The designs listed here form the factorial core of the central composite design. The default suggestion is the largest fraction that will produce a design under 1000 runs or maintain at least resolution V behavior. Smaller fractional cores can be used to save on the budget. For six or more factors we recommend using the Min Run Res V type to get both good estimates while keeping the number of runs under control. Small designs are the smallest recognized central composite design and are not recommended due to possible analysis issues.
Blocks: Central composite designs can be split into blocks. The factorial design is split into sub-fractions that support the two-factor interaction model and the axial points in the final block if needed.
Options button: Click on the options button to change the axial (alpha) distance which is how far the star points will be from the center in coded units. Use the help button on the options dialog for more information.
The Points area provides a preview of the number and type of runs that will be in the design.
