[R] Simple spectral analysis

From: Georg Hoermann <georg.hoermann_at_gmx.de>
Date: Mon 08 Jan 2007 - 16:12:34 GMT

Hello world,

I am actually trying to transfer a lecture from Statistica to R and I ran into problems with spectral analysis, I think I just don't get it 8-(
(The posting from "FFT, frequs, magnitudes, phases" from 2005 did not enlighten me)

As a starter for the students I have a 10year data set of air temperature with daily values and I try to get a periodogram where the annual period (365 days) should be clearly visible (in statistica I can get the frequencies and the period). I tried the spectrum() and pgram() functions, but did not find a way through... The final aim would be to get the periodogram (and the residuals from the reassembled data set...)

Thanks and greetings,

The data set:

air = read.csv("http://www.hydrology.uni-kiel.de/~schorsch/air_temp.csv") airtemp = ts(T_air, start=c(1989,1), freq = 365) plot(airtemp)

Georg Hoermann, Dep. of Hydrology, Ecology, Kiel University, Germany
+49/431/23761412, mo: +49/171/4995884, icq:348340729, skype: ghoermann

R-help@stat.math.ethz.ch mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Tue Jan 09 04:24:44 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 10 Jan 2007 - 04:30:27 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.