En primer lugar voy a preparar un ejemplo en base a la información que aportaste, te sugiero que cuando quieras compartir información uses las función dput()
por ejemplo dput(head(COM))
ya que esto genera un "script" para generar exactamente el objeto que quieres compartir, con su clase, con el tipo exacto para cada columna, con factores si los tiene, etc.
COM <- read.table(text="
Antant Apoimb Athpres Athspp Blenni Booboo Chelab
10 75 0 0 0 0 0
5 0 0 0 0 0 0
5 300 0 0 0 0 0
21 400 0 0 0 0 0
6 500 0 0 0 0 0
9 100 0 0 0 0 0 ", skip=1, header=TRUE)
TROPHIC <- read.table(text="
RN TG
Antant PL
Apoimb PL
Athpres PL
Athspp PL
Blenni CA
Booboo PL
Chelab DE
", skip=1, header=TRUE)
rownames(TROPHIC) <- TROPHIC[,1]
TROPHIC <- TROPHIC[2]
Por lo que finalmente tendríamos lo siguiente:
> COM
Antant Apoimb Athpres Athspp Blenni Booboo Chelab
1 10 75 0 0 0 0 0
2 5 0 0 0 0 0 0
3 5 300 0 0 0 0 0
4 21 400 0 0 0 0 0
5 6 500 0 0 0 0 0
6 9 100 0 0 0 0 0
> TROPHIC
RN TG
1 Antant PL
2 Apoimb PL
3 Athpres PL
4 Athspp PL
5 Blenni CA
6 Booboo PL
7 Chelab DE
Revisa por favor, que los datos se correspondan efectivamente con tus estructuras.
Solución:
Si en definitiva lo que buscas es llegar a un total por etiqueta, podemos hacer lo siguiente:
COM <- as.data.frame(cbind(Total=colSums(COM)))
FINAL <- as.data.frame(cbind(COM, TG=TROPHIC[match(rownames(COM), rownames(TROPHIC)),]))
aggregate(. ~ TG, FINAL, sum)
Resultado Final:
TG Total
1 CA 0
2 DE 0
3 PL 1431
Detalle:
En primer lugar, con COM <- as.data.frame(cbind(Total=colSums(COM)))
generamos un nuevo data.frame
con la suma para cada especie, algo así:
Total
Antant 56
Apoimb 1375
Athpres 0
Athspp 0
Blenni 0
Booboo 0
Chelab 0
Con as.data.frame(cbind(COM, TG=TROPHIC[match(rownames(COM), rownames(TROPHIC)),]))
le agregamos a cada especie la etiqueta que le corresponde, algo así:
Total TG
Antant 56 PL
Apoimb 1375 PL
Athpres 0 PL
Athspp 0 PL
Blenni 0 CA
Booboo 0 PL
Chelab 0 DE
Y por último, con aggregate(. ~ TG, FINAL, sum)
agrupamos por cada etiqueta y obtenemos la suma final:
TG Total
1 CA 0
2 DE 0
3 PL 1431
temp=merge(COM, TROPHIC, by = c("TG"), all.x=T, sort = FALSE) Error in fix.by(by.x, x) : 'by' must specify a uniquely valid column
COM es la matriz biológica, y TROPHIC la matriz de las etiquetasCOM
es undata.frame
dónde cada columna se corresponde con una fila deTROPHIC
. Esto efectivamente es así? Y en definitiva lo que buscas es tener las sumas de especies por cadaTG
?