data(causton)
## Find the 3000 most variable genes, according to sd/mean:
varMeas<-function(vec) sd(vec)/mean(vec)
variability<-apply(causton,1,varMeas)
rks<-rank(variability)
causton3000<-causton[rks>length(rownames(causton))-3000,]
Run the code above in your browser using DataLab