Fig5.3s <- function() { # postscript(file = "Fig5.3.ps", height = 5.46, width = 5.46, horiz = F) par(mar = c(0, 0, 0, 0)) x <- (1:1024)/1024 y <- sin(2 * pi * x) y[513:1024] <- y[513:1024] + sin(4 * pi * x[513:1024]) y[769:1024] <- y[769:1024] + ((1:256) * sin((10 * pi * (1:256))/256))/128 y <- y + (1:1024)/1024 y[201:280] <- y[201:280] + sin(((1:80) * 30 * pi)/80) plot(mrd(y, wavelet = "d10"), cex = 0.7) # This line is to clear out the x-axis. lines(c(0, 1050), c(-0.5, -0.5), lwd = 50, col = 0) # graphics.off() NULL }