Pavel Kúr
2009-Apr-03 20:29 UTC
[R] Discriminant Analysis - Obtaining Classification Functions
Hello! I need some help with the linear discriminant analysis in R. I have some plant samples (divided into several groups) on which I measured a few quantitative characteristics. Now, I need to infer some classification rules usable for identifying new samples. I have used the function lda from the MASS library in a usual fashion: lda.1 <- lda(groups~char1+char2+char3, data=xxx) I'd like to obtain the classification functions for the particular groups, with the aid of which I could classify unknown samples. I know I can use predict.lda to classify such samples, but I need to obtain some equations into which I could simply put the measured values of an unknown sample manually and the result would predict which group the sample most probably belongs to (like in eg. STATISTICA). I haven't found out how to extract these functions from the lda output. Could somebody give me some advice? Thank you in advance, Pavel Kur
Emmanuel Charpentier
2009-Apr-03 21:24 UTC
[R] Discriminant Analysis - Obtaining Classification Functions
reauire(MASS) ; ?predict.lda should enlighten you. Glancing at V&R4 might be a bit more illuminating... HTH Emmanuel Charpentier Le vendredi 03 avril 2009 ? 22:29 +0200, Pavel K?r a ?crit :> Hello! > > I need some help with the linear discriminant analysis in R. > I have some plant samples (divided into several groups) on which I > measured a few quantitative characteristics. Now, I need to infer some > classification rules usable for identifying new samples. > I have used the function lda from the MASS library in a usual fashion: > > lda.1 <- lda(groups~char1+char2+char3, data=xxx) > > I'd like to obtain the classification functions for the particular > groups, with the aid of which I could classify unknown samples. I know > I can use predict.lda to classify such samples, but I need to obtain > some equations into which I could simply put the measured values of an > unknown sample manually and the result would predict which group the > sample most probably belongs to (like in eg. STATISTICA). > I haven't found out how to extract these functions from the lda output. > Could somebody give me some advice? > > Thank you in advance, > > Pavel Kur > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >