similar to: [gently off topic] arima seasonal question

> It might clarify your thinking to note that a seasonal ARIMA model > is just an ``ordinary'' ARIMA model with some coefficients > constrained to be 0 in an efficient way. E.g. a seasonal AR(1) s = > 4 model is the same as an ordinary (nonseasonal) AR(4) model with > coefficients theta_1, theta_2, and theta_3 constrained to be 0. You > can get the same answer as from

Hi all, I just couldn't find a R function which can fit multiple seasonal patters... i.e. in the following code: *arima(x = data, order = c(p, d, q), seasonal = list(order = c(P, D, Q), period = S), ... *** * there can be only one "period", am I right? What if the data seem to have three different seasonality cycles, 5, 12, 21? Thanks a lot! * [[alternative HTML version

I am using AUTO ARIMA for forecasting. But it is not detecting 'seasonality term' of its own for any data. Is there any other method by which we can detect seasonality and its frequency for any data? Is there any method through which seasonality and its frequency can be automatically detected from ACF plot? -- Siddharth Arun, 4th Year Undergraduate student Industrial Engineering and

Dear Members, I have been searching for a package in R which can handle multiple seasonality suggested by taylor(2003). It will be great help if anybody has used this on R before (i.e. which package). Thanks in Advance. Best Regards Atul Malik [[alternative HTML version deleted]]

Hello, does anyone know about: a) simulating seasonal ARIMA models? arima out of package ts can fit it, but it does not look like it can simulates data from seasonal models b) fitting and simulating fractional seasonal ARIMA models? Hints will be appreciated, Henning -- Henning Rust Potsdam Institute for Climate Impact Research Dept. Integrated Systems Analysis Tel.: #49/331/288-2596

I'm trying to use the following command. arima (x, order = c(p,d,q), seasonal =list(order=c(P,D,Q), period=s) How can I write an estimated seasonal ARIMA model from the outputs. To be specifically, which sign to use? I know R uses a different signs from S plus. Is it correct that the model is: (1-ar1*B-ar2*B^2-...)(1-sar1*B^s-sar2*B^2s-....)(1-B)^d(1-B^s)^D

Hello, I am new to usin R - which is a great tool - and would like to know if R has a seasonal adjustment program for time series and/if it incorporates the Census Bureau's ARIMA x-12 seasonal adjustment program in any way? Thanks so much! Tony [[alternative HTML version deleted]]

Hi list, I would like to estimate and forecast the seasonal component of a series. My model which uses daily data would be something y t = alpha + beta x SeasComp t + gamma x OtherRegressors t. One approach to this would be use quarterly dummies, another to use a sine function. The first would cause a step change when we move from a season to another; the latter impose too much regularity in

Hello, I have daily wind speed data and need to fit seasonal ARIMA model, problem is that my period is 365. But when I use arima(...) function, with period 365, I?m getting error message: ?Error in makeARIMA(trarma[[1]], trarma[[2]], Delta, kappa) : maximum supported lag is 350?. Can someone help me with this problem? Thank you Sincerely yours, Laura Saltyte

Dear all, I stumbled upon what appears to be a troublesome issue when sampling from an ARIMA model (from Rob Hyndman's excellent 'forecast' package) that contains a seasonal AR component. Here's how to reproduce the issue. (I'm using R 2.9.2 with forecast 2.19; see sessionInfo() below). First some data: > x <- c( 0.132475, 0.143119, 0.108104, 0.247291, 0.029510,

Are any seasonal adjustment programs, like Tramo/Seats, Census X12 ARIMA or Berliner Verfahren implemented in R? I am doing a simulation study and I don't know how to adjust the series in R. The possibility to access external the exe.files of the seasonal adjustment programs seems to be quite difficult. Can anyone help me? Thanks, Ingo

hallo to all I've to calculate an arima model and I need only the first and 365 th parameter and also the sar1 and the intercept, so I'm traing with: arima(X,order=c(365,0,0),seasonal=list(order=c(1,0,0),..),fixed=c(NA,rep(0,363),NA,NA,NA),transform.pars=F) but the error answer is: Error in polyroot(z) : polynomial degree too high (49 max) also there are problems in allocating memory

Hello, I am very new to R and Time Series. I need some help including R codes about the following issues. I' ll really appreciate any number of answers... # I have a time series data composed of 24 values: myinput = c(n1,n2...,n24); # In order to make a forecasting a, I use the following codes result1 = arima(ts(myinput),order = c(p,d,q),seasonal = list(order=c(P,D,Q))) result2 =

Hi All, a) Is it possible to estimate the strength of seasonality in timeseries data. Say I have monthly mean prices of an ten different assets. I decompose the data using stl() and obtain the seasonal parameter for each month. Is it possible to order the assets based on the strength of seasonality? b) which gives a better estimate on seasonality stl() or a robust linear model like

Dear R community, I have a project with electricity load forecasting, and I got hourly data for system load. If you haven't worked with electricity before, seasonality comes in many flavors: a daily pattern, with a peak at around 7pm; a weekly pattern, in which we use more electricity on weekdays in comparison to weekends; a winter-summer pattern, with air conditioning and heaters playing an

Hi all, I was wondering how to construct a seasonal differenced time series variable. I used the following code to construct a 12 span seasonal difference seasonal<-diff(V2, lag=12, differences=1) is this correct? thank you in advance joe [[alternative HTML version deleted]]

Dear R-help, How could a cross-correlation plot be optimized such that the relationship between seasonal time-series can be studied? We are working with strong seasonal time-series and derived a cross-correlation plot to study the relationship between time-series. The seasonal variation however strongly influences the cross-correlation plot and the plot seems to be ?rather? symmetrical (max

Good evening! I have a question regarding forecast package and time series analysis. My syntax: x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279, 304, 275, 276, 313, 292, 302, 322, 281, 298, 305, 295, 286, 327, 286, 270, 289, 293, 287, 267, 267, 288, 304, 273, 264, 254, 263, 265, 278) library(forecast) arima(x, order=c(1,1,2),

Hello everybody, I need to calculate seasonal means with temperature data for my work. I have 70 files coming from weather stations, which looks like this for example: startdate <- as.POSIXct("01/01/2006", format = "%d/%m/%Y") enddate <- as.POSIXct("05/01/2006", format = "%d/%m/%Y") date <- seq(from = startdate, to = enddate, by =

Hello, In order to do a future forecast based on my past Time Series data sets (salespricesproduct1, salespricesproduct2, etc..), I used arima() functions with different parameter combinations which give the smallest AIC. I also used auto.arima() which finds the parameters with the smallest AICs. But unfortuanetly I could not get satisfactory forecast() results, even sometimes catastrophic