sciv 0.0.99__tar.gz → 0.0.101__tar.gz

{sciv-0.0.99 → sciv-0.0.101}/PKG-INFO

@@ -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 → sciv-0.0.101}/pyproject.toml

@@ -3,7 +3,7 @@ requires = ["hatchling"]
 build-backend = "hatchling.build"
 [project]
 name = "sciv"
-version = "0.0.99"
+version = "0.0.101"
 authors = [
     { name = "Zheng-Min Yu", email = "yuzmbio@163.com" },
 ]

{sciv-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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-0.0.99 → sciv-0.0.101}/src/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()