similar to: MAPE

I don't understand where is my problem. I tried the problem suggested by Prof. Brian Ripley data(USAccDeaths, package="ts") and it worked. So this is not a machine or configuration problem. This is my series (dataset1.ts) Jan Feb Mar Apr May Jun Jul Aug 1998 20957656 22280838 28048581 32286754 36456729 35575315 34127179 34363108 1999

Hi, I'm trying to go through the examples for function arima0() in ts package, i.e, >data(lh) >arima0(lh, order = c(1,0,0)) each time the call to arima0() causes a segmentation fault. I checked the earlier version (1.1.1) of R, the function arima0 works fine. Tracing the call indicates that the function "setup_starma" (in pacf.c under ts) interprets the addresses of the

Hi List, I am running R 2.8.0 on a Windows XP machine, running ggplot2 version 0.8.1 I want to label the bars in a faceted grid barplot. Reproducible R code is given below: #### reproducible facet barplot ##### library(ggplot2) # Dataset from which to create the barplot ml <- rep(1:10,2) vals <- rnorm(20,mean = 10, sd=1) type <- c(rep("MAPE",10),rep("AIC",10))

I did a regression with ARMA errors using arima0 with ari<-arima0(y,order=c(2,0,2),xreg=reg1,delta=-1) or ari<-arima0(y,order=c(2,0,2),xreg=reg1) where reg1 is the matrix of the regressors and when I see diag(ari$var.coef) I get negative terms. Do you know what this mean ? I try to change transform.pars to 0 or 1 but this crash R on Windows. Is it possible to test the significativity

I would like the text plotted with 'mtext' to be alighned like it is for printing on the console. Here is what I have: > print(emt) ME RMSE MAE MPE MAPE MASE original -1.034568e+07 1.097695e+08 2.433160e+07 -31.30554 37.47713 1.5100050 xreg 1.561235e+01 2.008599e+03 9.089473e+02 267.05490 280.66734

Dear R People: I'm trying to use the predict command on an arima0 object. I do the following: xm.arma <- arima0(xm2,order=c(1,0,1)) predict(xm.arma,n.ahead=2) and I get the message: Error in round(x, digits) : Non-numeric argument to mathematical function Any ideas what the problem might be, please? R version 1 4 1 on Windows. Thanks in advance! Sincerely, Erin Hodgess Associate

Dear all, this may be a stupid question but... The underlying model in procedure arima0 is X[t] = a[1]X[t-1] + ... + a[p]X[t-p] + e[t] + b[1]e[t-1] + ... +b[q]e[t-q] Is it possible to get an estimate of e for every point t, t-1 etc. or at least an estimate of the variance of e? Thanks a lot in advance for any hints Kai Arzheimer

What is the difference between arima and arima0?

[accidentally sent to owner-r-help -- please do NOT! it's "r-help" !] ------- start of forwarded message ------- From: Marcos_Sanches at gallup.com To: owner-r-help at stat.math.ethz.ch Subject: How to get residuasl with arima0? Date: Wed, 19 Sep 2001 15:19:07 -0300 I know this is a basic question, but I've never used the 'ts' package and I'm having some

Rob J Hyndman gives great explanation here (https://robjhyndman.com/hyndsight/estimation/) for reasons why results from R's arima may differ from other softwares. @iacobus, to cite one, 'Major discrepancies between R and Stata for ARIMA' (https://stackoverflow.com/questions/22443395/major-discrepancies-between-r-and-stata-for-arima), assign the, sometimes, big diferences from R

Full_Name: Arto Luoma Version: 1.1.0 OS: Windows 98 Submission from: (NULL) (153.1.53.119) In arima0 it is possible to specify whether the mean of the original series is included in the model or not. However, it is not possible to specify whether the mean of the differenced series is included. It seems that it is not included. However, if differencing is used to eliminate trend, the mean of the

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]

Good day, I'm using a VAR model to forecast sales with some extra variables (google trends data). I have divided my dataset into a trainingset (weekly sales + vars in 2006 and 2007) and a holdout set (2008). It is unclear to me how I should predict the out-of-sample data, because using the predict() function in the vars package seems to estimate my google trends vars as well. However, I want

Hi all, I want to fit the following model to my data: Y_t= a+bY_(t-1)+cY_(t-2) + Z_t +Z_(t-1) + Z_(t-2) + X_t + M_t i.e. it is an ARMA(2,2) with some additional regressors X and M. [Z_t's are the white noise variables] So, I run the following code: for (i in 1:rep) { index=sample(4,15,replace=T) final<-do.call(rbind,lapply(index,function(i)

Hello, I would like to study time series with ARIMA and GARCH models. I installed R-Plus and its libraries but when I try to execute the function arima0, It answers that the function does not exist. Could you help me or give me references of papers dealing with arima and garch in R-Plus? Thanks Beno?t, ___________________________________ Mr. Beno?t LACHERON Rue de l'industrie, 44, 1040

HI, there, When i use R, i tried to use function arima(), it complains: Error: couldn't find function "arima" But when I type "help.search("arima") ", I got arima() poped up.. arima(ts) ARIMA Modelling of Time Series arima.sim(ts) Simulate from an ARIMA Model arima0(ts) ARIMA Modelling of Time Series -- Preliminary

Hi, In my ongoing quest to track down the source of an error (see message "[R] optim error in arima" above), I find in the cource code for arima0 the following: arma0f <- function(p) { par <- as.double(fixed) par[mask] <- p .Call("arma0fa", G, par, PACKAGE = "stats") } I would like to know what the function

Dear R People: Suppose I want to estimate the parameters of the following AR model: (1 - phi_1 B - phi_2 B^2 - phi_9 B^9) x_t = a_t and I want to use the arima0 command from the ts library. How would I use the order subcommand, please? R Version 1.4.1 for Windows. Thanks! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston -

Hey, I have a data frame x which consists of say 10 vectors. I essentially want to find out the best fit exponential smoothing for each of the vectors. The problem while I'm getting results when i say > lapply(x,ets) I am getting an error when I say >> myprint function(x) { for(i in 1:length(x)) { ets(x[i],model="AZZ",opt.crit=c("amse")) } } The error message is

Hi, I would like to create variables in R containing parameters of summary(*Forecast Results*). Using the following code: library(forecast) data <- AirPassengers xets <- ets(data, model="ZZZ", damped=NULL) xfor <- forecast(xets,h=12, level=c(80,95)) summary(xfor) the output is: Forecast method: ETS(M,A,M) Model Information: ETS(M,A,M) Call: ets(y = data, model =