我正在尝试根据数据帧中连续的观察值计算一个新的值向量。我想从第一个值开始,将其乘以第二个值,然后将第三个值乘以我刚刚计算的值。我不确定如何在计算之前引用一个不存在的值。
如果我用手写的话,这就是我要做的事情:
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 但是,当您通过数据集索引引用时,我无法弄清楚如何使这项工作成功。
我确信必须有一个简单的方法来做到这一点,但我不知道从哪里开始。
谢谢你的帮助。
答案 0 :(得分: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)