In R,
lda.score<-lda(class ~., data=iris)
if we use it only with two classes, we can get LDA scores (the position for linear discriminant