rob.dwt(x, ..., span=5, cutoff1=2.5, cutoff2=3.75, max.resid.level=3)
All the information for the robust wavelet transform is stored in dictionary.
The details and algorithm are given in the S+WAVELETS User's Manual, in the section "Robust Smoother-Cleaner Wavelets".
xx <- make.signal("blocks") yy <- xx + rcauchy(length(xx)) # add Cauchy noisepar(mfrow=c(2,2)) plot(xx, type="l", xlab="Blocks", ylab="") plot(yy, type="l", xlab="Blocks with Cauchy Noise", ylab="") rw <- rob.dwt(yy, wavelet="d2", span=5) zz <- reconstruct(rw) plot(zz, type="l", xlab="Cleaned Noisy Blocks", ylab="") ss <- wavesmooth(zz, wavelet="d2") plot(ss, type="l", xlab="Smoothed Blocks", ylab="")