# load sda library
library("sda")
# load full Khan et al (2001) data set
data(khan2001)
dim(khan2001$x)
khan2001$y
# create data set containing only the SRBCT samples
idx = which( khan2001$y == "non-SRBCT" )
srbct.x = khan2001$x[-idx,]
srbct.y = factor(khan2001$y[-idx])
srbct.descr = khan2001$descr[-idx]
dim(srbct.x)
srbct.y
Run the code above in your browser using DataLab