[BioC] How to get and save normalized data from limma
Wang, Jixin
jixinwang at neo.tamu.edu
Thu Sep 11 04:56:24 CEST 2008
Hi Mark,
Thank you very much for your kind reply.
I want to get normalized cy3 and cy5 channel data and tend to look at the correlation (scatter plot) of normalized data between biological replicates and dye swaps. But I can’t figure out R code for this task. I will use MeV(TM4) and Cluster/TreeView software to do the clustering.
Here is my target file:
SlideNumber FileName Cy3 Cy5 Date
13845971 13845971.gpr R wild type 8/12/2008
13845972 13845972.gpr wild type R 8/12/2008
13845392 13845392.gpr R wild type 8/12/2008
13845393 13845393.gpr wild type R 8/12/2008
13845400 13845400.gpr R wild type 8/12/2008
13845401 13845401.gpr wild type R 8/12/2008
Regards,
Jixin
----- Original Message -----
From: "Mark Robinson" <mrobinson at wehi.EDU.AU>
To: "Jixin Wang" <jixinwang at neo.tamu.edu>
Cc: bioconductor at stat.math.ethz.ch
Sent: Wednesday, September 10, 2008 5:00:23 PM GMT -06:00 US/Canada Central
Subject: Re: [BioC] How to get and save normalized data from limma
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
_______________________________________________
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
More information about the Bioconductor
mailing list