Optimal transport for time-course single cell data
Uses time-course data to infer how the probability distribution of cells in gene-expression space evolves over time, by using the mathematical approach of Optimal Transport (OT)
- Prepare Gene Expression Matrix
- Optimal Transport
- Optimal Transport Transition Table
- Optimal Transport Validation
- Gene Set Scores
- File Formats
Python 3 is required.
Optionally install dependencies using conda:
conda install cython h5py flask gunicorn numexpr numpy pandas scikit-learn scipy simplejson psutil conda install -c conda-forge pot
pip install wot
WOT tools are run using the syntax wot tool. To see all tool options, type wot tool -h (e.g. wot optimal_transport -h)
wot optimal_transport calculates transport maps between consecutive time points and automatically learns cellular growth and death rates.
Common Options (required in bold)
|matrix||Normalized gene expression matrix.|
|cell_days||Assigns days to cells|
|gene_set_scores||Apoptosis and cell cycle scores used to compute growth rates. If not specified, a constant growth rate is used. The wot tool gene_set_scores can be used to compute gene set scores.|
|local_pca||Use principal component analysis to reduce the dimensionality of the expression matrix locally in the space of consecutive days. Thirty components are used by default.|
|day_pairs||Pairs of days to compute transport maps for|
|gene_filter||File with one gene id per line to use for computing cost matrices (e.g. variable genes)|
|out||Base name for output files|
wot optimal_transport_validation tests the performance of optimal transport by comparing interpolated distributions to held-out time points.
optimal_transport_validation has the same options as optimal_transport and
|t_interpolate||Interpolation fraction between two time points|
|covariate||Two column file with headers "id" and "covariate" indicating cell ids and covariate value|
wot force_layout generates a force-directed layout using the ForceAtlas2 algorithm. The tool can optionally perform dimensionality reduction using diffusion component embedding of the dataset. After you have generated the force layout coordinates, you can use wot wot_server to view trajectories and gene expression in the force layout.
Trajectory Trends plots the expression of a gene over time based on the transport maps
wot transition_table generate a transition table from one cell set to another cell set. Options
|dir||Directory of transport maps as produced by optimal_transport|
|cell_set||One or more gmt or gmx files containing cell sets. Each set id should end with _time (e.g. my_set_9).|
|start_time||The start time for the cell sets to compute the transitions to cell sets at end_time|
|end_time||The end time.|
wot gene_set_scores computes gene set scores for each cell given a gene expression matrix and gene sets.
|matrix||Normalized gene expression matrix to compute gene set scores (e.g. apoptosis and cell cycle scores) for each cell.|
Example Text File:
Two column file with header "id" and "day".
Two column file without header with pairs of days to compute transport maps for.
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
|Filename, size & hash SHA256 hash help||File type||Python version||Upload date|
|wot-0.1.5-py2.py3-none-any.whl (62.8 MB) Copy SHA256 hash SHA256||Wheel||py2.py3||Jun 20, 2018|
|wot-0.1.5.tar.gz (62.7 MB) Copy SHA256 hash SHA256||Source||None||Jun 20, 2018|