How to retrieve intensities using affy

From BITS wiki
Jump to: navigation, search
Go to parent Analyze your own microarray data in R/Bioconductor

Although AffyBatches have a slot assayData to hold the raw intensities, this slot is never directly accessed. Instead two methods give access to the raw intensities in an AffyBatch: exprs() and intensity(). These methods extract the intensities of all probes (both PM and MM probes) from the AffyBatch.

Since we will only work with PM probes, you might want to look at the intensities of the PM probes only using the pm() method.

It is of course much more useful to retrieve intensities based on probe set ID then on location in the CEL file. Often you'll want to retrieve the data of all the probes in the probe set that represents your favourite gene.

To plot the intensities of all the probes of a probe set, we will use ggplot(). It’s not the easiest way to go but ggplot() generates the nicest plots. To use ggplot() data has to be in the correct format (a data frame with informative column names) and categorical variables need to be transformed into factors.


BioC6.png


BioC7.png

More info on the ggplot() method.

Mortasecca.png Warning:
245027_at is a probe set from the Arabidopsis ATH1 array. This means that the code above will not work on other arrays. On other arrays you need to specify the ID of a probe set that is present on the array.