Quiero hacer un plot de mis muestras y asignar un símbolo distinto acorde con el tipo de especie que sea. Los paquetes que he usado son: library(vegan) library(reshape2) Los script que he realizado para la gráfica han sido: partiendo de un extracto del dataframe que tiene
Datraframe
Site Date Habitat Season Year Taxa
Q1F 08_09_2015 Oak Autumn 2015-2016 Adonis_flammea
E2F 08_09_2015 Oak Autumn 2015-2016 Agrimonia_eupatoria
Q4F 08_09_2015 Oak Autumn 2015-2016 Ajuga_chamaepitys
L1P 08_09_2015 Oak Autumn 2015-2016 Ajuga
la convierto en una matriz:
Matrix_Rawdata<-acast(Datos,Site~Taxa,fill = 0 )
Matrix_Rawdata<-data.frame(Matrix_Rawdata)
El dataframe Matrix_Rawdata es el siguiente: dput(head(Matrix_Rawdata))
Adonis_flammea Agrimonia_eupatoria Ajuga_chamaepitys Ajuga
Q1F 1 0 0 0
E2F 0 1 0 0
Q4F 0 0 1 0
L1P 0 0 0 1
y otro dataframe llamado Site, que tiene:
Plot Group
Q1F 1
E2F 2
Q4F 1
L1P 3
intento plotearlo usando estos scripts:
DCA.Rawdata <- decorana(Matrix_Rawdata)
plot(DCA.Rawdata,
choices = c(1, 2),
origin = TRUE,
display = "sites",
type = "n",
cex = 0.5,
main = "Raw Data",
xlab = "DCA1 (0.353)",
ylab = "DCA2 (0.247)")
points (DCA.Rawdata, col = Site$Group, pch = Site$Group)
pero la parte de pch del script points me da este error:
Error in plot.xy(xy.coords(x, y), type = type, ...) : símbolo de gráfico inválido
¿Podrían decirme dónde está mi fallo?
ggplot2
es una tarea bastante sencilladata.frame
para poder fabricar un ejemplo con tus datos. Para ello puedes usar el comandodput(head(Dataframe,20))
vegan
, por lo que me imagino que el tema pasa por los datos, trata de armar un dataframe más pequeño pero que produzca el error y publicalo tal como te indico Aldo,dput(head(Datos,20))
, esto no copia todo el dataframe solo las primeras 20 filas. Por otro lado, me llama la atención que el mensaje de error este en castellano, lo has traducido tu?