similar to: Time Series smoothing

Dear all, I am applying the StructTS function in ts-package. For some time series the program terminates and the following error appears: Error in optim(init[mask], getLike, method = "L-BFGS-B", lower = rep(0, : L-BFGS-B needs finite values of fn Do someone know what do I have to adjust in the original time series to avoid this error? It works fine for some subsets of the

Hi all, I have been experimenting with writing my own matrix column sum function. I want it to return a list. csum<-function(m) { a = data.frame(m) s = lapply(a,sum) return(s) } I wish to use the same code up until the return(s) that I have listed above. The problem is that s, I believe, is a data frame (looks like this:) $X1 [1] 148 $X2 [1] 156 $X3 [1] 164 $X4 [1] 172

Alright, so I am trying to write my own function to calculate column sums in a matrix. I want the result as a single list with the values. So far I have: csum<-function(m) { a = data.frame(m) s = lapply(a,sum) return(s) } What is the easiest way to have it return in a format such as [1] 6 15 24 ? Thanks. -- View this message in context:

Hi all, I've experienced this problem using StructTS(x) when the *first* element of x is a NA (R:R1.7.0, os: w2ksp4). Please look at the following code: a=rep(1:7,10) library(ts) #this works StructTS(a) #this works x=a x[2]=NA StructTS(x) #this doesn't work x=a x[1]=NA StructTS(x) The last command returns this error "Error in optim(init[mask], getLike,

I've just built R 2.9.2 from source on Slackware Linux 13.0 - 32-bit (will try 64-bit also next) - and seen: > Collecting examples for package 'stats' > Running examples in package 'stats' > Error: testing 'stats' failed > Execution halted > make[3]: *** [test-Examples-Base] Error 1 Looking at R-2.9.2/tests/Examples/stats-Ex.Rout.fail I see: ... >

Hello all, I am looking for some help in understanding how to change the way R optimizes the smoothing constant selection process for the HoltWinters function. I'm a SAS veteran but very new to R and still learning my way around. Here is some sample data and the current HoltWinters code I'm using: rawdata <- c(294, 316, 427, 487, 441, 395, 473, 423, 389, 422, 458, 411, 433, 454,

Hi, I've noted that not all `predict' methods are hidden in the namespace: > methods("predict") [1] predict.ar* predict.Arima* [3] predict.arima0* predict.glm [5] predict.HoltWinters* predict.lm [7] predict.loess* predict.mlm [9] predict.nls* predict.poly [11] predict.ppr* predict.prcomp* [13]

Hello everyone, I was wondering if anyone was aware of a way in which I could use ggplot's stat_smooth() function for add an exponential moving average. I was thinking that I could maybe use something like: >myggplot + stat_smooth (method = 'HoltWinters( data , .9 , 0, 0)') but my efforts were futile. Perhaps there is a way to write my own custom method to throw into

I want to use double-exponential smoothing to forecast time series datas,but I couldn't find it in the document,does R support this method?

Dear helpers I'm having a problem with function HoltWinters from package ts. I have a time series that I want to fit an Holt-Winters model and make predictions for the next values. I've already built an object of class ts to serve as input to HoltWinters. But then I get an error; I've used HoltWinters a lot of times and this never hapened > data.HW<-HoltWinters(data.ts) Error

Hi, I run R on a server via SSH, over a terminal. After loading a specific package, how do I know what functions are in that package? Is there any way to list or display them? tia. -- View this message in context: http://www.nabble.com/Check-functions-in-package-tp24707084p24707084.html Sent from the R help mailing list archive at Nabble.com.

Dear list members I want to analyze separately the months of a time series. In other words, I want to plot and fit models for each month separately. Taking the example of http://a-little-book-of-r-for-time-series.readthedocs.org/en/latest/src/timeseries.html births <- scan("http://robjhyndman.com/tsdldata/data/nybirths.dat") birthstimeseries <- ts(births, frequency=12,

Hi everyone, I'm using the HoltWinters() function to do a time series analysis. The function only returns the back fitted values ($fitted) after the first year of data, which is my case, is a little more than half. However, when I use the plot() function, it plots the back fit for almost the entire data set. Any ideas on how to extract the fitted values going all the way back to the start

Hi R-users, I have a question relating to the HoltWinters() function. I am trying to forecast a series using the Holt Winters methodology but I am getting some unusual results. I had previously been using R for Windows version 2.7.2 and have just started using R 2.9.1. While using version 2.7.2 I was getting reasonable results however upon changing versions I found I started to see unusual

Full_Name: John Bodley Version: 2.5.1 (2007-06-27) OS: Windows XP Submission from: (NULL) (12.144.182.66) I was fitting a number of time series in R using the stats::HoltWinters method to define a single exponential smoothing model, i.e., beta = gamma = 0. I came across an example where the fitted value of alpha was not defined in the [0, 1] interval which seems to violate the lower and upper

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)

Hey All, First time on these forums. Thanks in advance. Soooo... I have a process that was functioning well before the 2.14 update. Now the HoltWinters function is throwing an error whereby I get the following: Error in HoltWinters(sales.ts) : optimization failure I've been looking around to determine why this happens (see if I can test the data beforehand) but I haven't come

Hello, I would like to look at the compiled C code behind HoltWinters from the stats package. Is that possible? If so where do I find it? thanks, Spencer [[alternative HTML version deleted]]

Good Morning! May you help me? I need to understand the function predict. I need to understand the algorithm implemented, the calculations associated. Where can I find this information? Thank You!

Dear All I am a biginner of R. I have difficulty with reading the code of a method. I am using the vars package to estimate a VAR model and I want to view the code of "predict" method for objects with class attribute "varest". I thougt I could just type the name "predict" without anything to display the code of the method as I often do with generic function.