14.03.2014 Views

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 3 Fitting St<strong>and</strong>ard Least Squares Models 87<br />

Effect Screening<br />

The resulting posterior probabilities are listed <strong>and</strong> plotted with bars as shown in Figure 3.25. An overall<br />

posterior probability is also shown for the outcome that the sample is uncontaminated.<br />

Bayes Plot for Factor Activity<br />

JMP includes a script that helps you determine which factors are active in the design.<br />

1. Open the Reactor.jmp sample data table.<br />

1. Open the BayesPlotforFactors.jsl file in the Sample Scripts folder.<br />

2. Select Edit > Run Script.<br />

3. Select Y <strong>and</strong> click Y, Response.<br />

4. Select F, Ct, A, T, <strong>and</strong> Cn <strong>and</strong> click X, Factor.<br />

5. Click OK.<br />

Figure 3.26 Bayes Plot for Factor Activity<br />

In this case, we specified that the highest order interaction to consider is two. Therefore, all possible models<br />

that include (up to) second-order interactions are constructed. Based on the value assigned to Prior<br />

Probability (see the Controls section of the plot), a posterior probability is computed for each of the possible<br />

models. The probability for a factor is the sum of the probabilities for each of the models where it was<br />

involved.<br />

In this example, we see that Ct, T, <strong>and</strong> Cn are active <strong>and</strong> that A <strong>and</strong> F are not. These results match those of<br />

the Bayes Plot shown in Figure 3.25.<br />

Note: If the ridge parameter were zero (not allowed), all the models would be fit by least squares. As the<br />

ridge parameter gets large, the parameter estimates for any model shrink toward zero. Details on the ridge<br />

parameter (<strong>and</strong> why it cannot be zero) are explained in Box <strong>and</strong> Meyer (1993).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!