r - 是否可以对 R 中向量元素的顺序更新进行向量化？

Question

是否可以像下面这样矢量化代码？

length(x) <- 100;
x[1]      <- 1;
y         <- rnorm(100);

for(i in 2:100) {
    x[i] <- 2 * y[i] * x[i-1];
}

我很欣赏这是一个微不足道的例子，但它有助于说明这个想法。

我经常需要编写代码，其中向量中的第 i 个值取决于第 (i-1) 个值，如果可能，我想在不需要 for 的情况下编写此代码。循环，因为分析表明具有这种类型操作的函数是我代码中的主要瓶颈。

此操作是否可矢量化，因此我不需要使用 for()循环计算？

Answer 1

一般来说，如果你想要一个矢量化的解决方案，你需要解决 recurrence relation .