Question: Is it possible to convert Methylumi object into RGChannelSet
gravatar for poojitha.stemcell
22 months ago by
poojitha.stemcell0 wrote:

I am trying to integrate over 15 datasets into one RGChannelSet in order to normalise it using the package funtoonorm. I used GEOquery to get the datasets from GEO. However, since minfi only reads IDAT files, I could not get an RGset from GEO datasets that do not have IDAT files. Hence, I had to create a methylumi object from the series matrix text file from GEO. However, funtoonorm only works with minfi object. Would any of you know if there is a way of converting methylumi object into RGset? 

Thank you,


ADD COMMENTlink modified 11 months ago • written 22 months ago by poojitha.stemcell0
gravatar for James W. MacDonald
22 months ago by
United States
James W. MacDonald47k wrote:

Yes, there is. See the last few pages of the methylumi vignette.

ADD COMMENTlink written 22 months ago by James W. MacDonald47k

Thank you James:) I created a methylumiset object from series matrix in GEO and tried using the coercing function from methylumi vignette, but got the following error message:

> lumiMethy <- methylumiR("GSE50192_GPL13534_Matrix_raw_data.txt", sep="\t")

> rgSet450k <- as(mset450k, 'RGChannelSet') 

Error in methylumiToMinfi(lumiMethy) : 

  Cannot construct an RGChannelSet without full (OOB) intensities

Would you know how to resolve this? Is there a way of retaining OOB intensities while creating the methylumi object?







ADD REPLYlink written 22 months ago by poojitha.stemcell0

It depends on what data the submitter sent to GEO. In the case of this particular GSE, it appears they just sent the processed beta values, and you can't convert those data to an RGChannelSet. An RGChannelSet implies that you have the original Red and Green channel data, and you can't get those values from the betas alone.

ADD REPLYlink written 22 months ago by James W. MacDonald47k
gravatar for poojitha.stemcell
11 months ago by
poojitha.stemcell0 wrote:

Thank you James. The data i am now dealing with is a matrix of signal intensities from GSE35069. I still get the same error of missing OOB intensities. Would you know why this is happening?



ADD COMMENTlink written 11 months ago by poojitha.stemcell0

Its the same issue. When you get intensities only from GEO you do not get

  - control probes
  - out of band probes

To get this information, you generally need the original IDAT files. This means that any algorithm which makes use of either control probes or out of band probes will fail.

ADD REPLYlink written 11 months ago by Kasper Daniel Hansen6.4k
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.2.0
Traffic: 295 users visited in the last hour