EGSEA with Affymetrix data
Entering edit mode
cverschoor • 0
Last seen 18 days ago

Hi Everyone,

I have been using EGSEA recently on my RNA seq datasets with no issues. Recently, I suggested a student try it with their Affymetrix data and hit a snag. We are pre-processing with RMA through minfi, filtering and using that as the expression matrix. We then generate a probe.annotation data frame including the probe IDs and entrez IDs from the array annotation file. Design and contrast matrices and the group vector are specified as suggested by the authors. When running the function we get the following error...

Error in dimnames(x[[a]]) <- value : 
  length of 'dimnames' [1] not equal to array extent

I am failrly certain the problem is with the expression matrix or probe.annotation frame since the error comes up with this stripped down version..., group=t, probe.annot = probe.annot)

I can easily run this function with either of the examples the authors provide and can't for the life of me figure out how those objects differ from the ones I am trying to use. If anyone knows where in the source code this error is likely occuring or knows what the check is for, please let me know. I realize I haven't provided any reproducible code to work off of, so let me know what I can post to be more helplful.

Many thanks!


I figured out the issue.... I had too many NA and duplicated Entrez IDs in my probe.annotation object. Once I removed these and trimmed my expression matrix accordingly, everything ran fine. I should have thought of that in the first place!

EGSEA • 106 views
Entering edit mode
Last seen 1 day ago
United States

Just to clarify, you say you used minfi, but I assume that's a typo?

Login before adding your answer.

Traffic: 593 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6