Dear Users,

I follow Andreas idea to simulate an ar(1) model with a new kind of innovation process.
The new argument rand.gen, for the arima.sim function, I'm trying to generate as:

tGarchGen <- function(a, b, c) {
    # must return a vector of random deviates (eta(t))
  for (t in 1:100){
  z(t) <- c+a*(eta(t)^2)+b*z(t-1)
  eta(t) <-rt(100, 5)*sqrt(z(t))  #rt is the R random t-Student generator function rt(n,df)
  }
}

arModel <- list(ar=0, ma = 0, order = c(0, 1, 0))
arima.sim(arModel, n = 100, rand.gen = tGarchGen)

But, since I'm a newbie in R, course the loop doesn't work.
Finally, how I do to declare the parameters (a,b,c,n,df), including the loop's length in a way I can change it at the beginning of the program?
Any help, please?

Thanks in advance,

Rick
	[[alternative HTML version deleted]]


