每三行数据框相加,它们乘以新结果-R

时间:2019-06-12 02:30:21

标签: r matrix aggregate

我有一个数据框,我想每隔三行添加一次。然后,将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

高级问候。

3 个答案:

答案 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