假设我有一个矩阵(可以是可转换为矩阵的数据框,反之亦然)。
获取列总和很容易(使用colSums、apply 或者其他喜欢的方法)。但是,我有一个位置索引,比方说。
X=matrix(rnorm(45), nr=15 )
[,1] [,2] [,3]
[1,] -0.35659671 -0.5014418 -0.90051481
[2,] -2.84945987 -0.6558198 -1.59163998
[3,] -0.03612257 0.7118217 -0.76520379
[4,] 2.45922337 -1.4965018 -0.28132609
[5,] -0.46939481 -2.2576689 -0.91045703
[6,] 0.65121124 -0.4879071 -1.81893484
[7,] -0.28633954 0.3570085 -2.40230043
[8,] -1.32283279 -0.4904968 0.42214466
[9,] 0.55199355 -1.5181311 -0.62465292
[10,] -1.00114458 0.5140162 -0.01261708
[11,] -1.49307374 0.1734213 1.38231439
[12,] -1.71157588 -0.7030724 -1.05475532
[13,] 1.13823870 1.3380262 2.74988880
[14,] -0.41600520 -0.9527629 0.68447386
[15,] -1.25893418 -0.3633260 0.38870010
获取列总和很容易(使用colSums、apply 或者其他喜欢的方法)。但是,我有一个位置索引,比方说。
idx=c(3,7,13)
所以我需要将矩阵分割并对子部分的列进行求和(或者一般地执行某些函数),例如对于行1到3、4到7、8到13和14到15(最后一行),并获得4乘3矩阵。
[,1] [,2] [,3]
[1,] -3.2421792 -0.4454399 -3.2573586
[2,] 2.3547000 -3.8850690 -5.4130180
[3,] -3.8383947 -0.6862365 2.8623225
[4,] -1.6749390 -1.3160890 1.0731740
这是一个最小可重现的示例。实际矩阵要大得多,子部分的数量和位置是可变的,并且函数更加复杂 :)
Map
的方法提供了一个好的一行代码,但是@akrun的第一个解决方案更快。 - Colonel Beauvel