1

例如:我有一个矩阵列表,我想评估它们的差异,类似于 3-D 差异。所以如果我有:

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])

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

4

2 回答 2

1

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

mat.diff <- lapply( tail( seq_along(mat.list), -1 ), 
                    function(i) mat.list[[i]] - mat.list[[ i-1 ]] )
于 2013-02-06T17:46:16.950 回答
0

您可以使用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
于 2013-02-06T17:38:13.687 回答