Matrisi

Question

matrisinin kartezyen çarpımını üretmek için 10 satır ve 2 sütun içeren x matrisim var. x'dan her bir benzersiz satır çiftini içeren yeni bir matris M oluşturmak istiyorum - yani 55 satır ve 4 sütun içeren yeni bir matris.

 [,1] [,2] 
[1,] 1 11 
[2,] 2 12 
[3,] 3 13 
[4,] 4 14 
[5,] 5 15 
[6,] 6 16 
[7,] 7 17 
[8,] 8 18 
[9,] 9 19 
[10,] 10 20 

Sonra M 4 sütun, ilk iki bir x gelen satır ve sonraki 2 x gelen başka bir satır vardır şunlardır vardır:

Ör

x <- matrix (nrow=10, ncol=2, 1:20) 

M <- data.frame(matrix(ncol=4, nrow=55)) 
k <- 1 
for (i in 1:nrow(x)) 
for (j in i:nrow(x)) 
{ 
    M[k,] <- unlist(cbind (x[i,], x[j,])) 
    k <- k + 1 
} 

Yani, x geçerli:

> head(M,10) 
    X1 X2 X3 X4 
1 1 11 1 11 
2 1 11 2 12 
3 1 11 3 13 
4 1 11 4 14 
5 1 11 5 15 
6 1 11 6 16 
7 1 11 7 17 
8 1 11 8 18 
9 1 11 9 19 
10 1 11 10 20 

is R'de bunu daha hızlı mı yoksa daha basit (ya da her ikisi) mi?

Answer 1

expand.grid() fonksiyonu kullanışlı: Burada

iki sütun Kartezyen ürün olarak düşünmek ne

Şimdi

R> GG <- expand.grid(1:10,1:10) 
R> GG <- GG[GG[,1]>=GG[,2],]  # trim it to your 55 pairs 
R> dim(GG) 
[1] 55 2 
R> head(GG) 
    Var1 Var2 
1 1 1 
2 2 1 
3 3 1 
4 4 1 
5 5 1 
6 6 1 
R> 

Eğer 'n * var (n + 1)/2 'altkümeleri ve basitçe orijinal matrisinizi indeksleyebilirsiniz.

Answer 2

Yaptığınız şeyden çok emin değilim, bu yüzden size yardımcı olabilecek veya yardımcı olamayacak bir şeyler atacağım. Bunun için

expand.grid(x[,1],x[,2]) 

Answer 3

"İlişkiler" paketini de deneyebilirsiniz. Here is the vignette. Şöyle çalışması gerekir:

relation_table(x %><% x) 

Answer 4

Dirk en cevabını kullanma:

idx <- expand.grid(1:nrow(x), 1:nrow(x)) 
idx<-idx[idx[,1] >= idx[,2],] 
N <- cbind(x[idx[,2],], x[idx[,1],]) 

> all(M == N) 
[1] TRUE 

Teşekkür herkese!