# NOT RUN {
## calculate the MODWT of the sunspots series
W <- wavMODWT(sunspots)
## identify the boundary coefficients
z <- wavBoundary(W)
## plot the results
plot(wavShift(z))
## obtain a summary
summary(z)
# }
Run the code above in your browser using DataLab