如何使用 R語(yǔ)言模擬時(shí)間序列模型的數(shù)據(jù)。

如題，我必須在不使用R包(package)的情況下，在R環(huán)境模擬整數(shù)屬性GARCH （INGARCH）模型。與原始GARCH不同的是，這個(gè)模型不是常態(tài)分布（normal distribution） 而是泊松分佈 (poisson distribution).

模型如下。
INGARCH (1,1)

X_t ~ Poisson(lamb_t)

lamb_t = alpha_0 + alpha_1(X_t-1) + beta_1(lamb_t-1)

我寫(xiě)的程序無(wú)法成功的模擬出該模式的數(shù)據(jù)。

set.seed(100)
a0 <- 2
a1 <- 0.4
b1 <- 1
mu <- 5
n <-1000
w <- rpois(n,mu)
X <- lamb <- rep(0, n)

for (i in 2:n) {

X[i] <- w[i]
lamb[i] <- a0 + a1 (X[i-1]) + b1 lamb[i-1]
}

這個(gè)問(wèn)題困擾了我很久，請(qǐng)各方高人指點(diǎn)，先此致謝。