cellrank.pl.gene_trends

cellrank.pl.gene_trends(adata, model, genes, lineages=None, backward=False, data_key='X', time_key='latent_time', transpose=False, time_range=None, callback=None, conf_int=True, same_plot=False, hide_cells=False, perc=None, lineage_cmap=None, abs_prob_cmap=<matplotlib.colors.ListedColormap object>, cell_color='black', cell_alpha=0.6, lineage_alpha=0.2, size=15, lw=2, cbar=True, margins=0.015, sharex=None, sharey=None, gene_as_title=None, legend_loc='best', ncols=2, suptitle=None, return_models=False, n_jobs=1, backend='loky', show_progress_bar=True, figsize=None, dpi=None, save=None, plot_kwargs=mappingproxy({}), **kwargs)

Plot gene expression trends along lineages.

Each lineage is defined via it’s lineage weights which we compute using cellrank.tl.lineages(). This function accepts any model based off cellrank.ul.models.BaseModel to fit gene expression, where we take the lineage weights into account in the loss function.

Parameters

• adata (anndata.AnnData) – Annotated data object.

• model (Union[BaseModel, Mapping[str, Mapping[str, BaseModel]]]) –

  Model based on cellrank.ul.models.BaseModel to fit.

  If a dict, gene and lineage specific models can be specified. Use '*' to indicate all genes or lineages, for example {'Map2': {'*': ...}, 'Dcx': {'Alpha': ..., '*': ...}}.

• genes (Union[str, Sequence[str]]) – Genes in adata.var_names or in adata.raw.var_names, if use_raw=True.

• lineages (Union[str, Sequence[str], None]) – Names of the lineages to plot. If None, plot all lineages.

• backward (bool) – Direction of the process.

• data_key (str) – Key in adata.layers or ‘X’ for adata.X where the data is stored.

• time_key (str) – Key in adata.obs where the pseudotime is stored.

• time_range (Union[float, Tuple[Optional[float], Optional[float]], None, List[Union[float, Tuple[Optional[float], Optional[float]], None]]]) –

  Specify start and end times:

  • If a tuple, it specifies the minimum and maximum pseudotime. Both values can be None, in which case the minimum is the earliest pseudotime and the maximum is automatically determined.

  • If a float, it specifies the maximum pseudotime.

  This can also be specified on per-lineage basis.

• transpose (bool) – If same_plot=True, group the trends by lineages instead of genes. This enforces hide_cells=True. If same_plot=False, show lineages in rows and genes in columns.

• callback (Union[Callable, Mapping[str, Mapping[str, Callable]], None]) – Function which takes a cellrank.ul.models.BaseModel and some keyword arguments for cellrank.ul.models.BaseModel.prepare() and returns the prepared model. Can be specified in gene- and lineage-specific manner, similarly to model.

• conf_int (Union[bool, float]) – Whether to compute and show confidence interval. If the model is cellrank.ul.models.GAMR, it can also specify the confidence level, the default is 0.95.

• same_plot (bool) – Whether to plot all lineages for each gene in the same plot.

• hide_cells (bool) – If True, hide all cells.

• perc (Union[Tuple[float, float], Sequence[Tuple[float, float]], None]) – Percentile for colors. Valid values are in interval [0, 100]. This can improve visualization. Can be specified individually for each lineage.

• lineage_cmap (Optional[ListedColormap]) – Categorical colormap to use when coloring in the lineages. If None and same_plot, use the corresponding colors in adata.uns, otherwise use ‘black’.

• abs_prob_cmap (ListedColormap) – Continuous colormap to use when visualizing the absorption probabilities for each lineage. Only used when same_plot=False.

• cell_color (str) – Color of the cells when not visualizing absorption probabilities. Only used when same_plot=True.

• cell_alpha (float) – Alpha channel for cells.

• lineage_alpha (float) – Alpha channel for lineage confidence intervals.

• size (float) – Size of the points.

• lw (float) – Line width of the smoothed values.

• cbar (bool) – Whether to show colorbar. Always shown when percentiles for lineages differ. Only used when same_plot=False.

• margins (float) – Margins around the plot.

• sharex (Union[str, bool, None]) – Whether to share x-axis. Valid options are ‘row’, ‘col’ or ‘none’.

• sharey (Union[str, bool, None]) – Whether to share y-axis. Valid options are ‘row’, ‘col’ or ‘none’.

• gene_as_title (Optional[bool]) – Whether to show gene names as titles instead on y-axis.

• legend_loc (Optional[str]) – Location of the legend displaying lineages. Only used when same_plot=True.

• ncols (int) – Number of columns of the plot when plotting multiple genes. Only used when same_plot=True.

• suptitle (Optional[str]) – Suptitle of the figure.

• return_models (bool) – If True, return the fitted models for each gene in genes and lineage in lineages.

• show_progress_bar (bool) – Whether to show a progress bar. Disabling it may slightly improve performance.

• n_jobs (Optional[int]) – Number of parallel jobs. If -1, use all available cores. If None or 1, the execution is sequential.

• backend (str) – Which backend to use for parallelization. See joblib.Parallel for valid options.

• figsize (Optional[Tuple[float, float]]) – Size of the figure.

• dpi (Optional[int]) – Dots per inch.

• save (Union[Path, str, None]) – Filename where to save the plot.

• plot_kwargs (Mapping) – Keyword arguments for cellrank.ul.models.BaseModel.plot().

• **kwargs – Keyword arguments for cellrank.ul.models.BaseModel.prepare().

Return type

Optional[Mapping[str, Mapping[str, BaseModel]]]

Returns

• None – If return_models=False, just plots the figure and optionally saves it based on save.

• Dict[str, Dict[str, cellrank.ul.models.BaseModel]] – Otherwise returns the fitted models as {'gene_1': {'lineage_1': <model_11>, ...}, ...}. Models which have failed will be instances of cellrank.ul.models.FailedModel.

Examples

Plot gene trends