sciv 0.0.99__py3-none-any.whl → 0.0.101__py3-none-any.whl

sciv/plot/_bar_.py

@@ -76,6 +76,7 @@ def two_bar(
     bottom: float = 0,
     rotation: float = 65,
     text_left_move: float = 0.15,
+    y_limit: Tuple[float, float] = (0, 1),
     title: str = None,
     output: path = None,
     show: bool = True,
@@ -91,6 +92,9 @@ def two_bar(
 
     ax.legend()
 
+    ax.set_ylim(y_limit)
+
+    ax.set_xticks(range(len(ax_x)))
     ax.set_xticklabels(labels=list(ax_x), rotation=rotation)
 
     # Draw numerical values
@@ -103,6 +107,13 @@ def two_bar(
             color=text_color
         )
 
+    for spine in ["top", "left", "right", "bottom"]:
+        ax.spines[spine].set_linewidth(1)
+
+    ax.spines['bottom'].set_linewidth(1)
+    ax.grid(axis='y', ls='--', c='gray')
+    ax.set_axisbelow(True)
+
     plot_end(fig, title, x_name, y_name, output, show, close)
 
 
@@ -123,6 +134,7 @@ def class_bar(
     rotation: float = 65,
     title: str = None,
     text_left_move: float = 0.15,
+    y_limit: Tuple[float, float] = (0, 1),
     output: path = None,
     show: bool = True,
     close: bool = False,
@@ -159,6 +171,7 @@ def class_bar(
         bottom=bottom,
         rotation=rotation,
         text_left_move=text_left_move,
+        y_limit=y_limit,
         title=title,
         output=output,
         show=show,
@@ -185,6 +198,7 @@ def bar_trait(
     rotation: float = 65,
     title: str = None,
     text_left_move: float = 0.15,
+    y_limit: Tuple[float, float] = (0, 1),
     output: path = None,
     show: bool = True,
     close: bool = False,
@@ -216,6 +230,7 @@ def bar_trait(
             bottom=bottom,
             rotation=rotation,
             text_left_move=text_left_move,
+            y_limit=y_limit,
             text_color=text_color,
             output=os.path.join(output, f"cell_{trait_}_enrichment_bar.pdf") if output is not None else None,
             show=show,
@@ -224,6 +239,7 @@ def bar_trait(
         )
 
     trait_list = list(set(trait_df[trait_column_name]))
+
     # judge trait
     if trait_name != "All" and trait_name not in trait_list:
         ul.log(__name__).error(
@@ -237,8 +253,10 @@ def bar_trait(
 
     # plot
     if trait_name == "All":
+
         for trait in trait_list:
             trait_plot(trait_=trait, cell_df_=trait_df)
+
     else:
         trait_plot(trait_name, trait_df)
 
@@ -376,10 +394,9 @@ def bar_significance(
         y=y,
         hue=legend,
         hue_order=hue_order,
-        ci=ci,
+        errorbar=('ci', ci),
         capsize=capsize,
-        errwidth=line_width,
-        errcolor=errcolor,
+        err_kws={'color': errcolor, 'linewidth': line_width},
         ax=ax,
         palette=palette,
         edgecolor=errcolor,
@@ -437,7 +454,7 @@ def bar_significance(
     ax.set_axisbelow(True)
 
     if x_rotation != 0:
-        ax.set_xticklabels(labels=ax.get_xticklabels(), rotation=x_rotation)
+        ax.tick_params(axis='x', rotation=x_rotation)
 
     plt.legend(loc='upper left', bbox_to_anchor=(0.0, legend_gap), ncol=2)
 

sciv/plot/_box_.py

@@ -52,7 +52,6 @@ def box_base(
 
     if "color" in df_columns:
         new_df_color: DataFrame = df.groupby(group_columns, as_index=False)["color"].first()
-
         new_df = new_df.merge(new_df_color, how="left", on=clusters)
 
     colors: list = []
@@ -61,22 +60,31 @@ def box_base(
     if is_sort:
         new_df.sort_values([value], ascending=False, inplace=True)
         y_names: Union[list, None] = list(new_df[clusters])
+
         if "color" in df_columns:
-            colors = new_df["color"]
+            colors = list(new_df["color"])
+
     else:
         new_df.index = new_df[clusters]
+
         if order_names is not None:
             y_names: list = order_names
+
             if "color" in df_columns:
+
                 for i in order_names:
+
                     for j, c in zip(new_df[clusters], new_df["color"]):
+
                         if i == j:
                             colors.append(c)
                             break
+
         else:
             y_names = new_df[clusters]
+
             if "color" in df_columns:
-                colors = new_df["color"]
+                colors = list(new_df["color"])
 
     # scatter
     sns.boxplot(
@@ -103,6 +111,7 @@ def box_base(
         line.set_linewidth(line_width)
 
     # set coordinate
+    ax.set_xticks(range(len(y_names)))
     ax.set_xticklabels(labels=y_names, rotation=rotation)
     ax.spines['top'].set_linewidth(line_width)
     ax.spines['bottom'].set_linewidth(line_width)

sciv/plot/_core_.py

@@ -1,7 +1,7 @@
 # -*- coding: UTF-8 -*-
 
 import os
-from typing import Tuple, Literal, Optional, Union
+from typing import Tuple, Literal, Optional, Union, Any
 
 import pandas as pd
 from anndata import AnnData
@@ -314,9 +314,11 @@ def rate_bar_plot(
     rotation: float = 65,
     title: str = None,
     text_left_move: float = 0.15,
+    y_limit: Tuple[float, float] = (0, 1),
     plot_output: path = None,
     show: bool = True,
-    close: bool = False
+    close: bool = False,
+    **kwargs: Any
 ) -> None:
 
     if dir_name is not None:
@@ -349,9 +351,11 @@ def rate_bar_plot(
         bottom=bottom,
         rotation=rotation,
         text_left_move=text_left_move,
+        y_limit=y_limit,
         output=new_path if plot_output is not None else None,
         show=show,
-        close=close
+        close=close,
+        **kwargs
     )
 
 

sciv/plot/_graph_.py

@@ -143,8 +143,10 @@ def communities_graph(
     partition: list = [0 for _ in range(g.number_of_nodes())]
 
     for c_i, nodes in enumerate(communities):
+
         for i in nodes:
             partition[i] = type_colors[start_color_index + color_index * color_step_size + c_i]
+
         color_index += 1
 
     pos = nx.spring_layout(g)

sciv/plot/_kde_.py

@@ -42,7 +42,7 @@ def kde(
     # Random sampling
     if axis == -1:
         matrix = down_sampling_data(data.X, sample_number)
-        sns.kdeplot(matrix, shade=True, fill=True)
+        sns.kdeplot(matrix, fill=True)
     elif axis == 0:
         col_number = data.shape[1]
         if data.shape[0] * data.shape[1] > sample_number:
@@ -50,10 +50,10 @@ def kde(
 
             for i in tqdm(range(col_number)):
                 _vector_ = down_sampling_data(data.X[:, i], row_number)
-                sns.kdeplot(np.array(_vector_).flatten(), shade=True, fill=True, **kwargs)
+                sns.kdeplot(np.array(_vector_).flatten(), fill=True, **kwargs)
         else:
             for i in tqdm(range(col_number)):
-                sns.kdeplot(np.array(data.X[:, i]).flatten(), shade=True, fill=True, **kwargs)
+                sns.kdeplot(np.array(data.X[:, i]).flatten(), fill=True, **kwargs)
 
         if is_legend:
             ax.legend(list(adata.var.index))
@@ -65,10 +65,10 @@ def kde(
 
             for i in tqdm(range(row_number)):
                 _vector_ = down_sampling_data(data.X[i, :], col_number)
-                sns.kdeplot(np.array(_vector_).flatten(), shade=True, fill=True, **kwargs)
+                sns.kdeplot(np.array(_vector_).flatten(), fill=True, **kwargs)
         else:
             for i in tqdm(range(row_number)):
-                sns.kdeplot(np.array(data.X[i, :]).flatten(), shade=True, fill=True, **kwargs)
+                sns.kdeplot(np.array(data.X[i, :]).flatten(), fill=True, **kwargs)
 
         if is_legend:
             ax.legend(list(adata.obs.index))

sciv/plot/_line_.py

@@ -119,7 +119,7 @@ def stability_line(
         locator = mdates.DayLocator(interval=1)
         chart.xaxis.set_major_locator(locator)
 
-        ax.set_xticklabels(ax.get_xticklabels(), rotation=x_name_rotation)
+        ax.tick_params(axis='x', rotation=x_name_rotation)
     else:
         plt.xticks(x_ticks, rotation=x_name_rotation)
 

sciv/plot/_violin_.py

@@ -48,7 +48,8 @@ def violin_base(
 
     if hue is not None and hue not in df_columns:
         ul.log(__name__).error(
-            f"The `hue` ({hue}) parameter must be in the `df` parameter data column name ({df_columns})")
+            f"The `hue` ({hue}) parameter must be in the `df` parameter data column name ({df_columns})"
+        )
         raise ValueError(f"The `hue` ({hue}) parameter must be in the `df` parameter data column name ({df_columns})")
 
     fig, ax = plot_start(width, height, bottom, output, show)
@@ -67,20 +68,29 @@ def violin_base(
     if is_sort:
         new_df.sort_values([value], ascending=False, inplace=True)
         y_names: Union[list, None] = list(new_df[clusters])
+
         if "color" in df_columns:
             colors = list(new_df["color"])
+
     else:
         new_df.index = new_df[clusters]
+
         if order_names is not None:
             y_names: list = order_names
+
             if "color" in df_columns:
+
                 for i in order_names:
+
                     for j, c in zip(new_df[clusters], new_df["color"]):
+
                         if i == j:
                             colors.append(c)
                             break
+
         else:
             y_names = list(new_df[clusters])
+
             if "color" in df_columns:
                 colors = list(new_df["color"])
 

sciv/tool/_algorithm_.py

@@ -455,11 +455,11 @@ def semi_mutual_knn_weight(
         del data
         np.fill_diagonal(new_data, 0)
 
-    def _knn_k_(_mat: matrix_data, k: int):
+    def _knn_k_(_mat: matrix_data, k: int, info: str = "LOG"):
         n_rows = _mat.shape[0]
         adj = sparse.lil_matrix((n_rows, n_rows), dtype=np.int8)
 
-        ul.log(__name__).info("Calculate the k-nearest neighbors of each node.")
+        ul.log(__name__).info(f"Calculate the k-nearest neighbors of each node. ({info})")
 
         for i in tqdm(range(n_rows)):
             row = np.array(_mat[i]).ravel()
@@ -477,7 +477,7 @@ def semi_mutual_knn_weight(
 
         return adj.tocsr()
 
-    def _knn(_mat: matrix_data, k: int) -> matrix_data:
+    def _knn(_mat: matrix_data, k: int, info: str = "LOG") -> matrix_data:
         """
         Return k-nearest-neighbor 0/1 adjacency matrix (int8 to save memory).
         Supports both sparse and dense inputs.
@@ -486,11 +486,11 @@ def semi_mutual_knn_weight(
         if sparse.issparse(_mat):
             # Sparse path: sort each row's data to find the k-th largest
             _mat = _mat.tocsr(copy=False)
-            return _knn_k_(_mat, k)
+            return _knn_k_(_mat, k, info)
         else:
 
             if is_for:
-                return _knn_k_(_mat, k)
+                return _knn_k_(_mat, k, info)
             else:
                 # Dense path: vectorized thresholding
                 kth_val = np.sort(_mat, axis=1)[:, -(k + 1)]
@@ -499,12 +499,12 @@ def semi_mutual_knn_weight(
                 return adj
 
     # Compute adjacency matrices for AND/OR logic
-    adj_and = _knn(new_data, neighbors)
+    adj_and = _knn(new_data, neighbors, "AND")
 
     if neighbors == or_neighbors:
         adj_or = adj_and
     else:
-        adj_or = _knn(new_data, or_neighbors)
+        adj_or = _knn(new_data, or_neighbors, "OR")
 
     # Symmetrize
     if sparse.issparse(adj_and):
@@ -519,7 +519,12 @@ def semi_mutual_knn_weight(
 
     # Ensure full connectivity if required
     if is_mknn_fully_connected:
-        adj_1nn = _knn(new_data, 1)
+        adj_1nn = _knn(new_data, 1, "ONE")
+
+        if sparse.issparse(adj_1nn):
+            adj_1nn = adj_1nn.maximum(adj_1nn.T)
+        else:
+            adj_1nn = np.maximum(adj_1nn, adj_1nn.T)
 
         if sparse.issparse(adj_and):
             adj_and = adj_and.maximum(adj_1nn)

sciv/tool/_random_walk_.py

@@ -280,7 +280,7 @@ class RandomWalk:
             to `False`, `is_ablation` will only take effect;
         :return: Stable distribution score.
         """
-        ul.log(__name__).info("Random walk.")
+        ul.log(__name__).info("Random walk with weighted seed cells.")
 
         start_time = time.time()
 

{sciv-0.0.99.dist-info → sciv-0.0.101.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: sciv
-Version: 0.0.99
+Version: 0.0.101
 Summary: Unveiling the pivotal cell types involved in variant function regulation at a single-cell resolution
 Project-URL: github, https://github.com/YuZhengM/sciv
 Author-email: Zheng-Min Yu <yuzmbio@163.com>

{sciv-0.0.99.dist-info → sciv-0.0.101.dist-info}/RECORD

@@ -5,20 +5,20 @@ sciv/file/_write_.py,sha256=W3M9CmPi7BuKAffz1fdi-vA5DzAFZ7wmcggp33N9Xtg,7848
 sciv/model/__init__.py,sha256=k8SO9FpJaGn2ANqJyaz3HXMas7jH9toPVtpw703kOqg,149
 sciv/model/_core_.py,sha256=K4X1kGCdFW_5zxwMtACfjKeg_X34eYSPcSl5F2hMX_8,33690
 sciv/plot/__init__.py,sha256=2tRNT6TZNz9r38lnna712RGsH7OJ2QkGa37XKgzejHQ,1865
-sciv/plot/_bar_.py,sha256=xWpFbJTHgQMLuUSXa4uE69RGHXSCGinapxL-1imdDcU,14355
+sciv/plot/_bar_.py,sha256=iQoFsFFbPbFj3JX2BgG45dqG-hjOka90LDw0WfdYNY8,14823
 sciv/plot/_barcode_.py,sha256=RDOedQ8ZtXWFyJ2c772RDfqO4TMIpHMvcMZMAVqky90,5073
-sciv/plot/_box_.py,sha256=485YNmSBj2IB7JlRHaKHMVj0yr0Pz02B3-r57SJqtMk,6009
+sciv/plot/_box_.py,sha256=I_F6imlAEZNSM4ztdtFZnq5JSZpsDgITACKJ45R6M9g,6077
 sciv/plot/_bubble_.py,sha256=vwBs-voOdkBfyMho5w4Rt68b8_Dw1C3NClLdauonC4U,1066
-sciv/plot/_core_.py,sha256=kM3Ym0sZYl0AsfndJ6m9hOlvoMYvii-e8ukcF_USXiA,22222
-sciv/plot/_graph_.py,sha256=h5sI7p4Y1YKplR9e-bfKrvtaWN_LB0SoIndqH-2qZ7I,12613
+sciv/plot/_core_.py,sha256=9MRGhgcA4x1Gky2LFNFNu7ggBsCgVxt0fHceHPuL-Pw,22336
+sciv/plot/_graph_.py,sha256=-AasDYbFnRbWS4eTPDcZ4fzySqrLihUqiY_EBixrQUs,12617
 sciv/plot/_heat_map_.py,sha256=pne4DSydBg0sc4KT88N5A_7QkDnvvmGUwyqZ7-WWtcM,14896
-sciv/plot/_kde_.py,sha256=9VC6DKHbQAzKenz4DoiL5LgUzEJBtAWXXPCLxOV74kY,2459
-sciv/plot/_line_.py,sha256=P944YoGxSocHncjtme4n_iIPK9sghNW3-Xee6Tyy4HE,3850
+sciv/plot/_kde_.py,sha256=A7_NFoSD0ih9h0Hh_1Y8ZOOz06AyANfi3M6sdkz61ZE,2399
+sciv/plot/_line_.py,sha256=qOSf4MJSXCTh8WCcd_dFlTsqMp-MeR2fGqR9lzEVPwo,3834
 sciv/plot/_pie_.py,sha256=OVPv85MdjkjABh-uP5Y-KVk1VygFIzz87QT3OE1hZjU,6399
 sciv/plot/_radar_.py,sha256=g_abzmzibJIR6it59TendUI236inbzYl0IzzdoA3Uuc,6304
 sciv/plot/_scatter_.py,sha256=HPy1eFqDk4vksAiKjn_Qk00fYO4GyqRs1O3ZZh44Ezw,16998
 sciv/plot/_venn_.py,sha256=TfNTuxog2pT7sicKBEEMtleoHXwnenbl3CMWJu9c2vs,2675
-sciv/plot/_violin_.py,sha256=40LYeHFYyNL1XbKHTjmwLT0zzh7o5emqN9Gl2sbl-DA,6458
+sciv/plot/_violin_.py,sha256=IoeAv-ipGkBO0LouI2Y2huPP44u_-6O9CfBkgrPsK6s,6486
 sciv/preprocessing/__init__.py,sha256=56RgDai5I3sZ4hl3aaV80ogOeUscYsU3nJUWE80jZ-k,742
 sciv/preprocessing/_anndata_.py,sha256=3d1cHFs1YA9UZkIPE089nZHFi-DjK-c1fRyi2shfSh4,6302
 sciv/preprocessing/_gencode_.py,sha256=HKhRgK76khGepdv4FaKiOvTys1uJTbvIyrKUta5W0K8,2108
@@ -27,13 +27,13 @@ sciv/preprocessing/_scanpy_.py,sha256=tB8BD2wpLAU8_YxdqrgNtcjpNXNRo-JCdm2lxaKDBL
 sciv/preprocessing/_scvi_.py,sha256=7QxwPA2kR_g15X28aEak7AFA4kyQ-UbtpiLH-rc5Ksg,10780
 sciv/preprocessing/_snapatac_.py,sha256=Dq8CHF7Psl3CQszaEokQYO56Oe2uzyWOy_cGlaOywfc,27798
 sciv/tool/__init__.py,sha256=WXzHkWt6RgBC3qqD-98nR5wQmt6oC850ox_VpMrapSU,2468
-sciv/tool/_algorithm_.py,sha256=mYKfSuYGelLd2secwyqGPxBQYd3x2yDKw1z7HK8mqYE,53773
+sciv/tool/_algorithm_.py,sha256=82ATnZ9-6Cxl6UV9n6jzQawTt7s53crl0hLuZiW-JZ0,54012
 sciv/tool/_matrix_.py,sha256=SnC3sXic_ufuEXStcD_HncvYH6apBdNK6nhG6jFLmjA,24324
-sciv/tool/_random_walk_.py,sha256=JOB97XLxlZYHvlIST1wlXgA0mw6fybkWnJGq6X_kbsk,48871
+sciv/tool/_random_walk_.py,sha256=L7bpYDVNwGSpeB5ipDw81Y6KCJoAjTx1Cs9o4JWZej8,48896
 sciv/util/__init__.py,sha256=nOxZ8if27X7AUJ6hZwTwxOJwIBJb0obWlHjqCzjg_Gc,1964
 sciv/util/_constant_.py,sha256=w0wKQd8guLd1ZTW24_5aECrWsIWDiNQmEpLsWlHar1A,3000
 sciv/util/_core_.py,sha256=TUWfBNRJzWuoQ9ffew_DjnlkNydG-Rmujl_RH4Ln9io,14917
-sciv-0.0.99.dist-info/METADATA,sha256=mObiKCwRCM2mfoao7bZZF-4D_AfqDbClUAaZ0XEc-DI,3465
-sciv-0.0.99.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
-sciv-0.0.99.dist-info/licenses/LICENSE,sha256=4UvHVf3qCOZjHLs4LkYz8u96XRpXnZrpTKrkUQPs5_A,1075
-sciv-0.0.99.dist-info/RECORD,,
+sciv-0.0.101.dist-info/METADATA,sha256=AHOc-TQXHovRT5hgIsJ0tSRLdd6L5-axMiiKYF3_NGU,3466
+sciv-0.0.101.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
+sciv-0.0.101.dist-info/licenses/LICENSE,sha256=4UvHVf3qCOZjHLs4LkYz8u96XRpXnZrpTKrkUQPs5_A,1075
+sciv-0.0.101.dist-info/RECORD,,