[BioC] how to get normalized matrix
eyal avisar
enfer0666 at yahoo.com
Thu Oct 4 16:21:10 CEST 2012
hello,
my name is Eyal and i am new to R, cut and paste level.
i have created a DGEList object - d following case study 3.3.3.
i ran the calcNormFactors(d).
the output i need is a normalized table.
> library(edgeR)
> library(limma)
> raw.data <- read.delim("myFile.txt")
> names(raw.data)
d <- raw.data[, 2:64]
rownames(d) <- raw.data[, 1]
i did not group, irrelevant for me
d <- DGEList(counts = d)
dim(d)
result was [1] number... 3
> cpm.d <- cpm(d)
> d <- d[ rowSums(cpm.d > 1) >=3, ]
> d <- calcNormFactors(d)
thx,
Eyal
More information about the Bioconductor
mailing list