[R] Linear Discriminant Analysis
michael watson (IAH-C)
michael.watson at bbsrc.ac.uk
Wed Jun 6 19:57:10 CEST 2007
Region and Name are effectively the same variable
cor(olive[,4:11])
will also show you that there are strong correlations between some of the variables - this is something you might want to avoid....
________________________________
From: r-help-bounces at stat.math.ethz.ch on behalf of Soare Marcian-Alin
Sent: Wed 06/06/2007 4:45 PM
To: Uwe Ligges; R-help at stat.math.ethz.ch
Subject: Re: [R] Linear Discriminant Analysis
Thanks for explaining...
Im just sitting at the homework for 6 hours after taking for one week
antibiotica, because i had an amygdalitis...
I just wanted some tipps for solving this homework, but thanks, I will try
to get help on another way :)
I think i solved it, but I still get this Error :(
## Loading Data
library(MASS)
olive <- url("
http://www.statistik.tuwien.ac.at/public/filz/students/multi/ss07/olive.R")
print(load(olive))
dim(olive)
summary(olive)
index <- sample(nrow(olive), 286)
train <- olive[index,-11]
test <- olive[-index,-11]
summary(train)
summary(test)
table(train$Region)
table(test$Region)
# Linear Discriminant Analysis
z <- lda(Region ~ . , train)
zn <- predict(z, newdata=test)$class
mean(zn != test$Region)
2007/6/6, Uwe Ligges <ligges at statistik.uni-dortmund.de>:
>
>
> So what about asking your teacher (who seems to be Peter Filzmoser) and
> try to find out your homework yourself?
> You might want to think about some assumptions that must hold for LDA
> and look at the class of your explaining variables ...
>
> Uwe Ligges
>
>
>
> Soare Marcian-Alin wrote:
> > Hello,
> >
> > I want to make a linear discriminant analysis for the dataset olive, and
> I
> > get always this error:#
> > Warning message:
> > variables are collinear in: lda.default(x, grouping, ...)
> >
> > ## Loading Data
> > library(MASS)
> > olive <- url("
> >
> http://www.statistik.tuwien.ac.at/public/filz/students/multi/ss07/olive.R
> ")
> > print(load(olive))
> >
> > y <- 1:572
> > x <- sample(y)
> > y1 <- x[1:286]
> >
> > train <- olive[y1,-11]
> > test <- olive[-y1,-11]
> >
> > summary(train)
> > summary(test)
> >
> > table(train$Region)
> > table(test$Region)
> >
> > # Linear Discriminant Analysis
> > z <- lda(Region ~ . , train)
> > predict(z, train)
> >
> > z <- lda(Region ~ . , test)
> > predict(z, test)
> >
> > Thanks in advance!
> >
> >
> >
> > ------------------------------------------------------------------------
> >
> > ______________________________________________
> > R-help at stat.math.ethz.ch 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.
>
--
Mit freundlichen Grüssen / Best Regards
Soare Marcian-Alin
[[alternative HTML version deleted]]
More information about the R-help
mailing list