[R] ggplot and spec.pgram
Any Ideas to get an interactive periodogram? -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis __ R-help@r-project.org mailing list 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.
[R] Identifying x, y coordinates from periodogram generated by spec.pgram
I would like to identify the x,y coordinates in a plot of the periodogram spec.pgram(a, taper=0, log=no) I have tried the identify() function to no avail I have also tried to pass the periodogram over to rggobi e-plot.spec(spec.pgram(a, taper=0, log=no)) ggobi(e) to no avail I need to identify the predominant frequencies so that I can work back to the time of recurrence in minutes. I don't necessarily need a interactive graphical solution (although this may be easier). thanks -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis __ R-help@r-project.org mailing list 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.
[R] time series
I am trying to work with a time series of one observation every fifteen minutes for a length of two years. My question really lies in how to code the timeseries to represent every observation is a fifteen minute interval. Below is the code that I used any help would be appreciated. x-read.table(doconc.txt, header=TRUE) x-ts(x[,1], start=0, frequency=1/15) plot(stl(x, s.window=periodic)) Error in stl(x, s.window = periodic) : series is not periodic or has less than two periods Thanks Stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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.
[R] Power Spectral Sensity
I am working with a dissolved oxygen dataset. continuous readings are taken at 15 minute intervals and we have been recording these data at 12 stations along the savannah river for two years now. The longest set of readings that are continuous without interuption is 53 days. I would like to look at the power spectral density at each of these sites (most likely one day will be the overridding trend due to photosynthesis and respiration). I would also like to look at the cross spectrums between to contiguous site to look for advective transport from the upstream sections.I have used spectrum to get a periodogram. This I don't believe is the power spectral density. If it is I am unsure how to determine the # in the span argument to smooth the periodogram. I imagine that I am overlooking something- I appreciate all of your help. stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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.