Fix regions for bumphunter
2
0
Entering edit mode
@isbrorson-13349
Last seen 4.1 years ago

Hi,

I'm working with DNA methylation data obtained from the EPIC chip. I've previously used the R package Bumphunter for differentially methylated regions (DMRs). In my current project, I would like to specify a list of specific regions for bumphunter to analyse. In other words, I'm looking for a way to "fix" the regions in advance, and then tell Bumphunter to look at those specific regions - in stead of Bumphunter analysing DMRs genome wide.

Does anyone have any experience or know if it is possible to tweak bumphunter commands to restrict such regions?

I'm very grateful for any response and help to get me a step closer to solve this :)

Thanks for reading!

bumphunter DMR methylation minfi • 882 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 14 hours ago
United States

There is a 'cluster' argument to bumphunter that you can use to specify the regions. Not sure how that would work in practice if you wanted to ignore other regions. I would imagine you would need to subset your GenomicRatioSet to include only those regions as well. In other words, the cluster just tells bumphunter which CpGs are in a given cluster, not that certain clusters should be ignored.

Restricting to certain regions may adversely affect the inference provided by the p-values as well (which aren't necessarily easily interpreted as it stands). So caveat emptor?

ADD COMMENT
0
Entering edit mode
@isbrorson-13349
Last seen 4.1 years ago

I see what you're saying . Yup, seems like restricting such regions would not be optimal for what I have in mind. I'll give the 'cluster' argument a try, though not sure how that will work since I'd like to specify at least 1000 regions. I'll take a closer look at it, and maybe I need to rethink my line of reasoning .

Thank you for your input!

ADD COMMENT
0
Entering edit mode

If you want to make a comment, please use ADD COMMENT, rather than Add Answer (unless of course you are actually answering your own question).

ADD REPLY
0
Entering edit mode

Opps, will do next time! I'm new to this support site.

  • Ina
ADD REPLY

Login before adding your answer.

Traffic: 627 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