如何在 r 中绘制 LDA 的双标图？

我使用该函数进行了线性判别分析lda()来自 MASS 包。现在我会尝试绘制像 ade4 包（forLDA）中那样的双标图。你知道我该怎么做吗？

如果我尝试使用biplot()功能不起作用。例如，如果我使用 Iris 数据并制作 LDA：

dis2 <- lda(as.matrix(iris[, 1:4]), iris$Species)

然后我可以使用该函数绘制它plot()，但是如果我使用这个函数biplot()它不起作用：

biplot(dis2)
Error in nrow(y) : argument "y" is missing, with no default

如何绘制变量的箭头？

我编写了以下函数来执行此操作：

lda.arrows <- function(x, myscale = 1, tex = 0.75, choices = c(1,2), ...){
  ## adds `biplot` arrows to an lda using the discriminant function values
  heads <- coef(x)
  arrows(x0 = 0, y0 = 0, 
         x1 = myscale * heads[,choices[1]], 
         y1 = myscale * heads[,choices[2]], ...)
  text(myscale * heads[,choices], labels = row.names(heads), 
    cex = tex)
}

对于你的例子：

dis2 <- lda(as.matrix(iris[, 1:4]), iris$Species)
plot(dis2, asp = 1)
lda.arrows(dis2, col = 2, myscale = 2)

箭头的长度相对于 lda 图是任意的（当然，箭头之间的长度不是任意的！）。如果您想要更长或更短的箭头，请更改myscale因此。默认情况下，这会绘制第一轴和第二轴的箭头。如果要绘制其他轴，请更改choices来反映这一点。