#!/usr/bin/Rscript

n <- 5
data <- array(rnorm(21*n*100), c(20,n+1,100))
data[,1,] <- as.numeric(runif(n*100)>0.5)
dimnames(data) <- list(paste("sample", 1:20),c("label", paste("dim", as.character(1:n))),paste("data set", 1:100))
f <- data.frame(data[,,2])
z <- lda(label ~ ., f, CV=TRUE)