Fig5.14n <- function() { # postscript(file = "Fig5.14.ps", height = 3.9, width = 7.37, horiz = F) par(mfrow = c(1, 2), mar = c(1.5, 1.5, 1.5, 1.5), mgp = c(5, 0.4, 0)) rs <- c(57, 14, 55, 51, 30, 0, 53, 44, 34, 53, 49, 2) x <- (1:128)/128 f <- c(rep(0, 32), rep(4, 64), rep(0, 32)) .Random.seed <- rs y <- f + rnorm(128) plot(wd(y, filter.number = 5), main = "", sub = "") mtext("Raw data", side = 3, line = 0.1) yshrunk <- threshold(wd(y, filter.number = 5), type = "soft", policy = "manual", value = 1.669) plot(yshrunk, main = "", sub = "") mtext("Treated data", side = 3, line = 0.1) # graphics.off() NULL }