schist.inference._multi_nested
Module Contents
Functions
|
Cluster cells into subgroups using multiple modalities. |
- schist.inference._multi_nested.nested_model_multi(adatas: List[anndata.AnnData], deg_corr: bool = True, tolerance: float = 1e-06, n_sweep: int = 10, beta: float = np.inf, n_init: int = 100, collect_marginals: bool = True, n_jobs: int = -1, refine_model: bool = False, refine_iter: int = 100, overlap: bool = False, max_iter: int = 100000, *, random_seed: int | None = None, key_added: str = 'multi_nsbm', adjacency: List[scipy.sparse.spmatrix] | None = None, neighbors_key: List[str] | None = ['neighbors'], directed: bool = False, use_weights: bool = False, save_model: str | None = None, copy: bool = False, dispatch_backend: str | None = 'threads') List[anndata.AnnData] | None
Cluster cells into subgroups using multiple modalities.
Cluster cells using the nested Stochastic Block Model [Peixoto14], performing Bayesian inference on node groups. This function takes multiple experiments, possibly across different modalities, and perform joint clustering.
This requires having ran
neighbors()
orbbknn()
first. It also requires cells having the same names if coming from paired experimentsParameters
- adatas
A list of processed AnnData. Neighbors must have been already calculated.
- deg_corr
Whether to use degree correction in the minimization step. In many real world networks this is the case, although this doesn’t seem the case for KNN graphs used in scanpy.
- tolerance
Tolerance for fast model convergence.
- n_sweep
Number of iterations to be performed in the fast model MCMC greedy approach
- beta
Inverse temperature for MCMC greedy approach
- n_init
Number of initial minimizations to be performed. The one with smaller entropy is chosen
- refine_model
Wether to perform a further mcmc step to refine the model
- refine_iter
Number of refinement iterations.
- max_iter
Maximum number of iterations during minimization, set to infinite to stop minimization only on tolerance
- overlap
Whether the different layers are dependent (overlap=True) or not (overlap=False)
- n_jobs
Number of parallel computations used during model initialization
- key_added
adata.obs key under which to add the cluster labels.
- adjacency
Sparse adjacency matrix of the graph, defaults to adata.uns[‘neighbors’][‘connectivities’] in case of scanpy<=1.4.6 or adata.obsp[neighbors_key][connectivity_key] for scanpy>1.4.6
- neighbors_key
The key passed to sc.pp.neighbors. If all AnnData share the same key, one only has to be specified, otherwise the full tuple of all keys must be provided
- directed
Whether to treat the graph as directed or undirected.
- use_weights
If True, edge weights from the graph are used in the computation (placing more emphasis on stronger edges). Note that this increases computation times
- save_model
If provided, this will be the filename for the PartitionModeState to be saved
- copy
Whether to copy adata or modify it inplace.
- random_seed
Random number to be used as seed for graph-tool
Returns
- adata.obs[key_added]
Array of dim (number of cells) that stores the subgroup id (‘0’, ‘1’, …) for each cell.
- adata.uns[‘schist’][‘multi_level_params’]
A dict with the values for the parameters resolution, random_state, and n_iterations.
- adata.uns[‘schist’][‘multi_level_stats’]
A dict with the values returned by mcmc_sweep
- adata.obsm[‘CA_multi_nsbm_level_{n}’]
A np.ndarray with cell probability of belonging to a specific group
- adata.uns[‘schist’][‘multi_level_state’]
The NestedBlockModel state object