cellrank icon indicating copy to clipboard operation
cellrank copied to clipboard

How to get differential kinetics

Open roskcl opened this issue 6 months ago • 1 comments

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) scvelo_stream_single kinetics

With this code I get the corrected vectors: scv.tl.velocity(adata, diff_kinetics=True) scv.tl.velocity_graph(adata) scvelo_stream

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() cellrank_transition_stream

Is this expected? Do I need to add any additional argument for the differential kinetics to be considered?

Thank you for your help!

roskcl avatar Dec 04 '23 22:12 roskcl