workbench 0.8.162__py3-none-any.whl → 0.8.220__py3-none-any.whl

workbench/algorithms/dataframe/fingerprint_proximity.py

@@ -1,162 +1,498 @@
 import pandas as pd
 import numpy as np
 from sklearn.neighbors import NearestNeighbors
-from typing import Union, List
+from typing import Union, List, Optional
 import logging
 
 # Workbench Imports
-from workbench.algorithms.dataframe.proximity import Proximity, ProximityType
+from workbench.algorithms.dataframe.proximity import Proximity
+from workbench.algorithms.dataframe.projection_2d import Projection2D
+from workbench.utils.chem_utils.fingerprints import compute_morgan_fingerprints
 
 # Set up logging
 log = logging.getLogger("workbench")
 
 
 class FingerprintProximity(Proximity):
+    """Proximity computations for binary fingerprints using Tanimoto similarity.
+
+    Note: Tanimoto similarity is equivalent to Jaccard similarity for binary vectors.
+    Tanimoto(A, B) = |A ∩ B| / |A ∪ B|
+    """
+
     def __init__(
-        self, df: pd.DataFrame, id_column: Union[int, str], fingerprint_column: str, n_neighbors: int = 5
+        self,
+        df: pd.DataFrame,
+        id_column: str,
+        fingerprint_column: Optional[str] = None,
+        target: Optional[str] = None,
+        include_all_columns: bool = False,
+        radius: int = 2,
+        n_bits: int = 1024,
     ) -> None:
         """
         Initialize the FingerprintProximity class for binary fingerprint similarity.
 
         Args:
-            df (pd.DataFrame): DataFrame containing fingerprints.
-            id_column (Union[int, str]): Name of the column used as an identifier.
-            fingerprint_column (str): Name of the column containing fingerprints.
-            n_neighbors (int): Default number of neighbors to compute.
+            df: DataFrame containing fingerprints or SMILES.
+            id_column: Name of the column used as an identifier.
+            fingerprint_column: Name of the column containing fingerprints (bit strings).
+                If None, looks for existing "fingerprint" column or computes from SMILES.
+            target: Name of the target column. Defaults to None.
+            include_all_columns: Include all DataFrame columns in neighbor results. Defaults to False.
+            radius: Radius for Morgan fingerprint computation (default: 2).
+            n_bits: Number of bits for fingerprint (default: 1024).
+        """
+        # Store fingerprint computation parameters
+        self._fp_radius = radius
+        self._fp_n_bits = n_bits
+
+        # Store the requested fingerprint column (may be None)
+        self._fingerprint_column_arg = fingerprint_column
+
+        # Determine fingerprint column name (but don't compute yet - that happens in _prepare_data)
+        self.fingerprint_column = self._resolve_fingerprint_column_name(df, fingerprint_column)
+
+        # Call parent constructor with fingerprint_column as the only "feature"
+        super().__init__(
+            df,
+            id_column=id_column,
+            features=[self.fingerprint_column],
+            target=target,
+            include_all_columns=include_all_columns,
+        )
+
+    @staticmethod
+    def _resolve_fingerprint_column_name(df: pd.DataFrame, fingerprint_column: Optional[str]) -> str:
+        """
+        Determine the fingerprint column name, validating it exists or can be computed.
+
+        Args:
+            df: Input DataFrame.
+            fingerprint_column: Explicitly specified fingerprint column, or None.
+
+        Returns:
+            Name of the fingerprint column to use.
+
+        Raises:
+            ValueError: If no fingerprint column exists and no SMILES column found.
         """
-        self.fingerprint_column = fingerprint_column
+        # If explicitly provided, validate it exists
+        if fingerprint_column is not None:
+            if fingerprint_column not in df.columns:
+                raise ValueError(f"Fingerprint column '{fingerprint_column}' not found in DataFrame")
+            return fingerprint_column
+
+        # Check for existing "fingerprint" column
+        if "fingerprint" in df.columns:
+            log.info("Using existing 'fingerprint' column")
+            return "fingerprint"
+
+        # Will need to compute from SMILES - validate SMILES column exists
+        smiles_column = next((col for col in df.columns if col.lower() == "smiles"), None)
+        if smiles_column is None:
+            raise ValueError(
+                "No fingerprint column provided and no SMILES column found. "
+                "Either provide a fingerprint_column or include a 'smiles' column in the DataFrame."
+            )
 
