How to get probe names for raw intensity values?
2
0
Entering edit mode
vleblanc • 0
@vleblanc-7465
Last seen 7.5 years ago
Canada

I'm trying to get raw intensity values for a particular probeset from the Affy HuGene 2.0 array. I've tried getProbeInfo:

> subset(getProbeInfo(raw, target = 'probeset'), man_fsetid == "16862673")
           fid man_fsetid
378679  956480   16862673
450843 1140178   16862673
560361 1420274   16862673
602825 1535007   16862673

Is 'fid' equivalent to the rownames of the intensity matrix (exprs(raw))? The number of rows are not equivalent so I hesitate to assume this.

Thanks in advance for any help.

oligo • 838 views
ADD COMMENT
0
Entering edit mode
@eleni-christodoulou-2653
Last seen 5.6 years ago
Singapore

Hello,

I am not working with Affy Arrays but with Illumina. I guess the same things exist for both cases, though. For illumina, there is a manifest file for each array, which maps the  rownames of the array to symbols, entrez_ids etc. One of the manifest's columns is called 'array_address_id'. This is a big integer, similar to the fid number you described above. Could you find the manifest  for your case?

 

ADD COMMENT
0
Entering edit mode

Illumina and Affy arrays are completely different. There is no manifest file for Affy arrays.

ADD REPLY
0
Entering edit mode
@james-w-macdonald-5106
Last seen 8 minutes ago
United States

The fid IS equivalent to the row of the matrix you get from accessing the exprs slot from your GeneFeatureSet. So you could do

> z <- getProbeInfo(dat, target = "probeset")
> zlst <- split(z[,1], z[,2])
> zlst["16862873"]
$`16862873`
[1]  855406 1605173 1878834
ADD COMMENT
0
Entering edit mode

Great, thanks for your help!

ADD REPLY

Login before adding your answer.

Traffic: 874 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6