gwaslab 3.6.7__tar.gz → 3.6.9__tar.gz

{gwaslab-3.6.7/src/gwaslab.egg-info → gwaslab-3.6.9}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gwaslab
-Version: 3.6.7
+Version: 3.6.9
 Summary: A collection of handy tools for GWAS SumStats
 Author-email: Yunye <gwaslab@gmail.com>
 Project-URL: Homepage, https://cloufield.github.io/gwaslab/

{gwaslab-3.6.7 → gwaslab-3.6.9}/pyproject.toml

@@ -7,7 +7,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "gwaslab"
-version = "3.6.7"
+version = "3.6.9"
 authors = [
   { name="Yunye", email="gwaslab@gmail.com" },
 ]

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/g_Sumstats.py

@@ -102,7 +102,7 @@ from gwaslab.viz.viz_plot_mqqplot import mqqplot
 from gwaslab.viz.viz_plot_trumpetplot import plottrumpet
 from gwaslab.viz.viz_plot_compare_af import plotdaf
 from gwaslab.viz.viz_plot_credible_sets import _plot_cs
-
+from gwaslab.viz.viz_plot_associations import _plot_associations
 from gwaslab.io.io_read_pipcs import _read_pipcs
 from gwaslab.io.io_load_ld import tofinemapping_using_ld
 from gwaslab.io.io_preformat_input import preformat
@@ -811,6 +811,9 @@ class Sumstats():
         self.associations = associations_full
         
         return associations_summary
+    def plot_associations(self,**kwargs):
+
+        _plot_associations(self.associations, **kwargs)
 
     def check_cis(self, gls=False, **kwargs):
         if "SNPID" in self.data.columns:

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/g_version.py

@@ -17,8 +17,8 @@ def _get_version():
 def gwaslab_info():
     # version meta information
     dic={
-   "version":"3.6.7",
-   "release_date":"20250827"
+   "version":"3.6.9",
+   "release_date":"20251012"
     }
     return dic   
 

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/util/util_ex_phewwas.py

@@ -7,6 +7,11 @@ def _extract_associations(sumstats, rsid="rsID", log = Log(), verbose=True):
 
     assoc, traits, studies, variants = get_associations_from_gwascatalog(sumstats, rsid=rsid, log=log, verbose=verbose)
     
+    if len(assoc)==0:
+        # if no associations
+        log.write("No associations!")
+        return None, None
+    
     assoc = _fix_beta(assoc)
 
     traits_agg = traits.groupby("associationId")[["trait","shortForm"]].agg(lambda x: ",".join(x)).reset_index()
@@ -14,21 +19,25 @@ def _extract_associations(sumstats, rsid="rsID", log = Log(), verbose=True):
     assoc_traits_agg= pd.merge(assoc, traits_agg, on ="associationId",how="left")
     
     assoc_traits_agg= pd.merge(assoc_traits_agg, studies, on ="associationId", how="left")
-
+    
     assoc_traits_agg= pd.merge(assoc_traits_agg, variants, on ="associationId",how="left")
     
     assoc_traits_agg = assoc_traits_agg.rename(columns={"trait":"GWASCATALOG_TRAIT",
                                                         "riskFrequency":"RAF",
                                                         "betaNum":"Beta",
-                                                        "pvalue":"P-value"
+                                                        "pvalue":"P-value",
+                                                        "betaUnit":"Unit"
                                                         })
     
     summary_columns=['GWASCATALOG_TRAIT','associationId', 'rsID', "geneName", 
-                     'RA', 'RAF','Beta', 'P-value','cohort','initialSampleSize','publicationInfo.pubmedId',
+                     'RA', 'RAF','Beta',"Unit", 'P-value','cohort','initialSampleSize','publicationInfo.pubmedId',
                      "functionalClass","gene.geneName"]
     
     assoc_traits_agg_summary = assoc_traits_agg[summary_columns]
 
+    assoc_traits_agg_summary = _align_beta(sumstats, assoc_traits_agg_summary, log, verbose)
+    assoc_traits_agg = _align_beta(sumstats, assoc_traits_agg, log, verbose)
+
     return assoc_traits_agg, assoc_traits_agg_summary
 
 def get_associations_from_gwascatalog(sumstats, rsid="rsID", log=Log(), verbose=True):
