Question

例如：我有一个矩阵列表，我想评估它们的差异，有点像三维差异。如果我有：

m1 <- matrix(1:4, ncol=2)
m2 <- matrix(5:8, ncol=2)
m3 <- matrix(9:12, ncol=2)
mat.list <- list(m1,m2,m3)

我想获得

mat.diff <- list(m2-m1, m3-m2)

我找到的解决方案如下：

mat.diff <- mapply(function (A,B) B-A, mat.list[-length(mat.list)], mat.list[-1])

有没有更好/内置的方法来做到这一点？

Answer 1

您可以使用lapply或其他循环方式执行此操作：

mat.diff <- lapply( tail( seq_along(mat.list), -1 ), 
                    function(i) mat.list[[i]] - mat.list[[ i-1 ]] )

Answer 2

您可以使用combn生成矩阵索引并在每个组合上应用函数。

       combn(1:length(l),2,FUN=function(x) 
            if(diff(x) == 1)         ## apply just for consecutive index 
                  l[[x[2]]]-l[[x[1]]], 
                  simplify = FALSE)  ## to get a list

使用@Arun数据，我得到：

[[1]]
     [,1] [,2]
[1,]    4    4
[2,]    4    4

[[2]]
NULL

[[3]]
     [,1] [,2]
[1,]    4    4
[2,]    4    4

如何将后续列表元素组合到R中的新列表中？

2 个答案: