Output of BOXMEYER from input file PBExamp1.txt HEADING LINE in PBExamp1.txt: From Table 7.5 p287 of Box-Hunter-Hunter Values of variables: Nrows=12 Ncols=11 Nvars=5 Pi=0.25 Gamma=2.5 Alpha=2 Labels in PBExamp1.txt: A=Feed_Rate B=Catalyst C=Agitation_Rate D=Temperature E=Concentration Data in PBExamp1.txt: Rownum, Design matrix, Response variable Row A B C D E F G H J K L 1. 1 -1 1 -1 -1 -1 1 1 1 -1 1 56 2. 1 1 -1 1 -1 -1 -1 1 1 1 -1 93 3. -1 1 1 -1 1 -1 -1 -1 1 1 1 67 4. 1 -1 1 1 -1 1 -1 -1 -1 1 1 60 5. 1 1 -1 1 1 -1 1 -1 -1 -1 1 77 6. 1 1 1 -1 1 1 -1 1 -1 -1 -1 65 7. -1 1 1 1 -1 1 1 -1 1 -1 -1 95 8. -1 -1 1 1 1 -1 1 1 -1 1 -1 49 9. -1 -1 -1 1 1 1 -1 1 1 -1 1 44 10. 1 -1 -1 -1 1 1 1 -1 1 1 -1 63 11. -1 1 -1 -1 -1 1 1 1 -1 1 1 63 12. -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 61 Factor labels: ABCDEFGHJKL Factors analyzed: ABCDE Models and submodels (n=25): A B C D E AB AC AD AE BC BD BE CD CE DE ABC ABD ABE ACD ACE ADE BCD BCE BDE CDE Best 5 for each model size: Model AdjRsq Rmse Fstat* Likelihood Rmsec BPostProb ----------------------------------------------------------------------- B 0.45899 11.4054 10.3325 2.5e-020 11.4827 0.066 E 0.03756 15.2124 1.4292 7.8e-022 15.2267 0.003 D -0.03592 15.7824 0.6186 5e-022 15.7888 0.002 A -0.05754 15.9463 0.4015 4.4e-022 15.9505 0.0018 C -0.09719 16.2424 0.0256 3.6e-022 16.2427 0.0015 ------------------------------------------------------------------- BD 0.81146 6.7330 16.7812 5.3e-017 7.0062 0.23 BE 0.51479 10.8012 4.8902 1.8e-019 10.9308 0.0017 AB 0.38656 12.1450 3.3105 4.5e-020 12.2437 0.00049 BC 0.32764 12.7148 2.7868 2.6e-020 12.8019 0.0003 DE 0.25174 13.4133 2.2336 1.4e-020 13.4870 0.00017 -------------------------------------------------------------------- BDE 0.96985 2.6926 51.5452 2e-010 3.7616 0.7 ABD 0.82013 6.5765 8.1649 4.5e-015 7.1772 0.00057 BCD 0.79985 6.9372 7.2799 2.4e-015 7.5942 0.00031 ACE 0.71927 8.2158 5.0263 3.1e-016 8.8427 5.7e-005 ADE 0.41671 11.8427 2.1227 3.9e-018 12.1587 1.7e-006 All models sorted together by Fstat (best 15 out of 25): Model AdjRsq Rmse Fstat* Likelihood Rmsec BPostProb ----------------------------------------------------------------------- BDE 0.96985 2.6926 51.5452 2e-010 3.7616 0.7 BD 0.81146 6.7330 16.7812 5.3e-017 7.0062 0.23 B 0.45899 11.4054 10.3325 2.5e-020 11.4827 0.066 ABD 0.82013 6.5765 8.1649 4.5e-015 7.1772 0.00057 BCD 0.79985 6.9372 7.2799 2.4e-015 7.5942 0.00031 ACE 0.71927 8.2158 5.0263 3.1e-016 8.8427 5.7e-005 BE 0.51479 10.8012 4.8902 1.8e-019 10.9308 0.0017 AB 0.38656 12.1450 3.3105 4.5e-020 12.2437 0.00049 BC 0.32764 12.7148 2.7868 2.6e-020 12.8019 0.0003 DE 0.25174 13.4133 2.2336 1.4e-020 13.4870 0.00017 ADE 0.41671 11.8427 2.1227 3.9e-018 12.1587 1.7e-006 ABE 0.37148 12.2933 1.9288 2.5e-018 12.5684 1.2e-006 BCE 0.33405 12.6541 1.7883 1.7e-018 12.9665 8.5e-007 ABC 0.30078 12.9663 1.6760 1.3e-018 13.1786 7.1e-007 E 0.03756 15.2124 1.4292 7.8e-022 15.2267 0.003 All models sorted together by BayesPostProb (best 15 out of 25): Model AdjRsq Rmse Fstat Likelihood Rmsec BPostProb* ----------------------------------------------------------------------- BDE 0.96985 2.6926 51.5452 2e-010 3.7616 0.7 BD 0.81146 6.7330 16.7812 5.3e-017 7.0062 0.23 B 0.45899 11.4054 10.3325 2.5e-020 11.4827 0.066 E 0.03756 15.2124 1.4292 7.8e-022 15.2267 0.003 D -0.03592 15.7824 0.6186 5e-022 15.7888 0.002 A -0.05754 15.9463 0.4015 4.4e-022 15.9505 0.0018 BE 0.51479 10.8012 4.8902 1.8e-019 10.9308 0.0017 C -0.09719 16.2424 0.0256 3.6e-022 16.2427 0.0015 ABD 0.82013 6.5765 8.1649 4.5e-015 7.1772 0.00057 AB 0.38656 12.1450 3.3105 4.5e-020 12.2437 0.00049 BCD 0.79985 6.9372 7.2799 2.4e-015 7.5942 0.00031 BC 0.32764 12.7148 2.7868 2.6e-020 12.8019 0.0003 DE 0.25174 13.4133 2.2336 1.4e-020 13.4870 0.00017 ACE 0.71927 8.2158 5.0263 3.1e-016 8.8427 5.7e-005 AC 0.07256 14.9332 1.2869 3.7e-021 14.9805 5.3e-005 Posterior weights that factors are active: Weighting models by Fstat: A: 26.064 0.20510 ****** B: 110.484 0.86940 ************************** C: 21.8341 0.17181 ***** D: 90.8754 0.71510 ********************* E: 73.0175 0.57458 ***************** Weighting models by Bayesian post.prob. (pi=0.25, gam=2.5): A: 0.00299699 0.00300 B: 0.991472 0.99147 ****************************** C: 0.00219343 0.00219 D: 0.925193 0.92519 **************************** E: 0.700799 0.70080 *********************