similar to: ts and stl functions

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr <- ts(t(tkr), deltat=1/12) or tstkr <- ts(c(tkr), deltat=1/12) and use the stl function I get the following error: Error in

Hi, I have imported a csv file into R which contains one column (the rate er 100,000 population of a disease, by month over 11 years) I coerced into a time series using the following function, tstkr<-ts(tkr,deltat=1/12) This seems to work fine, and when I check for the class of the object using class(tstkr) I get "ts" as the response. When I try to use the stl function in stats

Hi, I am getting unexpected behaviour from aggregate.ts(). The 'ndeltat' argument is effectively being reduced by 1 in some cases, even when it is an integer, with the result that the blocks to be aggregated are not of the expected size, and also that the end() of the aggregated series is much later than the end() of the original series. rawts <- ts(rep(1:10, each = 5), start = 1) ##

Hello, Does anyone know of a way I can speed this up? Basically I'm attempting to get the data item on the same row as the report date for each report date available. In reality, I have over 11k of columns, not just A, B, C, D and I have to do that over 100 times. My solution is slow, but it works. The loop is slow because of merge. # create sample data z.dates =

I am running R 1.2.2 under Linux. When using the function stl in the ts library, how can I save the seasonal component? What I would like was something like: library(ts) data(nottem) data.stl <- stl(nottem, "per") x <- data.stl$sea This what I get: > x NULL I would, however, like to store in x the seasonal component. Thanks in advance. Francisco. -- Francisco

Hi R community! I'd like to create a temporal series with ts command for the interval 1st January 2002 - 31 December 2008. I have daily values, so this is a 2557 days temporal series. I'm using this command : ts(observations, start=2002, end=2009, freq=365) However, I don't get the correct temporal series since both frequency (365 OR 366 for leap years) and deltat (1/365 OR

Hi, I'm wondering why the spectrum() phase of quadrature couple isn't purely +/-pi. But mostly, I'm looking for a recommended way to take a 2-D spectrum and convert it into a single complex array. Kindly consider: # 10 Hz sine wave 10 seconds long sampled at 50 Hz deltaT = 1/50 t = seq(0, 10, deltaT) w = 2 * pi * 10 x = ts( sin( w * t ), deltat = deltaT ) y = ts( sin(

Dear R user, I have written a function which returns max,min and variation of a power (see below) Power is a given matrix(1,n) I call the function >Variation<-VAR(p,(n-deltat)) Now the problem is when I want plot(Results[1],Results[2]). Not possible! I become the following error (in english it means: Error in as.double.default(x) :Object cannot be transformed in double) >

Time Series functions in R ========================== I think a good basic S-like functionality for library(ts) in base R would include ts class, tsp, is.ts, as.ts plot methods start end window frequency cycle deltat lag diff aggregate filter spectrum, spec.pgram, spec.taper, cumulative periodogram, spec.ar? ar -- at least univariate by Yule-Walker arima -- sim, filter, mle, diag, forecast

This message is in MIME format --_=XFMail.1.3.p0.Linux:990702182137:16900=_ Content-Type: text/plain; charset=us-ascii There was some discussion a while back on R-devel between Ross Ihaka, Paul Gilbert and myself about row subsetting in time series. I think the consensus was that "[.ts" should not try to coerce its result back to a time series object (which is underlying the problem

>>>>> On Fri, 02 Jul 1999, Adrian Trapletti <Adrian.Trapletti@wu-wien.ac.at> said: Adrian> There seems to be a problem with plot.ts (R Version 0.64.2) > x<-cbind(1:10,2:11) > x<-as.ts(x) > plot(x) Adrian> Error: subscript (20) out of bounds, should be at most 10 This is definitely a bug --> CC: R-bugs ALL NOTE : This is *not* new

Hi. I have to forecast a time series of a Internet network traffic bitrate. The data are in file http://www.forumaltavilla.it/joomla/datitesi/dati.datand the sampling time is every 0.05 seconds. Now, i want to use HoltWinters forecasting. This is my script. dt=1.58443823e-9 #0.05 seconds in years dati.ts=ts(scan("dati.dat"),start=0,deltat=dt) model=HoltWinters(dati.ts)

OK. What if I have a time series which is collected every Monday, please? What is the proper way to use the start option within the ts command in order to indicate that this is Monday data, please? Thanks again! Sincerely, Erin

Hello R Users, I have a question regarding fitting a model with GAM{mgcv}. I have data from several predictor (X) variables I wish to use to develop a model to predict one Y variable. I am working with ecological data, so have data collected many times (about 20) over the course of two years. Plotting data independently for each date there appears to be relationships between Y (fish

Hi, I try to figure out what the stl-function exactly do. I was reading the paper by Cleveland et al. (1990) and tested some features of stl (the ability to decompose time series with missing values and the robustness feature). I tried the following: > data(co2) > co2.na <- co2 > is.na(co2.na[c(50, 100)]) <- TRUE > plot(stl(co2.na, s.window = 12, na.action = na.exclude)) With

Dear R users, I'm trying to write a function which returns minimum,maximum,mean of a vector(power) I've done the following : VAR<-function(power,length){ for(i in tml:length)){ tvar[i]<-i pmean[i]<-mean(power[i:i+deltat]) pmin[i]<-min(power[i:i+deltat]) pmax[i]<-max(power[i:i+deltat]) varmax[i]<-100*(pmax[i]-pmean[i])/pmean[i]

Hi there, My name is Renan X. Cortes, student of Statistics, from south of Brazil, and I'd like to ask you a few questions about decomposition of time series. In R, when I fit the decomposition using the "stl" function, an object is returned when ask the summary of the fit, called STL.seasonal (%), STL.trend (%) and STL.remainder (%). Once the decomposition is additive,

I am trying to understand the nuances of STL (seasonal trend decomposition with loess) based on William Cleveland's (and others?) original development. I do not understand the specification or use of "frequency components" or equivalent "low-pass filter" components in the stl() function. I have run the stl() function on a standard example data (co2) in both S-Plus and

Dear helpeRs, Is it possible to change the shading colour of the range bars in the plot generated by plot.stl? By default they are grey, but I would prefer them white (I am preparing some graphics for a powerpoint presentation so I'm inverting all colours). As far as I can see plot.stl allows you to turn off the range bars, but nothing about the shading colour. I tried to look at the

Hi everyone. I'm having some troubles with STL function to decompose some data. My issue is that I have monthly data from September 2005 up to August 2010 i. e. 60 observations. I define it in the following way: *u<-read.csv("C:/CELEBREX.csv",header = TRUE) u.ts<-ts(u, start=c(2005,9), frequency=12) * The issue is that when I try to use stl(u.ts, 'per') Then the