alors je vais essayer d'être la plus claire possible...
pour ma matrice 1 appelée "trait" :
Code : Tout sélectionner
trait<-read.csv("trait_ok.csv", header=T,sep=";", row.names=1)
trait[trait>=1]<-1
trait<-as.matrix(trait)
me donnant une matrice (ce n'est qu'une partie):
a1 a2 a3 a4 a5 a6 a7 b1 b2 c1 c2 c3 e1 e2 e3 e4
E1 0 0 1 0 0 0 0 1 0 0 1 0 0 0 0 1
E2 0 0 1 0 0 0 0 1 0 0 1 1 0 0 0 1
E3 0 0 1 1 0 0 0 1 0 0 1 1 0 0 1 1
E4 0 0 0 1 1 0 0 1 0 0 1 0 0 0 0 1
ma matrice 2 appelée "m" :
Code : Tout sélectionner
m<-read.csv("macro_hab-esp.csv",header=T,row.names=1,sep=";")
m<-data.frame(t(m))
m<-as.data.frame(lapply(m,as.numeric))
m<-as.matrix(m)
me donnant une matrice :
E1 E2 E3 E4
[1,] 0 0 4 0
[2,] 0 0 0 0
[3,] 0 0 0 7
[4,] 0 0 0 4
[5,] 0 0 7 0
[6,] 0 1 0 0
[7,] 0 0 0 0
[8,] 0 0 0 0
[9,] 0 0 2 4
J'effectue ensuite ma multiplication :
me donnant :
a1 a2 a3 a4 a5 a6 a7 b1 b2 c1 c2 c3 e1 e2 e3 e4
[1,] NA NA NA NA NA NA NA NA NA NA NA NA NA NA
[2,] NA NA NA NA NA NA NA NA NA NA NA NA NA NA
[3,] NA NA NA NA NA NA NA NA NA NA NA NA NA NA
[4,] NA NA NA NA NA NA NA NA NA NA NA NA NA NA