-        # Call the parent class constructor
-        super().__init__(df, id_column=id_column, features=[fingerprint_column], n_neighbors=n_neighbors)
+        # Fingerprints will be computed in _prepare_data
+        return "fingerprint"
 
-    # Override the build_proximity_model method
-    def build_proximity_model(self) -> None:
+    def _prepare_data(self) -> None:
+        """Compute fingerprints from SMILES if needed."""
+        # If fingerprint column doesn't exist yet, compute it
+        if self.fingerprint_column not in self.df.columns:
+            log.info(f"Computing Morgan fingerprints (radius={self._fp_radius}, n_bits={self._fp_n_bits})...")
+            self.df = compute_morgan_fingerprints(self.df, radius=self._fp_radius, n_bits=self._fp_n_bits)
+
+    def _build_model(self) -> None:
         """
-        Prepare the fingerprint data for nearest neighbor calculations.
-        Converts fingerprint strings to binary arrays and initializes NearestNeighbors.
+        Build the fingerprint proximity model for Tanimoto similarity.
+
+        For binary fingerprints: uses Jaccard distance (1 - Tanimoto)
+        For count fingerprints: uses weighted Tanimoto (Ruzicka) distance
         """
-        log.info("Converting fingerprints to binary feature matrix...")
-        self.proximity_type = ProximityType.SIMILARITY
+        # Convert fingerprint strings to matrix and detect format
+        self.X, self._is_count_fp = self._fingerprints_to_matrix(self.df)
 
-        # Convert fingerprint strings to binary arrays
+        if self._is_count_fp:
+            # Weighted Tanimoto (Ruzicka) for count vectors: 1 - Σmin(A,B)/Σmax(A,B)
+            log.info("Building NearestNeighbors model (weighted Tanimoto for count fingerprints)...")
 
-        fingerprint_bits = self.df[self.fingerprint_column].apply(
-            lambda fp: np.array([int(bit) for bit in fp], dtype=np.bool_)
-        )
-        self.X = np.vstack(fingerprint_bits)
+            def ruzicka_distance(a, b):
+                """Ruzicka distance = 1 - weighted Tanimoto similarity."""
+                min_sum = np.minimum(a, b).sum()
+                max_sum = np.maximum(a, b).sum()
+                if max_sum == 0:
+                    return 0.0
+                return 1.0 - (min_sum / max_sum)
 
-        # Use Jaccard similarity for binary fingerprints
-        log.info("Computing NearestNeighbors with Jaccard metric...")
-        self.nn = NearestNeighbors(metric="jaccard", n_neighbors=self.n_neighbors + 1).fit(self.X)
+            self.nn = NearestNeighbors(metric=ruzicka_distance, algorithm="ball_tree").fit(self.X)
+        else:
+            # Standard Jaccard for binary fingerprints
+            log.info("Building NearestNeighbors model (Jaccard/Tanimoto for binary fingerprints)...")
+            self.nn = NearestNeighbors(metric="jaccard", algorithm="ball_tree").fit(self.X)
 
-    # Override the prep_features_for_query method
-    def prep_features_for_query(self, query_df: pd.DataFrame) -> np.ndarray:
+    def _transform_features(self, df: pd.DataFrame) -> np.ndarray:
         """
-        Prepare the query DataFrame by converting fingerprints to binary arrays.
+        Transform fingerprints to matrix for querying.
 
         Args:
-            query_df (pd.DataFrame): DataFrame containing query fingerprints.
+            df: DataFrame containing fingerprints to transform.
 
         Returns:
-            np.ndarray: Binary feature matrix for the query fingerprints.
+            Feature matrix for the fingerprints (binary or count based on self._is_count_fp).
         """
-        fingerprint_bits = query_df[self.fingerprint_column].apply(
-            lambda fp: np.array([int(bit) for bit in fp], dtype=np.bool_)
-        )
-        return np.vstack(fingerprint_bits)
+        matrix, _ = self._fingerprints_to_matrix(df)
+        return matrix
 
