Example of a Noise Factor in the Prediction Profiler

For example, use the Tiretread.jmp sample data set. This data set shows the results of a tire manufacturer’s experiment whose objective is to match a target value of HARDNESS= 70 based on three factors: SILICA, SILANE, and SULFUR content. Suppose the SILANE and SULFUR content are easily (and precisely) controllable, but SILICA expresses variability that is worth considering.

For comparison, first optimize the factors for hardness without considering variation from the noise factor.

Select Graph > Profiler to launch the Prediction Profiler.

Assign Pred Formula HARDNESS to the Y, Prediction Formula role.

Click OK.

Select Optimization and Desirability > Desirability Functions in the Prediction Profiler menu.

Double-click in the Desirability plot to open the Response Goal window. Select Match Target from the list.

Select Optimization and Desirability > Maximize Desirability to find the optimum factor settings for our target value of HARDNESS.

We get the following Prediction Profiler display. Notice that the SILICA factor’s optimum value is on a sloped part of a profile curve. This means that variations in SILICA are transmitted to become variations in the response, HARDNESS.

Note: You might get different results from these because different combinations of factor values can all hit the target.

Figure 2.20 Maximizing Desirability for HARDNESS

Now, we would like to not just optimize for a specific target value of HARDNESS, but also be on a flat part of the curve with respect to Silica. So, repeat the process and add SILICA as a noise factor.

Select Graph > Profiler.

Select Pred Formula HARDNESS and click Y, Prediction Formula.

Select SILICA and click Noise Factors.

Click OK.

Change the Pred Formula Hardness desirability function as before.

The resulting profiler has the appropriate derivative of the fitted model with respect to the noise factor, set to be maximized at zero, its flattest point.

Figure 2.21 Derivative of the Prediction Formula with Respect to Silica

Select Optimization and Desirability > Maximize Desirability to find the optimum values for the process factor, balancing for the noise factor.

This time, we have also hit the targeted value of HARDNESS, but our value of SILICA is on its flatter region. This means variation in SILICA does not transmit as much variation to HARDNESS.

Figure 2.22 Maximize Desirability

You can see the effect this has on the variance of the predictions by following these steps for each profiler (one without the noise factor, and one with the noise factor):

Select Simulator from the platform menu.

Assign SILICA to have a random Normal distribution with a standard deviation of 0.05.

Figure 2.23 Setting a Random Normal Distribution

Click Simulate.

Click the Make Table button under the Simulate to Table node.

Doing these steps for both the original and noise-factor-optimal simulations results in two similar data tables, each holding a simulation. In order to make two comparable histograms of the predictions, we need the two prediction columns in a single data table.

Copy the Pred Formula HARDNESS column from one of the simulation tables into the other table. They must have different names, like Without Noise Factor and With Noise Factor.

Select Analyze > Distribution and assign both prediction columns as Y.

When the histograms appear, select Uniform Scaling from the Distribution main title bar.

Figure 2.24 Comparison of Distributions with and without Noise Factors

The histograms show that there is much more variation in Hardness when the noise factor was not included in the analysis.

It is also interesting to note the shape of the histogram when the noise factor was included. In the comparison histograms above, note that the With Noise Factor distribution has data trailing off in only one direction. The predictions are skewed because Hardness is at a minimum with respect to SILICA, as shown in Figure 2.25. Therefore, variation in SILICA can make only HARDNESS increase. When the non-robust solution is used, the variation could be transmitted either way.

Figure 2.25 Prediction Profiler Showing the Minima of HARDNESS by SILICA