自己顶个

以前有人遇到过，得到的答复是——
plot.lda does not accept type="n" precisely because it does support a panel function to plot within the axes it has set up.
所以这个暂时实现不了了

很简单吧，假设以iris 数据为例子
library(MASS)
model1=lda(Species~.,data=iris,prior=c(1,1,1)/3)
ld=predict(model1)$x

#画出ld1 与ld2
p=ggplot(cbind(iris,as.data.frame(ld))
,aes(x=LD1,y=LD2))
p+geom_point(aes(colour=Species),alpha=0.8,size=3)