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]]