矩阵的逐块和

Question

我想做与Blockwise sum of matrix elements 类似的事情，但使用任意给定的狭缝图案。例如，如果行被拆分为row_split = list(block_1=1, block_2=2:3, block_3=4)，列被拆分为col_split = list(block_1=1:3, block_2=4)，则输出为2x3 矩阵，其中每个元素都是按块求和。

有什么想法吗？

Answer 1

一个简单的 vapply 语句就可以解决问题：

vapply(row_split,function(x){
    a<-numeric(4)
    a[x]=1
    a
},numeric(4))

结果：

     block_1 block_2 block_3
[1,]       1       0       0
[2,]       0       1       0
[3,]       0       1       0
[4,]       0       0       1

您需要转置其中一个，但它应该可以解决问题。