Fig8.3n <- function() { # postscript(file = "Fig8.3.ps", horizontal = F, width = 6.5, height = 6.5) par(mfrow = c(2, 2), mar = c(1.5, 1.5, 1.5, 0.5), mgp = c(5, 0.4, 0)) rs <- c(57, 14, 55, 51, 30, 0, 53, 44, 34, 53, 49, 2) .Random.seed <- rs x <- (1:512)/512 f <- sin(2 * pi * x) ssig <- sqrt(var(f)) f <- ((f - mean(f)) * 2)/ssig fnoise <- f + rnorm(512) fwd <- wd(fnoise, fil = 5) plot(x, f, type = "l") mtext(side = 3, line = 0.1, "True function") plot(x, fnoise, type = "l") mtext(side = 3, line = 0.1, "Noisy observations") fwdsure <- threshsure(fwd) plot(x, wr(fwdsure), type = "l") mtext(side = 3, line = 0.1, "SURE estimator") fwdhyb <- threshhyb(fwd, lowlev = 2, seed = rs) plot(x, wr(fwdhyb), type = "l") mtext(side = 3, line = 0.1, "SURE hybrid estimator") # graphics.off() NULL }