scvi.external.Decipher#
- class scvi.external.Decipher(adata, **kwargs)[source]#
Decipher model for single-cell data analysis [Nazaret et al., 2024].
- Parameters:
adata (
AnnData) – AnnData object that has been registered viasetup_anndata().dim_v – Dimension of the interpretable latent space v.
dim_z – Dimension of the intermediate latent space z.
layers_v_to_z – Hidden layer sizes for the v to z decoder network.
layers_z_to_x – Hidden layer sizes for the z to x decoder network.
beta – Regularization parameter for the KL divergence.
Attributes table#
Data attached to model instance. |
|
Manager instance associated with self.adata. |
|
The current device that the module's params are on. |
|
What the get normalized functions name is |
|
Returns computed metrics during training. |
|
Whether the model has been trained. |
|
Data attached to model instance. |
|
Returns the run id of the model. |
|
Returns the run name of the model. |
|
Summary string of the model. |
|
Observations that are in test set. |
|
Observations that are in train set. |
|
Observations that are in validation set. |
Methods table#
|
Compute the decipher time for each cell, based on the inferred trajectories. |
|
Compute the gene patterns for a trajectory. |
|
Impute gene expression from the decipher model. |
|
Converts a legacy saved model (<v0.15.0) to the updated save format. |
|
Returns the object in AnnData associated with the key in the data registry. |
|
Deregisters the |
|
Not implemented for this model class. |
|
Retrieves the |
|
Returns the object in AnnData associated with the key in the data registry. |
|
Get the latent representation of the data. |
|
Not implemented for this model class. |
|
Returns the string provided to setup of a specific setup_arg. |
|
Returns the state registry for the AnnDataField registered with this instance. |
|
Variable names of input data. |
|
Instantiate a model from the saved output. |
|
Return the full registry saved with the model. |
|
Registers an |
|
Save the state of the model. |
|
Set up a model from on-disk AnnData files via the annbatch streaming loader. |
|
Sets up the |
|
Move the model to the device. |
|
Train the model. |
|
Transfer fields from a model to an AnnData object. |
|
Update setup method args. |
|
Print summary of the setup for the initial AnnData or a given AnnData object. |
|
Prints summary of the registry. |
|
Print args used to setup a saved model. |
Prints setup kwargs used to produce a given registry. |
Attributes#
Methods#
- static Decipher.compute_decipher_time(adata, cluster_obs_key, trajectory, n_neighbors=10)[source]#
Compute the decipher time for each cell, based on the inferred trajectories.
The decipher time is computed by KNN regression of the cells’ decipher v on the trajectories.
- Parameters:
adata (AnnData) – The annotated data matrix.
cluster_obs_key (str) – The key in adata.obs containing cluster assignments.
trajectory (Trajectory) – A
Trajectoryobject containing the trajectory information.n_neighbors (int) – The number of neighbors to use for the KNN regression.
- Return type:
- Returns:
The decipher time of each cell.
- Decipher.compute_gene_patterns(adata, trajectory, l_scale=10000, n_samples=100)[source]#
Compute the gene patterns for a trajectory.
The trajectory’s points are sent through the decoders, thus defining distributions over the gene expression. The gene patterns are computed by sampling from these distribution.
- Parameters:
adata (AnnData) – The annotated data matrix.
trajectory (Trajectory) – A
Trajectoryobject containing the trajectory information.l_scale (float) – The library size scaling factor.
n_samples (int) – The number of samples to draw from the decoder to compute the gene pattern statistics.
- Return type:
- Returns:
The gene patterns for the trajectory. Dictionary keys:
mean: the mean gene expression pattern
q25: the 25% quantile of the gene expression pattern
q75: the 75% quantile of the gene expression pattern
times: the times of the trajectory
- Decipher.compute_imputed_gene_expression(adata=None, indices=None, batch_size=None, compute_covariances=False, v_obsm_key=None, z_obsm_key=None)[source]#
Impute gene expression from the decipher model.
- Parameters:
adata (
AnnData|None(default:None)) – The annotated data matrix.indices (
Sequence[int] |None(default:None)) – Indices of the data to get the latent representation of.batch_size (
int|None(default:None)) – Batch size to use for the data loader.compute_covariances (
bool(default:False)) – Whether to compute the covariances between the Decipher v and each gene.v_obsm_key (
str|None(default:None)) – Key in adata.obsm to use for the Decipher v. Required if compute_covariances is True.z_obsm_key (
str|None(default:None)) – Key in adata.obsm to use for the Decipher z. Required if compute_covariances is True.
- Return type:
- Returns:
The imputed gene expression, and the covariances between the Decipher v and each gene if compute_covariances is True.
- classmethod Decipher.convert_legacy_save(dir_path, output_dir_path, overwrite=False, prefix=None, **save_kwargs)[source]#
Converts a legacy saved model (<v0.15.0) to the updated save format.
- Parameters:
dir_path (
str) – Path to the directory where the legacy model is saved.output_dir_path (
str) – Path to save converted save files.overwrite (
bool(default:False)) – Overwrite existing data or not. IfFalseand directory already exists atoutput_dir_path, an error will be raised.prefix (
str|None(default:None)) – Prefix of saved file names.**save_kwargs – Keyword arguments passed into
save().
- Return type:
- Decipher.data_registry(registry_key)[source]#
Returns the object in AnnData associated with the key in the data registry.
- Decipher.deregister_manager(adata=None)[source]#
Deregisters the
AnnDataManagerinstance associated with adata.If adata is None, deregisters all
AnnDataManagerinstances in both the class and instance-specific manager stores, except for the one associated with this model instance.
- Decipher.differential_abundance(*args, **kwargs)[source]#
Not implemented for this model class.
Available in models that inherit from
VAEMixin.- Raises:
- Decipher.get_anndata_manager(adata, required=False)[source]#
Retrieves the
AnnDataManagerfor a given AnnData object.Requires
self.idhas been set. Checks for anAnnDataManagerspecific to this model instance.- Parameters:
- Return type:
- Decipher.get_from_registry(adata, registry_key)[source]#
Returns the object in AnnData associated with the key in the data registry.
AnnData object should be registered with the model prior to calling this function via the
self._validate_anndatamethod.
- Decipher.get_latent_representation(adata=None, indices=None, batch_size=None, give_z=False)[source]#
Get the latent representation of the data.
- Parameters:
adata (
AnnData|None(default:None)) – AnnData object with the data to get the latent representation of.indices (
Sequence[int] |None(default:None)) – Indices of the data to get the latent representation of.batch_size (
int|None(default:None)) – Batch size to use for the data loader.give_z (
bool(default:False)) – Whether to return the intermediate latent space z or the top-level latent space v.
- Return type:
- Decipher.get_normalized_expression(*args, **kwargs)[source]#
Not implemented for this model class.
Available in RNA models that inherit from
RNASeqMixin.- Raises:
- Decipher.get_setup_arg(setup_arg)[source]#
Returns the string provided to setup of a specific setup_arg.
- Return type:
- Decipher.get_state_registry(registry_key)[source]#
Returns the state registry for the AnnDataField registered with this instance.
- Return type:
- Decipher.get_var_names(legacy_mudata_format=False)[source]#
Variable names of input data.
- Return type:
- classmethod Decipher.load(dir_path, adata=None, accelerator='auto', device='auto', prefix=None, backup_url=None, datamodule=None, allowed_classes_names_list=None)[source]#
Instantiate a model from the saved output.
- Parameters:
dir_path (
str) – Path to saved outputs.adata (
AnnData|MuData|None(default:None)) – AnnData organized in the same way as data used to train model. It is not necessary to run setup_anndata, as AnnData is validated against the saved scvi setup dictionary. If None, will check for and load anndata saved with the model. If False, will load the model without AnnData.accelerator (
str(default:'auto')) – Supports passing different accelerator types (“cpu”, “gpu”, “tpu”, “ipu”, “hpu”, “mps, “auto”) as well as custom accelerator instances.device (
int|str(default:'auto')) – The device to use. Can be set to a non-negative index (int or str) or “auto” for automatic selection based on the chosen accelerator. If set to “auto” and accelerator is not determined to be “cpu”, then device will be set to the first available device.prefix (
str|None(default:None)) – Prefix of saved file names.backup_url (
str|None(default:None)) – URL to retrieve saved outputs from if not present on disk.datamodule (
LightningDataModule|None(default:None)) –EXPERIMENTALALightningDataModuleinstance to use for training in place of the defaultDataSplitter. Can only be passed in if the model was not initialized withAnnData.allowed_classes_names_list (
list[str] |None(default:None)) – list of allowed classes names to be loaded (besides the original class name)
- Returns:
Model with loaded state dictionaries.
Examples
>>> model = ModelClass.load(save_path, adata) >>> model.get_....
- static Decipher.load_registry(dir_path, prefix=None)[source]#
Return the full registry saved with the model.
- classmethod Decipher.register_manager(adata_manager)[source]#
Registers an
AnnDataManagerinstance with this model class.Stores the
AnnDataManagerreference in a class-specific manager store. Intended for use in thesetup_anndata()class method followed up by retrieval of theAnnDataManagervia the_get_most_recent_anndata_manager()method in the model init method.Notes
Subsequent calls to this method with an
AnnDataManagerinstance referring to the same underlying AnnData object will overwrite the reference to previousAnnDataManager.
- Decipher.save(dir_path, prefix=None, overwrite=False, save_anndata=False, save_kwargs=None, legacy_mudata_format=False, datamodule=None, **anndata_write_kwargs)[source]#
Save the state of the model.
Neither the trainer optimizer state nor the trainer history are saved. Model files are not expected to be reproducibly saved and loaded across versions until we reach version 1.0.
- Parameters:
dir_path (
str) – Path to a directory.prefix (
str|None(default:None)) – Prefix to prepend to saved file names.overwrite (
bool(default:False)) – Overwrite existing data or not. If False and directory already exists at dir_path, an error will be raised.save_anndata (
bool(default:False)) – If True, also saves the anndatasave_kwargs (
dict|None(default:None)) – Keyword arguments passed intosave().legacy_mudata_format (
bool(default:False)) – IfTrue, saves the modelvar_namesin the legacy format if the model was trained with aMuDataobject. The legacy format is a flat array with variable names across all modalities concatenated, while the new format is a dictionary with keys corresponding to the modality names and values corresponding to the variable names for each modality.datamodule (
LightningDataModule|None(default:None)) –EXPERIMENTALALightningDataModuleinstance to use for training in place of the defaultDataSplitter. Can only be passed in if the model was not initialized withAnnData.anndata_write_kwargs – Kwargs for
write()
- classmethod Decipher.setup_annbatch(cls, collection_path=None, paths=None, batch_key=None, labels_key=None, sample_key=None, unlabeled_category='Unknown', layer=None, categorical_covariate_keys=None, continuous_covariate_keys=None, rebuild=True, batch_size=4096, chunk_size=256, preload_nchunks=32, preload_to_gpu=True, dataset_size='20GB', shuffle=False, var_subset=None, merge=None, adatas=None, use_class_sampler=False, class_sampler_key=None, class_weights=None)[source]#
Set up a model from on-disk AnnData files via the annbatch streaming loader.
Builds (or reuses) a zarr-backed
DatasetCollectionfrom the supplied h5ad file paths, then wraps it in aAnnbatchDataModuleready for training.- Parameters:
collection_path (
str|None(default:None)) – Directory where the zarr collection is written (or already exists). IfNone(default), a path is auto-generated as"./{ModelName}_annbatch.zarr"in the current working directory.paths (
list[str] |None(default:None)) – Paths to h5ad files that make up the training dataset. IfNone, an existing collection atcollection_pathis opened without rebuilding — useful when the zarr store was created by a previous call. Must be provided when the store does not exist yet.batch_key (
str|None(default:None)) – Column inobsto use as the batch variable.labels_key (
str|None(default:None)) – Column inobsto use as the cell-type / label variable.sample_key (
str|None(default:None)) – Column inobsto use as the sample variable. Used by models like MrVI.unlabeled_category (
str(default:'Unknown')) – Value used to mark unlabeled cells inlabels_key. Required by semi-supervised models such as SCANVI. Defaults to"Unknown".layer (
str|None(default:None)) – Layer in the h5ad files to use as the count matrix.Noneusesadata.X(falling back toadata.raw.Xwhen a raw slot exists).categorical_covariate_keys (
list[str] |None(default:None)) – Additional categorical covariate columns inobs.continuous_covariate_keys (
list[str] |None(default:None)) – Additional continuous covariate columns inobs.rebuild (
bool(default:True)) – IfTrue, always rebuild the zarr collection even if it already exists on disk. IfFalse(default), reuse an existing collection and skip the (potentially expensive)add_adatasstep.batch_size (
int(default:4096)) – Number of cells per batch yielded by theLoader.chunk_size (
int(default:256)) – Number of cells loaded from disk contiguously per read.preload_nchunks (
int(default:32)) – Number of chunks to preload and shuffle in memory.preload_to_gpu (
bool(default:True)) – Whether the loader should move data to GPU before yielding.dataset_size (
int|str(default:'20GB')) – Number of observations to load into memory for shuffling / pre-processing when building the collection, or annbatch’s human-readable size strings, e.g."20GB".shuffle (
bool(default:False)) – Whether to pre-shuffle cells when building the collection.var_subset (
list[str] |None(default:None)) – Optional list of gene names to restrict the collection to (passed asvar_subsettoadd_adatas).merge (
Optional[Literal['same','unique','first','only']] (default:None)) – How annbatch should mergevarmetadata across inputs. Passed through toannbatch.DatasetCollection.add_adatas().
- Returns:
AnnbatchDataModuleA configured datamodule whoseregistrycan be passed directly to the model constructor, e.g.model = SCVI(registry=dm.registry).
Notes
After training, saving and reloading the model requires reconstructing the datamodule manually and passing it to
load()via thedatamoduleargument — the same pattern used by all custom datamodules.
- classmethod Decipher.setup_anndata(adata, layer=None, **kwargs)[source]#
Sets up the
AnnDataobject for this model.A mapping will be created between data fields used by this model to their respective locations in adata. None of the data in adata are modified. Only adds fields to adata.
- Decipher.to_device(device)[source]#
Move the model to the device.
- Parameters:
device (
str|int|device) – Device to move model to. Options: ‘cpu’ for CPU, integer GPU index (e.g., 0), ‘cuda:X’ where X is the GPU index (e.g. ‘cuda:0’), or a torch.device object (including XLA devices for TPU). See torch.device for more info.
Examples
>>> adata = scvi.data.synthetic_iid() >>> model = scvi.model.SCVI(adata) >>> model.to_device("cpu") # moves model to CPU >>> model.to_device("cuda:0") # moves model to GPU 0 >>> model.to_device(0) # also moves model to GPU 0
- Decipher.train(max_epochs=None, accelerator='auto', device='auto', train_size=0.9, validation_size=None, shuffle_set_split=True, batch_size=128, early_stopping=False, training_plan=None, datasplitter_kwargs=None, plan_kwargs=None, **trainer_kwargs)[source]#
Train the model.
Wraps
train()with Decipher-specific defaults (early_stopping_monitor="nll_validation"anddrop_last=True).- Parameters:
max_epochs (
int|None(default:None)) – Number of passes through the dataset.accelerator (
str(default:'auto')) – Supports passing different accelerator types("cpu", "gpu", "tpu", "ipu", "hpu", "mps", "auto")as well as custom accelerator instances.device (
int|str(default:'auto')) – The device to use. Can be set to a non-negative index (int or str) or"auto"for automatic selection.train_size (
float(default:0.9)) – Size of training set in the range[0.0, 1.0].validation_size (
float|None(default:None)) – Size of the validation set. IfNone, defaults to1 - train_size.shuffle_set_split (
bool(default:True)) – Whether to shuffle indices before splitting.batch_size (
int(default:128)) – Minibatch size to use during training.early_stopping (
bool(default:False)) – Perform early stopping. Additional arguments can be passed in**trainer_kwargs.training_plan (
DecipherTrainingPlan|None(default:None)) – Training plan instance. IfNone, a defaultDecipherTrainingPlanis used.datasplitter_kwargs (
dict|None(default:None)) – Additional keyword arguments passed intoDataSplitter.plan_kwargs (
dict|None(default:None)) – Keyword arguments forDecipherTrainingPlan.**trainer_kwargs – Additional keyword arguments passed to
Trainer.
- Decipher.transfer_fields(adata, **kwargs)[source]#
Transfer fields from a model to an AnnData object.
- Return type:
- Decipher.update_setup_method_args(setup_method_args)[source]#
Update setup method args.
- Parameters:
setup_method_args (
dict) – This is a bit of a misnomer, this is a dict representing kwargs of the setup method that will be used to update the existing values in the registry of this instance.
- Decipher.view_anndata_setup(adata=None, hide_state_registries=False)[source]#
Print summary of the setup for the initial AnnData or a given AnnData object.
- Parameters:
adata (
AnnData|MuData|None(default:None)) – AnnData object setup withsetup_anndataortransfer_fields().hide_state_registries (
bool(default:False)) – If True, prints a shortened summary without details of each state registry.
- Return type: