similar to: bug in arima? (PR#9404)

Full_Name: Allan McRae Version: 1.6.0 OS: Win 2000 P Submission from: (NULL) (129.215.190.229) When using predict.Arima in library ts(), it appears differencing is only accounted for in the first step of prediction and so any trend is not apparent in the predictions. The example shows the difference between the predictions of an arima(1,1,1) model and the backtransformed predictions of an

There is currently a bug in the arima function. Namely, for arima models with differencing or seasonal differencing, the innovation variance estimator uses the wrong denominator whenever xreg is non-null. This is the case, for example, when fitting an ARIMA(p,1,q) model with a drift term (common in financial applications). I reported the bug (and a fix) at

On 21 May 2015, at 10:35 , Martin Maechler <maechler at lynne.stat.math.ethz.ch> wrote: >> >> I noticed that the 3.2.1 release cycle is about to start. Is there any >> chance that this fix will make it into the next version of R? >> >> This bug is fairly serious: getting the wrong variance estimate leads to >> the wrong log-likelihood and the wrong

I noticed that the 3.2.1 release cycle is about to start. Is there any chance that this fix will make it into the next version of R? This bug is fairly serious: getting the wrong variance estimate leads to the wrong log-likelihood and the wrong AIC, BIC etc, which can and does lead to suboptimal model selection. If it's not fixed, this issue will affect every student taking our time series

On 21 May 2015, at 12:49 , Martin Maechler <maechler at lynne.stat.math.ethz.ch> wrote: >>>>>> peter dalgaard <pdalgd at gmail.com> >>>>>> on Thu, 21 May 2015 11:03:05 +0200 writes: > >> On 21 May 2015, at 10:35 , Martin Maechler <maechler at lynne.stat.math.ethz.ch> wrote: > >>>> >>>> I noticed that

Hi, I have been using this website ( http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm ) to help me to fit ARIMA models to my data. At the moment I have two possible methods to use. Method 1 If I use arima(ts.data, order=c(1,2,0), xreg=1:length(ts.data)) then the wrong value for the intercept/mean is given (checked on SPSS and Minitab) and

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, All, This is my program ts1.sim <- arima.sim(list(order = c(1,1,0), ar = c(0.7)), n = 200) ts2.sim <- arima.sim(list(order = c(1,1,0), ar = c(0.5)), n = 200) tdata<-ts(c(ts1.sim[-1],ts2.sim[-1])) tre<-c(rep(0,200),rep(1,200)) gender<-rbinom(400,1,.5) x<-matrix(0,2,400) x[1,]<-tre x[2,]<-gender fit <- arima(tdata, c(1, 1, 0), method = "CSS",xreg=t(x))

I have the following model: m1.dis=arima(diff(diff(log(ts1),lag=12)),order=c(0,1,1),seasonal=list(order=c(0,1,1),period=12)) I would like to know how to plot the correct predictions in the original units because I am trying the following code but it is not working. I believe that there must be something to account for the differencing.

Dear R-helpers, I would appreicate if someone can help me on the transfer parameter in ARIMAX and also see what I am doing is correct. I am using ARIMAX with 2 Exogeneous Variables and 10 years data are as follows: DepVar Period, depVar, IndepVar1 Period, indepVar1, IndepVar2 Period, indepVar2 Jan 1998,708,Jan 1998,495,Jan 1998,245.490 Feb 1998,670,Feb 1998,421.25,Feb 1998,288.170 Mar

Hi, My query is related to ARIMA function in stats package. While looking for the time series literature I found following link which highlights discrepancy in "arima" function while dealing with differenced time series. Is there a substitute function similar to "sarima" mentioned in the following website implemened in R? Any pointers would be of great help.

Full_Name: Gang Liang Version: 1.7.1 OS: Debian/Woody Submission from: (NULL) (192.6.19.190) > print(arima.sim(list(ar=.3,order=c(1,1,1)), 30)) [1] 0.00000000 0.10734243 0.02907301 -1.23441659 -0.98819317 -2.82731975 [7] -2.69052512 -4.22884756 -5.02820635 -5.41514613 -6.20486350 -7.01040649 [13] -6.78121289 -5.41111810 -4.96338053 -5.42395408 -6.22741444 -5.75228153 [19] -6.07346580

Hi all, I have two time series that I would like to correlate but as they are autocorrelated, I am "pre-whitening" them first by fitting ARIMA models, then correlating their residuals....as described in https://onlinecourses.science.psu.edu/stat510/?q=node/75 However, http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm discusses some issues with ARIMA in R. In particular, for issue 2, if

Hi I'm using fracdiff package to estimate the parameters of a fractionally-differenced ARIMA (p,d,q) model, and it works fine, but I wanted to have also the filtered series and the series of residuals. I understand these are calculated in the subroutine fdfilt, in the program fdcore.f, but I can't manage to get them out. Any suggestion would be much appreciated Thanks Susana Barbosa

Hi, I have been looking at arima.sim to simulate the output from an ARMA model fed with a normal and uncorrelated input series but I cannot find a way to pass an intercept / constant into the model. In other words, the model input in the function allows only for the AR and MA components but I need to pass a constant. Can anyone help? Thanks Paolo [[alternative HTML version deleted]]

Hello everyone, I was trying to fit an arima model to a daily historical data, but, for some reason, havent been able to. I basically have 212 observations (from 12/1/2012 to 06/30/2013) containing the number of transits for a particular vessel. The following messages are produced by R: dailytrans.fit<-arima(dailytrans$transits, order=c(0,1,2), seasonal=list(order=c(0,1,2), period=365),

The bug repository is like an elephant: It doesn't forget, but the gestation period is long. In the present case, it is clear that something is not right, but someone needs to have sufficient recall and insight to check that your proposed fix is not unfixing a deliberate change. We should get to it eventually. (For some value of "we" not including "me"...) -pd On 20 Apr

>>>>> peter dalgaard <pdalgd at gmail.com> >>>>> on Thu, 21 May 2015 11:03:05 +0200 writes: > On 21 May 2015, at 10:35 , Martin Maechler <maechler at lynne.stat.math.ethz.ch> wrote: >>> >>> I noticed that the 3.2.1 release cycle is about to start. Is there any >>> chance that this fix will make it into

I have been trying to simulate from a time series with trend but I don't see how to include the trend in the arima.sim() call. The following code illustrates the problem: # Begin demonstration program x <- c(0.168766559, 0.186874000, 0.156710548, 0.151809531, 0.144638812, 0.142106888, 0.140961714, 0.134054659, 0.138722419, 0.134037018, 0.122829846, 0.120188714,

Good morning to you all, Sorry for taking your time from your research and teaching schedules.   If you have a non-stationary univariate time Series data that has the transformation: Say; l.dat<-log (series) d.ldat<-diff (l.dat, differences=1) and you fit say arima model. predit.arima<-predict (fit.series, n.ahead=10, xregnew= (n+1) :( n+10)) How could I re-transform