nobs <- 1000
nu <- 8
a <- c(0.1,0.2,0.7) # a <- c(a constant, ARCH parameter, GARCH parameter)
# with normal innovations
eps <- uni.vola.sim(a, nobs)
# with t innovations
eps.t <- uni.vola.sim(a, nobs, d.f = df)Run the code above in your browser using DataLab