missMethyl - get names of differentially methylated genes from gometh()
1
0
Entering edit mode
stewart999 ▴ 10
@stewart999-11769
Last seen 6 months ago
United Kingdom

I'm using the gometh() function of missMethyl and am trying to pull out the names of the differentially methylated genes from the output table, which looks like this:

go <- missMethyl::gometh(sig.cpg=sig.cpg, all.cpg=all.cpg, array.type="EPIC", collection="GO")
head(go)
           ID   ONT                                             TERM   N DE         P.DE          FDR
1: GO:0015748    BP                  organophosphate ester transport 116  5 4.907299e-08 0.0002769925
2: GO:0015914    BP                           phospholipid transport  90  5 1.505420e-08 0.0002769925
3: GO:0034382    BP                    chylomicron remnant clearance   9  3 3.849721e-08 0.0002769925
4: GO:0071830    BP triglyceride-rich lipoprotein particle clearance   9  3 3.849721e-08 0.0002769925
5: GO:0033700    BP                              phospholipid efflux  12  3 6.832103e-08 0.0003085104
6: GO:0008203    BP                    cholesterol metabolic process 146  5 1.082826e-07 0.0003628179

The DE column is the number of genes that are differentially methylated. How do I find out, for example, the names of the 5 differentially methylated genes for GO:0015748?

Thanks, Stewart

missMethyl • 254 views
ADD COMMENT
0
Entering edit mode
@jovana-maksimovic-4422
Last seen 10 months ago
Australia

Hi Stewart,

We recently implemented a parameter in gometh, sig.genes, allowing you to output an extra column in the results table showing the genes in the gene set that overlap the DM CpGs. It is set to FALSE by default.

Ensure that you have the latest version of missMethyl installed to be able to use it.

Cheers, Jovana

ADD COMMENT

Login before adding your answer.

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