迭代编程

Question

我想写一个关于时间迭代方差的代码。该函数应如下所示：

西格马（T）=阿尔法 X（t-1）+测试西格马（T-1）

不幸的是我无法弄清楚如何引入这个时间组件，例如在我的for循环中。有谁知道如何处理这类问题？

Answer 1

你的意思是??

alpha <- 2; beta <- .1
x     <- 1:5
sigma <- numeric(length(x))
sigma[1] <- 1
for (t in 2:length(x)) sigma[t] <- alpha*x[t-1] + beta*sigma[t-1]
sigma
# [1] 1.0000 2.1000 4.2100 6.4210 8.6421