j2s: Deriving stratified effects from joint models investigating Gene-Environment Interactions
======
The python3 script **j2s.py** allows for the estimation of genetic effect sizes in unexposed and exposed individuals separately from joint models investigating Gene-Environment interactions when the exposure is binary. A pre-print of the publication is available [here](https://www.biorxiv.org/content/10.1101/693218v1).
Prerequisite
------------
To execute the script, python3 must be installed as weel as the following Python packages:
+ pandas (version 0.22.0)
+ numpy (version 1.13.3)
+ scipy (version 0.19.1)
Execution
------------
To use the script, type the following command:
``` bash
python3 INFILE N N_EXPOSED OUTFILE
```
with:
* INFILE is the path to the input file describing the summary statitics in the joint model,
* N is the total sample size,
* N_EXPO is the number of exposed (E = 1) individuals,
* OUTFILE is the path to the outputfile.
Description of the input files
------------
The input file has 7 mandatory columns:
* the identifier of the variant (e.g rs number) labelled 'MarkerName'
* the main genetic effect size labelled 'Effect'
* the standard error of the main genetic effect size labelled 'StdErr'
* the interaction effect size labelled 'IntEffect'
* the standard error of the interaction effect size labelled 'IntStdErr'
* the covariance between the main genetic effect size and the interaction effect size labelled 'IntCov'
This format corresponds to the output of the METAL software performing the joint test ([Manning et al, 2011](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3312394/))
Short tutorial
--------------
Bug report / Help
-----------------
Please open an issue if you find a bug.
Code of conduct
---------------
Please note that this project is released with a [Contributor Code of Conduct](https://gitlab.pasteur.fr/statistical-genetics/j2s/blob/master/code-of-conduct.md). By participating in this project you agree to abide by its terms.
License
-------
This project is licensed under the MIT License - see the [LICENSE.md](https://gitlab.pasteur.fr/statistical-genetics/j2s/blob/master/LICENSE) file for details