Question: Reading Affymetrix data from using oligo for RNA degradation plot
gravatar for Vasu
3.0 years ago by
Vasu300 wrote:


I have some Affymetrix CEL files that use the HuGene 2.0 ST design.

I'm using oligo package to read CEL files. My question is how can I get RNA degradation plot using oligo? 


ADD COMMENTlink modified 3.0 years ago by Charles Warden6.1k • written 3.0 years ago by Vasu300
gravatar for Charles Warden
3.0 years ago by
Charles Warden6.1k
Duarte, CA
Charles Warden6.1k wrote:

You can get probe sequences and coordinates (usually) from the Affy website:

You'll want to extract the probe-level summarization and you can then plot intensity as a function of position with your genes of interest.

With respect to the probe-level summarization from the oligo package, the code should look something like this:


celFiles = list.celfiles("/path/to/folder", full.names=T)
rawData = read.celfiles(celFiles)
probesetData = rma(rawData, target="probeset")
probeset.expression = exprs(probesetData)
probesets = rownames(probesetData)

However, selecting the optimal genes may be tricky. I've tried using the validated genes from this paper, but they didn't seem to show much of a position degradation effect (in the array data that I have checked so far):

ADD COMMENTlink modified 5 months ago by RamRS20k • written 3.0 years ago by Charles Warden6.1k

Hi Charles,

Thanks for the reply. I have used "xps" package for this.

And I got the following plot. Help Needed to Interpret RNA degradation plot.

Thank you in Advance

ADD REPLYlink modified 5 months ago by RamRS20k • written 2.9 years ago by Vasu300

I'm not exactly sure how to interpret the figure - I would have expected more of a 5' degradation bias. At least the intensities are relatively consistent within each gene, so I would say the signal isn't too noisy.

ADD REPLYlink written 2.9 years ago by Charles Warden6.1k
Please log in to add an answer.


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