#train classifier on Singh's Prostate Cancer Data set,
# selecting genes by the Extended HC scheme
ldarule <- RFlda(SinghDS[,-1],SinghDS[,1])
# get in-sample classification results
predict(ldarule,SinghDS[,-1],grpcodes=levels(SinghDS[,1]))$class
# compare classifications with true assignments
cat("Original classes:
")
print(SinghDS[,1])
# show set of selected genes
cat("Genes kept in discrimination rule:
")
print(colnames(SinghDS)[ldarule$vkpt])
cat("Number of selected genes =",ldarule$nvkpt,"")Run the code above in your browser using DataLab