## calculate a MODWT of the sunpots series and
## verify the class
W <- wavMODWT(sunspots)
print(class(W))
## summarize the object
summary(W)
## reconstruct the MODWT of the sunspots series
## and compare to the original
sunup <- reconstruct(W)
vecnorm(sunup - sunspots)
Run the code above in your browser using DataLab