Meaning of the value column in bumphunter output (minfi package)
2
0
Entering edit mode
Antonio M. • 0
@antonio-m-13072
Last seen 9 months ago
Spain

Hi everyone.

In the documentation, the value column in minfi::bumphunter output is the average difference level of methylation in the bump. My question is, does this mean that the methylation levels are greater in for example, controls? (If we are comparing control vs case) .

How can I know the meaning of this measure in this case?

 

Thank you all for your attention

minfi bumphunter • 2.4k views
ADD COMMENT
0
Entering edit mode
@kasper-daniel-hansen-2979
Last seen 9 months ago
United States
That depends on your design matrix On Wed, May 17, 2017 at 10:28 AM, Antonio M. [bioc] < noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User Antonio M. <https: support.bioconductor.org="" u="" 13072=""/> wrote Question: > Meaning of the value column in bumphunter output (minfi package) > <https: support.bioconductor.org="" p="" 96046=""/>: > > Hi everyone. > > In the documentation, the value column in minfi::bumphunter output is the > average difference level of methylation in the bump. My question is, does > this mean that the methylation levels are greater in for example, controls? > (If we are comparing control vs case) . > > How can I know the meaning of this measure in this case? > > > > Thank you all for your attention > > ------------------------------ > > Post tags: minfi, bumphunter > > You may reply via email or visit Meaning of the value column in bumphunter output (minfi package) >
ADD COMMENT
0
Entering edit mode

I have this.

g <- as.factor(c(rep("c", 17), rep("t", 40))) # c for controls vs t for tumors
matriz <- model.matrix(~ g) # The design matrix

The output after bumphunter:

head(dmrs$table)
       chr     start       end     value     area cluster indexStart indexEnd  L clusterL p.value fwer p.valueArea fwerArea
2685 chr18  74961727  74962672 0.2516403 3.271323   88118     424245   424257 13       35       0    0           0        0

Thanks for your help

ADD REPLY
1
Entering edit mode

The way you specified your model matrix will by default set the second coefficient to be tumor - control. You can test that by looking at your factor levels (the first factor level is by default the baseline level) or by looking at the design matrix itself.
 

ADD REPLY
0
Entering edit mode

Thank you very much.

ADD REPLY
0
Entering edit mode
src3 • 0
@src3-15184
Last seen 6.1 years ago

Hi all,

I was wondering how average difference in methylation could be > 1 if beta values range from 0 to 1? 

Thank you!

ADD COMMENT

Login before adding your answer.

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