[R] Lomb periodograms

Josué Polanco jomopo at gmail.com
Tue Mar 31 22:45:09 CEST 2009

Hi All, hi Bernardo,

I think that the problem is when you remove 2*n data, then the problem
is when the number of the frequencies is computed. Look in the source
code, may be there is the key. If you have further  problems, you
could consider to use a Lomb-Scargle function that I wrote, send me an
email if you want.  Remember that the RAW Lomb-Scargle periodogram is
not a good  estimation  of the spectral power.

Cheers, and good luck,

Josue Polanco

On Tue, Mar 31, 2009 at 6:33 PM, Garcia Carreras, Bernardo
<bernardo.garcia-carreras08 at imperial.ac.uk> wrote:
> Hi,
> I have recently used the CTS package in order to use the Lomb-Scargle periodogram (spec.ls) function. I have noticed an issue that I hoped you may be able to explain. If a regularly spaced time series has two points removed, one at either side of a single data point (thus making an irregularly spaced time series), a spectrum with a very large peak at the highest frequencies is produced. An example of this is shown below:
> a <- runif(100)
> x <- 1:length(a)
> a <- a[-c(4,6)]
> x <- x[-c(4,6)]
> spec.ls(x=x, y=a, fast=FALSE, taper=0)
> Does anyone know why this is? Thank you very much in advance for any help you can give me with this issue.
> Best regards,
> Bernardo
> ______________________________________________
> R-help at 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.

Josué Mosés Polanco Martínez
Correo-e alternativo jomopo at linuxmail.org
It is a wasted day unless you have learned something new and made
someone smile -Mark Weingartz.

More information about the R-help mailing list