[BioC] how to get normalized matrix
Mark Robinson
mark.robinson at imls.uzh.ch
Thu Oct 4 23:12:21 CEST 2012
Hi Eyal,
You might read the posting guide:
http://www.bioconductor.org/help/mailing-list/posting-guide/
I don't really understand your question, but from your subject and your cryptic description, I think you've used the right function - cpm() - to get a normalized matrix. You might also consider using :
cpm.d <- cpm(d, normalized.lib.sizes=TRUE)
More description by typing:
?cpm
Best, Mark
On 04.10.2012, at 16:21, eyal avisar wrote:
> 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_______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list