2017-02-27 33 views
1

我試圖根據數據幀中的一行中的觀察值計算新的值向量。我想從第一個值開始,乘以第二個值,然後將第三個值乘以剛剛計算的值。我不確定如何引用一個不存在的值,直到我計算出來。迭代地乘以R中的一行的值

這是什麼,如果我在寫這一點,手寫我會做:

index1 <- rnorm(12) 
index2 <- rnorm(12) 
index3 <- rnorm(12) 
index4 <- rnorm(12) 
index5 <- rnorm(12) 
data <- rbind(index1, index2, index3, index4, index5) 

splice1=data[2,1] 
splice2=splice1*data[2,2] 
splice3=splice2*data[2,3] 
splice4=splice3*data[2,4] 
splice5=splice4*data[2,5] 

splices <- c(splice1, splice2, splice3, splice4, splice5) 

我需要爲25行數據的做到這一點。

我發現此文章 Multiply previous row value by constant R 但我不知道如何使這項工作,當你通過從數據集索引引用。

我確定必須有一個簡單的方法來做到這一點,但我不知道從哪裏開始。

謝謝你的幫助。

回答

3

?apply在每行數據的?cumprod功能:

t(apply(data, 1, cumprod)) 

#    [,1]  [,2]  [,3]  [,4]   [,5] 
#index1 0.37891415 1.02182094 0.62823228 0.03616281 -0.0262988311 
#index2 -0.06586518 -0.03332061 -0.02814185 0.01755643 -0.0238307051 
#index3 -0.09113830 0.13385726 -0.11245386 0.10917885 -0.1466176651 
#index4 -0.53727635 0.71947970 0.87007225 0.19779160 -0.2119473091 
#index5 -0.83982393 -0.08571211 -0.08986865 0.02597184 -0.0006549512 

第2行這裏看起來很熟悉:

splices 
#  index2  index2  index2  index2  index2 
#-0.06586518 -0.03332061 -0.02814185 0.01755643 -0.02383071 

cumprod()本質上只是乘以沿各供貨組數字,一個cum ulative prod uct:

cumprod(1:5) 
#[1] 1 2 6 24 120 

如果速度是一個問題,你可以看一下matrixStats包:

library(matrixStats) 
rowCumprods(data) 
+0

非常感謝這一點。非常感激。 – Laura