Displaying 20 results from an estimated 6000 matches similar to: "arima.sim"
2005 Oct 02
2
arima.sim bug?
Hi,
I am using the arima.sim function to generate some AR time series. However, the function does not seem to produce exactly the same time series when I specify the innov parameter. For example
> r <- rnorm(300)
> x <- arima.sim(300, model=list(order=c(1,0,0),ar=c(.96)), innov=r, n.start=10)
> y <- arima.sim(300, model=list(order=c(1,0,0),ar=c(.96)), innov=r, n.start=10)
>
2005 Oct 10
1
using innov in arima.sim
Hello,
I have used the arima.sim function to generate a lot of time series, but to day I got som results that I didn't quite understand. Generating two time series z0 and z1 as
eps <- rnorm(n, sd=0.03)
z0 <- arima.sim(list(ar=c(0.9)), n=n, innov=eps)
and
z1 <- arima.sim(list(ar=c(0.9)), n=n, sd=0.03),
I would expect z0 and z1 to be qualitatively similar. However, with n=10 the
2011 Nov 22
1
arima.sim: innov querry
Apologies for thickness - I'm sure that this operates as documented and with good reason. However...
My understanding of arima.sim() is obviously imperfect. In the example below I assume that x1 and x2 are similar white noise processes with a mean of 5 and a standard deviation of 1. I thought x3 should be an AR1 process but still have a mean of 5 and a sd of 1. Why does x3 have a mean of ~7?
2009 Jan 20
0
arima.sim help
I am trying to simulate time series data for an ar(1) and ma(1) process. I want the error term to have either a t distribution with 1 degree of freedom or a normal distribution with mean=0 and sd=1. Here is my code:
error.model=function(n){rnorm(n,mean=0, sd=1)}
data<-arima.sim(model=list(ar=c(0.1)), n=1000,
n.start=200, start.innov=rnorm(200,mean=0, sd=1),
rand.gen=error.model )
data
2009 Jul 21
0
Specifying initial values for arima.sim
Hi Everyone,
I'm having a problem with arima.sim. Namely specifying inital values
for the series.
If I generate a random walk
> vs = rnorm(100,0,1)
> xs = cumsum(vs)
and fit an ARIMA(1,0,0) to it
> xarima = arima(xs,order=c(1,0,0))
> xarima
Call:
arima(x = xs, order = c(1, 0, 0))
Coefficients:
ar1 intercept
0.9895 8.6341
s.e. 0.0106 6.1869
I should
2008 Jul 12
1
Help with arima.sim
Hey,
as a bloddy beginner in R I want to simulate a Arima (2,1,0) process with R.
My problem is, that I don't know how to specify the AR.
For a ARIMA(1,1,0) I use the following syntax: S <-
arima.sim(list(order=c(1,1,0), ar=0.9), n=100).
I think it is a stupid question with an easy answer. But when I google
the only results are arima.sim for an ARIMA (1,1,0).
Regards,
Chris
--
2004 May 24
1
Null model for arima.sim().
In some time series simulations I'm doing, I occasionally want the
model to be ``white noise'', i.e. no model at all. I thought it
would be nice if I could fit this into the arima.sim() context,
without making an exceptional case. I.e. one ***could*** do
something to the effect
if(length(model)==0) x <- rnorm(n) else x <- arima.sim(model,n)
but it would be more suave if one
2010 Aug 19
1
How to include trend (drift term) in arima.sim
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,
2009 Sep 29
0
Incoherence between arima.sim and auto.arima
Hello,
I have a question about function arima.sim
I tried to somulate a AR(1) process, with no innovation, no error term.
I used this code:
library(forecast)
e=rnorm(100,mean=0,sd=0)
series=arima.sim(model=list(ar=0.75),n=100,innov=e)+20
Then I tried to applicate ti this series auto.arima function:
mod1<-auto.arima(series,stepwise=FALSE,trace=TRUE,ic='aicc')
The best model returned
2003 Jul 16
1
arima.sim problems (PR#3495)
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
2000 Sep 22
2
arima.sim
Hi,
Before I re-invent the wheel, is there a function in R similar to S+'s
arima.sim, i.e., a function that simulates arima processes.
ts and tseries packages don't seem to have such function, but I may have
overlooked it.
Thank you for your time,
Alvaro Novo
R Version 1.1.1
SuSE 6.4 Linux
KDE 2.0 Beta 5
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
2005 Mar 31
2
how to simulate a time series
Dear useRs,
I want to simulate a time series (stationary; the distribution of
values is skewed to the right; quite a few ARMA absolute standardized
residuals above 2 - about 8% of them). Is this the right way to do it?
#--------------------------------
load("rdtb") #the time series
> summary(rdtb)
Min. 1st Qu. Median Mean 3rd Qu. Max.
-1.11800 -0.65010 -0.09091
2007 Jul 02
0
ARIMA prediction
Hi
This is my first post to this group, so apologies in advance if I get it wrong.
I would like to know how the prediction for arima models works in R. I
have a time series to which I fit an arima model, of varying AR and MA
orders. I then use the predict function to project it forward. I have
also written my own function to perform the prediction, but it gives
different answers to Arima.predict
2012 Jun 15
1
Replication of linear model/autoregressive model
Hi,
I would like to make a replication of 10 of a linear, first order
Autoregressive function, with respect to the replication of its innovation,
e. for example:
#where e is a random variables of innovation (from GEV distribution-that
explains the rgev)
#by using the arima.sim model from TSA package, I try to produce Y
replicates, with respect to every replicates of e,
#means for e[,1], I want
2009 Apr 05
0
Question about arima.sim()
Hi,
I tried to simulate an ARIMA model by using arima.sim(), say arima.sim(n=100,list(order=c(1,0,1),ar=0.6,ma=0.9,sd=1), but the acf and pacf of simulated data using acf() and pacf() are so much different from the theoritcal acf and pacf. For instance, in my case, ar=0.6 and ma=0.9, so the acf for all lags should be greater than 0 based on the theoritical calculation, but the acf of simulated
2007 Jan 01
0
arima.sim with a periodic model
Hi all.
I have a periodiv arma model and I want to simulate it. In S-plus, the
following works for me:
phi <- 0.9
theta <- 0
p <- 1 # period
model <- list(ar=phi, ma=theta, period=p)
Yt <- arima.sim(model, n=250)
How do I do something like this "period=12" in R? I read help(arima.sim)
but it doesn't tell.
Thanks! Philip.
2007 Jan 16
2
ARIMA xreg and factors
I am using arima to develop a time series regression model, I am using arima
b/c I have autocorrelated errors. Several of my independent variables are
categorical and I have coded them as factors . When I run ARIMA I don't
get any warning or error message, but I do not seem to get estimates for all
the levels of the factor. Can/how does ARIMA handle factors in xreg?
here is some example
2007 Aug 23
1
Estimate Intercept in ARIMA model
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))
2004 Mar 04
2
adding trend to an arima model
Hi,
Does anyone know a method for adding a linear/polynominal trend to a
simulated arima model using the arima.sim function?
Any help will be greatly appreciated.
Cheers,
Sam.
2003 Jan 09
2
using arima() function
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