R：このループをどのように変更して適用できますか？

Question

私は現在、Rプログラムに取り組んでいます。このプログラムの一部は、相互依存性のあるループの2つの値を計算します。私は10万回の反復をしなければならないので、それは非常に長い時間がかかります。

私はこのループを適用ループまたはより効率的な関数の代わりに使用したいと思いますが、その方法はわかりません。誰か助けてくれますか？

p <- c() 
for(i in 1:n) { 
    if(i == 1) { 
    x <- b[i] 
    } 
    else { 
    x <- c(x, max(h[i - 1], p[i])) 
    } 
    h <- c(h, x[i] + y[i]) 
} 

ありがとうございます！

Answer 1

ここでは完全な実例はないようですが、xとhのベクトルを構築するには、c()の機能が非常に遅いという大きな問題があります。

x <- numeric(n) # allocate vector of size n 
h <- numeric(n) 

を、あなたがx[i]とh[i]に割り当てることで行くようにし、それらを記入：それはそれらを事前に割り当てる方が良いでしょう。たとえば、次のループ：

x <- c(); for (i in 1:100000) x <- c(x,1) 

は私のラップトップ上で実行するのに約10秒かかりますが、このバージョン：ほぼ瞬時に実行している間

x <- numeric(100000); for (i in 1:100000) x[i] <- 1 

は同じことを行います。