Covert MethylSet to RGChannelSet object
1
0
Entering edit mode
metamaden ▴ 10
@madensean-8348
Last seen 13 months ago
United States

Hi there,

I'm interested in converting a MethylSet back to an RGChannelSet object for purposes of preprocessing (I am using unprocessed data obtained from GEO), but I'm unsure of how to do this.

My understanding is that preprocessRaw() uses manifest info to associate idat channels with the correct array probes. I'm not aware of a function that does this in reverse. I'm thinking it should be possible to do this manually with a few for loops, but the rownames for getGreen() and getRed() matrices from the RGChannelSet are numeric identifiers. I am unsure of how to go about identifying which channels go to which array probe.

I am running minfi 1.18.6 and Bioconductor 3.3.

Thanks!

Sean

minfi hm450 • 2.3k views
ADD COMMENT
2
Entering edit mode
@kasper-daniel-hansen-2979
Last seen 18 months ago
United States
You can't do this. For example, you have no knowledge of the so-called out of band probes, which are probes measured in the "other" color channel than they were designed for. You are also missing control probes measurements. Perhaps it is easier for you to look at the dimensions of the RGsetEx and MsetEx in the minfiData package to see that you are missing a great deal of numbers. Best, Kasper On Wed, Aug 31, 2016 at 2:53 PM, maden.sean [bioc] <noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User maden.sean <https: support.bioconductor.org="" u="" 8348=""/> wrote Question: > Covert MethylSet to RGChannelSet object > <https: support.bioconductor.org="" p="" 86673=""/>: > > Hi there, > > I'm interested in converting a MethylSet back to an RGChannelSet object > for purposes of preprocessing (I am using unprocessed data obtained from > GEO), but I'm unsure of how to do this. > > My understanding is that preprocessRaw() uses manifest info to associate > idat channels with the correct array probes. I'm not aware of a function > that does this in reverse. I'm thinking it should be possible to do this > manually with a few for loops, but the rownames for getGreen() and getRed() > matrices from the RGChannelSet are numeric identifiers. I am unsure of how > to go about identifying which channels go to which array probe. > > I am running minfi 1.18.6 and Bioconductor 3.3. > > Thanks! > > Sean > > ------------------------------ > > Post tags: minfi, hm450 > > You may reply via email or visit Covert MethylSet to RGChannelSet object >
ADD COMMENT
0
Entering edit mode

Kasper,

Thanks for your reply. I was hoping that, despite the missing data, I would be able to fill in the limited data I have in order to perform normalizations that don't require the other kinds of probes and intensities. 

Sean

ADD REPLY
1
Entering edit mode
So it is theoretically possible to fill in the matrices with some numbers and lots of NAs. In general, we have not assumed the presence of NAs in the code, so that could very well end badly. If you still want to try, you need to understand the data structures in the manifest package. On Wed, Aug 31, 2016 at 3:13 PM, maden.sean [bioc] <noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User maden.sean <https: support.bioconductor.org="" u="" 8348=""/> wrote Comment: > Covert MethylSet to RGChannelSet object > <https: support.bioconductor.org="" p="" 86673="" #86675="">: > > Kasper, > > Thanks for your reply. I was hoping that, despite the missing data, I > would be able to fill in the limited data I have in order to perform > normalizations that don't require the other kinds of probes and > intensities. > > Sean > > ------------------------------ > > Post tags: minfi, hm450 > > You may reply via email or visit https://support.bioconductor. > org/p/86673/#86675 >
ADD REPLY

Login before adding your answer.

Traffic: 944 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6