-    def all_neighbors(
-        self,
-        min_similarity: float = None,
-        include_self: bool = False,
-        add_columns: List[str] = None,
-    ) -> pd.DataFrame:
+    def _fingerprints_to_matrix(self, df: pd.DataFrame) -> tuple[np.ndarray, bool]:
         """
-        Find neighbors for all fingerprints in the dataset.
+        Convert fingerprint strings to a numpy matrix.
+
+        Supports two formats (auto-detected):
+            - Bitstrings: "10110010..." → binary matrix (bool), is_count=False
+            - Count vectors: "0,3,0,1,5,..." → count matrix (uint8), is_count=True
 
         Args:
-            min_similarity: Minimum similarity threshold (0-1)
-            include_self: Whether to include self in results
-            add_columns: Additional columns to include in results
+            df: DataFrame containing fingerprint column.
 
         Returns:
-            DataFrame containing neighbors and similarities
+            Tuple of (2D numpy array, is_count_fingerprint boolean)
         """
+        # Auto-detect format based on first fingerprint
+        sample = str(df[self.fingerprint_column].iloc[0])
+        if "," in sample:
+            # Count vector format: preserve counts for weighted Tanimoto
+            fingerprint_values = df[self.fingerprint_column].apply(
+                lambda fp: np.array([int(x) for x in fp.split(",")], dtype=np.uint8)
+            )
+            return np.vstack(fingerprint_values), True
+        else:
+            # Bitstring format: binary values
+            fingerprint_bits = df[self.fingerprint_column].apply(
+                lambda fp: np.array([int(bit) for bit in fp], dtype=np.bool_)
+            )
+            return np.vstack(fingerprint_bits), False
 
