Pathway analysis using rare-variant SKAT-O summary statistics
0
1
Entering edit mode
18 months ago
tacrolimus ▴ 140

Dear Biostars,

I have summary stats from a SAIGE-GENE analysis of a case-control cohort. This is a collapsing rare variant test and the output has a heading like this:

GENE      Pvalue
ABC         0.093 
LOL2       1.123e-78
KJD         0.192

etc

I would like to do a gene-set analysis looking at the enrichment of pathways and was wondering if there was a method that used summary stats like this only? I have the individual markers, their AFs,BETA's, Pvalues etc but I would rather not format them into something like MAGMA if I can avoid the extra effort.

I have found this: https://github.com/yaowuliu/ACAT/blob/master/doc/ACAT_manual.pdf and was wondering if anyone had any experience with similar tools?

Many thanks for your help

variant rare pathway analysis R • 553 views
ADD COMMENT
2
Entering edit mode

You can use fgsea with the p-value ranked list of genes, and a simple hypergeometric test using a p-value cutoff. Tons of R packages are out there which implement pathway enrichment analysis; my personal choice is piano.

ADD REPLY

Login before adding your answer.

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