cellrank
cellrank copied to clipboard
How to get differential kinetics
I normally run the differential kinetics analysis with the original cellrank. But when I tried cellrank2, adding the diff_kinetics=True argument did not make any difference.
I have then checked, what happens if I run:
scv.tl.velocity(adata, mode='dynamical')
scv.tl.velocity_graph(adata)
With this code I get the corrected vectors:
scv.tl.velocity(adata, diff_kinetics=True)
scv.tl.velocity_graph(adata)
But in both cases, if I then run the cellrank kernel, I get the same results:
vk = cr.kernels.VelocityKernel(adata)
vk.compute_transition_matrix()
ck = cr.kernels.ConnectivityKernel(adata)
ck.compute_transition_matrix()
combined_kernel = 0.8 * vk + 0.2 * ck
vk.plot_projection()
Is this expected? Do I need to add any additional argument for the differential kinetics to be considered?
Thank you for your help!