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.