@@ -101,7 +110,12 @@ def get_associations_from_gwascatalog(sumstats, rsid="rsID", log=Log(), verbose=
     return association, traits, studies, variants
 
 def _fix_beta(association):
-
+    if "betaNum" not in association:
+        association["betaNum"] = pd.NA
+    if "orPerCopyNum" not in association:
+        association["orPerCopyNum"] = pd.NA
+    if "range" not in association:
+        association["range"] = pd.NA
     is_or_available = (association["betaNum"].isna()) & (~association["orPerCopyNum"].isna())
     is_range_available = (association["betaNum"].isna()) & (association["orPerCopyNum"].isna()) & (~association["range"].isna())
 
@@ -114,4 +128,42 @@ def parse_range(x):
     high = np.log(range_list[1])
     low = np.log(range_list[0])
     beta = (high + low)/2
-    return beta
+    return beta
+
+def _align_beta(sumstats, assoc_traits_agg_summary, log, verbose):
+    # merge with sumstats
+    log.write("Merging GWAS Catalog associationms with Sumstats...",)
+    assoc_traits_agg_summary = pd.merge(sumstats[["rsID","CHR","POS","EA","NEA","BETA"]], assoc_traits_agg_summary, on = "rsID", how="inner")
+    # check direction 
+    assoc_traits_agg_summary["GC_Beta_GL_RA"] = np.nan
+    assoc_traits_agg_summary["GC_Beta_GL_EA"] = np.nan
+
+    # get sumstats RA
+    assoc_traits_agg_summary["_GL_RA_BETA"] = assoc_traits_agg_summary["BETA"].abs()
+    assoc_traits_agg_summary["_GL_RA"] = assoc_traits_agg_summary["EA"]
+    assoc_traits_agg_summary["_GL_NRA"] = assoc_traits_agg_summary["NEA"]
+    is_GL_RA_flipped = assoc_traits_agg_summary["BETA"]<0
+    assoc_traits_agg_summary.loc[is_GL_RA_flipped, "_GL_RA"] = assoc_traits_agg_summary.loc[is_GL_RA_flipped, "NEA"]
+    assoc_traits_agg_summary.loc[is_GL_RA_flipped, "_GL_NRA"] = assoc_traits_agg_summary.loc[is_GL_RA_flipped, "EA"]
+    
+    is_matched = assoc_traits_agg_summary["_GL_RA"] == assoc_traits_agg_summary["RA"]
+    log.write(" -GL_RA matched RA for {} variants ".format(sum(is_matched)), verbose=verbose)
+    assoc_traits_agg_summary.loc[is_matched, "GC_Beta_GL_RA"] = assoc_traits_agg_summary.loc[is_matched, "Beta"]
+    
+    is_flipped = assoc_traits_agg_summary["_GL_NRA"] == assoc_traits_agg_summary["RA"]
+    log.write(" -GL_NRA matched RA for {} variants... Aligning with Sumstats Risk allele (BETA>0) in GC_Beta_GL_RA... ".format(sum(is_flipped)), verbose=verbose)
+    assoc_traits_agg_summary.loc[is_flipped, "GC_Beta_GL_RA"] = - assoc_traits_agg_summary.loc[is_flipped, "Beta"]
+
+    is_matched = assoc_traits_agg_summary["EA"] == assoc_traits_agg_summary["RA"]
+    log.write(" -EA matched RA for {} variants ".format(sum(is_matched)), verbose=verbose)
+    assoc_traits_agg_summary.loc[is_matched, "GC_Beta_GL_EA"] = assoc_traits_agg_summary.loc[is_matched, "Beta"]
+    
+    is_flipped = assoc_traits_agg_summary["NEA"] == assoc_traits_agg_summary["RA"]
+    log.write(" -NEA matched RA for {} variants... Aligning with Sumstats Risk allele (BETA>0) in GC_Beta_GL_RA... ".format(sum(is_flipped)), verbose=verbose)
+    assoc_traits_agg_summary.loc[is_flipped, "GC_Beta_GL_EA"] = - assoc_traits_agg_summary.loc[is_flipped, "Beta"]
+
+    not_match = (assoc_traits_agg_summary["_GL_RA"] != assoc_traits_agg_summary["RA"])& (assoc_traits_agg_summary["_GL_NRA"] != assoc_traits_agg_summary["RA"])
+    log.write(" -NEA and EA not matching RA for {} variants... ".format(sum(not_match)), verbose=verbose)
+
+
+    return assoc_traits_agg_summary

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/viz/viz_aux_save_figure.py

@@ -53,7 +53,8 @@ def get_default_path(keyword,fmt="png"):
                         "afc":"allele_frequency_comparision",
                         "gwheatmap":"genome_wide_heatmap",
                         "scatter":"scatter",
-                        "forest":"forest"
+                        "forest":"forest",
+                        "gwaslab":"gwaslab"
                         }
     prefix = path_dictionary[keyword]
     count = 1

