Me gustaría saber como poder hacer rotar la imagen del DCA en el script de R
Para que este DCA quedara en la misma orientación que este otro DCA
Para hacer ambos DCA he seguido el mismo script pero diferentes bases de datos:
Load packages.
library(readr)
library(vegan)
library(reshape2)
Load Raw data.
colnames(ngs.data)[c(1)] <- c("habitat") ; head(ngs.data); dim(ngs.data)
ngs.data.coll.dca <- acast(ngs.data, collection+habitat ~ taxon , fill = 0); head(ngs.data.coll.dca); dim(ngs.data.coll.dca)
rownames(ngs.data.coll.dca)
Create plot parameters from rowname attributes.
raw.site.hab <- sapply(strsplit(rownames(ngs.data.coll.dca), "_"), "[[", 2); head(raw.site.hab)
site.des <- c("Wl1F", "Wl2F", "Wl4F", "Wl4P", "Ed1F", "Ed1P", "Ed1V", "Ed2F", "Ed2P",
"Ed2V", "Ed3F", "Ed3P", "Ed3V", "Ed4F", "Ed4P", "M1V", "M2V",
"M3V", "M4V", "M5V", "M7V", "M8V", "Oa1F", "Oa1P", "Oa2F", "Oa2P",
"Oa3F", "Oa4F", "Oa4P", "Z1V", "Z2V")
length(site.des)
Perform DCA
coll.dca <- decorana(ngs.data.coll.dca)
Variance contribtutions of axes.
coll.dca$evals/sum(coll.dca$evals)
plot(coll.dca,
type = "n",
cex.main = 0.75,
cex.lab = .75,
axes = TRUE,
cex.axis = 0.5,
yaxt = "n",
xlim = c(-8, 18),
ylim = c(-4, 10),
#main = "Raw data",
xlab = "DCA1 (0.334 total variance)",
ylab = "DCA2 (0.229 total variance)")
axis(2, at=c(-4, -2, 0, 2, 4, 6, 8, 10), tick = TRUE, cex.axis = 0.5)
points(coll.dca, col = as.integer(as.factor(raw.site.hab)),
pch = as.integer(as.factor(raw.site.hab)))
toda sugerencia es bien recibida, muchas gracias de antemano.