scikit-learn LDA produces incorrect results when used for dimensionality reduction. I went back and implemented an LDA step-by-step, and compared my results with Matlab's and R's LDA functions on the same dataset (Iris). There is definitely something weird about the scikit-function.