有谁知道创建矩阵的快速方法
就像 R 中的以下内容。

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

上面的矩阵是 4x4,我想创建类似 10000x10000 的东西。

最佳答案

你可以做:

N <- 4
m <- matrix(nrow = N, ncol = N)
m[] <- pmin.int(col(m), row(m))

或@dickoa 建议的较短版本:
m <- outer(1:N, 1:N, pmin.int)

这些也有效并且都更快:
m <- pmin.int(matrix(1:N, nrow = N, byrow = TRUE),
              matrix(1:N, nrow = N, byrow = FALSE))

m <- matrix(pmin.int(rep(1:N, each = N), 1:N), nrow = N)

最后,这是一个使用矩阵乘积的可爱方法,但速度相当慢:
x <- matrix(1, N, N)
m <- lower.tri(x, diag = TRUE) %*% upper.tri(x, diag = TRUE)

请注意,R 的 10k×10k 矩阵看起来很大,我希望您不会耗尽内存。

关于r - 在 R 中创建某个矩阵,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/17554779/

10-16 23:45