This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem in processing microarray data with HuEx-1_0-st-v1 platform

When I am trying to process a micro-array data with Array Express id E-GEOD-21298, I get the following error:

Code Run:

  >celFile<-list.celfiles(full.names=TRUE)
  >rawDat<- read.celfiles(celFile)
  >norDat<- oligo::rma(rawDat,background=TRUE,normalize=TRUE,target="core")

Platform design info loaded.
Reading in : ./GSM532260.CEL
Reading in : ./GSM532261.CEL
Reading in : ./GSM532262.CEL
Reading in : ./GSM532263.CEL
Reading in : ./GSM532264.CEL
Reading in : ./GSM532265.CEL
Reading in : ./GSM532274.CEL
Reading in : ./GSM532277.CEL
Reading in : ./GSM532278.CEL
Reading in : ./GSM532279.CEL
[1] "pd.huex.1.0.st.v2"
Background correcting
Normalizing
Calculating Expression
  **Error in (function (classes, fdef, mtable)  : 
      unable to find an inherited method for function ‘columns’ for signature ‘"AffyExonPDInfo"’**

I am unable to understand this error. Any help with this error ?

Also, this array belongs to HuEx-1_0-st-v1 while the current available annotation package in Bioconductor is HuEx-1_0-st-v2 . [Which might be the reason of above error] Running the command:

biocLite("pd.huex.1.0.st.v1")

gives error: 

package ‘pd.huex.1.0.st.v1’ is not available (for R version 3.4.4)

So how do I annotate the processed files for the older version ?

microarray r sequencing bioconductor next-gen
• 3,482 views • link
updated by Kevin Blighe Kevin Blighe • written by ashuein • ashuein

1 answer

There have indeed been issues with this array's annotation package. There is further discussion here: pd.huex.1.0.st.v2 package doubt.

What you should do is download the CDF from GenoMax's first link - you will require a login for Affymetrix and then, also, to unzip the downloaded file. Then set it up as a custom CDF by following this tutorial: How To Install A Specific Cdf Package For Use In Bioconductor

Once you have it set up, you can then explicitly specify it like this:

project <- read.celfiles(CELFiles, package="HuEx-1_0-st-v2-cdf")

If you don't specify 'package', here, then the function will just attempt to download the erronous package yet again. It even failed on my computer.

By the way, I'm not sure why you were trying to install via biocLite() 'pd.huex.1.0.st.v1' when these arrays are actually 'pd.huex.1.0.st.v2'.

Kevin

• 1 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this