################################################### ### chunk number 1: ################################################### library(ALL) data(ALL) show(ALL) ################################################### ### chunk number 2: ################################################### print(summary(pData(ALL))) ################################################### ### chunk number 3: ################################################### hist(cvv <- apply(exprs(ALL),1,function(x)sd(x)/mean(x))) ################################################### ### chunk number 4: ################################################### ok <- cvv > .08 & cvv < .18 fALL <- ALL[ok,] show(fALL) allx2 <- data.frame(t(exprs(fALL)), class=ALL$BT) ################################################### ### chunk number 5: ################################################### library(rpart) rp1 <- rpart(class~.,data=allx2) plot(rp1) text(rp1)