[R] Simple spectral analysis

2007-01-08 Thread Georg Hoermann
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 o

Re: [R] Simple spectral analysis

2007-01-08 Thread Earl F. Glynn
"Georg Hoermann" <[EMAIL PROTECTED]> wrote in message news:[EMAIL PROTECTED] > 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) Maybe this will get you started using fft or spectrum -- I

Re: [R] Simple spectral analysis

2007-01-08 Thread Peter Dalgaard
Earl F. Glynn wrote: > "Georg Hoermann" <[EMAIL PROTECTED]> wrote in message > news:[EMAIL PROTECTED] > >> 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) >> > > Maybe this w

Re: [R] Simple spectral analysis

2007-01-08 Thread Georg Hoermann
Peter Dalgaard wrote: > Earl F. Glynn wrote: >> > The defaults for detrending and tapering could be involved. Putting, > e.g., detrend=F gives me a spectrum with substantially higher > low-frequency components. > > But what was the problem in the first place? > understanding how this things

Re: [R] Simple spectral analysis

2007-01-08 Thread Petr Pikal
HTH Petr On 8 Jan 2007 at 17:12, Georg Hoermann wrote: Date sent: Mon, 08 Jan 2007 17:12:34 +0100 From: Georg Hoermann <[EMAIL PROTECTED]> To: r-help@stat.math.ethz.ch Subject: [R] Simple spectral analysis > Hello world,

Re: [R] Simple spectral analysis

2007-01-09 Thread Earl F. Glynn
"Georg Hoermann" <[EMAIL PROTECTED]> wrote in message news:[EMAIL PROTECTED] > Peter Dalgaard wrote: >> Earl F. Glynn wrote: > Thanks a lot for the help. I will post the script when its ready > (an introduction for our biology students to time series, just 8 hours) I've been working with one of