gwaslab 3.5.8__py3-none-any.whl → 3.6.1__py3-none-any.whl

gwaslab/data/reference.json

@@ -103,7 +103,9 @@
     "13to19":"https://s3-us-west-2.amazonaws.com/human-pangenomics/T2T/CHM13/assemblies/chain/v1_nflo/chm13v2-hg19.chain",
     "13to38":"https://s3-us-west-2.amazonaws.com/human-pangenomics/T2T/CHM13/assemblies/chain/v1_nflo/chm13v2-grch38.chain",
     "18to19":"https://hgdownload.soe.ucsc.edu/goldenPath/hg18/liftOver/hg18ToHg19.over.chain.gz",
-    "18to38":"https://hgdownload.soe.ucsc.edu/goldenPath/hg18/liftOver/hg18ToHg38.over.chain.gz"
+    "18to38":"https://hgdownload.soe.ucsc.edu/goldenPath/hg18/liftOver/hg18ToHg38.over.chain.gz",
+    "1kg_hm3_hg38_eaf":"https://www.dropbox.com/scl/fi/ymkqfsaec6mwjzlvxsm45/PAN.hapmap3.hg38.EAF.tsv.gz?rlkey=p1auef5y1kk7ui41k6j3s8b0z&dl=1",
+    "1kg_hm3_hg19_eaf":"https://www.dropbox.com/scl/fi/dmv9wtfchv6ahim86d49r/PAN.hapmap3.hg19.EAF.tsv.gz?rlkey=ywne2gj1rlm2nj42q9lt2d99n&dl=1"
 }
 
 

gwaslab/g_Sumstats.py

@@ -84,6 +84,8 @@ from gwaslab.util_ex_ldsc import _estimate_rg_by_ldsc
 from gwaslab.util_ex_ldsc import _estimate_h2_cts_by_ldsc
 from gwaslab.util_ex_ldsc import _estimate_partitioned_h2_by_ldsc 
 from gwaslab.util_ex_ldproxyfinder import _extract_ld_proxy
+from gwaslab.util_ex_run_magma import _run_magma
+from gwaslab.util_ex_infer_ancestry import _infer_ancestry
 from gwaslab.bd_get_hapmap3 import gethapmap3
 from gwaslab.util_abf_finemapping import abf_finemapping
 from gwaslab.util_abf_finemapping import make_cs
@@ -674,6 +676,9 @@ class Sumstats():
         fig,outliers = plotdaf(self.data, **kwargs)
         return fig, outliers
     
+    def infer_ancestry(self, **kwargs):
+        self.meta["gwaslab"]["inferred_ancestry"] = _infer_ancestry(self.data, **kwargs)
+
     def plot_gwheatmap(self, **kwargs):
         fig = _gwheatmap(self.data, **kwargs)
         return fig
@@ -882,8 +887,14 @@ class Sumstats():
         if build is None:
             build = self.meta["gwaslab"]["genome_build"]
         insumstats = gethapmap3(self.data.copy(), build=build, verbose=verbose , match_allele=match_allele, how=how )
-        _run_prscs(sst_file = insumstats[["rsID","CHR","POS","EA","NEA","BETA","SE"]], log=self.log, **kwargs)
-
+        _run_prscs(sst_file = insumstats[["rsID","CHR","POS","EA","NEA","BETA","SE"]], 
+                   log=self.log, 
+                   **kwargs)
+    
+    def run_magma(self, build=None, verbose=True, **kwargs):
+        _run_magma(self.data, 
+                   study=self.meta["gwaslab"]["study_name"], 
+                   build=build, verbose=verbose, log=self.log, **kwargs)
 ## LDSC ##############################################################################################
     def estimate_h2_by_ldsc(self, build=None, verbose=True, match_allele=True, how="right", **kwargs):
         if build is None:

gwaslab/g_meta.py

@@ -6,7 +6,7 @@ def _init_meta(object="Sumstats"):
                         "gwas_id":"Unknown", 
                         "samples":{
                             "sample_size":"Unknown", 
-                            "sample_ancestry":"European", 
+                            "sample_ancestry":"Unknown", 
                             "ancestry_method":"self-reported|genetically determined", 
                         } ,
                         "trait_description":"Unknown", 
@@ -37,6 +37,7 @@ def _init_meta(object="Sumstats"):
                     "species":"homo sapiens",
                     "genome_build":"99",
                     "sample_prevalence":"Unknown",
