“…MetalCutting[c(62,15,6,17,41,28,36,55), 1:7] R> y <- MetalCutting[c(62,15,6,17,41,28,36,55) R> y <-c(y, y_TOP_DES) R> es8.aug.OBsProb <-OBsProb(X = X, y = y, abeta = 1, bbeta = 1, blk = 1, + mFac = 6, mInt = 2, nTop = 64) R> print ( The update posterior probability of factors is displayed in Figure 5 on the right side. The result is comparable with that reported in Edwards et al (2014).…”