3

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 のようなものを作成したいと考えています。

4

1 に答える 1

13

できるよ:

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 x 10k の行列は大きく見えることに注意してください。メモリ不足にならないように願っています。

于 2013-07-09T17:44:50.300 に答える