+                    "inferred_ancestry":"Unknown",
                     "population_prevalence":"Unknown",
                     "variants":{
                         "variant_number":"Unknown",

gwaslab/g_version.py

@@ -15,8 +15,8 @@ def _get_version():
 def gwaslab_info():
     # version meta information
     dic={
-   "version":"3.5.8",
-   "release_date":"20250424"
+   "version":"3.6.1",
+   "release_date":"20250506"
     }
     return dic   
 

gwaslab/util_ex_infer_ancestry.py

@@ -0,0 +1,65 @@
+
+import pandas as pd
+from gwaslab.g_Log import Log
+
+def _infer_ancestry(sumstats, 
+                    ancestry_af=None,
+                    build="19",
+                    log=Log(),
+                    verbose=True):
+    log.write("Start to infer ancestry based on Fst...", verbose=verbose)
+    ref_af = pd.read_csv(ancestry_af, sep="\t")
+    
+    data_af = pd.merge(sumstats[["CHR","POS","EA","NEA","EAF"]] ,ref_af,on=["CHR","POS"],how="inner") 
+
+    log.write(f"  -Estimating Fst using {len(data_af)} variants...", verbose=verbose)
+
+    is_filp = data_af["EA"] == data_af["ALT"]
+    data_af.loc[is_filp, ["EA","NEA"]] = data_af.loc[is_filp, ["NEA","EA"]]
+    data_af.loc[is_filp, "EAF"] = 1 - data_af.loc[is_filp, "EAF"]
+    
+    headers = []
+    for i in ['GBR', 'FIN', 'CHS', 'PUR', 'CDX',
+        'CLM', 'IBS', 'PEL', 'PJL', 'KHV', 'ACB', 'GWD', 'ESN', 'BEB', 'MSL',
+        'STU', 'ITU', 'CEU', 'YRI', 'CHB', 'JPT', 'LWK', 'ASW', 'MXL', 'TSI',
+        'GIH', 'EUR', 'EAS', 'AMR', 'SAS', 'AFR']:
+        headers.append(f"FST_{i}")
+        data_af[f"FST_{i}"] = data_af.apply(lambda x: calculate_fst(x["EAF"], x[i]), axis=1)
+    
+    for i,value in data_af[headers].mean().sort_values().items():
+        log.write( f"  -{i} : {value}", verbose=verbose)
+        
+        closest_ancestry = data_af[headers].mean().sort_values().idxmin()
+
+    log.write(f"  -Closest Ancestry: {closest_ancestry.split('_')[1]}", verbose=verbose)
+    log.write("Finished inferring ancestry.", verbose=verbose)
+    return closest_ancestry.split("_")[1]
+
+def calculate_fst(p_1, p_2):
+    # https://bios1140.github.io/understanding-fst-the-fixation-index.html
+    # calculate q1 and q2
+    q_1 = 1 - p_1
+    q_2 = 1 - p_2
+
+    # calculate total allele frequency
+    p_t = (p_1 + p_2)/2
+    q_t = 1 - p_t
+
+    # calculate expected heterozygosity
+    # first calculate expected heterozygosity for the two populations
+    # pop1
+    hs_1 = 2*p_1*q_1
+    # pop2
+    hs_2 = 2*p_2*q_2
+    # then take the mean of this
+    hs = (hs_1 + hs_2)/2
+
+    # next calculate expected heterozygosity for the metapopulations
+    ht = 2*p_t*q_t
+
+    # calculate fst
+    fst = (ht - hs)/ht
+
+    # return output
+    return fst
+    

gwaslab/util_ex_run_magma.py

@@ -0,0 +1,74 @@
+import subprocess
+import os
+import gc
+import pandas as pd
+import numpy as np
+from gwaslab.g_Log import Log
+from gwaslab.util_in_filter_value import _exclude_hla
+
+def _run_magma(sumstats, 
+                magma="magma",
+                study="Study1",
+                exclude_hla=True,
+                window="35,10",
+                id_to_use="rsID",
+                ref=None,
+                ncbi=None,
+                set_annot=None,
+                out="./",
+                delete=True, 
+                ncol="N",
+                build="19",
+                log=Log(), 
+                verbose=True):
+    
+    log.write(" Start to run magma from command line:", verbose=verbose)
+
+    if exclude_hla==True:
+        sumstats = _exclude_hla(sumstats, build =build)
+
+    snploc="{}{}.rsid.chr.pos.tsv".format(out,study)
+    pval="{}{}.rsid.p.n.tsv".format(out, study)
+
+    log.write(f" -writing temp file for --snp-loc:{snploc}", verbose=verbose)
+    sumstats.dropna()[[id_to_use,"CHR","POS"]].rename(columns={id_to_use:"SNP"}).to_csv("{}{}.rsid.chr.pos.tsv".format(out,study),index=None, sep="\t")
+    
+    log.write(f" -writing temp file for --pval:{pval}", verbose=verbose)
+    sumstats.dropna()[[id_to_use,"P","N"]].rename(columns={id_to_use:"SNP"}).to_csv("{}{}.rsid.p.n.tsv".format(out,study),index=None, sep="\t")
+    
+    log.write(f" --annotate window: {window}", verbose=verbose)
+    log.write(f" --gene-loc: {ncbi}", verbose=verbose)
+    log.write(f" --bfile: {ref}", verbose=verbose)
+    log.write(f" Output prefix: {out}", verbose=verbose)
+
+    bash_script=f'''
+
+#!/bin/bash
+
+{magma} --annotate window={window}  --snp-loc {snploc} --gene-loc {ncbi} --out {study}
+
+{magma} --bfile {ref} --pval {pval} ncol={ncol} --gene-annot {study}.genes.annot --out {study}
+
+'''
+    
+    if set_annot is not None:
+        bash_script+=f'''
+{magma} --gene-results {study}.genes.raw --set-annot {set_annot} --out {study}
+'''
+    log.write(f"Script: {bash_script}")
+
+    try:
+        log.write(" Running magma from command line...", verbose=verbose)
+        output = subprocess.check_output(bash_script, stderr=subprocess.STDOUT, shell=True,text=True)
+        output =  output + "\n"
+        
+        if delete == True:
+            os.remove(snploc)
+            os.remove(pval)
+
+    except subprocess.CalledProcessError as e:
+        log.warning("ERROR!")
+        log.write(e.output)
+
+    log.write("Finished running magma.", verbose=verbose)
+

gwaslab/util_in_meta.py

@@ -7,7 +7,6 @@ from gwaslab.g_Log import Log
 from gwaslab.io_to_pickle import load_data_from_pickle
 from gwaslab.g_Sumstats import Sumstats
 import gc
-import statsmodels.api as sm
 
 def meta_analyze(sumstats_list,
                  random_effects=False, 

gwaslab/util_in_meta_polars.py

@@ -6,9 +6,7 @@ from scipy.stats import norm
 from gwaslab.g_Log import Log
 from gwaslab.io_to_pickle import load_data_from_pickle
 from gwaslab.g_Sumstats import Sumstats
-import gc
 import polars as pl
-import statsmodels.api as sm
 ########################################################################################################################################################################################################################################################################################################################################################
 ########################################################################################################################################################################################################################################################################################################################################################
 ########################################################################################################################################################################################################################################################################################################################################################

gwaslab/viz_aux_annotate_plot.py

@@ -26,6 +26,7 @@ def annotate_single(
     anno_alias,
     anno_style,
     anno_args,
+    anno_args_single,
     arm_scale,
     anno_max_iter,
     arm_scale_d,
@@ -216,13 +217,23 @@ def annotate_single(
             if anno_style == "tight" :
                 anno_default["rotation"] = 90
             ################################################################################################################################
-
+            # anno args for all
             for key,value in anno_args.items():
                 anno_default[key]=value
-            if len(highlight_i) >0 and highlight_chrpos==True:
+            
+            # anno args for highlight group
+            if len(highlight_i) >0:
                 if row["i"] in highlight_i:
                     for key,value in highlight_anno_args.items():
                         anno_default[key]=value
+            
+            # anno args for specifc
+            #try:
+            if row[snpid] in anno_args_single.keys():
+                for key,value in anno_args_single[row[snpid]].items():
+                    anno_default[key]=value
+            #except:
+            #    pass
             ################################################################################################################################
             if anno_adjust==True:
                 if  _invert==False:

gwaslab/viz_plot_mqqplot.py

@@ -154,6 +154,7 @@ def mqqplot(insumstats,
           anno_alias=None,
           anno_d=None,
           anno_args=None,
+          anno_args_single=None,
           anno_style="right",
           anno_fixed_arm_length=None,
           anno_source = "ensembl",
@@ -209,7 +210,9 @@ def mqqplot(insumstats,
           drop_chr_start=False,
           title =None,
           mtitle=None,
+          mtitle_pad=1.08,
           qtitle=None,
+          qtitle_pad=1.08,
           ylabel=None,
           xlabel=None,
           title_pad=1.08, 
@@ -256,6 +259,7 @@ def mqqplot(insumstats,
     anno_alias = _update_args(anno_alias, dict())
     anno_d = _update_args(anno_d,dict())
     anno_args = _update_args(anno_args,dict())
+    anno_args_single = _update_args(anno_args_single,dict())
     arrow_kwargs = _update_args(arrow_kwargs,dict())
 
     colors = _update_arg(colors, ["#597FBD","#74BAD3"])
@@ -1002,7 +1006,7 @@ def mqqplot(insumstats,
 
         
         if mtitle and anno and len(to_annotate)>0: 
-            pad=(ax1.transData.transform((skip, title_pad*maxy))[1]-ax1.transData.transform((skip, maxy)))[1]
+            pad=(ax1.transData.transform((skip, mtitle_pad*maxy))[1]-ax1.transData.transform((skip, maxy)))[1]
             ax1.set_title(mtitle,pad=pad,fontsize=title_fontsize,family=font_family)
         elif mtitle:
             ax1.set_title(mtitle,fontsize=title_fontsize,family=font_family)
@@ -1023,6 +1027,7 @@ def mqqplot(insumstats,
                                 anno_alias=anno_alias,
                                 anno_style=anno_style,
                                 anno_args=anno_args,
+                                anno_args_single=anno_args_single,
                                 arm_scale=arm_scale,
                                 anno_max_iter=anno_max_iter,
                                 arm_scale_d=arm_scale_d,
@@ -1074,6 +1079,7 @@ def mqqplot(insumstats,
                     fontsize=fontsize,
                     font_family=font_family,
                     qtitle=qtitle,
+                    qtitle_pad=qtitle_pad,
                     title_fontsize=title_fontsize,
                     include_chrXYMT=include_chrXYMT,
                     cut_line_color=cut_line_color,
@@ -1106,9 +1112,10 @@ def mqqplot(insumstats,
     # Titles 
     if title and anno and len(to_annotate)>0:
         # increase height if annotation 
-        fig.suptitle(title , fontsize = title_fontsize ,x=0.5, y=1.05)
+        fig.suptitle(title , fontsize = title_fontsize ,x=0.5, y=title_pad)
     else:
-        fig.suptitle(title , fontsize = title_fontsize, x=0.5,y=1)
+        title_pad = title_pad -0.05
+        fig.suptitle(title , fontsize = title_fontsize, x=0.5,y=title_pad)
         ## Add annotation arrows and texts
     
     # Saving figure

gwaslab/viz_plot_qqplot.py

@@ -30,6 +30,7 @@ def _plot_qq(
     fontsize,
     font_family,
     qtitle,
+    qtitle_pad,
     title_fontsize,
     include_chrXYMT,
     cut_line_color,
@@ -182,7 +183,8 @@ def _plot_qq(
     ax2.tick_params(axis='both', which='both', labelsize=fontsize,labelfontfamily=font_family)
     #
     if qtitle:
-        ax2.set_title(qtitle,fontsize=title_fontsize,pad=10,family=font_family)
+        pad=(ax2.transData.transform((skip, qtitle_pad*maxy))[1]-ax2.transData.transform((skip, maxy)))[1]
+        ax2.set_title(qtitle,fontsize=title_fontsize,pad=pad,family=font_family)
 
     log.write("Finished creating QQ plot successfully!",verbose=verbose)
     

{gwaslab-3.5.8.dist-info → gwaslab-3.6.1.dist-info}/METADATA

@@ -1,12 +1,12 @@
 Metadata-Version: 2.4
 Name: gwaslab
-Version: 3.5.8
+Version: 3.6.1
 Summary: A collection of handy tools for GWAS SumStats
 Author-email: Yunye <yunye@gwaslab.com>
 Project-URL: Homepage, https://cloufield.github.io/gwaslab/
 Project-URL: Github, https://github.com/Cloufield/gwaslab
 Classifier: Programming Language :: Python :: 3
-Classifier: License :: OSI Approved :: MIT License
+Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
 Classifier: Operating System :: OS Independent
 Requires-Python: <3.13,>=3.9
 Description-Content-Type: text/markdown
@@ -25,6 +25,7 @@ Requires-Dist: scikit-allel>=1.3.5
 Requires-Dist: pyensembl==2.2.3
 Requires-Dist: gtfparse==1.3.0
 Requires-Dist: h5py>=3.10.0
+Requires-Dist: pyarrow
 Dynamic: license-file
 
 # GWASLab
@@ -36,43 +37,63 @@ Dynamic: license-file
 ![badge_pip](https://img.shields.io/pypi/dm/gwaslab)
 ![badge_commit_m](https://img.shields.io/github/commit-activity/m/Cloufield/gwaslab)
 
-* A handy Python toolkit for handling GWAS summary statistics (sumstats).
+* A handy Python-based toolkit for handling GWAS summary statistics (sumstats).
 * Each process is modularized and can be customized to your needs.
 * Sumstats-specific manipulations are designed as methods of a Python object, `gwaslab.Sumstats`.
 
-Please check GWASLab documentation at [https://cloufield.github.io/gwaslab/](https://cloufield.github.io/gwaslab/)
+## Installation
 
-Note: GWASLab is being updated very frequently for now. We will release the first stable version soon! Please stay tuned.
+### install via pip
 
-Warning: Known issues of GWASLab are summarized in [https://cloufield.github.io/gwaslab/KnownIssues/](https://cloufield.github.io/gwaslab/KnownIssues/) .
+The latest version of GWASLab now supports Python 3.9, 3.10, 3.11, and 3.12.
 
-## Install
+```bash
+pip install gwaslab
+```
 
-### install via pip
+### install in conda environment
 
-The latest version of GWASLab now supports Python 3.9, 3.10, 3.11, and 3.12.
+Create a Python 3.9, 3.10, 3.11 or 3.12 environment and install gwaslab using pip:
+
+```bash
+conda env create -n gwaslab -c conda-forge python=3.12
+
+conda activate gwaslab
 
+pip install gwaslab
 ```
-pip install gwaslab==3.5.7
+
+or create a new environment using yml file [environment.yml](https://github.com/Cloufield/gwaslab/blob/main/environment.yml)
+
+```bash
+conda env create -n gwaslab -f environment.yml
 ```
 
+### install using docker (deprecated)
+
+A docker file is available [here](https://github.com/Cloufield/gwaslab/blob/main/docker/Dockerfile) for building local images.
+
+## Quick start
+
 ```python
+
 import gwaslab as gl
+
 # load plink2 output
-mysumstats = gl.Sumstats("t2d_bbj.txt.gz", fmt="plink2")
+mysumstats = gl.Sumstats("sumstats.txt.gz", fmt="plink2")
 
-# load sumstats with auto mode (auto-detecting common headers) 
+# or load sumstats with auto mode (auto-detecting commonly used headers) 
 # assuming ALT/A1 is EA, and frq is EAF
-mysumstats = gl.Sumstats("t2d_bbj.txt.gz", fmt="auto")
+mysumstats = gl.Sumstats("sumstats.txt.gz", fmt="auto")
 
 # or you can specify the columns:
-mysumstats = gl.Sumstats("t2d_bbj.txt.gz",
+mysumstats = gl.Sumstats("sumstats.txt.gz",
              snpid="SNP",
              chrom="CHR",
              pos="POS",
              ea="ALT",
              nea="REF",
-             neaf="Frq",
+             eaf="Frq",
              beta="BETA",
              se="SE",
              p="P",
@@ -85,26 +106,9 @@ mysumstats.plot_mqq()
 ...
 ```
 
-### install in conda environment
-
-Create a Python 3.9 environment and install gwaslab using pip:
-
-```
-conda env create -n gwaslab_test -c conda-forge python=3.9
-conda activate gwaslab
-pip install gwaslab==3.4.45
-```
-
-or create a new environment using yml file [environment_3.4.40.yml](https://github.com/Cloufield/gwaslab/blob/main/environment_3.4.40.yml)
-
-```
-conda env create -n gwaslab -f environment_3.4.40.yml
-```
-
-
-### install using docker
+## Documentation and tutorials
 
-A docker file is available [here](https://github.com/Cloufield/gwaslab/blob/main/docker/Dockerfile) for building local images.
+Documentation and tutorials for GWASLab are avaiable at [here](https://cloufield.github.io/gwaslab/).
 
 ## Functions
 
@@ -151,7 +155,7 @@ A docker file is available [here](https://github.com/Cloufield/gwaslab/blob/main
 - Scatter plot: allele frequency comparison 
 - Scatter plot: trumpet plot (plot of MAF and effect size with power lines)
 
-### Visualization Examples
+#### Visualization Examples
 
 <img width="600" alt="image" src="https://user-images.githubusercontent.com/40289485/233836639-34b03c47-5a59-4fd4-9677-5e13b02aab15.png">
 <img width="600" alt="image" src="https://user-images.githubusercontent.com/40289485/197393168-e3e7076f-2801-4d66-9526-80778d44f3da.png">
@@ -167,42 +171,15 @@ A docker file is available [here](https://github.com/Cloufield/gwaslab/blob/main
 - Sumstats summary: give you a quick overview of the sumstats. 
 - ...
 
-## Requirements (deprecated)
-
-environment.yml
+## Issues
 
-```
-name: gwaslab
-channels:
-  - conda-forge
-  - defaults
-dependencies:
-  - python=3.8.16=h7a1cb2a_3
-  - jupyter==1.0.0
-  - pip==23.1.2
-  - pip:
-      - adjusttext==0.8
-      - biopython==1.81
-      - gwaslab==3.4.16
-      - liftover==1.1.16
-      - matplotlib==3.7.1
-      - numpy==1.24.2
-      - pandas==1.4.4
-      - scikit-allel==1.3.5
-      - scikit-learn==1.2.2
-      - scipy==1.10.1
-      - seaborn==0.11.2
-      - statsmodels==0.13
-      - adjustText==0.8
-      - pysam==0.19
-      - pyensembl==2.2.3
-      - h5py==3.10.0
-```
+- GWASLab is currently under active development, with frequent updates.
+- Note: Known issues are documented at https://cloufield.github.io/gwaslab/KnownIssues/.
 
 ## How to cite
 - GWASLab preprint: He, Y., Koido, M., Shimmori, Y., Kamatani, Y. (2023). GWASLab: a Python package for processing and visualizing GWAS summary statistics. Preprint at Jxiv, 2023-5. https://doi.org/10.51094/jxiv.370
 
-## Sample Data
+## Sample data used for tutorial
 - Sample GWAS data used in GWASLab is obtained from: http://jenger.riken.jp/ (Suzuki, Ken, et al. "Identification of 28 new susceptibility loci for type 2 diabetes in the Japanese population." Nature genetics 51.3 (2019): 379-386.).
 
 ## Acknowledgement

{gwaslab-3.5.8.dist-info → gwaslab-3.6.1.dist-info}/RECORD

@@ -6,7 +6,7 @@ gwaslab/bd_get_hapmap3.py,sha256=qWTvIRZsd7F3nT9sN2NSXUsxZJRf5k4HLgJ6kN0qaUc,410
 gwaslab/cache_manager.py,sha256=HOTnSkCOyGEPLRl90WT8D_6pAdI8d8AzenMIDGuCeWc,28113
 gwaslab/g_Log.py,sha256=C3Zv-_6c3C9ms8bgQ-ytplz22sjk7euqXYkWr9zNeAs,1573
 gwaslab/g_Phenotypes.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-gwaslab/g_Sumstats.py,sha256=J7c8nbqgGu0Ih3rRZkFA_GvVFbsWlcSaDWzyvJGWqX8,42993
+gwaslab/g_Sumstats.py,sha256=l1zmEOifAZvnv7BoLzyLcTV6Yp1KX5n5HH930GQ9xlE,43502
 gwaslab/g_SumstatsMulti.py,sha256=2y8brnKAs5I9tUexQJiFwOsFLC4MyACn7IybvvXsRDU,13828
 gwaslab/g_SumstatsPair.py,sha256=PmVPGU32degDuU5RDG9igyT1sFAbEDh6alrRulUphyk,13181
 gwaslab/g_SumstatsSet.py,sha256=AiTISWPfmu8NTGa8j9Yuts8DNw1pEUENYyPoS0HXp5I,29866
@@ -14,11 +14,11 @@ gwaslab/g_SumstatsT.py,sha256=u_DighLMnMxwTLnqm-B58pA0G6WXRj6pudPyKMVKjSU,2133
 gwaslab/g_Sumstats_polars.py,sha256=F9d8y-6SJLRibrZrvWtwE2XpxZdG3A9DqqTT8TW4sS0,8528
 gwaslab/g_Sumstats_summary.py,sha256=FECvvFXJVKaCX5dggBvvk9YvJ6AbdbcLfjltysX7wEE,6380
 gwaslab/g_headers.py,sha256=B0bTVqzf2pdUwAX1U8MoYDO7nQEKAtgVnn-9bVXsJGE,6766
-gwaslab/g_meta.py,sha256=u0rLqPa3l0_C7b-HSOvuGb6EqrFQJwtI4ANu34xKaAo,6422
+gwaslab/g_meta.py,sha256=pXtSsQfFPv2UJmOxKMmLtcsp_Ku_H73YP7PnlfMd0sg,6472
 gwaslab/g_meta_update.py,sha256=dWgz4kcq9bvsXycCjfhoYtcJHlEP4HRdYoQ8KFqMyTQ,2490
 gwaslab/g_vchange_status.py,sha256=w3zsYYOcCaI3PTeboonvkQjudzUAfVIgATzRdiPViZs,1939
 gwaslab/g_vchange_status_polars.py,sha256=kxyGQCur0ibVFBCyZghA-XNf_kLDXKK-l7VC-Om2IdA,1839
-gwaslab/g_version.py,sha256=-t2SPXU9268hO_iB_dQHkssdawqapRZM9udTRo04Mrs,1889
+gwaslab/g_version.py,sha256=U5RKVaP711pSTX5OzNta51b7Np64kCdo7Gg09AFfD-o,1889
 gwaslab/hm_casting.py,sha256=xoq1E4Tp5VC4aLWfq9-_AfiQzb1WZAHrnZG33W4sCOE,14178
 gwaslab/hm_casting_polars.py,sha256=_3ZeMg3mRsa7vHpxObFErOledpaydhC_PxStcwlYO2c,8371
 gwaslab/hm_harmonize_sumstats.py,sha256=2MeWnWmMHpLWjUMTVqExpg_3mE1VIaBsh6Mz3ffJCMc,84761
@@ -51,6 +51,7 @@ gwaslab/util_abf_finemapping.py,sha256=LRcopjtkT-iXtKPAJIzR4qjPdhD7nrS_BGit4EW89
 gwaslab/util_ex_calculate_ldmatrix.py,sha256=TeXazNXInGYrBR58Q3-f51yrHz5-_A3tXNieYJA2-i8,16924
 gwaslab/util_ex_calculate_prs.py,sha256=9uJ588Sdj4V0vw3OZ9NeLECwOvW67f0IdLandVPS5RY,9442
 gwaslab/util_ex_gwascatalog.py,sha256=RZ_qIsfjONUbF3AGMp-0feTBubOEEZUufYGdvv0Sk7g,9211
+gwaslab/util_ex_infer_ancestry.py,sha256=AXLEhuJUBT6XIYjADB-envFRGIZ-misnegGdtK8SsFc,2192
 gwaslab/util_ex_ldproxyfinder.py,sha256=HadmnjEv5EvZCHR5SFSAbqxZWJryEErJftCS0WhZUbs,16898
 gwaslab/util_ex_ldsc.py,sha256=uEVrsGrT5A-ni7Vw7esvqQ7NwiR1P10Q_N1_uZNRzRY,21600
 gwaslab/util_ex_match_ldmatrix.py,sha256=8nFziC8LiUy83bG9vW_u_CzbZek-9R9eVIfFdRBDkrA,17495
@@ -62,6 +63,7 @@ gwaslab/util_ex_run_ccgwas.py,sha256=8uOAaJ44zAhlxj_A2vth5AzRLk2_BYlYBXEccFbB6T0
 gwaslab/util_ex_run_clumping.py,sha256=yvUhKi83KhXfE5yPy2i47B58BMACZ_r5gt3-uN8Znbo,7908
 gwaslab/util_ex_run_coloc.py,sha256=evLJ_ct_RV7UKLnYGJ_VfT6MeuGeYDtj8M8xqqk-X-U,6297
 gwaslab/util_ex_run_hyprcoloc.py,sha256=sElGLXsBKx3PpDTOET7ulvjBShtDOJ3OyDuOp3Yl_y0,4722
+gwaslab/util_ex_run_magma.py,sha256=0vUac7iVySnRx8QoYBncSfIEke5MwTfmYzDLcSmw1Us,2439
 gwaslab/util_ex_run_mesusie.py,sha256=g0x7-zwoN9fbHYJQAue3vAMIeFu3x8H9KjgiYI1u9Mk,5532
 gwaslab/util_ex_run_mtag.py,sha256=TxfDwjbGy5C1SiTVDMNopcSjl9p5mgCjJVBfeLq__bg,2952
 gwaslab/util_ex_run_prscs.py,sha256=_BOUqX3fgALdGx-xDyWlREXCafZZddhcK8ufoESUaac,3396
@@ -76,10 +78,10 @@ gwaslab/util_in_filter_value.py,sha256=dr1ZiTsMoYZKznQtucHBFW2mVYOv9OTESb24vXNcS
 gwaslab/util_in_get_density.py,sha256=kpKXH69acMkeYVG5vs-VbJC3COhmuLBfYco-wuOxgjc,3934
 gwaslab/util_in_get_sig.py,sha256=jOMdDu1iWq5-qs8crHLNvn4t12VHSGQp0kG-ce2l3pQ,42441
 gwaslab/util_in_merge.py,sha256=KB5VKRTUUZ1XGYUxgCP_l3cWYanjthMdwiZ_DincCZQ,2020
-gwaslab/util_in_meta.py,sha256=3OxpnigfaShV7995-1HCbGCwdIjSbsDhC2X7hnT28jo,18712
-gwaslab/util_in_meta_polars.py,sha256=i_2DGSDX8KlR3gtuDPXw88a_rbGYWIb58J1MGDMbJx8,8426
+gwaslab/util_in_meta.py,sha256=_oEWvOvj6DG28sljstS0SYl8sm7oP-5pspGRvqbnipA,18683
+gwaslab/util_in_meta_polars.py,sha256=WgAelsxtvheHA-Jq1sbxh02nDhyTFHxBpP3G-iSRREY,8387
 gwaslab/util_in_snphwe.py,sha256=-KpIDx6vn_nah6H55IkV2OyjXQVXV13XyBL069WE1wM,1751
-gwaslab/viz_aux_annotate_plot.py,sha256=3PyBioC-3ZBzqKT2JRv6vp-hZUXSqTTePIT5cblEvnQ,25334
+gwaslab/viz_aux_annotate_plot.py,sha256=IyGhqFQofrVW9A8UxJ5oNVH96czQ-2_i4s1lVvf2hww,25696
 gwaslab/viz_aux_chromatin.py,sha256=aWZaXOSvGyZY7wQcoFDaqHRYCSHZbi_K4Q70HruN9ts,4125
 gwaslab/viz_aux_property.py,sha256=UIaivghnLXYpTwkKnXRK0F28Jbn9L6OaICk3K73WZaU,33
 gwaslab/viz_aux_quickfix.py,sha256=cGX5i3WBmvKIiqck8V00caDg-pvKOO709Ux3DBXsUrM,18693
@@ -92,9 +94,9 @@ gwaslab/viz_plot_effect.py,sha256=qbM6c1IB2HlUlMNgFZlJ5G8ODQJ8-oSWD8t0Q8DDuz8,10
 gwaslab/viz_plot_forestplot.py,sha256=xgOnefh737CgdQxu5naVyRNBX1NQXPFKzf51fbh6afs,6771
 gwaslab/viz_plot_miamiplot.py,sha256=rCFEp7VNuVqeBBG3WRkmFAtFklbF79BvIQQYiSY70VY,31238
 gwaslab/viz_plot_miamiplot2.py,sha256=eVS2v1YFEO0xZ-2zs_DCCvEx-Hqbt0wQkbgCHiRfqeE,16247
-gwaslab/viz_plot_mqqplot.py,sha256=1pgbQuyuGCbwBiecv9eagpAC2GS3Ekr-nGHUTp_Wx20,70795
+gwaslab/viz_plot_mqqplot.py,sha256=mlcsA_wEnv7N9mP6CVL15WXWHJ0-8TwqHuznJXgPwCE,71103
 gwaslab/viz_plot_phe_heatmap.py,sha256=qoXVeFTIm-n8IinNbDdPFVBSz2yGCGK6QzTstXv6aj4,9532
-gwaslab/viz_plot_qqplot.py,sha256=2z2dtNRZJtuJYdAsOdusYr_frOrl-XWfzc0vSYIQblY,7359
+gwaslab/viz_plot_qqplot.py,sha256=cB4vRlFv69zWY9NMLfSkfAbirYp3_EEW2kQiBTEMDoc,7483
 gwaslab/viz_plot_regional2.py,sha256=BoL1V56ww9B2_vFkABgln_f6OrzobiFjUISI5X6XXMM,43146
 gwaslab/viz_plot_regionalplot.py,sha256=8u-5-yfy-UaXhaxVVz3Y5k2kBAoqzczUw1hyyD450iI,37983
 gwaslab/viz_plot_rg_heatmap.py,sha256=z-G4gxK5-H_e13jV8RQnNzXPrKSQ0c7q41-KpMrA-cs,13861
@@ -102,16 +104,16 @@ gwaslab/viz_plot_scatter_with_reg.py,sha256=PmUZDQl2q4Dme3HLPXEwf_TrMjwJADA-uFXN
 gwaslab/viz_plot_stackedregional.py,sha256=Jzu5NvX45L4iTLfzXjNkxeaA9SGggtGGxEJdGUzSsuU,19163
 gwaslab/viz_plot_trumpetplot.py,sha256=y4sAFjzMaSLuWrdr9_ao-wPYCK5DlP2ykiqulWsoN_k,42680
 gwaslab/data/formatbook.json,sha256=j6XRTlJvqex4HLBtGLhiaU8JV4_cD1LWl9uqB3Oou4s,40718
-gwaslab/data/reference.json,sha256=IrjwFnXjrpVUp3zYfcYClpibJE9Y-94gtrC1Aw8sXxg,12332
+gwaslab/data/reference.json,sha256=oaE_HLFeYuEpcZ_rCGY3GwiAJWsDxRtjJOAc188WG28,12620
 gwaslab/data/chrx_par/chrx_par_hg19.bed.gz,sha256=LocZg_ozhZjQiIpgWCO4EYCW9xgkEKpRy1m-YdIpzQs,83
 gwaslab/data/chrx_par/chrx_par_hg38.bed.gz,sha256=VFW11MnQVC-Iu-ZGvUDcEhVpb-HVRsVTg-W-GNJyxP4,82
 gwaslab/data/hapmap3_SNPs/hapmap3_db150_hg19.snplist.gz,sha256=qD9RsC5S2h6l-OdpW09DWXHIi2kcPebctMnhxt8mzfU,10282886
 gwaslab/data/hapmap3_SNPs/hapmap3_db151_hg38.snplist.gz,sha256=Y8ZT2FIAhbhlgCJdE9qQVAiwnV_fcsPt72usBa7RSBM,10225828
 gwaslab/data/high_ld/high_ld_hla_hg19.bed.gz,sha256=R7IkssKu0L4WwkU9SrS84xCMdrkkKL0gnTNO_OKbG0Y,219
 gwaslab/data/high_ld/high_ld_hla_hg38.bed.gz,sha256=76CIU0pibDJ72Y6UY-TbIKE9gEPwTELAaIbCXyjm80Q,470
-gwaslab-3.5.8.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-gwaslab-3.5.8.dist-info/licenses/LICENSE_before_v3.4.39,sha256=GhLOU_1UDEKeOacYhsRN_m9u-eIuVTazSndZPeNcTZA,1066
-gwaslab-3.5.8.dist-info/METADATA,sha256=_sOvj4V9krrRC2-3GUqRbzr4fGlkOlyvVSfAbd254H8,7595
-gwaslab-3.5.8.dist-info/WHEEL,sha256=SmOxYU7pzNKBqASvQJ7DjX3XGUF92lrGhMb3R6_iiqI,91
-gwaslab-3.5.8.dist-info/top_level.txt,sha256=PyY6hWtrALpv2MAN3kjkIAzJNmmBTH5a2risz9KwH08,8
-gwaslab-3.5.8.dist-info/RECORD,,
+gwaslab-3.6.1.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+gwaslab-3.6.1.dist-info/licenses/LICENSE_before_v3.4.39,sha256=GhLOU_1UDEKeOacYhsRN_m9u-eIuVTazSndZPeNcTZA,1066
+gwaslab-3.6.1.dist-info/METADATA,sha256=HE25IdF4FRLHwtbRzZoaERNZHJxhGLR1wUzTu2DSNwU,7043
+gwaslab-3.6.1.dist-info/WHEEL,sha256=0CuiUZ_p9E4cD6NyLD6UG80LBXYyiSYZOKDm5lp32xk,91
+gwaslab-3.6.1.dist-info/top_level.txt,sha256=PyY6hWtrALpv2MAN3kjkIAzJNmmBTH5a2risz9KwH08,8
+gwaslab-3.6.1.dist-info/RECORD,,

{gwaslab-3.5.8.dist-info → gwaslab-3.6.1.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (79.0.1)
+Generator: setuptools (80.3.1)
 Root-Is-Purelib: true
 Tag: py3-none-any