diff --git a/pertpy/tools/_mixscape.py b/pertpy/tools/_mixscape.py index 0b384c44..2184714c 100644 --- a/pertpy/tools/_mixscape.py +++ b/pertpy/tools/_mixscape.py @@ -227,7 +227,7 @@ def mixscape( X = adata_comp.layers["X_pert"] except KeyError: raise KeyError( - "No 'X_pert' found in .layers! Please run pert_sign first to calculate perturbation signature!" + "No 'X_pert' found in .layers! Please run perturbation_signature first to calculate perturbation signature!" ) from None # initialize return variables adata.obs[f"{new_class_name}_p_{perturbation_type.lower()}"] = 0 @@ -461,7 +461,13 @@ def _get_perturbation_markers( adata_split = adata[split_mask].copy() # find top DE genes between cells with targeting and non-targeting gRNAs sc.tl.rank_genes_groups( - adata_split, layer=layer, groupby=labels, groups=genes, reference=control, method="t-test" + adata_split, + layer=layer, + groupby=labels, + groups=genes, + reference=control, + method="t-test", + use_raw=False, ) # get DE genes for each gene for gene in genes: