Question

我有一个关于R的简单问题。

假设我有一个数组x，其长度是整数h的倍数。我想构造一个长度为（x）/ h的新数组，其元素是前h个元素的总和，然后是下一个h元素的总和等。

我可以使用错误的C相关编码来做到这一点，但我希望看到使用R进行此操作的最有效方法。

感谢您的帮助。

Answer 1

说出你的价值

x = 1:100
h = 10

您可以通过创建一个包含h行的矩阵并获取列总和来实现此目的：

colSums(matrix(x, nrow=h))
# [1]  55 155 255 355 455 555 655 755 855 955

这是有效的，因为矩阵按列填充值：

matrix(x, nrow=h)
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
 [1,]    1   11   21   31   41   51   61   71   81    91
 [2,]    2   12   22   32   42   52   62   72   82    92
 [3,]    3   13   23   33   43   53   63   73   83    93
 [4,]    4   14   24   34   44   54   64   74   84    94
 [5,]    5   15   25   35   45   55   65   75   85    95
 [6,]    6   16   26   36   46   56   66   76   86    96
 [7,]    7   17   27   37   47   57   67   77   87    97
 [8,]    8   18   28   38   48   58   68   78   88    98
 [9,]    9   19   29   39   49   59   69   79   89    99
[10,]   10   20   30   40   50   60   70   80   90   100

Answer 2

或者很难，

y <- unlist(sapply(0:((length(x)/h)-1) , function(k) sum(x[((k)*h + 1):((k+1)*h)] ) ) )

R中的数组操作

2 个答案: