ExploreGolub.R

Code snipits and programs written in R, S or S-PLUS library(multtest) data(golub)

nbreaks=100

x <- rowMeans(golub[,golub.cl==0]) y <- rowMeans(golub[,golub.cl==1])

mdplot(cbind(x,y))

hist(y-x, breaks=nbreaks) hist(2^(y-x), breaks=nbreaks)

1. Abundance apears to be a slight mixture (hump more prevelant with background subtraction)
2. Mixture could be of exponential and low expressing distribution on raw scale

hist((y+x)/2, breaks=nbreaks)

1. Abundance, observed exponential like decay on raw scale

hist(2^((y+x)/2), breaks=nbreaks) qqplot(rexp(10000,1), 2^((y+x)/2)) abline(0,1)