[FIX] orthogonal LDA transform (conform with R)

Removed an unnecessary dot product from the `LDA.transform()`. Now it gives the same results as R on the setosa data set. See discussion in #3500 for details.

Regression test included.
can we not call the output of this operation `X` and rather `X_new` as stated in the docstring?