search for: sndgene

...;) library("klaR") td <- read.table("dgdata1.txt", header=TRUE, sep=",") dgenes <- subset(td, dg == 1, select = dg:eg) ndgenes <- subset(td, dg == 0, select = dg:eg) n1 <- nrow(dgenes) n2 <- nrow(ndgenes) ndgrows <- 1:n2 selrows <- sample(ndgrows) sndgenes <- ndgenes[selrows[1:n1],] train <- rbind(dgenes, sndgenes) attach(train) traind <- subset(train, select = -dg) trainr <- factor(dg) detach(train) sc_res <- stepclass(traind, trainr, "svm", direction = "forward", criterion = "AC", fold = 10)