similar to: arima vs arima0

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

Hello, I'm fitting a series of ARIMA models to a data set to compare fits. After taking the logs of the data and then differencing them to induce stationarity, I execute arima( y, order=c( p, 0, q ), seasonal=list( order=c( P, 0, Q ), period=7 ) ) for various values of p, q, P and Q. For one set of these values, I get Error in optim(init[mask], armafn, method = "BFGS", hessian

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

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

[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

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

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 -

does the arima/arima0 function use the state space form of the model equation even when fitting with the "CSS"-method? regards Christoph [[alternative HTML version deleted]]

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

Hi, I work on order estimation for autoregressive processes and after some inconsistencies cropped up I implemented the AIC criterion myself. Its results do not match the implementation in R and there are a few things I can not understand even after reading the source code of R. I used the function called arima0 (with empty "ma" coeficient vector), and I do not understand how some of

Here is the summary on predict when using an arima0 object: The arima0 object must be based on a time series vector. That is; x <- ts(xm1, frequency=12, start=c(1975,1)) x.ar <- arima0(x,order=c(1,1,1)) predict(x.ar,n.ahead=3) Thanks so much to Prof. Brian Ripley and David Brahm and other! Sincerely, Erin Hodgess

Full_Name: Ahmad Abu Hammour Version: rw0651 OS: windows 95 Submission from: (NULL) (63.23.128.44) Although I know that "ts package" is preliminary, I wanted to compare the results from R and SPSS. I ran ARIMA(2,1,2) in both softwares. I got NaN in standard errors of coefficients from R and real figures from SPSS. I changed "delta" in R to match that used by SPSS, I received

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 again With arima0 the problem was solved but what are the diferences between arima and arima0? I have another question. I fit the model to the data and I make some predictions. But I also want to calculate MAPE based in the last 3 observations available. Is it possible? Can I obtain the fitted values from the model? thanks~ luis -- SAPO ADSL.PT, apanhe j? o comboio da Banda Larga. Kit

Hi, I have a suggestion for the fonction arima and arima0. I think you should not call the constant an intercept because it creates confusion. It is not really an intercept but a mean. For an AR(1) the intercept mu should be defined as: X(t)=mu + phi X(t-1) + e(t) What you call intercept mu is rather defined as (X(t)-mu) = phi (X(t-1)-mu)) + e(t) which is not a common way to define an

Could someone please give a brief explanation, or pointer to an explanation, of the following error: > arima(ts.growth, order = c(1,0,0),include.mean=T) Error in arima(ts.growth, order = c(1, 0, 0), include.mean = T) : non-stationary AR part from CSS and why it does not arise with > arima0(ts.growth, order = c(1,0,0)) Many thanks ____________________________ Dr. Daniel P. Bebber

Dear All, R 1.9.1, Windows When copying and pasting a few lines from the 'predict.Arima' help, I get an error message: > data(lh) > predict(arima(lh, order = c(3,0,0)), n.ahead = 12) Error in eval(expr, envir, enclos) : Object "xreg" not found On the other hand, the following is OK: > data(lh) > predict(arima0(lh, order = c(3,0,0)), n.ahead = 12) $pred Time

Hello R people, I am need some help using the try() function. Currently I am running a loop which uses arima() for some values of p and q, which sometimes crashes. When it crashes, I want the program to just ignore it and move on to the next values to loop through. I currently have this, looping through a range of values for p and q: lo = try( arima1 <- arima( y, order=c( p, 0, q ) ) if(