Me gustaría poder modificar los ejes de coordenadas en un plot que obtengo tras realizar un análisis estadístico (DCA). Este es un ejemplo (a pequeña escala) de los datos a analizar en el plot:
library(vegan)
library(reshape2)
Animal = sample(0:30,3800,replace = T)
Numero = sample(0:500,3800,replace = T)
df<-as.data.frame(cbind(Animal,Numero))
Matrix_df<-acast(df,Animal~Numero,fill = 0 )
Este es el análisis estadístico (DCA) donde ordena los datos dependiendo de los valores compartidos entre ambos ejes de la matriz:
DCA.df <- decorana(Matrix_df)
plot(DCA.df, origin = TRUE,
display = "sites",
type = "n",
cex.main = 0.75,
cex.axis=0.75,
cex.lab=0.75,
axes =FALSE,
xlim=c(-4,4),
ylim=c(-4,4),
main = "Raw data",
xlab = "DCA1 (0.353)",
ylab = "DCA2 (0.247)")
points (DCA.df)
La única manera factible de cambiarlos es usando este script
axis(1, at=c(-4,-2,0,2,4))
axis(2, at=c(-4,-2,0,2,4))
Pero cuando lo ejecuto obtengo un espacio muy grande entre el eje x e y
¿Cómo puedo juntarlo? ¿Cuál es mi fallo? ¿Existe otra manera más fácil de modificar los ejes dependiendo de mis necesidades?
Además me gustaría que el título del xlab y del ylab estuviera más junto al plot ¿Cómo podría juntarlo?