[R] spec.pgram returns different spectra when fast=TRUE and the number of samples is odd

Duncan Murdoch murdoch@dunc@n @end|ng |rom gm@||@com
Mon Oct 19 22:11:19 CEST 2020


Just curious:  did you read the help page to see what the "fast" 
argument does?  If you did, why are you surprised at this result?

Duncan Murdoch

On 19/10/2020 1:53 p.m., Jade Mosley via R-help wrote:
> Dear all,
> 
> This is potentially a bug in spec.pgram, when the number of samples is odd,spec.pgramreturns a different result withfast = TRUE, the example below contains the two varieties with a reference spectrum calculated manually. the number of returned spectra is also larger (50 compared to 49) whenfast = TRUE
> 
> x <- rnorm(
> 
> 99
> 
> )
> plot(spec.pgram(x, taper =
> 
> 0
> 
> , detrend =
> 
> FALSE
> 
> , plot =
> 
> FALSE
> 
> , fast =
> 
> FALSE
> 
> )$spec, ty =
> 
> 'l'
> 
> )
> lines(spec.pgram(x, taper =
> 
> 0
> 
> , detrend =
> 
> FALSE
> 
> ,plot =
> 
> FALSE
> 
> , fast =
> 
> TRUE
> 
> )$spec, col =
> 
> 'red'
> 
> )
> lines(((abs(fft(x))^
> 
> 2
> 
> )[
> 
> 2
> 
> :
> 
> 50
> 
> ])/
> 
> 99
> 
> , col =
> 
> 'blue'
> 
> , lty =
> 
> 2
> 
> )
> 
> Best regards
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list