Question

我是初学R用户，我需要编写一个函数，该函数在固定的时间间隔（每4行）对数据帧的行进行求和。我试过以下代码

camp<-function(X){
  i<-1
  n<-nrow(X)
  xc<-matrix(nrow=36,ncol=m)
  for (i in 1:n){
    xc<-apply(X[i:(i+4),],2,sum)
    rownames(xc[i])<-rownames(X[i])
    i<-i+5
  }
  return(xc)
}

结果是“X中的错误[i：（i + 4），]：索引超出范围”。我怎么解决？有什么建议吗？

感谢。

Answer 1

zoo包有rollapply，这对于像这样的东西非常方便......

#  Make some data
set.seed(1)
m <- matrix( sample( 10 , 32 , repl = TRUE ) , 8 )
#     [,1] [,2] [,3] [,4]
#[1,]    3    7    8    3
#[2,]    4    1   10    4
#[3,]    6    3    4    1
#[4,]   10    2    8    4
#[5,]    3    7   10    9
#[6,]    9    4    3    4
#[7,]   10    8    7    5
#[8,]    7    5    2    6

#  Sum every 4 rows
require( zoo )
tmp <- rollapply( m , width = 4 , by = 4 , align = "left" , FUN = sum )
#     [,1] [,2] [,3] [,4]
#[1,]   23   13   30   12
#[2,]   29   24   22   24

如果您确实希望将列聚合为4行中的每一行的单个值，也可以对结果使用rowSums() ...

rowSums( tmp )
#[1] 78 99

Answer 2

这是一种方法：

## Sample data
m <- matrix(1:36, nrow=12)
## Create a "group" index
fac <- (seq_len(nrow(m))-1) %/% 4
## Apply sum 
apply(m, 2, function(v) tapply(v, fac, sum))

示例数据：

      [,1] [,2] [,3]
 [1,]    1   13   25
 [2,]    2   14   26
 [3,]    3   15   27
 [4,]    4   16   28
 [5,]    5   17   29
 [6,]    6   18   30
 [7,]    7   19   31
 [8,]    8   20   32
 [9,]    9   21   33
[10,]   10   22   34
[11,]   11   23   35
[12,]   12   24   36

结果：

  [,1] [,2] [,3]
0   10   58  106
1   26   74  122
2   42   90  138

R：对行间隔求和

2 个答案: