Difference between revisions of "ExploreGolub.R"
From Organic Design wiki
(Adding golub ''Affymetrix'' prelim analysis) |
m (# {{R}}) |
||
| Line 1: | Line 1: | ||
| + | # {{R}} | ||
library(multtest) | library(multtest) | ||
data(golub) | data(golub) | ||
Revision as of 00:37, 6 June 2007
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)
- Abundance apears to be a slight mixture (hump more prevelant with background subtraction)
- Mixture could be of exponential and low expressing distribution on raw scale
hist((y+x)/2, breaks=nbreaks)
- 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)
