У меня есть данные с обработками, сделанными за два разных года (столбец "Treat.year") и некоторыми переменными среды.Я хочу создать диаграмму рассеяния для PCA, но не могу понять, как изменить цвета точек.
Я хочу, чтобы точки в данной обработке (C или Ca) были одного цвета воба периода выборки (например, все точки C синие и все точки Ca, зеленые независимо от года) и «кадры» для соединения точек в соответствии с годами обработки И (т. е. кадр, соединяющий только точки C.2009, другойтолько баллы C.2015 и т. д.).
ex <- data.frame(treat.year = c("C.2015", "C.2015", "C.2015", "C.2015", "Ca.2015", "Ca.2015", "Ca.2015", "Ca.2015", "C.2009", "C.2009", "C.2009", "C.2009", "Ca.2009", "Ca.2009", "Ca.2009", "Ca.2009"),
var.b = c(33, 18, 34, 28, 12, 11, 15, 13, 45, 52, 35, 39, 29, 27, 30, 27),
var.e = c(1.1, 0.5, 1.3, 1.2, 0, 0, 0, 0, 1.06, 0.96, 1.06, 0.58, 0.1, 0, 0, 0),
var.g = c(83.2, 53.6, 94.9, 80.8, 25.6, 29.2, 31.9, 29.8, 102.3, 83.9, 158.7, 126.6, 78.5, 186.9, 82.4, 62))
PCA1 <- prcomp(ex[ , 2:4], scale = TRUE)
library(ggfortify)
autoplot(PCA1, data = ex, colour = "treat.year", loadings = TRUE,
loadings.label = TRUE, frame = TRUE)
