Re: [R] lda()

From: Roger D. Peng <rpeng_at_jhsph.edu>
Date: Tue 13 Jul 2004 - 06:49:55 EST

I haven't done this in years but I think the `scaling' element in the list returned by lda is the original data matrix multiplied by the rotation matrix from the SVD. Taking a look at

getAnywhere("lda.default")

will probably answer your question.

-roger

marzban wrote:
>>Perhaps some reading would be helpful. I suggest you look first at the
>>help file for lda(). Second, I suggest you read Venables and Ripley, MASS,
>>4th Edition, where lda() is discussed extensively. Third, I suggest you
>>read Ripley's Pattern Recognition and Neural Networks, where the theory is
>>laid out clearly. Both of these latter books are referenced in lda's help
>>file.
>
>
> I am quite familiar with all of these references. The theory behind LDA
> is not where the problem is - I'm comfortable with that. The problem is
> that I do not know what R is computing when it prints "Coefficients of
> linear discriminants". (According to the source code (lda.R), it's
> x$scaling, but I don't know what that is either.)
>
>
>>Finally, you might want to tell us what version of lda() you're using, what
>>version of R you're using, and what platform you're running on. For all
>>we know, you're using a 2-year old version of R and lda, both long
>>superceded by vastly improved programs and packages.
>
>
> R-1.9.1 on linux
>
> Caren
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list

>
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>



R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Jul 13 06:59:55 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:54:52 EST