This macro calculates the scaled prediction variance for up to five design matrices given a specific model. You can assess the scaled prediction variance of a design, evaluate the rotatability of a design, and compare designs to determine how the prediction variance behaves across the design space. This will let you compare designs in coded units and their performance prior to collecting any data for a DOE. The designs are limited to 2 or more factors.
Be sure that Minitab knows where to find your downloaded macro. Choose . Under Macro location browse to the location where you save macro files.
If you use an older web browser, when you click the Download button, the file may open in Quicktime, which shares the .mac file extension with Minitab macros. To save the macro, right-click the Download button and choose Save target as.
You need up to five design matrices.
Suppose the factors of a CCD have been copied to a matrix M1, and the factors of a Box-Behnken design has been stored in M2. To run the macro, choose and type:
%VDG M1 M2;
MODEL "quadratic";
RESOLUTION 150.
Click Run.
Variance Dispersion Graphs and Design of Experiments
Variance dispersion graphs are graphical tools used in the process of selecting a design of experiment (DOE). These graphs show the scaled prediction variance of a DOE across the design space and are typically used to compare response surface designs. However, variance dispersion graphs can also be used to compare the performance of multiple designs for a specific model, such as a linear model, a linear model with interaction terms, a linear model with quadratic terms, or a full quadratic model.
History and Background of VDG
The VDG was originally developed by Giovannitti-Jensen and Myers (1989) and Myers et al. (1992) to display a snapshot of the stability of the scaled prediction variance across the design space.
Mathematically, the scaled prediction variance can be estimated with the equation below:
where v(x) is the scaled prediction variance at the point x and, as shown above, does not depend on the response data, only on the design matrix X. In other words, the scaled prediction variance is a function of the number of runs in your experiment, the covariance matrix, and the point at which you are evaluating the prediction variance. Because the scale prediction variance displayed by the VDG is independent of the response data, the VDG allows you to assess the goodness of a design a priori.
In general, a VDG is constructed with four basic components:
The VDG shows the average, the maximum, and the minimum scaled prediction variances at different values of r, which represents the distance of any given point to the origin or center of the design.
Using a VDG to Evaluate Other Aspects of a Design
Using a VDG, you can evaluate the performance of a single design in terms of its prediction variance. You can also use the VDG to compare two or more designs. Another application of a VDG is to evaluate the effect of adding center points to a design.
You can also use a VDG to assess whether a design is rotatable. In a rotatable design, all three lines (minimum, average and maximum) in the graph will coincide. If the minimum, average, and maximum lines do not coincide, the face centered design is not rotatable.
Monte Carlo Simulation to Estimate Prediction Variance
This advanced section describes the mathematical methods used to estimate the minimum, maximum, and average prediction variances displayed in a VDG. If you are interested only in how to run a macro to display and interpret a VDG you may want to skip this section.
As mentioned previously, for a fixed value of the radius, a search needs to be performed to find the minimum, maximum, and integrated scaled prediction variance. The method selected for the Minitab macro is based on Monte Carlo simulation and the ability to search uniformly at equidistant points from the center x = (0,0,…,0)T of the design.
Three quantities need to be calculated at different distances from the origin: 1) the minimum prediction variance, 2) the maximum prediction variance, and 3) the integrated prediction variance. There are multiple ways to find these quantities. One method is by using a nonlinear programming algorithm to find the minimum and maximum prediction variance. However, the method we have chosen to estimate these quantities is by uniformly sampling and finding a large set of distinct points that are the same distance from the center of the design for different values of r. From the space of equidistant points, one can guarantee the proper estimation of the minimum and maximum prediction variances by having two conditions satisfied: First, the sampling of equidistant points needs to be uniform. In other words, there can‟t be regions where some equidistant points are more likely to be chosen than others. Second, the number of points needs to be large enough to cover the whole space of possible design points.
Because all the points that need to be sampled are equidistant, the problem of sampling is equivalent to finding a good way of obtaining points from an n-ball or hypersphere for a given radius.
The procedure is described in Marsaglia (1972) and summarized below:
We have illustrated how to sample points in a hypersphere of constant radius r. We can use this sampling technique to estimate the maximum or minimum prediction variance. Using these samples, from the Weak Law of Large Numbers, the average prediction variance of the sampled design points can be used to estimate the true integrated prediction variance as long as n is large.
References
[1] Marsaglia G. "Choosing a Point from the Surface of a Sphere". The Annals of Mathematical Statistics, Vol. 43, No. 2 (1972), pp. 645-646.
[2] Myers R, Montgomery D. (2002). Response Surface Methodology - Process and Product Optimization Using Designed Experiments, Second Edition. John Wiley & Sons, Inc.