Question: Error while running affylmGUI
0
gravatar for devyani66
6 months ago by
devyani660
devyani660 wrote:

Hello, I am trying to run the affylmGUI R package and I get this error on CDF env during the normalization step. Below is the link: CDF env Error

Any idea what might be the reason? I googled it but didnt find much help. Thanks in advance.

microarray affy limma affylmgui • 184 views
ADD COMMENTlink modified 6 months ago by Gordon Smyth38k • written 6 months ago by devyani660
Answer: Error while running affylmGUI
1
gravatar for Gordon Smyth
6 months ago by
Gordon Smyth38k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth38k wrote:

The error message tells you what the problem is. You need install the rmhu01aa520485cdf CDF package so that affylmGUI knows how your arrays are designed.

The problem is that you are trying to analyse data from a custom Rosetta-Affymetrix platform and Bioconductor doesn't have a standard CDF package for this array: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL14189

A similar question was asked about 5 years ago and James MacDonald explained how to get the package: https://support.bioconductor.org/p/59723

This question and response from 8 years ago is also relevant: https://support.bioconductor.org/p/40455/

I wonder what you Googled? I searched for "rmhu01aa520485cdf" and James' answer came up straight away.

ADD COMMENTlink modified 6 months ago • written 6 months ago by Gordon Smyth38k

Thank you for your reply. My bad I didnt mention that I did go through the above link where MacDonald explained how to get the package. I had tried everything mentioned in and still, it used to give the same error.

Now going through again I realized the problem was with the name of the cdf file. If you see the discussion, in many places rmhu01 is mentioned as rmhu10. That simple change of 01 to 10 was the issue as even though my env was created but it was rmhu10aa520485cdf and not rmhu01aa520485cdf- the one its searching for. I had just copied the command from there and changed the file name (.gz) and kept the rest same...so this silly mistake. But to be honest that discussion thread is a big confusing one...with many things repeated and hard to understand which response is for which query.

But since you pointed again, I went back and realized this difference. Thank you.

ADD REPLYlink modified 6 months ago • written 6 months ago by devyani660

Hi Mr.Gordon, I have one more query regarding this. So, I ran these steps before carrying out normalization using affylmGUI.

> library(makecdfenv)
> make.cdf.package("GPL10687_RM-HU01Aa520485_custom_MMPM.cdf.gz",
> "rmhu01aa520485cdf", compress=TRUE, species="Homo_sapiens")
> BiocManager::install("AnnotationDbi")
> install.packages("rmhu01aa520485cdf", repos=NULL, type="source")

I exported the normalized values but the file has just probe ids and no annotation. Is there something else too that I need to run to include annotation? Thanks in advance.

ADD REPLYlink modified 6 months ago by Gordon Smyth38k • written 6 months ago by devyani660

I checked you answer to a similar question 8 yrs ago- https://support.bioconductor.org/p/40823/

I am not sure if I understood well, but I suppose I have to install the proper annotation package. For me, while I tried to install "rmhu01aa520485cdf" file, it wrote I did not have "AnnotationDbi" package and so I downloaded it.

ADD REPLYlink written 6 months ago by devyani660

You need to create CDF and annotation packages for your custom arrays. I have no experience with that. I can only answer questions about affylmGUI.

I suggest you ask a new question with an informative title and tags so that the right people will see your question.

But beware that affyLMGUI was never intended to be used with custom arrays, so you might have to consider using the affy and limma packages directly instead of through the GUI.

ADD REPLYlink modified 6 months ago • written 6 months ago by Gordon Smyth38k

Thank you for your response. And so sorry for a super late reply. To keep it simple, I just mapped the probe ids from the normalized file to the annotation file using simple R command. It worked. Regards

ADD REPLYlink written 6 months ago by devyani660
Please log in to add an answer.

Help
Access

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