Comments (3)
Hi Anaïs,
There is actually nothing in the analysis pipeline specific to quantitative versus disease traits.
This part of the config file references trait types:
trait_files:
BC: "clumped.genome_wide_ash.bc.max_r2.tsv"
trait_types:
BC: "disease"
trait_type_abbrevs:
disease: "Disease"
The entries here about trait types don't affecting anything except some plots.
If you are comfortable then following the pipeline in example_trait to process the GWAS summary statistics then nothing else in the config file would need to be altered aside from adding the new traits.
Thank you for trying out our software. Please let me know if that solves things or if there is anything else I can help with or improve in the code.
Evan
from smilenfer.
Hi Evan,
Thanks for your reply! I ran the software for the continuous trait and have a couple more questions:
- Where can I find the p values that go with the figure 3 plot in the preprint? I can find the plot itself and the excel sheet that goes with it, but am unsure about what all the headers (ll_neut I2_stab ll_stab I1_dir ll_dir I1_full I2_full ll_full Ip_plei ll_plei I2_nplei nn_nplei ll_nplei) mean and how to transform the values to p-values.
- In what direction is the directional selection tested? Is it both ways or only trait increasing? Should I change the sign for the beta's for a trait where I want to know if there is directional selection lowering the trait?
from smilenfer.
Hi Anaïs,
The columns with ll_*
in the title are log likelihoods. From these you can can obtain p-values for a likelihood ratio test where the neutral model is treated a the null hypothesis. For instance: 2*(ll_stab-ll_neut)
has a Chi2 distribution with one degree of freedom. The script I used to generate figure 3 is here and you may use the plotting function. Let me know if you have any issues with it. Also please ignore the nplei
columns for now as there are some issues currently with how this is calculated.
Directional selection is tested in both ways so you don't need to change the sign and re-run anything. If I1_dir is negative this means the fitted model has negative selection against trait-increasing alleles and if it is positive there is negative selection against trait-decreasing alleles.
from smilenfer.
Related Issues (2)
- Changing threshold HOT 8
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 smilenfer.