dynamo.tl.cell_growth_rate
- dynamo.tl.cell_growth_rate(adata, group, source, target, L0=0.3, L=1.2, k=0.001, birth_genes=None, death_genes=None, clone_column=None, **kwargs)[source]
Estimate the growth rate via clone information or logistic equation of population dynamics.
Growth rate is calculated as 1) number_of_cell_at_source_time_in_the_clone / number_of_cell_at_end_time_in_the_clone when there is clone information ([clone_column, time_column, source_time, target_time] are all not None); 2) estimate via logistic equation of population growth and death.
- Parameters:
adata (
AnnData
) – AnnData object that contains the reconstructed vector field function in the uns attribute.group (str or None (default: None)) – The column key in .obs points to the collection time of each cell, required for calculating growth rate with clone information.
source (str or None (default: None)) – The column key in .obs points to the starting point from collection time of each cell, required for calculating growth rate with clone information.
target (str or None (default: None)) – The column key in .obs points to the end point from collection time of each cell, required for calculating growth rate with clone information.
L0 (float (default: 0.3)) – The base growth/death rate.
L (float (default: 1.2)) – The maximum growth/death rate.
k (float (default: `0.001)) – The steepness of the curve.
birth_genes (list or None (default: None)) – The gene list associated with the cell cycle process. If None, GSEA’s KEGG_CELL_CYCLE will be used.
death_genes (list or None (default: None)) – The gene list associated with the cell cycle process. If None, GSEA’s KEGG_APOPTOSIS will be used.
clone_column (str or None (default: None)) – The column key in .obs points to the clone id if there is any. If a cell doesn’t belong to any clone, the clone id of that cell should be assigned as np.nan
kwargs – Additional arguments that will be passed to score_cells function.
- Returns:
An updated adata object that includes growth_rate column or growth_rate, birth_score, death_score in its .obs attribute when the clone based or purely expression based growth rate was calculated.