# NOT RUN {
#REX
library(psd)
##
## Multitaper PSD estimation
##
set.seed(1234)
X <- rnorm(1e3)
# use the defaults, and appeal to plot.spec
# sampling assumed to be 1
plot(psdcore(X))
# use more tapers, compare to stats::spectrum, and clear
# env data from the previous calculation
psdcore(X, ntaper=10, plot=TRUE, refresh=TRUE)
# change the sampling frequency to 20
psdcore(X, X.frq=20, ntaper=10, plot=TRUE, refresh=TRUE)
# }
# NOT RUN {
#REX
# }
Run the code above in your browser using DataLab