suggest = function(X,level=0.05){
spec = spectrum(X,log='no',plot=FALSE)
spx = spec$freq
spy = 2*spec$spec
z.score = scale(spy,scale=TRUE)
p.value = 2*(1-pnorm(abs(z.score)))
spx = spx[p.value<=level]
spy = spy[p.value<=level]
S = round(1/spx)
return(sort(S))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.