This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SNP‐based enrichment analysis for WGCNA modules using an hypergeometric test

Hello everyone

I have 13 clusters of co-expressed genes, obtained form an expression dataset of 6800 genes per 27 samples. Among the 6800 genes used for the netowrk analysis, 168 have a single SNPs (compared to the WT strain). I was wondering if an hypergeometric test can applied to find modules enriched with genes having a SNP:

Number of genes with SNP in the turquoise module 17
Total number of genes in the turquoise module 889
Total number of genes in the network 6800
Total number of genes with SNP 168

phyper(17-1,889,6800,168, lower.tail = TRUE, log.p = FALSE) #edited

Does it make any sense?

Thank you

hypergeometrictest snps wgcna
• 1,785 views • link
written by andres.firrincieli andres.firrincieli

1 answer

I think that it is [almost] correct. Take a look at THIS same question on StackOverflow where it would end up being:

phyper(16, 889, 6800 - 889, 168, lower.tail = TRUE, log.p = FALSE)

Note, also, whether it is 16 or 17? Please re-check your question.

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this