Difference between revisions of "ExploreSwirl.R"
From Organic Design wiki
m |
|||
| Line 13: | Line 13: | ||
MA <- MA.RG(RG) | MA <- MA.RG(RG) | ||
| − | # [[[[http://en. | + | # [[[[http://en.wikipedia.org/wiki/histogram|Histogram]]]] not centred |
hist(rowMeans(MA$M), breaks=nbreaks) | hist(rowMeans(MA$M), breaks=nbreaks) | ||
hist(rowMeans(MA$A), breaks=nbreaks) | hist(rowMeans(MA$A), breaks=nbreaks) | ||
Revision as of 04:09, 22 August 2006
library(marray) library(limma)
bc.method <- "subtract" nbreaks <- 100
dataDir <- system.file("swirldata", package = "marray")
RG <- read.maimages(path=dataDir, source="spot", ext="spot")
- Right skew
plotDensities(RG, log=FALSE)
MA <- MA.RG(RG)
- [[[[1]]]] not centred
hist(rowMeans(MA$M), breaks=nbreaks) hist(rowMeans(MA$A), breaks=nbreaks) for(i in seq(ncol(RG))) {
plotMA(RG, array=i)
}
MA <- normalizeWithinArrays(RG, method="loess", bc.method = bc.method)
for(i in seq(ncol(RG))) {
plotMA(MA, array=i) Sys.sleep(0.5)
}
- CLT: After normalization the distribution is centred around 0 (right long tail)
hist(rowMeans(MA$M), breaks=nbreaks) hist(2^rowMeans(MA$M), breaks=nbreaks)
- Abundance apears to be a slight mixture (hump more prevelant with background subtraction)
hist(rowMeans(MA$A), breaks=nbreaks)
- Abundance, observed exponential decay on raw scale
hist(2^rowMeans(MA$A), breaks=nbreaks) qqplot(rexp(10000,1/6000), 2^rowMeans(MA$A)) abline(0,1)
