小编Joh*_*Sal的帖子

我正在尝试使用 MASS 和 ggplot2 包绘制 Iris 数据集二次判别分析 (QDA) 的结果。脚本的第一部分显示了线性判别分析 (LDA)，但我不知道要继续为 QDA 执行此操作。“qda”类的对象与“lda”类对象有点不同，例如：我找不到解释的组间方差/判别分量的迹线比例/X%，无法将它们添加到图中轴。任何帮助或想法如何使用 ggplot2 编码此图？

代码：

require(MASS)
require(ggplot2)
require(scales)
 

irislda <- lda(Species ~ ., iris)
prop.lda = irislda$svd^2/sum(irislda$svd^2)
plda <- predict(irislda,   iris)

datasetLDA = data.frame(species = iris[,"Species"], irislda = plda$x)
ggplot(datasetLDA) + geom_point(aes(irislda.LD1, irislda.LD2, colour = species, shape = species), size = 2.5) + 
    labs(x = paste("LD1 (", percent(prop.lda[1]), ")", sep=""),
       y = paste("LD2 (", percent(prop.lda[2]), ")", sep=""))

 
irisqda <- qda(Species ~ ., iris)
pqda <- predict(irisqda,   iris)
datasetQDA = data.frame(species …