risk-network 0.0.6b10__py3-none-any.whl → 0.0.7b0__py3-none-any.whl

risk/__init__.py

@@ -7,4 +7,4 @@ RISK: RISK Infers Spatial Kinships
 
 from risk.risk import RISK
 
-__version__ = "0.0.6-beta.10"
+__version__ = "0.0.7-beta.0"

risk/annotations/annotations.py

@@ -132,27 +132,32 @@ def define_top_annotations(
             nx.connected_components(enriched_network), key=len, reverse=True
         )
         size_connected_components = np.array([len(c) for c in connected_components])
+
+        # Filter the size of connected components by min_cluster_size and max_cluster_size
+        filtered_size_connected_components = size_connected_components[
+            (size_connected_components >= min_cluster_size)
+            & (size_connected_components <= max_cluster_size)
+        ]
+        # Calculate the number of connected components and large connected components
         num_connected_components = len(connected_components)
-        num_large_connected_components = np.sum(
-            np.logical_and(
-                size_connected_components >= min_cluster_size,
-                size_connected_components <= max_cluster_size,
-            )
-        )
+        num_large_connected_components = len(filtered_size_connected_components)
+
+        # Assign the number of connected components
         annotations_enrichment_matrix.loc[attribute, "num connected components"] = (
             num_connected_components
         )
-        annotations_enrichment_matrix.at[attribute, "size connected components"] = (
-            size_connected_components
-        )
+        # Filter out attributes with more than one connected component
+        annotations_enrichment_matrix.loc[
+            annotations_enrichment_matrix["num connected components"] > 1, "top attributes"
+        ] = False
+        # Assign the number of large connected components
         annotations_enrichment_matrix.loc[attribute, "num large connected components"] = (
             num_large_connected_components
         )
-
-    # Filter out attributes with more than one connected component
-    annotations_enrichment_matrix.loc[
-        annotations_enrichment_matrix["num connected components"] > 1, "top attributes"
-    ] = False
+        # Assign the size of connected components, ensuring it is always a list
+        annotations_enrichment_matrix.at[attribute, "size connected components"] = (
+            filtered_size_connected_components.tolist()
+        )
 
     return annotations_enrichment_matrix
 

risk/neighborhoods/domains.py

@@ -23,7 +23,8 @@ def define_domains(
     linkage_method: str,
     linkage_metric: str,
 ) -> pd.DataFrame:
-    """Define domains and assign nodes to these domains based on their enrichment scores and clustering.
+    """Define domains and assign nodes to these domains based on their enrichment scores and clustering,
+    handling errors by assigning unique domains when clustering fails.
 
     Args:
         top_annotations (pd.DataFrame): DataFrame of top annotations data for the network nodes.
