problem with custom cdf
1
0
Entering edit mode
tadoktor • 0
@tadoktor-7140
Last seen 9.9 years ago
Belgium

Hi,

I am new in R and i trying to re-annotate some Affy cell files to Ensemle ID and then afterwards normalize the data with RMA. The problem  is that as an end result in stead of having  a reduced number of IDs (from the normal 54 000 Affy IDs to something around 20-15 000 ENSE IDs) which is the logic step, i end up with a list of around 77 000 IDs/genes. Does anyone has an idea what i am doing wrong?

Here is the code:

library(hgu133plus2hsenseprobe)
library(hgu133plus2hsensecdf)
CELs.start <- ReadAffy(celfile.path = "/Users/dimiter/genomark/NewTest/Working",cdfname="hgu133plus2hsensecdf")
eset.rma <- rma(CELs.start)
matrix <- exprs(eset.rma)

Thanks in advance!

 

 

microarray affy hgu133plus2 annotation cdf • 1.5k views
ADD COMMENT
2
Entering edit mode
@james-w-macdonald-5106
Last seen 1 day ago
United States

You haven't done anything wrong, except perhaps for being confused as to what ENSE means. Please note that ENSE is Ensembl Exon IDs, whereas ENSG is Ensembl Gene IDs. So if you want to summarize at the gene level, based on Ensembl genes, then you need the hgu133plus2hsensgcdf, which has 20,049 genes as compared to the ENSE version, which has 77,898 exons.

ADD COMMENT
0
Entering edit mode

Many thanks for the answer!  :)

ADD REPLY

Login before adding your answer.

Traffic: 503 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