genal-python 0.6__tar.gz → 0.8__tar.gz

{genal_python-0.6 → genal_python-0.8}/.gitignore

@@ -2,5 +2,4 @@ __pycache__/
 dist/
 .ipynb_checkpoints/
 ipynb_checkpoints/
-genal/.ipynb_checkpoints/
-docs/
+genal/.ipynb_checkpoints/

genal_python-0.8/.readthedocs.yaml

@@ -0,0 +1,22 @@
+# .readthedocs.yaml
+# Read the Docs configuration file
+# See https://docs.readthedocs.io/en/stable/config-file/v2.html for details
+
+# Required
+version: 2
+
+# Set the version of Python and other tools you might need
+build:
+  os: ubuntu-22.04
+  tools:
+    python: "3.11"
+
+# Build documentation in the docs/ directory with Sphinx
+sphinx:
+  configuration: docs/source/conf.py
+
+# We recommend specifying your dependencies to enable reproducible builds:
+# https://docs.readthedocs.io/en/stable/guides/reproducible-builds.html
+python:
+  install:
+  - requirements: docs/requirements.txt

genal_python-0.6/README.md → genal_python-0.8/PKG-INFO

@@ -1,3 +1,27 @@
+Metadata-Version: 2.1
+Name: genal-python
+Version: 0.8
+Summary: A python toolkit for polygenic risk scoring and mendelian randomization.
+Author-email: Cyprien Rivier <riviercyprien@gmail.com>
+Requires-Python: >=3.7
+Description-Content-Type: text/markdown
+Classifier: Programming Language :: Python :: 3
+Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)
+Classifier: Operating System :: OS Independent
+Requires-Dist: aiohttp==3.9.5
+Requires-Dist: nest_asyncio==1.5.5
+Requires-Dist: numpy>=1.24.4, <2.0
+Requires-Dist: pandas>=2.0.3
+Requires-Dist: plotnine==0.12.3
+Requires-Dist: psutil==5.9.1
+Requires-Dist: pyliftover==0.4
+Requires-Dist: scikit_learn>=1.3.0
+Requires-Dist: scipy>=1.11.4
+Requires-Dist: statsmodels==0.14.0
+Requires-Dist: tqdm==4.66.1
+Requires-Dist: wget==3.2
+Project-URL: Home, https://github.com/CypRiv/genal
+
 <center><h1> genal: A Python Toolkit for Genetic Risk Scoring and Mendelian Randomization </h1></center>
 
 
