[BioC] difference between hgu95av2probe library and pm matrix

Laurent Gautier laurent at cbs.dtu.dk
Fri Sep 19 13:15:18 MEST 2003


On Fri, Sep 19, 2003 at 09:55:25AM +0200, M.Ruschhaupt at dkfz-heidelberg.de wrote:
>  
> > Quantifying the differences would be a cute example to show the 
> > flexibility 
> > of the 'affy' package... it is not a too difficult thing do 
> > (hint: make use of the corresponding cdfenv packs). 
> > Do not hesitate to share your code.   
> >  
> Here is the code I used to see, which probes are not listed in the hgu95av2probe object. It is 
> far away from being efficient and the difference is only checked for the Dilution set of the 
> affycomp package but that should not matter. 
>  
> > load("Dilution.rda") 
> > library(hgu95av2probe) 
> >  p.index <- xy2indices(hgu95av2probe$x+1,hgu95av2probe$y+1,ab=Dilution) 
> > length(p.index) 
> [1] 199084 
> >  ind <- unlist(indexProbes(Dilution, "pm")) 
> > length(ind) 
> [1] 201800 
> > ind2 <- ind[! ind %in% p.index] 
> > length(ind2) 
> [1] 2716 
> > ind2[1:5] 
> 1142_at1 1142_at2 1142_at3 1142_at4 1142_at5 
>   387228    17915   239242   334935    48584 
>  
> Maybe this can help to find out why these probes are not in the hgu95av2probe object. 

I think one has to go to the annotation for these probes. For the particular
example you give, the information available are:
(copy/paste from http://www.medfac.leidenuniv.nl/lgtc/db/Affy/HG-U95A.txt)
HG-U95A	1142_at	HG3432-HT3618	"Fibroblast Growth Factor Receptor K-Sam, Alt. Splice 1"	 	individual sequence
HG-U95A	1143_s_at	HG3432-HT3620	"Fibroblast Growth Factor Receptor K-Sam, Alt. Splice 3, K-Sam III"	 	group members not included
HG-U95A	1144_at	HG3432-HT3621	"Fibroblast Growth Factor Receptor K-Sam, Alt. Splice 4, K-Sam IV"	 	individual sequence
HG-U95A	1145_g_at	HG3432-HT3621	"Fibroblast Growth Factor Receptor K-Sam, Alt. Splice 4, K-Sam IV"	 	group members not included


I shall be able mail a little more about alternative splicing and
Affymetrix probes within the very near future (hopefully).


L.



>  
> Markus 

-- 
--------------------------------------------------------------
Laurent Gautier			CBS, Building 208, DTU
PhD. Student			DK-2800 Lyngby,Denmark	
tel: +45 45 25 24 89		http://www.cbs.dtu.dk/laurent



More information about the Bioconductor mailing list