Re: [R] help with lda function from MASS package

2009-09-29 Thread David Winsemius
If you want to see how Venables and Ripley get their lda distances, then this is a quick path to the (uncommented) source: 1) > methods(lda) [1] lda.data.frame* lda.default*lda.formula*lda.matrix* 2) since you call involved a formula I looked first at: > getAnywhere(lda.formula) A singl

Re: [R] help with lda function from MASS package

2009-09-29 Thread Pete Shepard
Thanks David, Yes, I am talking about the MASS package.Thank you for pointing out that these scale the same. My question is, how do I get from the V1 data: V1 1 164.4283 2 166.2492 3 170.5232 4 156.5622 5 127.7540 6 136.7704 7 136.3436 to the other set of data: + 1 -2.3769280 + 2 -2.704

Re: [R] help with lda function

2009-09-28 Thread David Winsemius
Your results are the same (after scaling and sign reversal) out to the 4th decimal place as those from lda (which by the way is almost certainly from the MASS package and not from an impossible to find "lda package".) > read.table(textConnection(txt)) V1 1 164.4283 2 166.2492 3 170.

[R] help with lda function

2009-09-28 Thread Pete Shepard
I am having a problem understanding the lda package. I have a dataset here: [,1] [,2] [,3] [1,] 2.95 6.630 [2,] 2.53 7.790 [3,] 3.57 5.650 [4,] 3.16 5.470 [5,] 2.58 4.461 [6,] 2.16 6.221 [7,] 3.27 3.521 If I do the following; "names(d)<-c("y","x1","x2") d$x1 = d$x