Plotting of Posterior Probabilities from Bayesian Screening
Method function for plotting marginal factor posterior probabilities for Bayesian screening.
list. List of class
Magnification used for the axis annotation.
additional graphical parameters passed to
A spike plot, similar to barplots, is produced with a spike for each factor.
Marginal posterior probabilities are used for the vertical axis.
X2, ... are used to label the factors
otherwise the original factor names are used.
print.BsProb function is called
and the posterior probabilities are displayed.
BsProb is called for more than one value of gamma (
the spikes for each factor probability are overlapped to show the
resulting range of each marginal probability.
The function is called for its side effects. It returns an invisible
Box, G. E. P and R. D. Meyer (1986). "An Analysis for Unreplicated Fractional Factorials". Technometrics. Vol. 28. No. 1. pp. 11–18.
Box, G. E. P and R. D. Meyer (1993). "Finding the Active Factors in Fractionated Screening Experiments". Journal of Quality Technology. Vol. 25. No. 2. pp. 94–105.
1 2 3 4 5 6 7 8 9 10 11 12 13 14
library(BsMD) data(BM86.data,package="BsMD") X <- as.matrix(BM86.data[,1:15]) y <- BM86.data["y1"] # Using prior probability of p = 0.20, and k = 10 (gamma = 2.49) drillAdvance.BsProb <- BsProb(X = X, y = y, blk = 0, mFac = 15, mInt = 1, p = 0.20, g = 2.49, ng = 1, nMod = 10) plot(drillAdvance.BsProb) summary(drillAdvance.BsProb) # Using prior probability of p = 0.20, and a 5 <= k <= 15 (1.22 <= gamma <= 3.74) drillAdvance.BsProbG <- BsProb(X = X, y = y, blk = 0, mFac = 15, mInt = 1, p = 0.25, g = c(1.22, 3.74), ng = 3, nMod = 10) plot(drillAdvance.BsProbG, code = FALSE, prt = TRUE)