Bayes plot for effect screening

Ben_Sagalovsky

In the ''old" BayesPlotForFactors.jsl script there was an option to indicate the highest order interaction to consider when determining the posterior probability of a factor being active. Now that the Bayes plot has been included in Standard Least Squares under Effect Screening (hooray!) there is no place to indicate the highest order interactions to include in the calculation. Does anybody know what is the value used "by default", and why there is no longer that flexibility? 

 

Also, there is a change from Ridge Parameter Gamma to K contamination; what would the correspondence be among those two?

Thanks!

Mark_Bailey

Confusion Alert! The Bayes plot provided in the Fit Least Squares platform and the one provided by the script are not the same method, although they are related.

 

First, the built-in plot was added to JMP a long time ago. (JMP 3?) It is based on the 1986 method by Box et. al. that only considers main effects. You cannot change that computation.

 

Second, the scripted plot was added to JMP sample files 10-15 years ago. It is based on the 1993 method by Box et. al. that was adapted to consider higher order terms.

 

So you have the chronology backwards.

 

The computations are completely different. The older method used K to define the ratio of the variance of the estimates with and without contamination (presence of non-null effects). You could specify individual prior probabilities for each parameter. The newer method uses a ridge parameter for shrinkage. It also uses a prior probability of 0.25 for all parameters. They both consider all potential terms within their respective scopes.

 

The papers behind each of these methods is very accessible. Box was a wonderful teacher and writer.