我有一个数据框,我想每隔三行添加一次。然后,将cumprod应用于行,以便使用新的data.frame以及生成的新行。
最后,我将只有行数的三分之一。
下面您可以找到我编写的一些代码。我试图寻找对象的类,并复制适用于矢量而不是矩阵的代码。
XYZ<-read.xlsx2("XYZ.xlsx",1)
XYZ.CUT<-aggregate(XYZ~gl(nrow(XYZ)/3, 3), data.frame(XYZ), sum)
F.XYZ<-apply(t(XYZ.CUT+1),1,cumprod)
这就是我所拥有的:
X Y Z
-0,01% 0,32% 0,11%
-0,04% 0,01% 0,45%
-0,11% -0,06% 0,03%
0,03% -0,04% 0,45%
0,02% 0,04% 0,30%
-0,07% -0,11% 0,11%
-0,12% -0,13% 0,30%
-0,01% -0,07% 0,04%
-0,37% 0,08% 0,05%
首先我要
X Y Z
-0,16% 0,25% 0,59%
-0,02% -0,11% 0,86%
-0,50% -0,12% 0,39%
,然后在每个元素上加1:
X Y Z
(1-0,16%) (1+0,25%) (1+0,59%)
(1-0,02%) (1-0,11%) (1+0,86%)
(1-0,50%) (1-0,12%) (1+0,39%)
我想用行制作一个cumprod:
X Y Z
x1 y1 z1
x1*x2 y1*y2 z1*z2
x1*x2*x3 y1*y2*y3 z1*z2*z3
高级问候。
答案 0 :(得分:2)
我们可以使用tidyverse
。使用parse_number
中的readr
,从列(mutate_all
)中提取数字部分,并按gl
创建的索引进行分组,summarise
将所有列获取sum
library(tidyverse)
library(readr)
out <- XYZ %>%
mutate_all(parse_number) %>%
group_by(grp = as.integer(gl(n(), 3, n()))) %>%
summarise_all(sum)
然后,我们使用rowCumprods
中的matrixStats
来获取每一行的累积乘积
library(matrixStats)
rowCumprods(as.matrix(out[-1]) + 1)
# [,1] [,2] [,3]
#[1,] 0.84 1.0668 1.696212
#[2,] 0.98 0.8722 1.622292
#[3,] 0.50 0.4400 0.611600
XYZ <- structure(list(X = c("-0.01%", "-0.04%", "-0.11%", "0.03%", "0.02%",
"-0.07%", "-0.12%", "-0.01%", "-0.37%"), Y = c("0.32%", "0.01%",
"-0.06%", "-0.04%", "0.04%", "-0.11%", "-0.13%", "-0.07%", "0.08%"
), Z = c("0.11%", "0.45%", "0.03%", "0.45%", "0.30%", "0.11%",
"0.30%", "0.04%", "0.05%")), row.names = c(NA, -9L), class = "data.frame")
答案 1 :(得分:1)
我们可以仅使用基数R来通过以下方式进行操作:
#First remove the % symbol from the columns and convert the values to numeric
XYZ[] <- lapply(XYZ, function(x) as.numeric(sub("%", "", x)))
#Sum every 3 rows
XYZ.CUT <- aggregate(.~ gl(nrow(XYZ)/3, 3),XYZ, sum)[-1]
#Add 1 and take cumulative product for each row
t(apply(XYZ.CUT + 1, 1, cumprod))
#Or if you need it columnwise use
#apply(XYZ.CUT + 1, 2, cumprod)
# X Y Z
#[1,] 0.84 1.0668 1.696212
#[2,] 0.98 0.8722 1.622292
#[3,] 0.50 0.4400 0.611600
答案 2 :(得分:1)
使用rowsum
获取汇总和的另一个基本R版本(使用@Akrun的XYZ
数据集):
XYZ[] <- lapply(XYZ, sub, pat="%$", rep="")
XYZ[] <- lapply(XYZ, as.numeric)
out <- rowsum(XYZ, (seq_len(nrow(XYZ)) + 2) %/% 3) + 1
然后根据希望累积产品工作的方式选择自己的冒险方式:
## column-wise cumprod
out[] <- lapply(out, cumprod)
out
# X Y Z
#1 0.8400 1.270000 1.590000
#2 0.8232 1.130300 2.957400
#3 0.4116 0.994664 4.110786
## row-wise cumprod
out[] <- Reduce(`*`, out, accumulate=TRUE)
out
# X Y Z
#1 0.84 1.0668 1.696212
#2 0.98 0.8722 1.622292
#3 0.50 0.4400 0.611600