affy MAplot NA error
1
0
Entering edit mode
3.7 years ago
ciemanek ▴ 140

Hi all!

I am trying to draw MA plots for my Affymetrix microarrays using affy package and function MAplot. I am trying to take a look at my raw data in AffyBatch object (1st microarray) as follows:

MAplot(data.raw, which=1)

However, I am keep on gettin an error stating:

Error in quantile.default(as.numeric(x), c(0.25, 0.75), na.rm = na.rm,  : 
  missing values and NaN's not allowed if 'na.rm' is FALSE

I don't know how could I make it work? Any ideas and suggestions on what might be the issue will be appreciated.

Cheers! Agata

affymetrix maplot microarrays affy quantile • 1.4k views
ADD COMMENT
0
Entering edit mode
3.7 years ago

Hi Agata,

You are not supposed to generate a MA plot from raw data. You first have to normalise your raw data and then pass the normalised object to MAplot().

Please take a look at the first few pages of the tutorial affyPLM: Advanced use of the MAplot function

Kevin

ADD COMMENT
1
Entering edit mode

Thanks a lot Kevin!

However, if I understand correctly from the vignette you recommended, it should be possible to make an MA plot from AffyBatch object, including object with the raw data - the function should by default transform it to logarithms?

All the best, Agata

* EDIT *

I managed to plot M vs. A on the raw data by using only PM probes - otherwise there were missing values in my MM probes which were causing issues with MAplot() function:

MAplot(data.raw, which = 1, log = TRUE, type="pm")
ADD REPLY

Login before adding your answer.

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