Diciamo che ho una matrice x
che contiene 10 righe e 2 colonne. Voglio generare una nuova matrice M
che contenga ciascuna coppia di righe univoche da x
, ovvero una nuova matrice con 55 righe e 4 colonne.Genera rapidamente il prodotto cartesiano di una matrice
Es,
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
}
Quindi, 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 poi M
ha 4 colonne, i primi due sono una riga dalla x
e il successivo 2 sono un'altra riga da 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
È c'è un modo più veloce o più semplice (o entrambi) di farlo in R?
Mai saputo di expand.grid(). La risposta di Dirk mette tutto insieme (come sempre ...) –