[BioC] Re: a bug? non agreement btw. indexProbes and affy CDF file

ulas karaoz ulask at bu.edu
Mon Nov 1 17:03:30 CET 2004


thanks for all the help, i have tried to do some comparison with the 
values you have provided:
In the AFFY CDF file for U133A chip, the following line is for the 
x=467, y=181:
Cell1=467       181     N       control 1007_s_at       0       13      
G       C       G       0       129339  -1      -1      99
and to x=467, y=182:
Cell2=467       182     N       control 1007_s_at       0       13      
G       G       G       0       130051  -1      -1      99


Why these indices are offset by 1 compared to yours, has this anything 
to do with the package wide xy.offset = 1 option?
I am very puzzled about what are the correct indices, can I plug in 
simply the indices from the CDF file, which seem to be computed as: x+ 
(712 * y) for U133A chip. Clearly I am missing something.
thanks.




On Nov 1, 2004, at 9:21 AM, James W. MacDonald wrote:

> ulas karaoz wrote:
>> thanks. That is not exactly what i meant, what i meant was given a 
>> probe's coordinates from the AFFY files, i.e. HG-U133A_probe_tab, how 
>> do I figure out the coordinates for the corresponding mismatch probe.
>> Is this info stored somewhere if i read the cdf file with 
>> read.cdffile in makecdfenv package?
>> thanks.
>
> Although indexProbes only gets you the index for a given probe, this 
> information can be used to get the xy coordinates and show you how 
> things are arranged on the chip.
>
> > dat <- ReadAffy()
> > indexProbes(dat, "pm")["1007_s_at"]
> $"1007_s_at"
>  [1] 129340 213420 396671  82246 430968 427082 432610  72465 432865 
> 99501 504952 443862 341432 198778
> [15] 463575  10989
>
> > indexProbes(dat, "mm")["1007_s_at"]
> $"1007_s_at"
>  [1] 130052 214132 397383  82958 431680 427794 433322  73177 433577 
> 100213 505664 444574 342144 199490
> [15] 464287  11701
>
> > i2xy(129340)
>        x   y
> [1,] 467 181
> > i2xy(130052)
>        x   y
> [1,] 467 182
>
> So the mm probes are located directly under the pm probes. Therefore, 
> the coordinates of the mm probe are x and y+1, where (x,y) are the pm 
> coordinates.
>
> Jim
>
>
>
> -- 
> James W. MacDonald
> Affymetrix and cDNA Microarray Core
> University of Michigan Cancer Center
> 1500 E. Medical Center Drive
> 7410 CCGC
> Ann Arbor MI 48109
>
--------------------------------------------
Ulas Karaoz
Boston University
Ph.D.Student in Bioinformatics
zlab.bu.edu



More information about the Bioconductor mailing list