#install.packages("astsa") require(astsa) #help(package = astsa) par(mfrow=c(2,1)) k = kernel("daniell", 4) soi.ave = spec.pgram(soi, k, taper=0, log="no") abline(v=c(.25,1,2,3), lty=2) soi.ave = spec.pgram(soi, k, taper=0, log="yes") abline(v=c(.25,1,2,3), lty=2) out = SigExtract(soi, L=9, M=64, max.freq=.05)