-        # Call the parent class method to find neighbors
-        return self.neighbors(
-            query_df=self.df,
-            min_similarity=min_similarity,
-            include_self=include_self,
-            add_columns=add_columns,
+    def _precompute_metrics(self) -> None:
+        """Precompute metrics, adding Tanimoto similarity alongside distance."""
+        # Call parent to compute nn_distance (Jaccard), nn_id, nn_target, nn_target_diff
+        super()._precompute_metrics()
+
+        # Add Tanimoto similarity (keep nn_distance for internal use by target_gradients)
+        self.df["nn_similarity"] = 1 - self.df["nn_distance"]
+
+    def _set_core_columns(self) -> None:
+        """Set core columns using nn_similarity instead of nn_distance."""
+        self.core_columns = [self.id_column, "nn_similarity", "nn_id"]
+        if self.target:
+            self.core_columns.extend([self.target, "nn_target", "nn_target_diff"])
+
+    def _project_2d(self) -> None:
+        """Project the fingerprint matrix to 2D for visualization using UMAP."""
+        if self._is_count_fp:
+            # For count fingerprints, convert to binary for UMAP projection (Jaccard needs binary)
+            X_binary = (self.X > 0).astype(np.bool_)
+            self.df = Projection2D().fit_transform(self.df, feature_matrix=X_binary, metric="jaccard")
+        else:
+            self.df = Projection2D().fit_transform(self.df, feature_matrix=self.X, metric="jaccard")
+
+    def isolated(self, top_percent: float = 1.0) -> pd.DataFrame:
+        """
+        Find isolated data points based on Tanimoto similarity to nearest neighbor.
+
+        Args:
+            top_percent: Percentage of most isolated data points to return (e.g., 1.0 returns top 1%)
+
+        Returns:
+            DataFrame of observations with lowest Tanimoto similarity, sorted ascending
+        """
+        # For Tanimoto similarity, isolated means LOW similarity to nearest neighbor
+        percentile = top_percent
+        threshold = np.percentile(self.df["nn_similarity"], percentile)
+        isolated = self.df[self.df["nn_similarity"] <= threshold].copy()
+        isolated = isolated.sort_values("nn_similarity", ascending=True).reset_index(drop=True)
+        return isolated if self.include_all_columns else isolated[self.core_columns]
+
+    def proximity_stats(self) -> pd.DataFrame:
+        """
+        Return distribution statistics for nearest neighbor Tanimoto similarity.
+
+        Returns:
+            DataFrame with similarity distribution statistics (count, mean, std, percentiles)
+        """
+        return (
+            self.df["nn_similarity"]
+            .describe(percentiles=[0.01, 0.05, 0.1, 0.25, 0.5, 0.75, 0.9, 0.95, 0.99])
+            .to_frame()
         )
 
     def neighbors(
         self,
-        query_df: pd.DataFrame,
-        min_similarity: float = None,
-        include_self: bool = False,
-        add_columns: List[str] = None,
+        id_or_ids: Union[str, int, List[Union[str, int]]],
+        n_neighbors: Optional[int] = 5,
+        min_similarity: Optional[float] = None,
+        include_self: bool = True,
     ) -> pd.DataFrame:
         """
-        Find neighbors for each row in the query DataFrame.
+        Return neighbors for ID(s) from the existing dataset.
 
         Args:
-            query_df: DataFrame containing query fingerprints
-            min_similarity: Minimum similarity threshold (0-1)
-            include_self: Whether to include self in results (if present)
-            add_columns: Additional columns to include in results
+            id_or_ids: Single ID or list of IDs to look up
+            n_neighbors: Number of neighbors to return (default: 5, ignored if min_similarity is set)
+            min_similarity: If provided, find all neighbors with Tanimoto similarity >= this value (0-1)
+            include_self: Whether to include self in results (default: True)
 
         Returns:
-            DataFrame containing neighbors and similarities
-
-        Note: The query DataFrame must include the feature columns. The id_column is optional.
+            DataFrame containing neighbors with Tanimoto similarity scores
         """
-
-        # Calculate radius from similarity if provided
+        # Convert min_similarity to radius (Jaccard distance = 1 - Tanimoto similarity)
         radius = 1 - min_similarity if min_similarity is not None else None
 
-        # Call the parent class method to find neighbors
+        # Call parent method (returns Jaccard distance)
         neighbors_df = super().neighbors(
-            query_df=query_df,
+            id_or_ids=id_or_ids,
+            n_neighbors=n_neighbors,
             radius=radius,
             include_self=include_self,
-            add_columns=add_columns,
         )
 
-        # Convert distances to similarity
+        # Convert Jaccard distance to Tanimoto similarity
         neighbors_df["similarity"] = 1 - neighbors_df["distance"]
         neighbors_df.drop(columns=["distance"], inplace=True)
+
         return neighbors_df
 
+    def neighbors_from_smiles(
+        self,
+        smiles: Union[str, List[str]],
+        n_neighbors: int = 5,
+        min_similarity: Optional[float] = None,
+    ) -> pd.DataFrame:
+        """
+        Find neighbors for SMILES strings not in the reference dataset.
+
+        Args:
+            smiles: Single SMILES string or list of SMILES to query
+            n_neighbors: Number of neighbors to return (default: 5, ignored if min_similarity is set)
+            min_similarity: If provided, find all neighbors with Tanimoto similarity >= this value (0-1)
+
+        Returns:
+            DataFrame containing neighbors with Tanimoto similarity scores.
+            The 'query_id' column contains the SMILES string (or index if list).
+        """
+        # Normalize to list
+        smiles_list = [smiles] if isinstance(smiles, str) else smiles
+
+        # Build a temporary DataFrame with the query SMILES
+        query_df = pd.DataFrame({"smiles": smiles_list})
+
+        # Compute fingerprints using same parameters as the reference dataset
+        query_df = compute_morgan_fingerprints(query_df, radius=self._fp_radius, n_bits=self._fp_n_bits)
+
+        # Transform to matrix (use same format detection as reference)
+        X_query, _ = self._fingerprints_to_matrix(query_df)
+
+        # Query the model
+        if min_similarity is not None:
+            radius = 1 - min_similarity
+            distances, indices = self.nn.radius_neighbors(X_query, radius=radius)
+        else:
+            distances, indices = self.nn.kneighbors(X_query, n_neighbors=n_neighbors)
+
+        # Build results
+        results = []
+        for i, (dists, nbrs) in enumerate(zip(distances, indices)):
+            query_id = smiles_list[i]
+
+            for neighbor_idx, dist in zip(nbrs, dists):
+                neighbor_row = self.df.iloc[neighbor_idx]
+                neighbor_id = neighbor_row[self.id_column]
+                similarity = 1.0 - dist if dist > 1e-6 else 1.0
+
+                result = {
+                    "query_id": query_id,
+                    "neighbor_id": neighbor_id,
+                    "similarity": similarity,
+                }
+
+                # Add target if present
+                if self.target and self.target in self.df.columns:
+                    result[self.target] = neighbor_row[self.target]
+
+                # Include all columns if requested
+                if self.include_all_columns:
+                    for col in self.df.columns:
+                        if col not in [self.id_column, "query_id", "neighbor_id", "similarity"]:
+                            result[f"neighbor_{col}"] = neighbor_row[col]
+
+                results.append(result)
+
+        df_results = pd.DataFrame(results)
+
+        # Sort by query_id then similarity descending
+        if len(df_results) > 0:
+            df_results = df_results.sort_values(["query_id", "similarity"], ascending=[True, False]).reset_index(
+                drop=True
+            )
+
+        return df_results
+
 
 # Testing the FingerprintProximity class
 if __name__ == "__main__":
     pd.set_option("display.max_columns", None)
     pd.set_option("display.width", 1000)
 
-    # Example DataFrame
+    # Create an Example DataFrame with fingerprints
     data = {
-        "id": ["a", "b", "c", "d"],
-        "fingerprint": ["101010", "111010", "101110", "011100"],
+        "id": ["a", "b", "c", "d", "e"],
+        "fingerprint": ["101010", "111010", "101110", "011100", "000111"],
+        "Feature1": [0.1, 0.2, 0.3, 0.4, 0.5],
+        "Feature2": [0.5, 0.4, 0.3, 0.2, 0.1],
+        "target": [1, 0, 1, 0, 5],
     }
     df = pd.DataFrame(data)
 
-    # Initialize the FingerprintProximity class
-    proximity = FingerprintProximity(df, fingerprint_column="fingerprint", id_column="id", n_neighbors=3)
+    # Test basic FingerprintProximity with explicit fingerprint column
+    prox = FingerprintProximity(df, fingerprint_column="fingerprint", id_column="id", target="target")
+    print(prox.neighbors("a", n_neighbors=3))
+
+    # Test neighbors with similarity threshold
+    print(prox.neighbors("a", min_similarity=0.5))
+
+    # Test with include_all_columns=True
+    prox = FingerprintProximity(
+        df,
+        fingerprint_column="fingerprint",
+        id_column="id",
+        target="target",
+        include_all_columns=True,
+    )
+    print(prox.neighbors(["a", "b"]))
+
+    # Regression test: include_all_columns should not break neighbor sorting
+    print("\n" + "=" * 80)
+    print("Regression test: include_all_columns neighbor sorting...")
+    print("=" * 80)
+    neighbors_all_cols = prox.neighbors("a", n_neighbors=4)
+    # Verify neighbors are sorted by similarity (descending), not alphabetically by neighbor_id
+    similarities = neighbors_all_cols["similarity"].tolist()
+    assert similarities == sorted(
+        similarities, reverse=True
+    ), f"Neighbors not sorted by similarity! Got: {similarities}"
+    # Verify query_id column has correct value (the query, not the neighbor)
+    assert all(
+        neighbors_all_cols["id"] == "a"
+    ), f"Query ID column corrupted! Expected all 'a', got: {neighbors_all_cols['id'].tolist()}"
+    print("PASSED: Neighbors correctly sorted by similarity with include_all_columns=True")
+
+    # Test neighbors_from_smiles with synthetic data
+    print("\n" + "=" * 80)
+    print("Testing neighbors_from_smiles...")
+    print("=" * 80)
+
+    # Create reference dataset with known SMILES
+    ref_data = {
+        "id": ["aspirin", "ibuprofen", "naproxen", "caffeine", "ethanol"],
+        "smiles": [
+            "CC(=O)OC1=CC=CC=C1C(=O)O",  # aspirin
+            "CC(C)CC1=CC=C(C=C1)C(C)C(=O)O",  # ibuprofen
+            "COC1=CC2=CC(C(C)C(O)=O)=CC=C2C=C1",  # naproxen
+            "CN1C=NC2=C1C(=O)N(C(=O)N2C)C",  # caffeine
+            "CCO",  # ethanol
+        ],
+        "activity": [1.0, 2.0, 2.5, 3.0, 0.5],
+    }
+    ref_df = pd.DataFrame(ref_data)
+
+    prox_ref = FingerprintProximity(ref_df, id_column="id", target="activity", radius=2, n_bits=1024)
+
+    # Query with a single SMILES (acetaminophen - similar to aspirin)
+    query_smiles = "CC(=O)NC1=CC=C(C=C1)O"  # acetaminophen
+    print(f"\nQuery: acetaminophen ({query_smiles})")
+    neighbors = prox_ref.neighbors_from_smiles(query_smiles, n_neighbors=3)
+    print(neighbors)
+
+    # Query with multiple SMILES
+    print("\nQuery: multiple SMILES (theophylline, methanol)")
+    multi_query = [
+        "CN1C=NC2=C1C(=O)NC(=O)N2",  # theophylline - similar to caffeine
+        "CO",  # methanol - similar to ethanol
+    ]
+    neighbors_multi = prox_ref.neighbors_from_smiles(multi_query, n_neighbors=2)
+    print(neighbors_multi)
+
+    # Test with min_similarity threshold
+    print("\nQuery with min_similarity=0.3:")
+    neighbors_thresh = prox_ref.neighbors_from_smiles(query_smiles, min_similarity=0.3)
+    print(neighbors_thresh)
+
+    print("PASSED: neighbors_from_smiles working correctly")
+
+    # Test on real data from Workbench
+    from workbench.api import FeatureSet, Model
+
+    fs = FeatureSet("aqsol_features")
+    model = Model("aqsol-regression")
+    df = fs.pull_dataframe()[:1000]  # Limit to 1000 for testing
+    prox = FingerprintProximity(df, id_column=fs.id_column, target=model.target())
+
+    print("\n" + "=" * 80)
+    print("Testing Neighbors...")
+    print("=" * 80)
+    test_id = df[fs.id_column].tolist()[0]
+    print(f"\nNeighbors for ID {test_id}:")
+    print(prox.neighbors(test_id))
+
+    print("\n" + "=" * 80)
+    print("Testing isolated compounds...")
+    print("=" * 80)
+
+    # Test isolated data in the top 1%
+    isolated_1pct = prox.isolated(top_percent=1.0)
+    print(f"\nTop 1% most isolated compounds (n={len(isolated_1pct)}):")
+    print(isolated_1pct)
+
+    # Test isolated data in the top 5%
+    isolated_5pct = prox.isolated(top_percent=5.0)
+    print(f"\nTop 5% most isolated compounds (n={len(isolated_5pct)}):")
+    print(isolated_5pct)
+
+    print("\n" + "=" * 80)
+    print("Testing target_gradients...")
+    print("=" * 80)
+
+    # Test with different parameters
+    gradients_1pct = prox.target_gradients(top_percent=1.0, min_delta=1.0)
+    print(f"\nTop 1% target gradients (min_delta=1.0) (n={len(gradients_1pct)}):")
+    print(gradients_1pct)
+
+    gradients_5pct = prox.target_gradients(top_percent=5.0, min_delta=5.0)
+    print(f"\nTop 5% target gradients (min_delta=5.0) (n={len(gradients_5pct)}):")
+    print(gradients_5pct)
+
+    # Test proximity_stats
+    print("\n" + "=" * 80)
+    print("Testing proximity_stats...")
+    print("=" * 80)
+    stats = prox.proximity_stats()
+    print(stats)
 
-    # Test 1: All neighbors
-    print("\n--- Test 1: All Neighbors ---")
-    all_neighbors_df = proximity.all_neighbors()
-    print(all_neighbors_df)
+    # Plot the similarity distribution using pandas
+    print("\n" + "=" * 80)
+    print("Plotting similarity distribution...")
+    print("=" * 80)
+    prox.df["nn_similarity"].hist(bins=50, figsize=(10, 6), edgecolor="black")
 
-    # Test 2: Neighbors for a specific query
-    print("\n--- Test 2: Neighbors for Query ---")
-    query_df = pd.DataFrame({"id": ["a"], "fingerprint": ["101010"]})
-    query_neighbors_df = proximity.neighbors(query_df=query_df)
-    print(query_neighbors_df)
+    # Visualize the 2D projection
+    print("\n" + "=" * 80)
+    print("Visualizing 2D Projection...")
+    print("=" * 80)
+    from workbench.web_interface.components.plugin_unit_test import PluginUnitTest
+    from workbench.web_interface.components.plugins.scatter_plot import ScatterPlot
 
-    # Test 3: Neighbors with similarity threshold
-    print("\n--- Test 3: Neighbors with Minimum Similarity 0.5 ---")
-    query_neighbors_sim_df = proximity.neighbors(query_df=query_df, min_similarity=0.5)
-    print(query_neighbors_sim_df)
+    unit_test = PluginUnitTest(ScatterPlot, input_data=prox.df[:1000], x="x", y="y", color=model.target())
+    unit_test.run()