Displaying 1 result from an estimated 1 matches for "ts_lomb".
Did you mean:
ts_load
2010 Jan 21
0
Using spec.ls to interpolate very long time series
....) I want to use spec.ls as an
interpolator and then use the output to reconstruct the time series via
inverse fft. But so far I've been having difficulty doing this.
ts<-read.csv("timeseries.csv",header=TRUE) #file contains over 30000
irregularly spaced observations
attach(ts)
ts_lomb<-spec.ls(tsspec,y=NULL,demean=TRUE,detrend=FALSE)
>length(ts_lomb$freq)
[1]17496
#Now to see what some of the frequencies look like.
>ts_lomb$freq[1:10]
[1] 2.857796e-05 5.715592e-05 8.573388e-05 1.143118e-04 1.428898e-04
[6] 1.714678e-04 2.000457e-04 2.286237e-04 2.572016e-04 2.8577...