suite2p
suite2p copied to clipboard
Python deconvolution implementation seems very different from OASIS and old Matlab implementation
Sorry if I'm missing something obvious, but the deconvolution implementation appears to be quite different from the original OASIS implementation and the original C/MATLAB version of Suite2p. As far as I can tell, the Suite2p preprint and JNeuro paper are written about the MATLAB implementation, but I'm wondering if the newer, seemingly more minimalistic Python implementation is written about anywhere, specifically how it differs from the original OASIS method (aside from not inferring spike trains or imposing L1 penalties)?
Is the deconvolved output supposed to show:
- Underlying calcium activity after removing the influence of the sensor and frame rate?
- Underlying calcium transient amplitudes
- Underlying action potential amplitudes
- Something else?
Any clarification is greatly appreciated.
piling on @RandallJEllis question I open discussion post #1001
We re-implemented OASIS (this function here: https://github.com/j-friedrich/OASIS/blob/6dae97381b295f764530d952a3b2bf5033a009b1/oasis/oasis_methods.pyx#L25) in numba without L0 and L1 penalties. The OASIS paper discusses the interpretations of spike deconvolved traces, and we discuss this in our review paper, it is correlated to the firing rate of the neuron: https://www.sciencedirect.com/science/article/pii/S0959438818300977