Załóżmy, że mam macierz x
, która zawiera 10 wierszy i 2 kolumny. Chcę wygenerować nową macierz M
, która zawiera każdą unikalną parę wierszy od x
- czyli nową macierz z 55 wierszami i 4 kolumnami.Szybkie generowanie kartezjańskiego produktu macierzy
Np
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
}
Więc x
jest:
[,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
A potem M
ma 4 kolumny, dwa pierwsze są jednym rzędzie od x
i następny 2 to kolejny wiersz z x
:
> 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
Czy istnieje szybszy lub prostszy (lub oba) sposób robienia tego w R?
Nigdy nie wiedziałem o expand.grid(). Odpowiedź Dirka przynosi wszystko razem (jak zawsze ...) –