Question: normalization of a miRNA microarray
0
gravatar for alptaciroglu
3.7 years ago by
alptaciroglu50
Turkey
alptaciroglu50 wrote:

Dear All,

 

I am trying to process an affymetrix miRNA microarray that is shown to have 7815 rownames (miRNA identifiers). When I read the CEL files with ReadAffy() function :

 

RAW <- ReadAffy()

dim(exprs(RAW)) yields 52900 probes where there should be 7815. Affybatch object also says that there are 7815 genes but why expression file of the Affybatch object shows up 52900 rows ? This is dealt with in rma() normalization and there stays correctly 7815 rows but other normalization packages (vsn, ExiMiR) cannot deal with this and there remains 52900 probes after normalization. I would appreciate if you could let me know what is going on. Thanks 


With kind regards,


Alperen Taciroglu

 

ADD COMMENTlink modified 3.7 years ago by Wolfgang Huber13k • written 3.7 years ago by alptaciroglu50
Answer: normalization of a miRNA microarray
1
gravatar for James W. MacDonald
3.7 years ago by
United States
James W. MacDonald50k wrote:

Affy arrays have multiple 25-mers that together are intended to measure a single transcript. This is silly in the case of the miRNA arrays, as the target is shorter than the oligo that is used to measure (and there are multiple identical 25-mers for each miRNA).

Anyway, that's why there are more probes before you normalize. Also note that the vsn package has a function called vsnrma() that does understand what to do with Affy arrays.

ADD COMMENTlink written 3.7 years ago by James W. MacDonald50k

It seems vsnrma() correctly returns 7815 rows. Thanks very much

ADD REPLYlink written 3.7 years ago by alptaciroglu50
Answer: normalization of a miRNA microarray
0
gravatar for Wolfgang Huber
3.7 years ago by
EMBL European Molecular Biology Laboratory
Wolfgang Huber13k wrote:

What do you get from

cdfName(RAW)

Affymetrix arrays usually have multiple features (probes) per probe set (targeted gene). Could it be that 7815 is the number of probe sets and 52900 the number of probes?

 

Wolfgang

 

ADD COMMENTlink written 3.7 years ago by Wolfgang Huber13k

miRNA-1_0_2Xgain I got from cdfName(RAW). I do not know the answer to your question but I think James' reply might be answer to both of us. Anyway, thanks so much for reply

ADD REPLYlink written 3.7 years ago by alptaciroglu50

miRNA-1_0_2Xgain I got from cdfName(RAW). I do not know the answer to your question but I think James' reply might be answer to both of us. Anyway, thanks so much for reply

ADD REPLYlink written 3.7 years ago by alptaciroglu50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 136 users visited in the last hour