[BioC] Limma question_Intra-Spot Correlation Question
jgomez at uni-potsdam.de
jgomez at uni-potsdam.de
Sun Feb 5 11:45:12 CET 2012
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
More information about the Bioconductor
mailing list