スプリットマトリックスと再結合

Question

これは私の最初の投稿です。私の質問が愚かであれば、事前に謝罪してください。私はプログラミングが初めてです。

私はに152720行と2列のRを持っています。 それぞれ15272行の10個の別々の行列に分割したいと思います。

私はこれを試してみました：

> newmx <-split(as.data.frame(eBpvalues), rep(1:10, each = 15272))) 

> summary(newmx) 

    Length Class  Mode 
1 2  data.frame list  
2 2  data.frame list 
3 2  data.frame list 
4 2  data.frame list 
5 2  data.frame list 
6 2  data.frame list 
7 2  data.frame list 
8 2  data.frame list 
9 2  data.frame list 
10 2  data.frame list 

にはどうすればいいので、私は20列と15272行を持つ新しい行列を持つこれらの行列のサイド・バイ・サイドへの参加については行くだろうか？

乾杯、あなたがほとんど存在し

ニール

Answer 1

。これらの状況でよく使用される関数は、適用する関数とそれを適用するデータのリストを取るdo.callです。適用する関数は、cbind〜の列は、10データフレーム/行列を結合します。

文字通り、適切なサイズのmat（あなたのQにeBpvalues）の行列から始めます。データフレームに変換：

mat <- matrix(rnorm(152720 * 2), ncol = 2) 
df <- data.frame(mat) 

指標ファクターを製造する簡単な方法はgl()機能によってである：

newMat <- split(df, ind) 

最後：

ind <- gl(10, 15272) 

その後、我々はあなたのsplit()電話を持っていますステップはこれです。ここではdo.call()cbind()をnewMatのデータフレームのセットに適用する：

res <- do.call(cbind, newMat) 

これは、あなたが望むものを提供します（列名を整理する必要があるかもしれませんが）。インデックスを使用して

> str(res) 
'data.frame': 15272 obs. of 20 variables: 
$ 1.X1 : num -0.268 -0.8568 -0.0267 1.0894 1.5847 ... 
$ 1.X2 : num 0.71 -0.298 0.359 0.97 -2.158 ... 
$ 2.X1 : num -0.987 -0.222 2.991 0.443 0.228 ... 
$ 2.X2 : num -2.343 -1.023 -1.48 1.47 0.758 ... 
$ 3.X1 : num -0.305 -0.761 0.817 1.347 0.694 ... 
$ 3.X2 : num -0.0915 0.4816 1.4662 -1.2668 -1.3523 ... 
$ 4.X1 : num -0.678 -1.056 1.029 -0.468 0.836 ... 
$ 4.X2 : num -0.656 -0.459 -0.965 -1.666 0.877 ... 
$ 5.X1 : num -0.295 -1.255 1.395 -1.985 -1.71 ... 
$ 5.X2 : num 1.141 1.177 -1.003 -0.29 -0.234 ... 
$ 6.X1 : num -0.0548 1.8673 -1.5388 -1.1063 0.3923 ... 
$ 6.X2 : num -1.399 0.57 0.367 -0.811 -2.434 ... 
$ 7.X1 : num 0.389 -1.058 0.61 1.102 -0.063 ... 
$ 7.X2 : num 0.854 1.251 1.095 -0.485 0.451 ... 
$ 8.X1 : num -2.018 0.849 0.3 0.988 -1.993 ... 
$ 8.X2 : num -1.23 -1.025 -0.546 1.674 0.588 ... 
$ 9.X1 : num 0.814 0.726 1.04 0.985 1.781 ... 
$ 9.X2 : num -1.094 -1.051 0.749 1.426 0.402 ... 
$ 10.X1: num 0.3786 1.6131 -0.4149 0.0684 -0.815 ... 
$ 10.X2: num 0.383 -0.136 -0.751 -0.164 0.434 ... 

Answer 2

、あなたは簡単にRは、行列の列方向を満たすことを利用して、splitコマンドなしでこれを行うことができます。

おもちゃ例：

> eBpvalues <- matrix(1:100,ncol=2) # a matrix with 2 cols and 50 values 
> id <- c(1,11)+rep(0:9,each=2) 
> id 
[1] 1 11 2 12 3 13 4 14 5 15 6 16 7 17 8 18 9 19 10 20 
> tt <- matrix(eBpvalues,ncol=20)[,id] 
> tt # a matrix with the 2 cols split up and binded rowwise 
    [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] 
[1,] 1 51 6 56 11 61 16 66 21 71 26 76 31 81 36 86 41 91 46 96 
[2,] 2 52 7 57 12 62 17 67 22 72 27 77 32 82 37 87 42 92 47 97 
[3,] 3 53 8 58 13 63 18 68 23 73 28 78 33 83 38 88 43 93 48 98 
[4,] 4 54 9 59 14 64 19 69 24 74 29 79 34 84 39 89 44 94 49 99 
[5,] 5 55 10 60 15 65 20 70 25 75 30 80 35 85 40 90 45 95 50 100 

これは再びあなたに行列を与えます。分割を使用する場合は、Gavinのソリューションを使用できます。これにより、すべてが不要な副作用を伴う可能性のあるデータフレームに変換されます。開始のために大きな行列ではかなり遅くなります。