Gere rapidamente o produto cartesiano de uma matriz
Digamos que eu tenha uma matrizx
que contém 10 linhas e 2 colunas. Eu quero gerar uma nova matrizM
que contém cada par único de linhas dex
- isto é, uma nova matriz com 55 linhas e 4 colunas.
Por exemplo.,
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
}
Assim,x
é:
[,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
E depoisM
tem 4 colunas, as duas primeiras são uma linha dex
e os próximos 2 são outra linha dex
:
> 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
Existe uma maneira mais rápida ou mais simples (ou ambas) de fazer isso em R?