Re: [R] LDA and confidence ellipse
Thanks again -- View this message in context: http://r.789695.n4.nabble.com/LDA-and-confidence-ellipse-tp4671308p4671427.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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.
[R] LDA and confidence ellipse
Hi, I wish to add confidence ellipse on my LDA result of the iris data set. Therefore: Is there statistical logic to do that as I only wish it to make the species separation more visable? How can I add it to the script below (ggplot): require(MASS) require(ggplot2) iris.lda-lda(Species ~ Sepal.Length + Sepal.Width + Petal.Length + Petal.Width, data = iris) LD1-predict(iris.lda)$x[,1] LD2-predict(iris.lda)$x[,2] ggplot(iris, aes(x=LD1, y=LD2, col=iris$Species) ) + geom_point( size = 4, aes(color = iris$Species))+theme_bw() Could someone please help me. Thank you very much. -- View this message in context: http://r.789695.n4.nabble.com/LDA-and-confidence-ellipse-tp4671308.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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.
Re: [R] LDA and confidence ellipse
Hi,
May be this helps:
require(MASS)
require(ggplot2)
iris.lda-lda(Species ~ Sepal.Length + Sepal.Width + Petal.Length +
Petal.Width, data = iris)
datPred-data.frame(Species=predict(iris.lda)$class,predict(iris.lda)$x)
library(ellipse)
dat_ell - data.frame()
for(g in levels(datPred$Species)){
dat_ell - rbind(dat_ell,
cbind(as.data.frame(with(datPred[datPred$Species==g,], ellipse(cor(LD1, LD2),
scale=c(sd(LD1),sd(LD2)),
centre=c(mean(LD1),mean(LD2),Species=g))
}
ggplot(datPred, aes(x=LD1, y=LD2, col=Species) ) + geom_point( size = 4,
aes(color =
Species))+theme_bw()+geom_path(data=dat_ell,aes(x=x,y=y,color=Species),size=1,linetype=2)
A.K.
Hi,
I wish to add confidence ellipse on my LDA result of the iris data set.
Therefore:
Is there statistical logic to do that as I only wish it to make the species
separation more visable?
How can I add it to the script below (ggplot):
require(MASS)
require(ggplot2)
iris.lda-lda(Species ~ Sepal.Length + Sepal.Width + Petal.Length +
Petal.Width, data = iris)
LD1-predict(iris.lda)$x[,1]
LD2-predict(iris.lda)$x[,2]
ggplot(iris, aes(x=LD1, y=LD2, col=iris$Species) ) + geom_point( size = 4,
aes(color = iris$Species))+theme_bw()
Could someone please help me. Thank you very much.
__
R-help@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.
Re: [R] LDA and confidence ellipse
Hi, Thanks works like magic. BTW What is the confidence ellipses probability used? -- View this message in context: http://r.789695.n4.nabble.com/LDA-and-confidence-ellipse-tp4671308p4671357.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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.
Re: [R] LDA and confidence ellipse
Hi, No problem. The default should be 0.95 ?ellipse() level: The confidence level of a pairwise confidence region. The default is 0.95, for a 95% region. This is used to control the size of the ellipse being plotted. A vector of levels may be used. A.K. - Original Message - From: Lluis l.r...@gmx.us To: r-help@r-project.org Cc: Sent: Thursday, July 11, 2013 3:15 PM Subject: Re: [R] LDA and confidence ellipse Hi, Thanks works like magic. BTW What is the confidence ellipses probability used? -- View this message in context: http://r.789695.n4.nabble.com/LDA-and-confidence-ellipse-tp4671308p4671357.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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. __ R-help@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.
