search for: stl

Hi all, How can the label of the x-axes in the plot() of a stl.object be adapted? e.g., When plotting: plot(stl(nottem, "per")) In the labels of the x-axes is “time”. How can this be changed to e.g., “Time (dekade) “? It does not work with xlab or others anymore… Thanks, Jan ______________________________________________________________...

Hi Jack, I'm almost there with Android..... I've actually got the static and dynamic libraries built. I'm choking on the test suite. Do you want to take a shot at working around Android and [embedded, lame] STLport? I would try removing the call to rend() in std::map<_Key, _Tp, _Compare, _Alloc>::rend(). But I really should not modify program code. Fiddling with a Makefile is one thing, but modifying source code is another. Jeff arm-linux-androideabi-g++ --sysroot=/opt/android-ndk-r8e//platforms/...

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() funct...

Hi, I have a monthly time series with missing values and I would use stl function to identify seasonality. I tried all settings of na.action but the result is the same: stl(tm245,s.window=11, na.action=na.pass) Error in stl(tm245, s.window = 11, na.action = na.pass) : NA/NaN/Inf in foreign function call (arg 1) Can you help me? Thanks Andrea Toreti [[al...

I don't understand the output of stl. As a simple example: y <- numeric(1:365) y[250] = 1 stl <- stl(ts(y, frequency=7), s.window="periodic") This returns without error but the results are puzzling to me. If you plot the results it is probably easiest to visualize what I mean. plot(stl) This shows the original dat...

...ree series of monthly data from 2007 to 2009. I would like to show the season decomposition of these two series side by side on the one device, however using plot doesn't seem to respect any use of layout(matrix(1:3, ncol=3)) or par(mfcol=c(1,3)). I'm guessing that this means that the plot(stl) perhaps uses them, but I can't find anywhere the / a plot.stl() - ie, I can't work out where the plot() call is going? I've attached a small example of data and some probably overly verbose code. There are only two sets of example data, not three as mentioned above. # load the data...

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...

...m not sure how to handle. I have 15 years of quarterly data and I typically create a TS object via something like... data.ts <- ts(mydata, start = 2002, frequency = 4) this create a matric as opposed to a vector object as I receive a univariate error when I try to decompose the data using the STL function data.stl <- stl(data.ts, "periodic") Error in stl(data.ts, "periodic") : only univariate series are allowed ok so is.vector(data.ts) [1] FALSE so to convert to a vector I'll use data.ts <- as.vector(data.ts) but then I lose the frequency as the periods as...

On 1/21/2013 12:35 AM, Chris Lattner wrote: > > I'm confused here. You're acting as though we don't use the STL. In fact, we do use std::string, std::vector, std::map etc when they are the right solution for the job. I'm trying to understand the reasoning behind the decisions made at the beginning of LLVM. My working assumption is that ADT didn't exist when LLVM started (whereas STL did). In su...

On Jan 20, 2013, at 7:46 AM, Krzysztof Parzyszek <kparzysz at codeaurora.org> wrote: >> It is one of the motivations. > > The reason I ask is that STL comes all ready, with containers and algorithms. They may not be optimal for every task, but they do their job and they are part of the standard. There may be some price to pay in terms of performance/memory usage/etc. for a specific application, but overall it may be worth it. Evidently, in cas...

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...

Hello On Wed, 2008-01-23 at 00:10 +0300, Nikita Danilov wrote: > Hello, > > below is a tentative list of tasks into which WBC effort can be > sub-divided. I also provided a less exact list for the EPOCH component, > and an incomplete list for the STL component. > > WBC tasks are estimated in lines-of-code with the total of (9100 + 3000) > LOC, where LOC is a non-comment, non-whitespace line of the source > file. I believe it is too early to estimate all EPOCH tasks, hopefully > there will be more clear understanding of the situa...

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...

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 >...

Greetings: I'm using the time series decomposition routine "stl" from the package "stats". But how do I get the results into a vector to work with them? example: data(AirPassengers) m<-stl(AirPassengers,"per") print(m) This lists the output but can't figure out how to extract the individual series like seasonal, trend, irregu...

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 function by...

> Hello, > > I have lots of data in zoo format and would like to do some time > series analysis. (using library(zoo), library(ts) ) > > My data is usually from one year, and I try for example stl() to find > some seasonalities or trends. > > I have now accepted, that I might have to convert my series into ts() > but still I am not able to execute the comand since stl() is not > satisfied > > x<-zoo(rnorm(365), as.Date("2005-01-01"):as.Date("2005-12-3...

Hi All, I am in the process of writing an R extension in c++ and am using several STL containers (e.g., vector<double>, map<int, double>, multimap<int, double>). I make sure to clear all these containers at the end of the .Call. Everything compiles and runs just fine, but I'm a bit worried since I haven't found any other packages that use STL. So, my q...

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 thought that this would be some kind of decompositon of the variability of the time serie in terms of seasonal, trend and r...

Hi, There is a discussion thread on llvm-commits list about a possibility of using custom memory allocators for STL to improve the performance and reduce the memory pressure of STL containers, e.g. std::set. I thought that this discussion may be interesting for a wider audience and therefore I re-post the messages on llvm-dev as well. It would be interesting to hear what others think about - custom allocators...