[BioC] Problem getting the exact ProbeNames

Karsten Voigt karsten.voigt at biologie.uni-freiburg.de
Tue Jan 11 18:56:32 CET 2011


Dear all,

I am currently working on a project where I need to get the exact IDs of 
probes of a custom Affymetrix Chip in order to merge it with another 
list containing the sequence.

I am using this small R script for creating the list:

mitdata <- ReadAffy();
stddata <- apply(pm(mitdata), 2, bg.adjust);
nrmdata <- normalize.quantiles(stddata);
namedata <- probeNames(mitdata);
enddata <- cbind(namedata, nrmdata);
write.table(enddata, file="probesdata.txt",sep="\t");

This is an output example

...
145 TZG_ARR_0001_x_at    135.115780787133 ...
146 TZG_ARR_0001_x_at    147.346049115501 ...
147 TZG_ARR_0001_x_at    203.840215898533 ...
148 TZG_ARR_0003_x_at    48.7635207480323 ...
...

As you can see, a number of probes have the same name but refer to 
different oligos. The number in front of the row is just added by me, 
therefore you can ignore it.

I received a list containing the probe name, a couple of other 
information AND the sequence.

This is a part of it:

15    ggagattgtttgtaatcaaaatgaa    TGZ_ARR_0001_x    !    2398    0    
176    200    +    1
103  gcaaatttacttctaacagctgatc    TGZ_ARR_0001_x    !    2398    1    
264    288    +    1
188  ttgatgcaactgtaaacaaaagtgg    TGZ_ARR_0001_x    !    2398    2    
349    373    +    1
15    gatagattcttcaagtaacaatact    TGZ_ARR_0003_x    !    2400    0    
2046    2070    +    1

This should be the same area.

In this received list, I can identify the unique probes using the 2 
numbers right after the exclamation mark, which are referring to the 
position on the chip, I guess. How can I extract those coordinates for 
my own list? I tried it with indices2xy, however I failed to get it 
running since I don't understand how to use this function correctly.

Thanks in advance for all answers,

Karsten Voigt





-- 

_________________________________________________
Karsten Voigt, Msc.
Experimentelle Bioinformatik, Hess Group
University of Freiburg, BIO III
t: 0761-2032708
m: 0176-61110420
e: karsten.voigt at biologie.uni-freiburg.de



More information about the Bioconductor mailing list