Como puedo modificar mi codigo de estos gráficos triangulares para poder cambiarles los colores para poder obtener el gráfico en mejores condiciones, donde dice habitat es una data cualquiera de filas x columnas:
library (adespatial)
library (ade4)
library (betapart)
if (require(ade4, quietly = TRUE))
data(habitat)
habitat.pod.J <- beta.div.comp(habitat, coef = "J", quant = FALSE)
# Sorensen
habitat.pod.S <- beta.div.comp(habitat, coef = "S", quant = FALSE)
# Ruzicka
habitat.pod.qJ <- beta.div.comp(habitat, coef = "J", quant = TRUE)
# Percentage difference
habitat.pod.qS <- beta.div.comp(habitat, coef = "S", quant = TRUE)
# Data frames for the triangular plots
# Data frames for the triangular plots
habitat.pod.J.3 <- cbind((1-habitat.pod.J$D),
habitat.pod.J$repl,
habitat.pod.J$rich)
colnames(habitat.pod.J.3) <- c("Similarity", "Repl", "RichDiff")
habitat.pod.S.3 <- cbind((1-habitat.pod.S$D),
habitat.pod.S$repl,
habitat.pod.S$rich)
colnames(habitat.pod.S.3) <- c("Similarity", "Repl", "RichDiff")
habitat.pod.qJ.3 <- cbind((1-habitat.pod.qJ$D),
habitat.pod.qJ$repl,
habitat.pod.qJ$rich)
colnames(habitat.pod.qJ.3) <- c("Similarity", "Repl", "AbDiff")
habitat.pod.qS.3 <- cbind((1-habitat.pod.qS$D),
habitat.pod.qS$repl,
habitat.pod.qS$rich)
colnames(habitat.pod.qS.3) <- c("Similarity", "Repl", "AbDiff")
par(mfrow = c(2, 2))
triangle.plot(as.data.frame(habitat.pod.J.3[, c(3, 1, 2)]),
show = FALSE,
labeltriangle = FALSE,
addmean = TRUE
)
text(-0.45, 0.5, "RichDiff", cex = 1.5)
text(0.4, 0.5, "Repl", cex = 1.5)
text(0, -0.6, "Jaccard similarity", cex = 1.5)
triangle.plot(as.data.frame(habitat.pod.S.3[, c(3 ,1 ,2)])
Por favor ayudenme como puedo modificarle los colores para que no salga negro y si es posible hacer que salgan con la letra Time New Romann