search for: lda1

All Is there an R command to produce error bars around a plotted point. Crawley's book uses the command error.bar() but my version of R rejects it. Must be an S+ command(?). Thanks REX [[alternative HTML version deleted]]

Dear R-users How can I get the eigenvalues out of an lda analysis? thanks a lot christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

...also produces different (though rather similar) results for p=10 (no longer collinear). See here: library(MASS) set.seed(12345) n <- 50 p <- 200 # or p<- 10 testdata <- matrix(ncol=p,nrow=n) for (i in 1:p) testdata[,i] <- rnorm(n) class <- as.factor(c(rep(1,25),rep(2,25))) lda1 <- lda(x=testdata,grouping=class,CV=TRUE) table1 <- table(lda1$class,class) y.lda <- rep(NA, n) for(i in 1:n){ testset <- testdata[i,,drop=FALSE] trainset <- testdata[-i,] model.lda <- lda(x=trainset,grouping=class[-i]) y.lda[i] <- predict(model.lda, testset)$clas...

...oss-validation routine that I run (of course mine may be wrong, but I don't think so). See here: library(MASS) set.seed(12345) n <- 50 p <- 10 # or p<- 200 testdata <- matrix(ncol=p,nrow=n) for (i in 1:p) testdata[,i] <- rnorm(n) class <- as.factor(c(rep(1,25),rep(2,25))) lda1 <- lda(x=testdata,grouping=class,CV=TRUE) table1 <- table(lda1$class,class) y.lda <- rep(NA, n) for(i in 1:n){ testset <- testdata[i,,drop=FALSE] trainset <- testdata[-i,] model.lda <- lda(x=trainset,grouping=class[-i]) y.lda[i] <- predict(model.lda, testset)$clas...

Hi R experts. I performed an Linear Discriminant Analysis (lda) and now I want to plot the first two axes (LDA1 and LDA2). Well MASS package have the plot.lda and pairs.lda to do that. But, they don’t let me personalize them, once they don’t accept the type=’n’ plot. So I start looking at the lda object properties, trying to found out my groups (n=4) axis coordinates, curiously I couldn’t. There is only the...