0

我是初学者 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), ] 中的错误：索引超出范围”。我该如何解决？有什么建议吗？

谢谢。

4

2 回答 2

2

该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

rowSums()如果您实际上想要将列聚合为 4 行中的每一行的单个值，您也可以使用结果...

rowSums( tmp )
#[1] 78 99
于 2013-10-26T13:48:31.283 回答
1

这是一种方法：

## 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
于 2013-10-26T13:44:35.673 回答