gwaslab 3.4.45__py3-none-any.whl → 3.4.47__py3-none-any.whl

gwaslab/viz_plot_regionalplot.py

@@ -365,6 +365,10 @@ def  _plot_recombination_rate(sumstats,pos, region, ax1, rr_path, rr_chr_dict, r
 
     rc = rc.loc[(rc["Position(bp)"]<region[2]) & (rc["Position(bp)"]>region[1]),:]
     ax4.plot(rc_track_offset+rc["Position(bp)"],rc["Rate(cM/Mb)"],color="#5858FF",zorder=1)
+    
+    ax1.set_zorder(ax4.get_zorder()+1)
+    ax1.patch.set_visible(False)
+    
     if rr_ylabel:
         ax4.set_ylabel("Recombination rate(cM/Mb)")
     if rr_lim!="max":

gwaslab/viz_plot_stackedregional.py

@@ -22,6 +22,7 @@ from gwaslab.bd_common_data import get_number_to_chr
 from gwaslab.bd_common_data import get_recombination_rate
 from gwaslab.bd_common_data import get_gtf
 from gwaslab.viz_aux_reposition_text import adjust_text_position
+from gwaslab.viz_aux_chromatin import _plot_chromatin_state
 from gwaslab.viz_aux_quickfix import _quick_fix
 from gwaslab.viz_aux_quickfix import _get_largenumber
 from gwaslab.viz_aux_quickfix import _quick_add_tchrpos
@@ -37,21 +38,28 @@ from gwaslab.io_to_pickle import load_data_from_pickle
 from gwaslab.g_Sumstats import Sumstats
 from gwaslab.viz_aux_save_figure import save_figure
 from gwaslab.viz_plot_mqqplot import mqqplot
+import matplotlib.patches as patches
 
 def plot_stacked_mqq(objects,
                         vcfs=None,
                         mode="r",
                         mqqratio=3,
                         region=None,
+                        region_chromatin_height=0.1,
+                        region_chromatin_files = None,
+                        region_chromatin_labels= None,
                         titles= None,
                         title_pos=None,
                         title_args=None,
+                        #title_box = None,
                         gtf=None,
                         gene_track_height=0.5,
                         fig_args=None,
                         region_hspace=0.05,
                         subplot_height=4,
                         region_lead_grid_line=None,
+                        fontsize=9,
+                        font_family="Arial",
                         build="99", 
                         save=None,
                         save_args=None,
@@ -72,11 +80,16 @@ def plot_stacked_mqq(objects,
         fig_args = {"dpi":200}
     if region_lead_grid_line is None:
         region_lead_grid_line = {"alpha":0.5,"linewidth" : 2,"linestyle":"--","color":"#FF0000"}
-    if title_pos is None:
-        title_pos = [0.01,0.97]
+    if region_chromatin_files is None:
+        region_chromatin_files = []
+        region_chromatin_height = len(region_chromatin_files) * region_chromatin_height
+    if region_chromatin_labels is None:
+        region_chromatin_labels = []
     if title_args is None:
-        title_args = {}
-    
+        title_args = {"family":"Arial"}
+    else:
+        if "family" not in title_args.keys():
+            title_args["family"] = "Arial"
     # create figure and axes ##################################################################################################################
     if mode=="r":
         if len(vcfs)==1:
@@ -84,9 +97,15 @@ def plot_stacked_mqq(objects,
         n_plot = len(sumstats_list)
         n_plot_plus_gene_track = n_plot + 1
 
+        if len(region_chromatin_files)>0 and mode=="r":
+            height_ratios = [1 for i in range(n_plot_plus_gene_track-1)]+[region_chromatin_height]+[gene_track_height]
+            n_plot_plus_gene_track +=1
+        else:
+            height_ratios = [1 for i in range(n_plot_plus_gene_track-1)]+[gene_track_height]
+        
         fig_args["figsize"] = [16,subplot_height*n_plot_plus_gene_track]
         fig, axes = plt.subplots(n_plot_plus_gene_track, 1, sharex=True, 
-                             gridspec_kw={'height_ratios': [1 for i in range(n_plot_plus_gene_track-1)]+[gene_track_height]},
+                             gridspec_kw={'height_ratios': height_ratios},
                              **fig_args)
         plt.subplots_adjust(hspace=region_hspace)
     elif mode=="m":
@@ -119,6 +138,7 @@ def plot_stacked_mqq(objects,
     ##########################################################################################################################################
     # a dict to store lead variants of each plot
     lead_variants_is={}
+    lead_variants_is_color ={}
 
     ##########################################################################################################################################
     # plot manhattan plot
@@ -132,7 +152,7 @@ def plot_stacked_mqq(objects,
         #################################################################
         if index==0:
             # plot last m and gene track 
-            fig,log,lead_i,lead_i2 = mqqplot(sumstats,
+            fig,log,lead_snp_is,lead_snp_is_color = mqqplot(sumstats,
                             chrom="CHR",
                             pos="POS",
                             p="P",
@@ -141,6 +161,8 @@ def plot_stacked_mqq(objects,
                             snpid="SNPID",
                             vcf_path=vcfs[index],
                             mode=mode,
+                            fontsize=fontsize,
+                            font_family=font_family,
                             region_lead_grid=False,
                             gtf_path="default",
                             rr_ylabel=False,
@@ -153,10 +175,11 @@ def plot_stacked_mqq(objects,
                             log=log,
                             **mqq_args_for_each_plot[index]
                             )  
-            lead_variants_is[index] = (lead_i,lead_i2)
+            lead_variants_is[index] = lead_snp_is
+            lead_variants_is_color[index] = lead_snp_is_color
         else:
             # plot only the scatter plot
-            fig,log,lead_i,lead_i2 = mqqplot(sumstats,
+            fig,log,lead_snp_is,lead_snp_is_color = mqqplot(sumstats,
                             chrom="CHR",
                             pos="POS",
                             p="P",
@@ -165,6 +188,8 @@ def plot_stacked_mqq(objects,
                             snpid="SNPID",
                             vcf_path=vcfs[index],
                             region_lead_grid=False,
+                            fontsize=fontsize,
+                            font_family=font_family,
                             mode=mode,
                             rr_ylabel=False,
                             region_ld_legend=False,
@@ -178,24 +203,66 @@ def plot_stacked_mqq(objects,
                             log=log,
                             **mqq_args_for_each_plot[index]
                             )
-            lead_variants_is[index] = (lead_i,lead_i2)
-    
+            lead_variants_is[index] = lead_snp_is
+            lead_variants_is_color[index] = lead_snp_is_color
+    if len(region_chromatin_files)>0 and mode=="r":
+        xlim_i = axes[-1].get_xlim()
+        fig = _plot_chromatin_state(     region_chromatin_files = region_chromatin_files, 
+                                         region_chromatin_labels = region_chromatin_labels,
+                                         region = region, 
+                                         fig = fig, 
+                                         ax = axes[-2],
+                                         xlim_i=xlim_i,
+                                         log=log,
+                                         verbose=verbose,
+                                         fontsize = fontsize,
+                                         font_family = font_family)    
     # adjust labels
     # drop labels for each plot 
     # set a common laebl for all plots
 
+    #if title_box is None:
+    #    title_box = dict(boxstyle='square', facecolor='white', alpha=1.0, edgecolor="black")
+    #    title_box = {}
+
+    #if title_args is None:
+    #    title_args = {}   
+    #if titles is not None and mode=="r":    
+    #    if title_pos is None:
+    #        title_pos = [0.01,0.99]
+    #    for index,title in enumerate(titles):
+    #        
+    #        current_text = axes[index].text(title_pos[0], title_pos[1] , title, transform=axes[index].transAxes,ha="left", va='top',zorder=999999, **title_args)
+    #        r = fig.canvas.get_renderer()
+    #        bb = current_text.get_window_extent(renderer=r).transformed(axes[index].transAxes.inverted())
+    #        width = bb.width
+    #        height = bb.height
+#
+    #        rect = patches.Rectangle((0.0,1.0 - height),
+    #                        height=height + 0.02*2,
+    #                        width=width + 0.01*2, 
+    #                        transform=axes[index].transAxes,
+    #                        linewidth=1, 
+    #                        edgecolor='black', 
+    #                        facecolor='white',
+    #                        alpha=1.0,
+    #                        zorder=99998)
+    #        axes[index].add_patch(rect)
+    #        rect.set(zorder=99998) 
+    #else:
+    if title_pos is None:
+        title_pos = [0.01,0.97]
+    for index,title in enumerate(titles):
+        axes[index].text(title_pos[0], title_pos[1] , title, transform=axes[index].transAxes,ha="left", va='top',zorder=999999, **title_args)
     
-    if titles is not None:
-        for index,title in enumerate(titles):
-            axes[index].text(title_pos[0], title_pos[1] , title, transform=axes[index].transAxes,ha="left", va='top',**title_args)
     ##########################################################################################################################################
     # draw the line for lead variants
-    _draw_grid_line_for_lead_variants(mode, lead_variants_is, n_plot, axes, region_lead_grid_line)
+    _draw_grid_line_for_lead_variants(mode, lead_variants_is,lead_variants_is_color, n_plot, axes, region_lead_grid_line,region_chromatin_files)
     
     ##########################################################################################################################################  
     _drop_old_y_labels(axes, n_plot)
     
-    _add_new_y_label(mode, fig, gene_track_height,n_plot,subplot_height )
+    _add_new_y_label(mode, fig, gene_track_height,n_plot,subplot_height ,fontsize,font_family)
     
     ##########################################################################################################################################
     save_figure(fig = fig, save = save, keyword= "stacked_" + mode, save_args=save_args, log = log, verbose=verbose)
@@ -208,22 +275,30 @@ def _drop_old_y_labels(axes, n_plot):
     for index in range(n_plot):
         axes[index].set_ylabel("")
 
-def _draw_grid_line_for_lead_variants(mode, lead_variants_is, n_plot, axes, region_lead_grid_line):
+def _draw_grid_line_for_lead_variants(mode, lead_variants_is,lead_variants_is_color, n_plot, axes, region_lead_grid_line,region_chromatin_files):
+    if len(region_chromatin_files)>0:
+        n_plot_and_track = n_plot+2
+    else:
+        n_plot_and_track = n_plot+1
     if mode=="r":
         for index, sig_is in lead_variants_is.items():
-            for sig_i in sig_is:
+            for j, sig_i in enumerate(sig_is):
+                try:
+                    region_lead_grid_line["color"]=lead_variants_is_color[index][j]
+                except:
+                    pass
                 if sig_i is not None:
-                    for each_axis_index in range(n_plot + 1):    
+                    for each_axis_index in range(n_plot_and_track):    
                         axes[each_axis_index].axvline(x=sig_i, zorder=2,**region_lead_grid_line)
 
-def _add_new_y_label(mode, fig, gene_track_height,n_plot,subplot_height ):
+def _add_new_y_label(mode, fig, gene_track_height,n_plot,subplot_height ,fontsize,font_family):
     gene_track_height_ratio = gene_track_height/(gene_track_height + n_plot*subplot_height)
     ylabel_height = (1 - gene_track_height_ratio)*0.5 + gene_track_height_ratio
     if mode=="r":
-        fig.text(0.08, ylabel_height , "$-log_{10}(P)$", va='center', rotation='vertical')
-        fig.text(0.93, ylabel_height, "Recombination rate(cM/Mb)", va='center', rotation=-90)
+        fig.text(0.08, ylabel_height , "$-log_{10}(P)$", va='center', rotation='vertical',fontsize=fontsize,family=font_family)
+        fig.text(0.93, ylabel_height, "Recombination rate(cM/Mb)", va='center', rotation=-90,fontsize=fontsize,family=font_family)
     elif mode=="m":
-        fig.text(0.08, ylabel_height , "$-log_{10}(P)$", va='center', rotation='vertical')    
+        fig.text(0.08, ylabel_height , "$-log_{10}(P)$", va='center', rotation='vertical',fontsize=fontsize,family=font_family)    
 
 def _sort_args(mqq_args, n_plot):
     mqq_args_for_each_plot={i:{} for i in range(n_plot)}

{gwaslab-3.4.45.dist-info → gwaslab-3.4.47.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: gwaslab
-Version: 3.4.45
+Version: 3.4.47
 Summary: A collection of handy tools for GWAS SumStats
 Author-email: Yunye <yunye@gwaslab.com>
 Project-URL: Homepage, https://cloufield.github.io/gwaslab/
@@ -13,8 +13,8 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 License-File: LICENSE_before_v3.4.39
 Requires-Dist: pandas !=1.5,>=1.3
-Requires-Dist: numpy >=1.21.2
-Requires-Dist: matplotlib !=3.7.2,>=3.5
+Requires-Dist: numpy <2,>=1.21.2
+Requires-Dist: matplotlib !=3.7.2,<3.9,>=3.5
 Requires-Dist: seaborn >=0.12
 Requires-Dist: scipy >=1.12
 Requires-Dist: pySAM ==0.22.1
@@ -51,7 +51,7 @@ Warning: Known issues of GWASLab are summarized in [https://cloufield.github.io/
 ### install via pip
 
 ```
-pip install gwaslab==3.4.43
+pip install gwaslab==3.4.45
 ```
 
 ```python
@@ -90,7 +90,7 @@ 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.43
+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)

{gwaslab-3.4.45.dist-info → gwaslab-3.4.47.dist-info}/RECORD

@@ -1,25 +1,25 @@
-gwaslab/__init__.py,sha256=dFnrh4L620F5JirsSF98SmkuligA-fybIGdBF6r9Ims,2386
-gwaslab/bd_common_data.py,sha256=v98X3tdRNOVE2gCiSHkfyBb0pSIjTk5IFG8A725Oj3o,12639
+gwaslab/__init__.py,sha256=7TKJaODdpeuQKibL7gIEa4MtyQ0pmrU-vIHQ-Et27lQ,2433
+gwaslab/bd_common_data.py,sha256=qr6OMbBaTH2Smfu8347SO9NmF410tn8dq8pRGF5-OpY,13751
 gwaslab/bd_config.py,sha256=TP-r-DPhJD3XnRYZbw9bQHXaDIkiRgK8bG9HCt-UaLc,580
 gwaslab/bd_download.py,sha256=cDDk2C5IvjeAzvPvVYGTkI4Ss33DUtEDjGo8eAbQRvY,15663
 gwaslab/bd_get_hapmap3.py,sha256=asNjQYeGfQi8u3jnfenRvDdKMs5ptql5wpcUzqMlwUI,3937
 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=NOEQd00guGch_GIt5bHv1wcrAvETfChqzmtgm-nIx_I,35298
+gwaslab/g_Sumstats.py,sha256=TUcFQFyODS_-FYMdXDvrBijG4Qtfi1igIWM-eEgb0nc,35352
 gwaslab/g_SumstatsPair.py,sha256=20snPb4SlI6ftMGVjgxAuyxsxYRQF-GzzlBSnoB-3Lo,8851
 gwaslab/g_SumstatsT.py,sha256=u_DighLMnMxwTLnqm-B58pA0G6WXRj6pudPyKMVKjSU,2133
 gwaslab/g_Sumstats_summary.py,sha256=FECvvFXJVKaCX5dggBvvk9YvJ6AbdbcLfjltysX7wEE,6380
 gwaslab/g_meta.py,sha256=htWlgURWclm9R6UqFcX1a93WN27xny7lGUeyJZOtszQ,2583
 gwaslab/g_vchange_status.py,sha256=jLoVzMJFhB5k_cJKzHuBNc2HZGBWydAunCNa0n_d54g,1923
-gwaslab/g_version.py,sha256=49_gR8lEQ_jgmfO9XJszEzuzDIESj5dHj6gta3Ilkmw,1818
+gwaslab/g_version.py,sha256=_mujOAftzn4ic9p92turA93xxg-RqxoWIqVfJLiNaSg,1886
 gwaslab/hm_casting.py,sha256=FqP4EQl83Q2OKLw004OgLIvUH795TVCGwziLk5jsHqY,11368
-gwaslab/hm_harmonize_sumstats.py,sha256=ympk2MZkbb0MnZ1n2ajkV36L8EAm7nBEaYhjqjI38tU,78548
+gwaslab/hm_harmonize_sumstats.py,sha256=1hjUdle2DSKHGBp2BktfFqf-QHU_q2xWl_mPhiYc_ZA,78616
 gwaslab/hm_rsid_to_chrpos.py,sha256=ODWREO0jPN0RAfNzL5fRzSRANfhiksOvUVPuEsFZQqA,6552
-gwaslab/io_preformat_input.py,sha256=w62JLAr16Ru0EgUtBCEV2eXRO89OqhidQxwf2IPAM38,20014
+gwaslab/io_preformat_input.py,sha256=AZ43WGqVTzbo3XtClWhjRjsj6pBR9stw6JBL_TZ461U,20673
 gwaslab/io_read_ldsc.py,sha256=8S9n4imgl4d0WPms_GYld-6uUM5z7iWGiCA-M814kzY,12123
 gwaslab/io_read_tabular.py,sha256=EG-C6KhCutt4J4LlOMgXnqzJvU-EZXzVhMvaDFnHrMM,2380
-gwaslab/io_to_formats.py,sha256=QuGWdvnAamaZAuhymj-0SuNBaKz1maTTyH396gvVaO8,29229
+gwaslab/io_to_formats.py,sha256=m57dGoqmHzAE1E27j9YxYKVCA12_lKd1qCnZtp0WZLw,29401
 gwaslab/io_to_pickle.py,sha256=HhePU0VcaGni0HTNU0BqoRaOnrr0NOxotgY6ISdx3Ck,1833
 gwaslab/ldsc_irwls.py,sha256=83JbAMAhD0KOfpv4IJa6LgUDfQjp4XSJveTjnhCBJYQ,6142
 gwaslab/ldsc_jackknife.py,sha256=XrWHoKS_Xn9StG1I83S2vUMTertsb-GH-_gOFYUhLeU,17715
@@ -28,16 +28,16 @@ gwaslab/ldsc_parse.py,sha256=MBnfgcWlV4oHp9MoDRh1mpilaHhAR15Af77hMFn4-5k,10564
 gwaslab/ldsc_regressions.py,sha256=yzbGjgNV7u-SWXNPsh9S8y9mK97Bim_Nmad9G9V18ZU,30078
 gwaslab/ldsc_sumstats.py,sha256=O0olsDxKlh1MJ1gAuEN1t40rxhajOEwOQ20ak7xoDrI,26245
 gwaslab/qc_check_datatype.py,sha256=kW68uk4dTLOU2b1dHoVat6n0loundDysAjIqxsXW28Q,3379
-gwaslab/qc_fix_sumstats.py,sha256=cpJibJ_77p4cg39R4zRunhOK2deIK4PfQA9wmYZgyqk,92745
+gwaslab/qc_fix_sumstats.py,sha256=-DQz5dPW6YXXVP-LV2Txa4lJrpZHhqAoKNny6IYAW18,93100
 gwaslab/run_script.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-gwaslab/util_ex_calculate_ldmatrix.py,sha256=LpE__LoYRHLgVKlCHo6lYWlz9LEUVUDqYPEAP-Svbm0,14598
-gwaslab/util_ex_calculate_prs.py,sha256=5l1eiZs8YwIpEgp7i3IurP8n5KwQM5awbG9fWSm4iT4,9053
+gwaslab/util_ex_calculate_ldmatrix.py,sha256=Z_spxbq6SHDS0v84I59YTTF40iyLQIOZbt0dmEcNJjw,15417
+gwaslab/util_ex_calculate_prs.py,sha256=9uJ588Sdj4V0vw3OZ9NeLECwOvW67f0IdLandVPS5RY,9442
 gwaslab/util_ex_gwascatalog.py,sha256=--Gde9HrsikfYTeFqSaYz0gUODr9wmv_gV6GZGNlElg,7688
 gwaslab/util_ex_ldproxyfinder.py,sha256=wWNW9wITWozj23gT41LR00WxU-rrHpGKbxs2H_3jEyM,9431
 gwaslab/util_ex_ldsc.py,sha256=a-OwyWmDcyinTPRfCNSeNYm6beVsydFeYmRAXAPcRnI,18350
 gwaslab/util_ex_plink_filter.py,sha256=pK1Yxtv9-J4rMOdVAG7VU9PktvI6-y4FxBiVEH0QuRs,1673
 gwaslab/util_ex_process_h5.py,sha256=ynFvo3zxgvOxWYL565v2IQf8P6iEuq7UlKQ_ULxrd6Y,2831
-gwaslab/util_ex_process_ref.py,sha256=BwS4hmzIWdgmRY_yiG6RJQJHjYcOR-2yZX7kU3s8lPw,16552
+gwaslab/util_ex_process_ref.py,sha256=GQ0ZEWLxGpHLdBs3tqnAqKn3Pqx1A1YvNbYrBLBvXeg,17126
 gwaslab/util_ex_run_2samplemr.py,sha256=5c0DGF694T9j0Y58L2I7pr1_Z1hfpaatIgix7P5oPA8,9127
 gwaslab/util_ex_run_clumping.py,sha256=Y8KeUKuSpGc1a8iI_VKUgVQpAVBd7O9F11p85jY4oW4,7868
 gwaslab/util_ex_run_coloc.py,sha256=u57h8wPbTCOf6aY5u5DpzK1gv7inuDT8a15UGo-1ras,6288
@@ -47,23 +47,27 @@ gwaslab/util_in_calculate_power.py,sha256=JfHJFg3tNF0f4NHgWlzVW2mSxCiP07mAHIyEfV
 gwaslab/util_in_convert_h2.py,sha256=a8Cbudt3xn9WP2bPc-7ysuowB-LYub8j8GeDXl7Lk7Q,6483
 gwaslab/util_in_correct_winnerscurse.py,sha256=Gp--yAQ8MMzdkWIvXP9C1BHVjZc-YzqHfYWhAj19w9w,2110
 gwaslab/util_in_fill_data.py,sha256=gdTwYA6FvBMnrtxAeL0lEj_Z0aGIoRNPScWDlJvZWeQ,14021
-gwaslab/util_in_filter_value.py,sha256=6yz3omukfqhmkfGZwGtr2BPQ6FcSTj4l6o7EhPMXRz0,22100
+gwaslab/util_in_filter_value.py,sha256=dY4X66N9A4MHCRHjPqLYFufMM91ggLRwUBf_nJYh8Lg,23605
 gwaslab/util_in_get_density.py,sha256=kpKXH69acMkeYVG5vs-VbJC3COhmuLBfYco-wuOxgjc,3934
-gwaslab/util_in_get_sig.py,sha256=atyBJZCWGUSgy-nvIR8_a_isseq1nKhzTaRVG2LbKQk,37762
+gwaslab/util_in_get_sig.py,sha256=9kq1GXacknO2YnVmsTli1GlPA728ASweTZ3UKm3Wszo,38783
+gwaslab/util_in_meta.py,sha256=5K9lIZcIgUy0AERqHy1GvMN2X6dp45JUUgopuDLgt4o,11284
+gwaslab/util_in_snphwe.py,sha256=-KpIDx6vn_nah6H55IkV2OyjXQVXV13XyBL069WE1wM,1751
 gwaslab/viz_aux_annotate_plot.py,sha256=R-1GT89E4NEBAMNTYzNawdi9rjQV5LCnODgnYOOKsys,32184
+gwaslab/viz_aux_chromatin.py,sha256=7cGmej5EkKO7fxR1b5w8r1oRRl9ofVzFRG52SCYWtz0,4109
 gwaslab/viz_aux_quickfix.py,sha256=Z6ZNEAUFuWVDTzH-qGreNGxPxJLCmqhXtBrvDOgo4g8,18308
 gwaslab/viz_aux_reposition_text.py,sha256=iRIP-Rkltlei068HekJcVubiqPrunBqvAoSQ1eHk04M,4304
 gwaslab/viz_aux_save_figure.py,sha256=nL-aoE8Kg06h7FgleGRBIZjhI-6w5gpn3E1HWMwBig8,2664
 gwaslab/viz_plot_compare_af.py,sha256=qtXW45-Sq_ugK8ZfqBYMpmf58SKi3lB3YyHnzn_akcE,5344
-gwaslab/viz_plot_compare_effect.py,sha256=8om3y6YQfnOk4FfkKSpKr2KqJcsMeCwQ6FRRKbDrm3U,49366
+gwaslab/viz_plot_compare_effect.py,sha256=iA74jMzh-G65U6BeXyQro08tPlJWpNyvtrjFsYHLvFM,49505
 gwaslab/viz_plot_forestplot.py,sha256=xgOnefh737CgdQxu5naVyRNBX1NQXPFKzf51fbh6afs,6771
 gwaslab/viz_plot_miamiplot.py,sha256=rCFEp7VNuVqeBBG3WRkmFAtFklbF79BvIQQYiSY70VY,31238
 gwaslab/viz_plot_miamiplot2.py,sha256=SWv82D8UBbREKsk8EoKth-2w68l6FbXyVLsb_E1hh8o,15882
-gwaslab/viz_plot_mqqplot.py,sha256=PzRWnm11whxww7ut-bzFkj1sbPc_c0OP7yRpIgYo2iQ,61739
+gwaslab/viz_plot_mqqplot.py,sha256=udOTBKSX239KLFNauhHGDJHp0Viu18NkGwFn5GtAHuM,63297
 gwaslab/viz_plot_qqplot.py,sha256=psQgVpP29686CEZkzQz0iRbApzqy7aE3GGiBcazVvNw,7247
-gwaslab/viz_plot_regionalplot.py,sha256=PBIWkNj2fj-dRLKQJNpM8wor5jya2anqix0-UYLE0Is,37901
+gwaslab/viz_plot_regional2.py,sha256=6Dfbq5vFm5B63nLUpaTrlq96x4CruIgJu0TXXq3a_Ck,34546
+gwaslab/viz_plot_regionalplot.py,sha256=8u-5-yfy-UaXhaxVVz3Y5k2kBAoqzczUw1hyyD450iI,37983
 gwaslab/viz_plot_rg_heatmap.py,sha256=PidUsgOiEVt6MfBPCF3_yDhOEytZ-I1q-ZD6_0pFrV4,13713
-gwaslab/viz_plot_stackedregional.py,sha256=EAHz5SZGengZ_pxcOg62ZRimGRz6_goQlp9MLCwmeZc,11890
+gwaslab/viz_plot_stackedregional.py,sha256=HISwtPDV84mT3PS41e2sXq77GPBujgSKVX1mQ4bEeKs,15820
 gwaslab/viz_plot_trumpetplot.py,sha256=ZHdc6WcVx0-oKoj88yglRkmB4bS9pOiEMcuwKW35Yvo,42672
 gwaslab/data/formatbook.json,sha256=N2nJs80HH98Rsu9FxaSvIQO9J5yIV97WEtAKjRqYwiY,38207
 gwaslab/data/reference.json,sha256=k8AvvgDsuLxzv-NCJHWvTUZ5q_DLAFxs1Th3jtL313k,11441
@@ -73,9 +77,9 @@ gwaslab/data/hapmap3_SNPs/hapmap3_db150_hg19.snplist.gz,sha256=qD9RsC5S2h6l-OdpW
 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.4.45.dist-info/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-gwaslab-3.4.45.dist-info/LICENSE_before_v3.4.39,sha256=GhLOU_1UDEKeOacYhsRN_m9u-eIuVTazSndZPeNcTZA,1066
-gwaslab-3.4.45.dist-info/METADATA,sha256=5FN5dbVypNPET635Eooi01_1NDFD1dNr1T9Jv0JXmLc,7757
-gwaslab-3.4.45.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-gwaslab-3.4.45.dist-info/top_level.txt,sha256=PyY6hWtrALpv2MAN3kjkIAzJNmmBTH5a2risz9KwH08,8
-gwaslab-3.4.45.dist-info/RECORD,,
+gwaslab-3.4.47.dist-info/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+gwaslab-3.4.47.dist-info/LICENSE_before_v3.4.39,sha256=GhLOU_1UDEKeOacYhsRN_m9u-eIuVTazSndZPeNcTZA,1066
+gwaslab-3.4.47.dist-info/METADATA,sha256=G_RLrYWZe0z6Wbvyl_g7rXnJKoZ6HFKUHM66ogtfyug,7765
+gwaslab-3.4.47.dist-info/WHEEL,sha256=y4mX-SOX4fYIkonsAGA5N0Oy-8_gI4FXw5HNI1xqvWg,91
+gwaslab-3.4.47.dist-info/top_level.txt,sha256=PyY6hWtrALpv2MAN3kjkIAzJNmmBTH5a2risz9KwH08,8
+gwaslab-3.4.47.dist-info/RECORD,,

{gwaslab-3.4.45.dist-info → gwaslab-3.4.47.dist-info}/WHEEL

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