[BioC] How to get and save normalized data from limma

[Mark Robinson]

Jixin.

I believe you want the normalized log-ratios, which are stored in MA$M.

For example:

library(gplots)
heatmap.2(MA$M)

... of course, depending on the number of rows in your data, you may  
want to take a subset ...

HTH,
Mark

On 11/09/2008, at 6:27 AM, Wang, Jixin wrote:

> Dear All,
>
> I want to use the normalized signal intensities data to do  
> hierarchical clustering. But I don’t know how to save the normalized  
> data from LIMMA.
>
>
> Below is the code:
>
> library(limma)
> setwd("C:/Documents and Settings/JWang/Desktop/analysis with R/RAO")
> targets <- readTargets("equi.txt")
> targets
> RG <- read.maimages(targets$FileName, source="genepix")
> RG
>
> MA <- normalizeWithinArrays(RG)
>
> MA <- normalizeBetweenArrays(MA,method="scale")
>
> design <- c(-1,1,-1,1,-1,1)
> fit <- lmFit(MA,design)
>
>
> fit <- eBayes(fit)
> topTable(fit,number=500,adjust="BH")
>
> volcanoplot(fit)
>
> Thanks a lot,
>
> Jixin
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor

------------------------------
Mark Robinson
Epigenetics Laboratory, Garvan
Bioinformatics Division, WEHI
e: m.robinson at garvan.org.au
e: mrobinson at wehi.edu.au
p: +61 (0)3 9345 2628
f: +61 (0)3 9347 0852