gwaslab-3.6.9/src/gwaslab/viz/viz_plot_associations.py

@@ -0,0 +1,61 @@
+import seaborn as sns
+import matplotlib.pyplot as plt
+from gwaslab.io.io_process_args import _update_args
+from gwaslab.io.io_process_args import _update_arg
+from gwaslab.g_Log import Log
+
+def _plot_associations(associations, 
+                       values="Beta",
+                        fig_args=None,
+                        log=Log(),
+                        verbose=True,
+                        sort="P-value",
+                        fontsize=12,
+                        font_family="Arial",
+                        cmap=None,
+                        ylabel="Traits",
+                        xlabel="rsID - Gene Name",
+                        **args               
+                       ):
+    
+    log.write("Start to create heatmap for associations...", verbose=verbose)
+    
+    cmap = _update_arg(cmap, "RdBu")
+
+    log.write(f" -Value to plot: {values}", verbose=verbose)
+    log.write(f" -Total number of associations :{len(associations)} ", verbose=verbose)
+    log.write(f" -Sorting associations by :{sort} ", verbose=verbose)
+
+    associations = associations.sort_values(by=[sort]).drop_duplicates(subset=["rsID","GWASCATALOG_TRAIT"])
+
+    log.write(f" -Keeping unique rsID-traits pairs :{len(associations)} ", verbose=verbose)
+
+    associations = associations.dropna(subset=[values])
+    log.write(f" -Keeping associations without NA in {values} :{len(associations)} ", verbose=verbose)
+
+    log.write(f" -Total number of unique variants for plotting :{associations["rsID"].nunique()} ", verbose=verbose)
+    log.write(f" -Total number of unique traits for plotting:{associations["GWASCATALOG_TRAIT"].nunique()} ", verbose=verbose)
+
+    matrix_beta = associations.pivot_table(index=['rsID','gene.geneName'], 
+                                    columns='GWASCATALOG_TRAIT', 
+                                    values=[values])
+
+    matrix_beta = matrix_beta.astype("float64")
+    
+    matrix_beta.columns = matrix_beta.columns.droplevel(0)
+
+    height = max(2, len(matrix_beta)//2)
+    width = max(2, len(matrix_beta.columns)//2)
+
+    fig_args = _update_args(fig_args, dict(figsize=(width,height)))
+    fig,ax = plt.subplots(**fig_args)
+    
+    sns.heatmap(matrix_beta.T, annot=True, fmt=".2f",cmap="RdBu",ax=ax)
+    
+
+    ax.set_xlabel(xlabel,fontsize=fontsize, fontfamily=font_family)
+    ax.set_ylabel(ylabel,fontsize=fontsize, fontfamily=font_family)
+
+
+    log.write("Finished creating heatmap for associations.", verbose=verbose)
+    return fig, ax

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/viz/viz_plot_mqqplot.py

@@ -100,6 +100,8 @@ def mqqplot(insumstats,
           qq_scatter_args=None,
           qqscatterargs=None,
           qq_line_color = "grey",
+          qq_xlabels = None,
+          qq_xlim = None,
           # region
           region = None,
           region_title=None,
@@ -204,6 +206,8 @@ def mqqplot(insumstats,
           highlight_color="#CB132D",
           highlight_windowkb = 500,
           highlight_anno_args = None,
+          highlight_lim = None,
+          highlight_lim_mode = "absolute",
           pinpoint= None,
           pinpoint_color ="red",
           stratified=False,
@@ -531,6 +535,8 @@ def mqqplot(insumstats,
                                                     highlight=highlight, 
                                                     highlight_chrpos=highlight_chrpos, 
                                                     highlight_windowkb=highlight_windowkb, 
+                                                    highlight_lim = highlight_lim, 
+                                                    highlight_lim_mode = highlight_lim_mode,
                                                     snpid=snpid, 
                                                     chrom=chrom, 
                                                     pos=pos)
@@ -1096,6 +1102,8 @@ def mqqplot(insumstats,
                     linewidth=sc_linewidth,
                     ytick3 = ytick3,
                     ylabels = ylabels,
+                    xlabels = qq_xlabels,
+                    xlim = qq_xlim,
                     ylabels_converted = ylabels_converted,
                     verbose=verbose,
                     qq_scatter_args=qq_scatter_args,
@@ -1335,19 +1343,38 @@ def _process_p_value(sumstats, mode,p, mlog10p, scaled, log, verbose ):
 
     return sumstats
 
-def _process_highlight(sumstats, highlight, highlight_chrpos, highlight_windowkb, snpid, chrom, pos):
+def _process_highlight(sumstats, highlight, highlight_chrpos, highlight_windowkb,  highlight_lim,  highlight_lim_mode, snpid, chrom, pos):
         if pd.api.types.is_list_like(highlight[0]):
             if highlight_chrpos == False:
                 # highlight for multiple sets
                 for i, highlight_set in enumerate(highlight):
                     to_highlight = sumstats.loc[sumstats[snpid].isin(highlight_set),:]
                     #assign colors: 0 is hightlight color
-                    for index,row in to_highlight.iterrows():
+                    #assign colors: i is highlight color
+                    for j, (index, row) in enumerate(to_highlight.iterrows()):
                         target_chr = int(row[chrom])
                         target_pos = int(row[pos])
                         right_chr=sumstats[chrom]==target_chr
-                        up_pos=sumstats[pos]>target_pos-highlight_windowkb*1000
-                        low_pos=sumstats[pos]<target_pos+highlight_windowkb*1000
+                        # Check if highlight_lim is provided for this group and SNP
+                        if (highlight_lim is not None and 
+                            i < len(highlight_lim) and 
+                            j < len(highlight_lim[i]) and 
+                            highlight_lim[i][j] is not None):
+                            # Use custom limits for this SNP
+                            if highlight_lim_mode == "absolute":
+                                # Absolute positions: (start_pos, end_pos)
+                                start_pos, end_pos = highlight_lim[i][j]
+                                up_pos = sumstats[pos] >= start_pos
+                                low_pos = sumstats[pos] <= end_pos
+                            else:
+                                # Offset mode: (lower_kb, upper_kb) relative to SNP position
+                                lower_kb, upper_kb = highlight_lim[i][j]
+                                up_pos = sumstats[pos] > target_pos + lower_kb * 1000
+                                low_pos = sumstats[pos] < target_pos + upper_kb * 1000
+                        else:
+                            # Use default highlight_windowkb
+                            up_pos = sumstats[pos] > target_pos - highlight_windowkb * 1000
+                            low_pos = sumstats[pos] < target_pos + highlight_windowkb * 1000
                         sumstats.loc[right_chr&up_pos&low_pos,"HUE"]=i
             else:
                 for i, highlight_chrpos_tuple in enumerate(highlight):
@@ -1371,13 +1398,30 @@ def _process_highlight(sumstats, highlight, highlight_chrpos, highlight_windowkb
         else:
             # highlight for one set
             to_highlight = sumstats.loc[sumstats[snpid].isin(highlight),:]
-            #assign colors: 0 is hightlight color
-            for index,row in to_highlight.iterrows():
+            #assign colors: 0 is highlight color
+            for j, (index, row) in enumerate(to_highlight.iterrows()):
                 target_chr = int(row[chrom])
                 target_pos = int(row[pos])
                 right_chr=sumstats[chrom]==target_chr
-                up_pos=sumstats[pos]>target_pos-highlight_windowkb*1000
-                low_pos=sumstats[pos]<target_pos+highlight_windowkb*1000
+                # Check if highlight_lim is provided for this SNP
+                if (highlight_lim is not None and 
+                    j < len(highlight_lim) and 
+                    highlight_lim[j] is not None):
+                    # Use custom limits for this SNP
+                    if highlight_lim_mode == "absolute":
+                        # Absolute positions: (start_pos, end_pos)
+                        start_pos, end_pos = highlight_lim[j]
+                        up_pos = sumstats[pos] >= start_pos
+                        low_pos = sumstats[pos] <= end_pos
+                    else:
+                        # Offset mode: (lower_kb, upper_kb) relative to SNP position
+                        lower_kb, upper_kb = highlight_lim[j]
+                        up_pos = sumstats[pos] > target_pos + lower_kb * 1000
+                        low_pos = sumstats[pos] < target_pos + upper_kb * 1000
+                else:
+                    # Use default highlight_windowkb
+                    up_pos = sumstats[pos] > target_pos - highlight_windowkb * 1000
+                    low_pos = sumstats[pos] < target_pos + highlight_windowkb * 1000
                 sumstats.loc[right_chr&up_pos&low_pos,"HUE"]=0
         return sumstats
 

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/viz/viz_plot_qqplot.py

@@ -39,6 +39,8 @@ def _plot_qq(
     linewidth,
     ytick3,
     ylabels,
+    xlabels,
+    xlim,
     ylabels_converted,
     qq_scatter_args,
     expected_min_mlog10p,
@@ -155,32 +157,12 @@ def _plot_qq(
                         log=log,
                         verbose=verbose
                         )
-
-    #if cut == 0: 
-    #    ax2.set_ylim(skip, ceil(maxy*1.2) )
-    #    
-    #if cut:
-    #    qcutline = ax2.axhline(y=cut, linewidth = linewidth, linestyle="--",color=cut_line_color,zorder=1)
-    #    step=2
-    #    if ((maxticker-cut)/cutfactor + cut) > cut:
-    #        if ystep == 0:
-    #            if (cut - skip ) // step > 10:
-    #                step = (cut - skip ) // 10
-    #        else:
-    #            step = ystep
-#
-    #        ax2.set_yticks([x for x in range(skip,cut-step,step)]+[cut]+[(maxticker-cut)/cutfactor + cut])
-    #        ax2.set_yticklabels([x for x in range(skip,cut-step,step)]+[cut]+[maxticker],fontsize=fontsize,family=font_family)
-    #    else:
-    #        if ystep == 0:
-    #            if (cut - skip ) // step > 10:
-    #                step = (cut - skip ) // 10
-    #        else:
-    #            step = ystep
-#
-    #        ax2.set_yticks([x for x in range(skip,cut-step,step)]+[cut])
-    #        ax2.set_yticklabels([x for x in range(skip,cut-step,step)]+[cut],fontsize=fontsize,family=font_family)
-    #    ax2.set_ylim(bottom = skip)
+    if xlim is not None:
+        ax2.set_xlim(xlim)
+        
+    if xlabels is not None:
+        ax2.set_xticks(xlabels)
+        ax2.set_xticklabels(xlabels,fontsize=fontsize,family=font_family)
 
     ax2.tick_params(axis='both', which='both', labelsize=fontsize,labelfontfamily=font_family)
     #

gwaslab-3.6.9/src/gwaslab/viz/viz_plot_template.py

@@ -0,0 +1,98 @@
+import seaborn as sns
+import matplotlib.pyplot as plt
+from gwaslab.io.io_process_args import _update_args
+from gwaslab.io.io_process_args import _update_arg
+from gwaslab.g_Log import Log
+from gwaslab.viz.viz_aux_save_figure import save_figure
+
+def _plot( associations, 
+            values="Beta",
+            sort="P-value",
+            mode="gwaslab",
+            fontsize=12,
+            font_family="Arial",
+            cmap=None,
+            ylabel="Y",
+            xlabel="X",
+            xlim=None,
+            ylim=None,
+            yticks=None,
+            xticks=None,
+            title="Title",
+            title_pad=1.08, 
+            title_fontsize=13,
+            linewidth=None,
+            linestyle=None,
+            linecolor=None,
+            dpi=200,
+            fig_args= None,
+            scatter_args=None,
+            save=None,
+            save_args=None,
+            log=Log(),
+            verbose=True,
+            **args               
+            ):
+    
+    # update args
+    
+    cmap = _update_arg(cmap, "RdBu")
+
+    fig_args = _update_args(fig_args, dict(figsize=(10,10)))
+
+    # create fig and ax    
+    fig,ax = plt.subplots(**fig_args)
+
+    # draw lines
+    horizontal_line = ax.axhline(y=1, linewidth = linewidth,
+                            linestyle="--",
+                            color=linecolor,zorder=1000)
+    vertical_line   = ax.axvline(x=1, linewidth = linewidth,
+                            linestyle="--",
+                            color=linecolor,zorder=1000)
+
+    # ticks
+    if xticks is not None:
+        ax.set_xticks(xticks)
+        ax.set_xticklabels(xticks,
+                            fontsize=fontsize,
+                            family=font_family)    
+    
+    if yticks is not None:
+        ax.set_xticks(yticks)
+        ax.set_xticklabels(yticks,
+                            fontsize=fontsize,
+                            family=font_family)    
+
+    # labels
+    ax.set_xlabel(xlabel,
+                  fontsize=fontsize, 
+                  fontfamily=font_family)
+    ax.set_ylabel(ylabel,
+                  fontsize=fontsize, 
+                  fontfamily=font_family)
+
+    ax.tick_params(axis='x', 
+                    labelsize=fontsize,
+                    labelfontfamily=font_family)
+    ax.tick_params(axis='y', 
+                    labelsize=fontsize,
+                    labelfontfamily=font_family)
+    
+    # title
+    title_pad = title_pad -0.05
+    fig.suptitle(title , 
+                 fontsize = title_fontsize, 
+                 x=0.5,
+                 y=title_pad)
+    
+    # spines
+    ax.spines["top"].set_visible(False)
+    ax.spines["right"].set_visible(False)
+    ax.spines["left"].set_visible(True)
+    ax.spines["bottom"].set_visible(True)
+
+    save_figure(fig = fig, save = save, keyword=mode, save_args=save_args, log = log, verbose=verbose)
+
+    log.write("Finished creating plots.", verbose=verbose)
+    return fig, ax

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab/viz/viz_plot_trumpetplot.py

@@ -57,6 +57,8 @@ def plottrumpet(mysumstats,
                 highlight_color="#CB132D",
                 highlight_windowkb = 500,
                 highlight_anno_args = None,
+                highlight_lim = None,
+                highlight_lim_mode = "absolute",
                 pinpoint= None,
                 pinpoint_color ="red",
                 scatter_args=None,
@@ -235,6 +237,8 @@ def plottrumpet(mysumstats,
                                                     highlight=highlight, 
                                                     highlight_chrpos=highlight_chrpos, 
                                                     highlight_windowkb=highlight_windowkb, 
+                                                    highlight_lim=highlight_lim,
+                                                    highlight_lim_mode=highlight_lim_mode,
                                                     snpid=snpid, 
                                                     chrom=chrom, 
                                                     pos=pos)

{gwaslab-3.6.7 → gwaslab-3.6.9/src/gwaslab.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gwaslab
-Version: 3.6.7
+Version: 3.6.9
 Summary: A collection of handy tools for GWAS SumStats
 Author-email: Yunye <gwaslab@gmail.com>
 Project-URL: Homepage, https://cloufield.github.io/gwaslab/

{gwaslab-3.6.7 → gwaslab-3.6.9}/src/gwaslab.egg-info/SOURCES.txt

@@ -19,7 +19,6 @@ src/gwaslab/g_meta_update.py
 src/gwaslab/g_vchange_status.py
 src/gwaslab/g_vchange_status_polars.py
 src/gwaslab/g_version.py
-src/gwaslab/hm_harmonize_sumstats.py
 src/gwaslab/run_script.py
 src/gwaslab.egg-info/PKG-INFO
 src/gwaslab.egg-info/SOURCES.txt
@@ -108,6 +107,7 @@ src/gwaslab/viz/viz_aux_property.py
 src/gwaslab/viz/viz_aux_quickfix.py
 src/gwaslab/viz/viz_aux_reposition_text.py
 src/gwaslab/viz/viz_aux_save_figure.py
+src/gwaslab/viz/viz_plot_associations.py
 src/gwaslab/viz/viz_plot_compare_af.py
 src/gwaslab/viz/viz_plot_compare_effect.py
 src/gwaslab/viz/viz_plot_credible_sets.py
@@ -123,4 +123,5 @@ src/gwaslab/viz/viz_plot_regionalplot.py
 src/gwaslab/viz/viz_plot_rg_heatmap.py
 src/gwaslab/viz/viz_plot_scatter_with_reg.py
 src/gwaslab/viz/viz_plot_stackedregional.py
+src/gwaslab/viz/viz_plot_template.py
 src/gwaslab/viz/viz_plot_trumpetplot.py