У меня есть довольно простой набор данных, содержащий расположение (DCA) релевантных данных о растительности и сбалансированный набор данных с переменными среды.
Растительные релевантности структурированы по годам и блокам.Есть три раунда лет (2001, 2008 и 2018), и каждый год имеет 8 блоков (A, B, C, D, E, F, I и J).В настоящее время я использую этот код для построения группы:
plot(ord.dataveg, type="n", main="DCA of vegetation plots sorted block", xlab="DCA1 Eigenvalue = 0.45",
ylab="DCA2 Eigenvalue = 0.44")
points(ord.dataveg, display = "sites", cex = 0.4, pch = 19,
col = Envir.vars$Block)
В результате на следующей диаграмме: Расположение DCA растительных участков
Я хочу, чтобы все блоки имели разныецвета, но разные ПК для каждого года.Я не могу понять, как использовать базовый синтаксис R, чтобы сделать это.
Ниже приведен минимальный пример моих данных.
dataveg <- structure(list(ï.. = 1:10, DCA1 = c(0.48453, 0.791106, 0.688445,
0.364132, 0.21456, 1.3365, 0.9256, 2.3301, 1.3, 1.8), DCA2 = c(-0.4399,
-1.810391, -2.018287, -1.2409, -1.2365, 0.2564, 1.216, -2.31564,
-1.3, -1)), row.names = c(NA, 10L), class = "data.frame", na.action =
structure(11:19, .Names = c("11",
"12", "13", "14", "15", "16", "17", "18", "19"), class = "omit"))
Envir.vars<- structure(list(ï.. = 1:10, Year = c(2001L, 2001L, 2001L,
2008L,
2008L, 2008L, 2018L, 2018L, 2018L, 2018L), Block = structure(c(1L,
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 4L), .Label = c("A", "B", "C",
"D"), class = "factor")), class = "data.frame", row.names = c(NA,
-10L))
Надеюсь, этого достаточно.