6

私はベクトルを言ってvec = c(1,1)おり、それを列ごとに10回複製(cbind)したいので、のようなものを得ることができますmatrix(1, 10, 2)。このレプリケーションを実行できる機能はありvecますか? つまり、rep(vec, 10)?

ありがとう!

4

3 に答える 3

9
vec <- c(1,2)
rep(1,10) %*% t.default(vec)
      [,1] [,2]
 [1,]    1    2
 [2,]    1    2
 [3,]    1    2
 [4,]    1    2
 [5,]    1    2
 [6,]    1    2
 [7,]    1    2
 [8,]    1    2
 [9,]    1    2
[10,]    1    2

または@Joshuaが指摘するように:

tcrossprod(rep(1,10),vec)

いくつかのベンチマーク:

library(microbenchmark)

microbenchmark(rep(1,10) %*% t.default(vec),
               matrix(rep(vec, each=10), ncol=2),
               t.default(replicate(10, vec)),
               tcrossprod(rep(1,10),vec),times=1e5)

Unit: microseconds
                                   expr    min     lq  median      uq       max
1 matrix(rep(vec, each = 10), ncol = 2)  2.819  3.699  4.3970  5.3700  2132.240
2         rep(1, 10) %*% t.default(vec)  2.456  3.185  3.6750  5.5370  2121.746
3         t.default(replicate(10, vec)) 57.741 62.987 64.3740 65.9590 26654.678
4           tcrossprod(rep(1, 10), vec)  2.192  2.924  3.3745  5.2465  2145.709
于 2013-02-13T16:53:31.040 に答える
5

1 つのオプションは次のとおりです。

vec <- c(1,1)
t(replicate(10, vec))

> t(replicate(10, vec))
      [,1] [,2]
 [1,]    1    1
 [2,]    1    1
 [3,]    1    1
 [4,]    1    1
 [5,]    1    1
 [6,]    1    1
 [7,]    1    1
 [8,]    1    1
 [9,]    1    1
[10,]    1    1
于 2013-02-13T16:46:53.720 に答える
3

どうですか:

R> vec = c(1,2)
R> matrix(rep(vec, each=10), ncol=2)
     [,1] [,2]
 [1,]    1    2
 [2,]    1    2
 [3,]    1    2
 [4,]    1    2
 ....
于 2013-02-13T16:45:59.207 に答える