Comments (4)
I've done a GWAS. And so far only have one hit. Now I want to fine map it. Would it work with this package considering the above? What is your advice. Do you know of alternatives?
from fgwas.
I don't think that any fine mapping / annotation enrichment methods are likely to work well when you have such a limited number of associations. If you think there are lots of sub-threshold associations you could still try. A danger is that if you test fgwas with lots of different annotations then you are likely to overfit your data, and the "best" annotation might not truly reflect enrichment of causal variants.
For fine mapping just one association, I would suggest instead looking at your credible set of SNPs (e.g. those that comprise 95% or 99% of the statistical probability of association without any annotations), and seeing what overlaps you find in HaploReg (http://archive.broadinstitute.org/mammals/haploreg/haploreg.php). (If your SNP positions are in hg37 you need to use an older version of HaploReg.)
This doesn't estimate a posterior probability using the annotation data, but such an estimate would be unreliable in an underpowered GWAS anyway.
from fgwas.
Thanks for the reply. To be clear, with one hit, I meant: one locus with a leadSNP, obvious there are more variants for that locus at genome-wide significance. Does it then work? Or still not?
from fgwas.
Yes, one associated locus is very little data from which to estimate annotation enrichments. It might give you an answer - I just would be skeptical how reliable the answer is. You would still want to use your genome-wide data, as sub-threshold associations might contribute something, and would make the estimate more robust (the prevalence of non-associated SNPs in an annotation is important for estimating enrichment of causal SNPs in the annotation).
from fgwas.
Related Issues (10)
- Dealing with polymorphism HOT 2
- Annotation enrichments from fgwas - what is the base of the logartihm? HOT 2
- 'std::bad_alloc' HOT 2
- current segment is 70755198 73174848, position is 7 HOT 2
- Interpretation of log likelihood
- NCASE and NCONTROL are not available in the summary statistics of a case-control study
- Input bayes factors HOT 1
- multithreading
- Average over SNP priors HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from fgwas.