dynamo.pl.sensitivity_heatmap
- dynamo.pl.sensitivity_heatmap(adata, cell_idx, skey='sensitivity', basis='pca', regulators=None, effectors=None, figsize=(7, 5), ncols=1, cmap='bwr', save_show_or_return='show', save_kwargs={}, **kwargs)[source]
Plot the Jacobian matrix for each cell as a heatmap.
Note that Jacobian matrix can be understood as a regulatory activity matrix between genes directly computed from the reconstructed vector fields.
- Parameters:
adata (
AnnData
) – an Annodata object with Jacobian matrix estimated.cell_idx (int or list) – The numeric indices of the cells that you want to draw the sensitivity matrix to reveal the regulatory activity.
skey (str (default: sensitivity)) – The key to the sensitivity dictionary in .uns.
basis (str) – The reduced dimension basis.
regulators (list or None (default: None)) – The list of genes that will be used as regulators for plotting the Jacobian heatmap, only limited to genes that have already performed Jacobian analysis.
effectors (List or None (default: None)) – The list of genes that will be used as targets for plotting the Jacobian heatmap, only limited to genes that have already performed Jacobian analysis.
figsize (None or [float, float] (default: None)) – The width and height of each panel in the figure.
ncols (int (default: 1)) – The number of columns for drawing the heatmaps.
cmap (str (default: bwr)) – The mapping from data values to color space. If not provided, the default will depend on whether center is set.
save_show_or_return (str {‘save’, ‘show’, ‘return’} (default: show)) – Whether to save, show or return the figure.
save_kwargs (dict (default: {})) – A dictionary that will passed to the save_fig function. By default it is an empty dictionary and the save_fig function will use the {“path”: None, “prefix”: ‘scatter’, “dpi”: None, “ext”: ‘pdf’, “transparent”: True, “close”: True, “verbose”: True} as its parameters. Otherwise you can provide a dictionary that properly modify those keys according to your needs.
kwargs – Additional arguments passed to sns.heatmap.
- Returns:
Nothing but plots the n_cell_idx heatmaps of the corresponding Jacobian matrix for each selected cell.
Examples
>>> import dynamo as dyn >>> adata = dyn.sample_data.hgForebrainGlutamatergic() >>> dyn.pp.recipe_monocle(adata) >>> dyn.tl.dynamics(adata) >>> dyn.tl.reduceDimension(adata) >>> dyn.tl.cell_velocities(adata, basis='pca') >>> dyn.vf.VectorField(adata, basis='pca') >>> valid_gene_list = adata[:, adata.var.use_for_transition].var.index[:2] >>> dyn.vf.sensitivity(adata, regulators=valid_gene_list[0], effectors=valid_gene_list[1]) >>> dyn.pl.sensitivity_heatmap(adata)