ベクトル要素の再帰的操作

Question

私はベクトルaを持っていて、ループを使わずに各要素に再帰的にbを掛けたいと思っています。

a <- rep(0, 10) 
a[1] <- 1 
b <- 2 

# with a loop 
for (i in 2:length(a)) a[i] <- a[i-1] * b 

ループを使用せずにこれに対処する方法についてのヒントに感謝します。

Answer 1

一般に、明示的なループなしでこれを行うことはできません。非常に単純に

a <- rep(2, 10) 
a[1] <- 1 
cumprod(a) 
# [1] 1 2 4 8 16 32 64 128 256 512 

Answer 2

指数関数^がベクトル化されているので、：：この特定のケースでは、暗黙のcumprodが提供するループを使用することができますあなたも

を書きたいかもしれません

2^(0:9) 
# [1] 1 2 4 8 16 32 64 128 256 512 

2^seq(from=0, to=9) 

長いベクトルの場合、@ JoshuaUlrichの方法ははるかに高速ですが、確かに非常にコンパクトです。あなたは、スピードについて特に心配していないとも述べました。

y[i] = x[i] + f[1]*y[i-1] + ... + f[p]*y[i-p] 

あなたがfilter機能を使用することができます。フォームの一般的な再帰的なシリーズで

Answer 3

。あなたの場合、x[i] = 0,f[1] = 2およびf[i] = 0はi > 1です。あなたはそれを使用する方法を学習した後

filter(rep(0,10), 2, method="recursive", init=1/2) 
# Time Series: 
# Start = 1 
# End = 10 
# Frequency = 1 
# [1] 1 2 4 8 16 32 64 128 256 512 

、常に初めて明らかにされていない、filterは非常に強力かつ効率的である：これはに変換されます。しかし、おそらくあなたの幾何学的な事例には残念です。