4

私は2つのことを成し遂げようとしています。まず、ベクトルがある場合1:5、同じ数の2倍を含み、繰り返しを除く、これらの要素の一意の組み合わせを示す行列（または2つのベクトル）を取得します。

今、私はマトリックスを使用してこれを行うことができます：

foo <- matrix(1:5,5,5)
cbind(foo[upper.tri(foo,diag=TRUE)],foo[lower.tri(foo,diag=TRUE)])
      [,1] [,2]
 [1,]    1    1
 [2,]    1    2
 [3,]    2    3
 [4,]    1    4
 [5,]    2    5
 [6,]    3    2
 [7,]    1    3
 [8,]    2    4
 [9,]    3    5
[10,]    4    3
[11,]    1    4
[12,]    2    5
[13,]    3    4
[14,]    4    5
[15,]    5    5

しかし、もっと簡単な方法が必要です。で使用しようとしVectorizeましseqたが、エラーが発生します：

cbind(Vectorize(seq,"from")(1:5,5),Vectorize(seq,"to")(5,1:5))
    Error in Vectorize(seq, "from") : 
      must specify formal argument names to vectorize

私がしたい2番目のことは、ベクトルを含むリストがある場合bar、その要素の要素の数に等しい繰り返しリストの要素を含むベクトルを取得することです。私はこれを行うことができます：

unlist(apply(rbind(1:length(bar),sapply(bar,length)),2,function(x)rep(x[1],x[2])))
 [1] 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3

しかし、やはりもっと簡単な方法があるはずです。ここで再試行Vectorizeしましたが、同じエラーが発生しました：

Vectorize(rep,"each")(1:length(bar),each=sapply(bar,length))
 in Vectorize(rep, "each") : 
  must specify formal argument names to vectorize
4

3 に答える 3

5 
> unlist(lapply(1:5, seq, from=1))
 [1] 1 1 2 1 2 3 1 2 3 4 1 2 3 4 5
> unlist(lapply(1:5, seq, 5))
 [1] 1 2 3 4 5 2 3 4 5 3 4 5 4 5 5

と

> bar = lapply(1:5, seq, from=1)
> rep(seq_along(bar), sapply(bar, length))
 [1] 1 2 2 3 3 3 4 4 4 4 5 5 5 5 5
于 2011-02-20T23:40:19.023 に答える
5

あなたの最初の質問に：combn()ベースの単純な関数はどうですか？

> combn(1:5,2)
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]    1    1    1    1    2    2    2    3    3     4
[2,]    2    3    4    5    3    4    5    4    5     5

作成したマトリックスを配置する必要がある場合は、次のt()ように転置するだけです。t(combn(1:5,2))

注：これでは、seqの繰り返し要素の組み合わせは返されませんが、それらをマトリックスに簡単に追加できます。

于 2011-02-20T23:48:22.307 に答える
4

最初の部分に対するMartinMorganのソリューションのより高速なバリエーション：

rep(1:5,5:1)
 [1] 1 1 1 1 1 2 2 2 2 3 3 3 4 4 5
unlist(lapply(1:5,function(x) x:5))
 [1] 1 2 3 4 5 2 3 4 5 3 4 5 4 5 5

それぞれ約7倍と3倍速くなります。

第二部であなたが何を意味するのか私が従うかどうかはわかりませんが、以下はあなたの説明に合うようです：

lapply(bar,function(x) rep(x,length(x)))
于 2011-02-21T12:56:24.903 に答える