1

我有一个像这样的矩阵:

myarray=cov(matrix(rexp(200),50,10))

我想生成所有可能的列组合并计算每个组合的相关矩阵,如果可能的话,使用列号而不是名称。在第二步中,我想计算每个矩阵的行列式,所以也许有一种有效的方法来做到这一点。

4

1 回答 1

2

这是一种方法:

list.of.matrices <- apply(expand.grid(rep(list(c(FALSE, TRUE)), ncol(myarray))),
                          1, function(j)myarray[, j, drop = FALSE])

length(list.of.matrices)
# [1] 1024

然后执行以下操作:

result <- sapply(list.of.matrices, function_of_your_choice)

但请注意,det只能应用于方阵......请澄清。

于 2013-02-02T00:56:50.653 回答