Код ниже
require(MASS)
LDA.nursery <- lda(formula = V9 ~., data = nursery)
LDA.nursery$scaling
LDA.nursery$svd
(prop = LDA.nursery$svd^2/sum(LDA.nursery$svd^2))
prop <- scales::percent(prop)
pred <- predict(LDA.nursery, newdata = nursery)
scores <- data.frame(V9 = nursery$V9, pred$x)
require(ggplot2)
ggplot() + geom_point(data = scores,
aes(x = LD1, y = LD2, shape = V9,
colour = V9), size = 3) +
scale_colour_manual(values = c('purple', 'green', 'blue', 'yellow', 'pink')) +
labs(x = paste("LD1 (", prop[1], ")", sep = ""),
y = paste("LD2 (", prop[2], ")", sep = "")) + theme_bw()
lda(scale(nursery[, 1:9]), gr = nursery$V9)$scaling