@@ -7,10 +31,11 @@
 
 # Table of contents
 1. [Introduction](#introduction)
-2. [Citation] (#citation)
+2. [Citation](#citation)
 3. [Requirements for the genal module](#paragraph1)
 4. [Installation and how to use genal](#paragraph2)
     1. [Installation](#paragraph2.1)
+    2. [Documentation](#paragraph2.2)
 5. [Tutorial and presentation of the main tools](#paragraph3)
     1. [Data loading](#paragraph3.1)
     2. [Data preprocessing](#paragraph3.2)
@@ -19,6 +44,7 @@
     5. [Mendelian Randomization](#paragraph3.5)
     6. [SNP-association testing](#paragraph3.6)
     7. [Lifting](#paragraph3.7)
+    8. [GWAS Catalog](#paragraph3.8)
 
 
 ## Introduction <a name="introduction"></a>
@@ -55,6 +81,16 @@ Once downloaded, the path to the plink executable can be set with:
 ```
 genal.set_plink(path="/path/to/plink/executable/file")
 ```
+### Documentation <a name="paragraph2.2"></a>
+
+For detailed information on how to use the functionalities of Genal, please refer to the documentation: https://genal.rtfd.io
+
+The documentation covers:
+- Installation 
+- This tutorial
+- The list of the main functions with complete description of their arguments
+- An exhaustive API reference
+
 
 ## Tutorial <a name="paragraph3"></a>
 For this tutorial, we will obtain genetic instruments for systolic blood pressure (SBP), compute a Polygenic Risk Score (PRS), and run a Mendelian Randomization analysis to investigate the genetically-determined effect of SBP on the risk of stroke. We will utilize summary statistics from Genome-Wide Association Studies (GWAS) and individual-level data from the UK Biobank. The steps include:
@@ -75,11 +111,11 @@ For this tutorial, we will obtain genetic instruments for systolic blood pressur
   - Data lifting to another genomic build
     - In pure Python
     - Using LiftOver
-  - Phenoscanner (to be added)
+  - Querying the GWAS Catalog
 
 ### Data loading <a name="paragraph3.1"></a>
 
-We begin with publicly available summary statistics from a large GWAS study of systolic blood pressure. [Link to study](https://www.nature.com/articles/s41588-018-0205-x). After downloading and unzipping the summary statistics, we load them into a pandas DataFrame:
+We start this tutorial with publicly available summary statistics from a large GWAS study of systolic blood pressure. [Link to study](https://www.nature.com/articles/s41588-018-0205-x). After downloading and unzipping the summary statistics, we load them into a pandas DataFrame:
 
 ```python
 import pandas as pd
@@ -108,6 +144,10 @@ The `genal.Geno` takes as input a pandas dataframe where each row corresponds to
 - **P**: Column name for effect p-value. Defaults to `'P'`.
 - **EAF**: Column name for effect allele frequency. Defaults to `'EAF'`.
 
+> **Note:**
+> 
+> You do not need all columns to move forward, as not all columns are required by every function. Additionally, some columns can be imputed as we will see in the next paragraph.
+
 After inspecting the dataframe, we first need to extract the chromosome and position information from the `MarkerName` column into two new columns `CHR` and `POS`:
 
 ```python
@@ -133,7 +173,7 @@ The last argument (`keep_columns = False`) indicates that we do not wish to keep
 
 > **Note:**
 > 
-> Make sure to read the readme file usually provided with the summary statistics to identify the correct columns. It is particularly important to correctly identify the allele that represents the effect allele. Also, you do not need all columns to move forward, as some can be inputted as we will see next.
+> Make sure to read the readme file usually provided with the summary statistics to identify the correct columns. It is particularly important to correctly identify the allele that represents the effect allele.
 
 ### Data preprocessing <a name="paragraph3.2"></a>
 
@@ -312,7 +352,7 @@ and the output is:
     The PRS computation was successful and used 1330/1538 (86.476%) SNPs.
     PRS data saved to SBP_prs.csv
 
-In our case, we have been able to find proxies for 571 of the 786 SNPs that were missing in the population genetic data (7 potential proxies have been removed because they were identical to SNPs already present in our data).
+In our case, we have been able to find proxies for 578 of the 786 SNPs that were missing in the population genetic data (7 potential proxies have been removed because they were identical to SNPs already present in our data).
 
 You can customize how the proxies are chosen with the following arguments:
 - `reference_panel`: The reference population used to derive linkage disequilibrium values and find proxies. Defaults to `eur`.
@@ -322,7 +362,7 @@ You can customize how the proxies are chosen with the following arguments:
 
 > **Note:**
 > 
-> You can call the `genal.Geno.prs` method on any `Geno` instance (containing at least the EA, BETA, and either SNP or CHR/POS columns). The data does not need to be clumped, and there is no limit to the number of instruments used to compute the scores.
+> You can call the `genal.Geno.prs` method on any `Geno` instance (containing at least the EA, BETA, and either SNP or CHR/POS columns). The data does not need to be clumped, and there is no limit to the number of SNPs used to compute the scores.
 
 
 ### Mendelian Randomization <a name="paragraph3.5"></a>
@@ -437,7 +477,7 @@ If you want to visualize the obtained MR results, you can use the `genal.Geno.MR
 SBP_clumped.MR_plot(filename="MR_plot_SBP_AS")
 ```
 
-![MR plot](docs/Images/MR_plot_SBP_AS.png)
+![MR plot](docs/build/_images/MR_plot_SBP_AS.png)
 You can select which MR methods you wish to plot with the `methods` argument. Note that for an MR method to be plotted, they must be included in the latest `genal.Geno.MR` call of this `genal.Geno` instance.
 
 If you wish to include the heterogeneity values (Cochran's Q) in the results, you can use the heterogeneity argument in the `genal.Geno.MR` call. Here, the heterogeneity for the inverse-variance weighted method:
@@ -457,7 +497,7 @@ As expected, many MR methods indicate that SBP is strongly associated with strok
 To investigate horizontal pleiotropy in more details, a very useful method is Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO). MR-PRESSO is a method designed to detect and correct for horizontal pleiotropy. It will identify which instruments are likely to be pleiotropic on their effect on the outcome, and it will rerun an inverse-variance weighted MR after excluding them. It can be run using the `genal.Geno.MRpresso` method:
 
 ```python
-SBP_clumped.MRpresso(action = 2, n_iterations = 30000)
+mod_table, GlobalTest, OutlierTest, BiasTest = SBP_clumped.MRpresso(action = 2, n_iterations = 30000)
 ```
 
 As with the `genal.Geno.MR` method, the `action` argument determines how the pleiotropic SNPs will be treated. The output is a list containing:
@@ -482,7 +522,7 @@ df_pheno = pd.read_csv("path/to/trait/data")
 
 > **Note:**
 > 
-> One important detail is to make sure that the individual IDs are identical between the phenotypic data and the genetic data for the target population.
+>    One important detail is to make sure that the IDs of the participants are identical in the phenotypic data and in the genetic data.
 
 Then, it is advised to make a copy of the `genal.Geno` instance containing our instruments as we are going to update their coefficients and to avoid any confusion:
 
@@ -560,5 +600,54 @@ You can specify the path of the LiftOver executable to the `liftover_path` argum
 SBP_Geno.lift(start = "hg19", end = "hg38", replace = False, liftover_path = "path/to/liftover/exec")
 ```
 
+### GWAS Catalog <a name="paragraph3.8"></a>
 
+It is sometimes interesting to determine the traits associated with our SNPs. In Mendelian Randomization, for instance, we may want to exclude instruments that are associated with traits likely causing horizontal pleiotropy. For this purpose, we can use the `genal.Geno.query_gwas_catalog` method. This method will query the GWAS Catalog API to determine the list of traits associated with each of our SNPs and store the results in a list in the `ASSOC` column of the `.data` attribute:
 
+```python
+SBP_clumped.query_gwas_catalog(p_threshold=5e-8)
+```
+Which will output:
+
+    Querying the GWAS Catalog and creating the ASSOC column. 
+    Only associations with a p-value <= 5e-08 are reported. Use the p_threshold argument to change the threshold.
+    To report the p-value for each association, use return_p=True.
+    To report the study ID for each association, use return_study=True.
+    The .data attribute will be modified. Use replace=False to leave it as is.
+    100%|██████████| 1545/1545 [00:34<00:00, 44.86it/s]
+    The ASSOC column has been successfully created.
+    701 (45.37%) SNPs failed to query (not found in GWAS Catalog) and 7 (0.5%) SNPs timed out after 34.33 seconds. You can increase the timeout value with the timeout argument.
+| EA  | NEA | EAF   | BETA   | SE     | CHR | POS        | SNP        | ASSOC                                                                 |
+|-----|-----|-------|--------|--------|-----|------------|------------|------------------------------------------------------------------------|
+| A   | G   | 0.1784| 0.2330 | 0.0402 | 10  | 102075479  | rs603424   | [eicosanoids measurement, decadienedioic acid (...]                     |
+| A   | G   | 0.0706| -0.3873| 0.0626 | 10  | 102403682  | rs2996303  | FAILED_QUERY                                                           |
+| T   | G   | 0.8872| 0.6846 | 0.0480 | 10  | 102553647  | rs1006545  | [diastolic blood pressure, systolic blood pressure...]                  |
+| T   | G   | 0.6652| -0.2098| 0.0340 | 10  | 102558506  | rs12570050 | FAILED_QUERY                                                           |
+| T   | C   | 0.3057| -0.2448| 0.0334 | 10  | 102603924  | rs4919502  | FAILED_QUERY                                                           |
+| ... | ... | ...   | ...    | ...    | ... | ...        | ...        | ...                                                                    |                                          |
+| T   | C   | 0.3514| 0.2203 | 0.0314 | 9   | 9350706    | rs1332813  | [diastolic blood pressure, systolic blood pressure...]                  |
+| T   | C   | 0.6880| -0.1897| 0.0332 | 9   | 94201341   | rs10820855 | FAILED_QUERY                                                           |
+| A   | T   | 0.3669| -0.1862| 0.0313 | 9   | 95201540   | rs7045409  | [protein measurement, pulse pressure measurement...]                   |
+
+If you are also interested in the p-values of each SNP-trait association, or the ID of the study from which the association was reported, you can use the `return_p = True` and `return_study = True` arguments. Then, the `ASSOC` column will contain a list of tuples, where each tuple contains the trait name, the p-value, and the study ID:
+
+```python
+SBP_clumped.query_gwas_catalog(p_threshold=5e-8, return_p=True, return_study=True)
+```
+
+| EA  | NEA | EAF   | BETA   | SE     | CHR | POS        | SNP        | ASSOC                                                                 |
+|-----|-----|-------|--------|--------|-----|------------|------------|------------------------------------------------------------------------|
+| A   | G   | 0.1784| 0.2330 | 0.0402 | 10  | 102075479  | rs603424   | TIMEOUT                                                                |
+| A   | G   | 0.0706| -0.3873| 0.0626 | 10  | 102403682  | rs2996303  | FAILED_QUERY                                                           |
+| T   | G   | 0.8872| 0.6846 | 0.0480 | 10  | 102553647  | rs1006545  | [(heart rate response to exercise, 6e-12, GCST...                      |
+| T   | G   | 0.6652| -0.2098| 0.0340 | 10  | 102558506  | rs12570050 | FAILED_QUERY                                                           |
+| T   | C   | 0.3057| -0.2448| 0.0334 | 10  | 102603924  | rs4919502  | FAILED_QUERY                                                           |
+| ... | ... | ...   | ...    | ...    | ... | ...        | ...        | ...                                                                    |                                                         |
+| T   | C   | 0.3514| 0.2203 | 0.0314 | 9   | 9350706    | rs1332813  | [(diastolic blood pressure, 1e-12, GCST9031029...                      |
+| T   | C   | 0.6880| -0.1897| 0.0332 | 9   | 94201341   | rs10820855 | FAILED_QUERY                                                           |
+| A   | T   | 0.3669| -0.1862| 0.0313 | 9   | 95201540   | rs7045409  | [(systolic blood pressure, 9e-13, GCST006624),...                      |
+
+
+> **Note:**
+> 
+> As you can see, many SNPs failed to be queried. This is normal as the GWAS Catalog is not exhaustive.

genal_python-0.6/PKG-INFO → genal_python-0.8/README.md

@@ -1,25 +1,3 @@
-Metadata-Version: 2.1
-Name: genal-python
-Version: 0.6
-Summary: A python toolkit for polygenic risk scoring and mendelian randomization.
-Author-email: Cyprien Rivier <riviercyprien@gmail.com>
-Requires-Python: >=3.7
-Description-Content-Type: text/markdown
-Classifier: Programming Language :: Python :: 3
-Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)
-Classifier: Operating System :: OS Independent
-Requires-Dist: numpy>=1.26.2
-Requires-Dist: pandas>=2.0.3
-Requires-Dist: plotnine>=0.12.3
-Requires-Dist: psutil>=5.9.1
-Requires-Dist: pyliftover>=0.4
-Requires-Dist: scikit_learn>=1.3.0
-Requires-Dist: scipy>=1.11.3
-Requires-Dist: statsmodels>=0.14.0
-Requires-Dist: tqdm>=4.66.1
-Requires-Dist: wget>=3.2
-Project-URL: Home, https://github.com/CypRiv/genal
-
 <center><h1> genal: A Python Toolkit for Genetic Risk Scoring and Mendelian Randomization </h1></center>
 
 
@@ -29,10 +7,11 @@ Project-URL: Home, https://github.com/CypRiv/genal
 
 # Table of contents
 1. [Introduction](#introduction)
-2. [Citation] (#citation)
+2. [Citation](#citation)
 3. [Requirements for the genal module](#paragraph1)
 4. [Installation and how to use genal](#paragraph2)
     1. [Installation](#paragraph2.1)
+    2. [Documentation](#paragraph2.2)
 5. [Tutorial and presentation of the main tools](#paragraph3)
     1. [Data loading](#paragraph3.1)
     2. [Data preprocessing](#paragraph3.2)
@@ -41,6 +20,7 @@ Project-URL: Home, https://github.com/CypRiv/genal
     5. [Mendelian Randomization](#paragraph3.5)
     6. [SNP-association testing](#paragraph3.6)
     7. [Lifting](#paragraph3.7)
+    8. [GWAS Catalog](#paragraph3.8)
 
 
 ## Introduction <a name="introduction"></a>
@@ -77,6 +57,16 @@ Once downloaded, the path to the plink executable can be set with:
 ```
 genal.set_plink(path="/path/to/plink/executable/file")
 ```
+### Documentation <a name="paragraph2.2"></a>
+
+For detailed information on how to use the functionalities of Genal, please refer to the documentation: https://genal.rtfd.io
+
+The documentation covers:
+- Installation 
+- This tutorial
+- The list of the main functions with complete description of their arguments
+- An exhaustive API reference
+
 
 ## Tutorial <a name="paragraph3"></a>
 For this tutorial, we will obtain genetic instruments for systolic blood pressure (SBP), compute a Polygenic Risk Score (PRS), and run a Mendelian Randomization analysis to investigate the genetically-determined effect of SBP on the risk of stroke. We will utilize summary statistics from Genome-Wide Association Studies (GWAS) and individual-level data from the UK Biobank. The steps include:
@@ -97,11 +87,11 @@ For this tutorial, we will obtain genetic instruments for systolic blood pressur
   - Data lifting to another genomic build
     - In pure Python
     - Using LiftOver
-  - Phenoscanner (to be added)
+  - Querying the GWAS Catalog
 
 ### Data loading <a name="paragraph3.1"></a>
 
-We begin with publicly available summary statistics from a large GWAS study of systolic blood pressure. [Link to study](https://www.nature.com/articles/s41588-018-0205-x). After downloading and unzipping the summary statistics, we load them into a pandas DataFrame:
+We start this tutorial with publicly available summary statistics from a large GWAS study of systolic blood pressure. [Link to study](https://www.nature.com/articles/s41588-018-0205-x). After downloading and unzipping the summary statistics, we load them into a pandas DataFrame:
 
 ```python
 import pandas as pd
@@ -130,6 +120,10 @@ The `genal.Geno` takes as input a pandas dataframe where each row corresponds to
 - **P**: Column name for effect p-value. Defaults to `'P'`.
 - **EAF**: Column name for effect allele frequency. Defaults to `'EAF'`.
 
+> **Note:**
+> 
+> You do not need all columns to move forward, as not all columns are required by every function. Additionally, some columns can be imputed as we will see in the next paragraph.
+
 After inspecting the dataframe, we first need to extract the chromosome and position information from the `MarkerName` column into two new columns `CHR` and `POS`:
 
 ```python
@@ -155,7 +149,7 @@ The last argument (`keep_columns = False`) indicates that we do not wish to keep
 
 > **Note:**
 > 
-> Make sure to read the readme file usually provided with the summary statistics to identify the correct columns. It is particularly important to correctly identify the allele that represents the effect allele. Also, you do not need all columns to move forward, as some can be inputted as we will see next.
+> Make sure to read the readme file usually provided with the summary statistics to identify the correct columns. It is particularly important to correctly identify the allele that represents the effect allele.
 
 ### Data preprocessing <a name="paragraph3.2"></a>
 
@@ -334,7 +328,7 @@ and the output is:
     The PRS computation was successful and used 1330/1538 (86.476%) SNPs.
     PRS data saved to SBP_prs.csv
 
-In our case, we have been able to find proxies for 571 of the 786 SNPs that were missing in the population genetic data (7 potential proxies have been removed because they were identical to SNPs already present in our data).
+In our case, we have been able to find proxies for 578 of the 786 SNPs that were missing in the population genetic data (7 potential proxies have been removed because they were identical to SNPs already present in our data).
 
 You can customize how the proxies are chosen with the following arguments:
 - `reference_panel`: The reference population used to derive linkage disequilibrium values and find proxies. Defaults to `eur`.
@@ -344,7 +338,7 @@ You can customize how the proxies are chosen with the following arguments:
 
 > **Note:**
 > 
-> You can call the `genal.Geno.prs` method on any `Geno` instance (containing at least the EA, BETA, and either SNP or CHR/POS columns). The data does not need to be clumped, and there is no limit to the number of instruments used to compute the scores.
+> You can call the `genal.Geno.prs` method on any `Geno` instance (containing at least the EA, BETA, and either SNP or CHR/POS columns). The data does not need to be clumped, and there is no limit to the number of SNPs used to compute the scores.
 
 
 ### Mendelian Randomization <a name="paragraph3.5"></a>
@@ -459,7 +453,7 @@ If you want to visualize the obtained MR results, you can use the `genal.Geno.MR
 SBP_clumped.MR_plot(filename="MR_plot_SBP_AS")
 ```
 
-![MR plot](docs/Images/MR_plot_SBP_AS.png)
+![MR plot](docs/build/_images/MR_plot_SBP_AS.png)
 You can select which MR methods you wish to plot with the `methods` argument. Note that for an MR method to be plotted, they must be included in the latest `genal.Geno.MR` call of this `genal.Geno` instance.
 
 If you wish to include the heterogeneity values (Cochran's Q) in the results, you can use the heterogeneity argument in the `genal.Geno.MR` call. Here, the heterogeneity for the inverse-variance weighted method:
@@ -479,7 +473,7 @@ As expected, many MR methods indicate that SBP is strongly associated with strok
 To investigate horizontal pleiotropy in more details, a very useful method is Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO). MR-PRESSO is a method designed to detect and correct for horizontal pleiotropy. It will identify which instruments are likely to be pleiotropic on their effect on the outcome, and it will rerun an inverse-variance weighted MR after excluding them. It can be run using the `genal.Geno.MRpresso` method:
 
 ```python
-SBP_clumped.MRpresso(action = 2, n_iterations = 30000)
+mod_table, GlobalTest, OutlierTest, BiasTest = SBP_clumped.MRpresso(action = 2, n_iterations = 30000)
 ```
 
 As with the `genal.Geno.MR` method, the `action` argument determines how the pleiotropic SNPs will be treated. The output is a list containing:
@@ -504,7 +498,7 @@ df_pheno = pd.read_csv("path/to/trait/data")
 
 > **Note:**
 > 
-> One important detail is to make sure that the individual IDs are identical between the phenotypic data and the genetic data for the target population.
+>    One important detail is to make sure that the IDs of the participants are identical in the phenotypic data and in the genetic data.
 
 Then, it is advised to make a copy of the `genal.Geno` instance containing our instruments as we are going to update their coefficients and to avoid any confusion:
 
@@ -582,6 +576,54 @@ You can specify the path of the LiftOver executable to the `liftover_path` argum
 SBP_Geno.lift(start = "hg19", end = "hg38", replace = False, liftover_path = "path/to/liftover/exec")
 ```
 
+### GWAS Catalog <a name="paragraph3.8"></a>
+
+It is sometimes interesting to determine the traits associated with our SNPs. In Mendelian Randomization, for instance, we may want to exclude instruments that are associated with traits likely causing horizontal pleiotropy. For this purpose, we can use the `genal.Geno.query_gwas_catalog` method. This method will query the GWAS Catalog API to determine the list of traits associated with each of our SNPs and store the results in a list in the `ASSOC` column of the `.data` attribute:
+
+```python
+SBP_clumped.query_gwas_catalog(p_threshold=5e-8)
+```
+Which will output:
+
+    Querying the GWAS Catalog and creating the ASSOC column. 
+    Only associations with a p-value <= 5e-08 are reported. Use the p_threshold argument to change the threshold.
+    To report the p-value for each association, use return_p=True.
+    To report the study ID for each association, use return_study=True.
+    The .data attribute will be modified. Use replace=False to leave it as is.
+    100%|██████████| 1545/1545 [00:34<00:00, 44.86it/s]
+    The ASSOC column has been successfully created.
+    701 (45.37%) SNPs failed to query (not found in GWAS Catalog) and 7 (0.5%) SNPs timed out after 34.33 seconds. You can increase the timeout value with the timeout argument.
+| EA  | NEA | EAF   | BETA   | SE     | CHR | POS        | SNP        | ASSOC                                                                 |
+|-----|-----|-------|--------|--------|-----|------------|------------|------------------------------------------------------------------------|
+| A   | G   | 0.1784| 0.2330 | 0.0402 | 10  | 102075479  | rs603424   | [eicosanoids measurement, decadienedioic acid (...]                     |
+| A   | G   | 0.0706| -0.3873| 0.0626 | 10  | 102403682  | rs2996303  | FAILED_QUERY                                                           |
+| T   | G   | 0.8872| 0.6846 | 0.0480 | 10  | 102553647  | rs1006545  | [diastolic blood pressure, systolic blood pressure...]                  |
+| T   | G   | 0.6652| -0.2098| 0.0340 | 10  | 102558506  | rs12570050 | FAILED_QUERY                                                           |
+| T   | C   | 0.3057| -0.2448| 0.0334 | 10  | 102603924  | rs4919502  | FAILED_QUERY                                                           |
+| ... | ... | ...   | ...    | ...    | ... | ...        | ...        | ...                                                                    |                                          |
+| T   | C   | 0.3514| 0.2203 | 0.0314 | 9   | 9350706    | rs1332813  | [diastolic blood pressure, systolic blood pressure...]                  |
+| T   | C   | 0.6880| -0.1897| 0.0332 | 9   | 94201341   | rs10820855 | FAILED_QUERY                                                           |
+| A   | T   | 0.3669| -0.1862| 0.0313 | 9   | 95201540   | rs7045409  | [protein measurement, pulse pressure measurement...]                   |
+
+If you are also interested in the p-values of each SNP-trait association, or the ID of the study from which the association was reported, you can use the `return_p = True` and `return_study = True` arguments. Then, the `ASSOC` column will contain a list of tuples, where each tuple contains the trait name, the p-value, and the study ID:
+
+```python
+SBP_clumped.query_gwas_catalog(p_threshold=5e-8, return_p=True, return_study=True)
+```
 
+| EA  | NEA | EAF   | BETA   | SE     | CHR | POS        | SNP        | ASSOC                                                                 |
+|-----|-----|-------|--------|--------|-----|------------|------------|------------------------------------------------------------------------|
+| A   | G   | 0.1784| 0.2330 | 0.0402 | 10  | 102075479  | rs603424   | TIMEOUT                                                                |
+| A   | G   | 0.0706| -0.3873| 0.0626 | 10  | 102403682  | rs2996303  | FAILED_QUERY                                                           |
+| T   | G   | 0.8872| 0.6846 | 0.0480 | 10  | 102553647  | rs1006545  | [(heart rate response to exercise, 6e-12, GCST...                      |
+| T   | G   | 0.6652| -0.2098| 0.0340 | 10  | 102558506  | rs12570050 | FAILED_QUERY                                                           |
+| T   | C   | 0.3057| -0.2448| 0.0334 | 10  | 102603924  | rs4919502  | FAILED_QUERY                                                           |
+| ... | ... | ...   | ...    | ...    | ... | ...        | ...        | ...                                                                    |                                                         |
+| T   | C   | 0.3514| 0.2203 | 0.0314 | 9   | 9350706    | rs1332813  | [(diastolic blood pressure, 1e-12, GCST9031029...                      |
+| T   | C   | 0.6880| -0.1897| 0.0332 | 9   | 94201341   | rs10820855 | FAILED_QUERY                                                           |
+| A   | T   | 0.3669| -0.1862| 0.0313 | 9   | 95201540   | rs7045409  | [(systolic blood pressure, 9e-13, GCST006624),...                      |
 
 
+> **Note:**
+> 
+> As you can see, many SNPs failed to be queried. This is normal as the GWAS Catalog is not exhaustive.

genal_python-0.8/docs/build/.buildinfo

@@ -0,0 +1,4 @@
+# Sphinx build info version 1
+# This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done.
+config: 1a3c03fa317dbf0f46b6f7567774d6c5
+tags: 645f666f9bcd5a90fca523b33c5a78b7