Hi, I was wondering if you can show me how to plot the discriminant boundary lines for an LD analysis in R. I am curious as to perform this on a larger scale, so I was wondering if you can provide me an example on the infamous 'iris' data. Here's what I have so far: iris.lda=lda(Species ~., data=iris) iris.lda plot(iris.lda) col <- rep(c("red", "black", "blue"), each=50) pch <- rep(1:3, each=50) plot(iris.lda, col=col, pch=pch) % This will not work to plot the coefficients. abline(lda(Species~.,data=iris)$scaling, col="black") Thanks, Lee _________________________________________________________________ [[alternative HTML version deleted]]