[BioC] Limma question_Intra-Spot Correlation Question
Gordon K Smyth
smyth at wehi.EDU.AU
Sun Feb 5 22:52:26 CET 2012
Dear Judith,
Getting negative intensities has nothing to do with normalization, but
everything to do with background correction. Please see page 24 of the
limma User's Guide:
"we often find
RG <- backgroundCorrect(RG, method="normexp", offset=50)
to be preferable to the simple background subtraction when using output
from most image analysis programs. This method adjusts the foreground
adaptively for the background intensities and results in strictly positive
adjusted intensities, i.e., negative or zero corrected intensities are
avoided."
Best wishes
Gordon
> Message: 10
> Date: Sun, 05 Feb 2012 11:45:12 +0100
> From: jgomez at uni-potsdam.de
> To: Bioconductor mailing list <bioconductor at r-project.org>
> Subject: [BioC] Limma question_Intra-Spot Correlation Question
>
> Dear List,
> I have a set of Agilent chips that I would like to analyse following
> the procedure for separate channel analysis.
> I just want to add that following the brief procedure outline in Limma
> User Guide wasn't useful.
> I have made the targets file as usual
> FileName Cy3 Cy5
> 1 US83800208_252412610022_1_4.txt WT_4 OX_4
> 2 US83800208_252412610019_1_1.txt KD_4 WT_4
> 3 US83800208_252412610019_1_2.txt OX_4 KD_4
> 4 US83800208_252412610019_1_3.txt OX_21 OX_4
> 5 US83800208_252412610019_1_4.txt WT_21 WT_4
> 6 US83800208_252412610020_1_1.txt KD_4 KD_21
> 7 US83800208_252412610020_1_2.txt WT_21 OX_21
> 8 US83800208_252412610021_2_1.txt KD_21 WT_21
> 9 US83800208_252412610020_1_3.txt OX_21 KD_21
> I read the chips using read.maimages and creating an RG object along
> the way. Then I normalized the arrays using Aquantile (Normalization
> Between Arrays). I convert the targets file and then looks like this:
> channel.col FileName Target
> 1.1 1 US83800208_252412610022_1_4.txt WT_4
> 1.2 2 US83800208_252412610022_1_4.txt OX_4
> 2.1 1 US83800208_252412610019_1_1.txt KD_4
> 2.2 2 US83800208_252412610019_1_1.txt WT_4
> 3.1 1 US83800208_252412610019_1_2.txt OX_4
> 3.2 2 US83800208_252412610019_1_2.txt KD_4
> 4.1 1 US83800208_252412610019_1_3.txt OX_21
> 4.2 2 US83800208_252412610019_1_3.txt OX_4
> 5.1 1 US83800208_252412610019_1_4.txt WT_21
> 5.2 2 US83800208_252412610019_1_4.txt WT_4
> 6.1 1 US83800208_252412610020_1_1.txt KD_4
> 6.2 2 US83800208_252412610020_1_1.txt KD_21
> 7.1 1 US83800208_252412610020_1_2.txt WT_21
> 7.2 2 US83800208_252412610020_1_2.txt OX_21
> 8.1 1 US83800208_252412610021_2_1.txt KD_21
> 8.2 2 US83800208_252412610021_2_1.txt WT_21
> 9.1 1 US83800208_252412610020_1_3.txt OX_21
> 9.2 2 US83800208_252412610020_1_3.txt KD_21
> When using the function intraspotCorrelation I got an error regarding
> " Missing or infinite values found in M or A".
> Checking older post regarding the erro, it may be becaouse few probes
> have after normalization negative intensities.
> I would like to know if the process I started is the right one for
> this kind of analysis and second If there is a kind of filter I can
> use in limma to get rid of those neg intensities to proceed to the
> next step.
> Thanks in advance for any help.
> Cheers,
> Judy
>
> --
> Judith Lucia Gomez, PhD
> Centre for Plant Biotechnology and Genomics - CBGP
> 28223 Pozuelo de Alarc?n (Madrid)
> Spain
______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}
More information about the Bioconductor
mailing list