@@ -35,31 +36,29 @@ def define_domains(
     Returns:
         pd.DataFrame: DataFrame with the primary domain for each node.
     """
-    # Check if there's more than one column in significant_neighborhoods_enrichment
-    if significant_neighborhoods_enrichment.shape[1] == 1:
-        print("Single annotation detected. Skipping clustering.")
-        top_annotations["domain"] = 1  # Assign a default domain or handle appropriately
-    else:
-        # Perform hierarchical clustering on the binary enrichment matrix
+    try:
+        # Transpose the matrix to cluster annotations
         m = significant_neighborhoods_enrichment[:, top_annotations["top attributes"]].T
         best_linkage, best_metric, best_threshold = _optimize_silhouette_across_linkage_and_metrics(
             m, linkage_criterion, linkage_method, linkage_metric
         )
-        try:
-            Z = linkage(m, method=best_linkage, metric=best_metric)
-        except ValueError as e:
-            raise ValueError("No significant annotations found.") from e
-
+        # Perform hierarchical clustering
+        Z = linkage(m, method=best_linkage, metric=best_metric)
         print(
             f"Linkage criterion: '{linkage_criterion}'\nLinkage method: '{best_linkage}'\nLinkage metric: '{best_metric}'"
         )
         print(f"Optimal linkage threshold: {round(best_threshold, 3)}")
-
+        # Calculate the optimal threshold for clustering
         max_d_optimal = np.max(Z[:, 2]) * best_threshold
-        domains = fcluster(Z, max_d_optimal, criterion=linkage_criterion)
         # Assign domains to the annotations matrix
+        domains = fcluster(Z, max_d_optimal, criterion=linkage_criterion)
         top_annotations["domain"] = 0
         top_annotations.loc[top_annotations["top attributes"], "domain"] = domains
+    except ValueError:
+        # If a ValueError is encountered, handle it by assigning unique domains
+        n_rows = len(top_annotations)
+        print(f"Error encountered. Skipping clustering and assigning {n_rows} unique domains.")
+        top_annotations["domain"] = range(1, n_rows + 1)  # Assign unique domains
 
     # Create DataFrames to store domain information
     node_to_enrichment = pd.DataFrame(

risk/network/plot.py

@@ -1123,7 +1123,9 @@ def _to_rgba(
     """
     # Handle single color case (string, RGB, or RGBA)
     if isinstance(color, str) or (
-        isinstance(color, (list, tuple, np.ndarray)) and len(color) in [3, 4]
+        isinstance(color, (list, tuple, np.ndarray))
+        and len(color) in [3, 4]
+        and not any(isinstance(c, (list, tuple, np.ndarray)) for c in color)
     ):
         rgba_color = np.array(mcolors.to_rgba(color))
         # Only set alpha if the input is an RGB color or a string (not RGBA)

{risk_network-0.0.6b10.dist-info → risk_network-0.0.7b0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: risk-network
-Version: 0.0.6b10
+Version: 0.0.7b0
 Summary: A Python package for biological network analysis
 Author: Ira Horecka
 Author-email: Ira Horecka <ira89@icloud.com>

{risk_network-0.0.6b10.dist-info → risk_network-0.0.7b0.dist-info}/RECORD

@@ -1,21 +1,21 @@
-risk/__init__.py,sha256=nUkWz8VqnztwzCLGVbN6G8bFixZFkDyLsEOsb5mLAGc,113
+risk/__init__.py,sha256=Qyktssx5ZswUpqcPtSMq9Zn-zzJXl2fka6MqbHS-JxQ,112
 risk/constants.py,sha256=XInRaH78Slnw_sWgAsBFbUHkyA0h0jL0DKGuQNbOvjM,550
 risk/risk.py,sha256=CKDIzVo9Jvl-fgzIlk5ZtJL9pIBMma24WK6EYdVu5po,20648
 risk/annotations/__init__.py,sha256=vUpVvMRE5if01Ic8QY6M2Ae3EFGJHdugEe9PdEkAW4Y,138
-risk/annotations/annotations.py,sha256=DRUTdGzMdqo62NWSapBUksbvPr9CrzD76qtOcxeNKmo,10554
+risk/annotations/annotations.py,sha256=K7cUA6vYTKYAvj0xHqrAwNEYtmPq4H7LDYENAOVQdQ0,11014
 risk/annotations/io.py,sha256=lo7NKqOVkeeBp58JBxWJHtA0xjL5Yoxqe9Ox0daKlZk,9457
 risk/log/__init__.py,sha256=xuLImfxFlKpnVhzi_gDYlr2_c9cLkrw2c_3iEsXb1as,107
 risk/log/console.py,sha256=im9DRExwf6wHlcn9fewoDcKIpo3vPcorZIaNAl-0csY,355
 risk/log/params.py,sha256=Rfdg5UcGCrG80m6V79FyORERWUqIzHFO7tGiY4zAImM,6347
 risk/neighborhoods/__init__.py,sha256=tKKEg4lsbqFukpgYlUGxU_v_9FOqK7V0uvM9T2QzoL0,206
 risk/neighborhoods/community.py,sha256=7ebo1Q5KokSQISnxZIh2SQxsKXdXm8aVkp-h_DiQ3K0,6818
-risk/neighborhoods/domains.py,sha256=5V--Nj-TrSdubhD_2PI57ffcn_PMSEgpX_iY5OjT6R8,10626
+risk/neighborhoods/domains.py,sha256=bxJUxqFTynzX0mf3E8-AA4_Rfccje1reeVVhfzb1-pE,10672
 risk/neighborhoods/neighborhoods.py,sha256=sHmjFFl2U5qV9YbQCRbpbI36j7dS7IFfFwwRb1_-AuM,13945
 risk/network/__init__.py,sha256=iEPeJdZfqp0toxtbElryB8jbz9_t_k4QQ3iDvKE8C_0,126
 risk/network/geometry.py,sha256=H1yGVVqgbfpzBzJwEheDLfvGLSA284jGQQTn612L4Vc,6759
 risk/network/graph.py,sha256=7haHu4M3fleqbrIzs6HC9jnKizSERzmmAYSmUwdoSXA,13953
 risk/network/io.py,sha256=gG50kOknO-D3HkW1HsbHMkTMvjUtn3l4W4Jwd-rXNr8,21202
-risk/network/plot.py,sha256=_g5xHolMTAfZCBvYYEX1CYME4s4zA2hTHtN-utaMPik,61978
+risk/network/plot.py,sha256=F6KPjmBYWrThKZScHs9SuzoKQiytBvzrmGhGberHjwo,62063
 risk/stats/__init__.py,sha256=e-BE_Dr_jgiK6hKM-T-tlG4yvHnId8e5qjnM0pdwNVc,230
 risk/stats/fisher_exact.py,sha256=-bPwzu76-ob0HzrTV20mXUTot7v-MLuqFaAoab-QxPg,4966
 risk/stats/hypergeom.py,sha256=lrIFdhCWRjvM4apYw1MlOKqT_IY5OjtCwrjdtJdt6Tg,4954
@@ -23,8 +23,8 @@ risk/stats/stats.py,sha256=kvShov-94W6ffgDUTb522vB9hDJQSyTsYif_UIaFfSM,7059
 risk/stats/permutation/__init__.py,sha256=neJp7FENC-zg_CGOXqv-iIvz1r5XUKI9Ruxhmq7kDOI,105
 risk/stats/permutation/permutation.py,sha256=qLWdwxEY6nmkYPxpM8HLDcd2mbqYv9Qr7CKtJvhLqIM,9220
 risk/stats/permutation/test_functions.py,sha256=HuDIM-V1jkkfE1rlaIqrWWBSKZt3dQ1f-YEDjWpnLSE,2343
-risk_network-0.0.6b10.dist-info/LICENSE,sha256=jOtLnuWt7d5Hsx6XXB2QxzrSe2sWWh3NgMfFRetluQM,35147
-risk_network-0.0.6b10.dist-info/METADATA,sha256=aOz9JrsPIpByvzMCwsbBBcBHHMLu3JzLF9FZMp9-IuM,43143
-risk_network-0.0.6b10.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
-risk_network-0.0.6b10.dist-info/top_level.txt,sha256=NX7C2PFKTvC1JhVKv14DFlFAIFnKc6Lpsu1ZfxvQwVw,5
-risk_network-0.0.6b10.dist-info/RECORD,,
+risk_network-0.0.7b0.dist-info/LICENSE,sha256=jOtLnuWt7d5Hsx6XXB2QxzrSe2sWWh3NgMfFRetluQM,35147
+risk_network-0.0.7b0.dist-info/METADATA,sha256=Yokjvu7qlqWV6F_qJQ9O6TIwKw_9XpD_2qgwQHyimRY,43142
+risk_network-0.0.7b0.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
+risk_network-0.0.7b0.dist-info/top_level.txt,sha256=NX7C2PFKTvC1JhVKv14DFlFAIFnKc6Lpsu1ZfxvQwVw,5
+risk_network-0.0.7b0.dist-info/RECORD,,