Actual code for LDA

Hi All,
         I've been trying to write my own code for LDA (linear discrim) so I
can modify it to be weighted LDA since some of my groups are outliers.

However, the code I write for standard LDA gives me slightly different
results to those from R (slightly different LDAs.....and not just scalar
differences.

Does R simply invert the pooled within groups scatter matrix e.g. solve(Sw)
and multiply it by the between groups scatter Sb, then take the eigen
vectors?

This is the approach that I use, yet get slightly different
results..........does R by default use a slightly more robust inverse or
estimate of Sw??

Actual code for LDA

Thread (2 messages)