Continuous SSP projection about mne-python HOT 7 CLOSED

we could add this but maybe a better way is indeed to make epochs from raw with fixed length such as 1s. I would allow you do reject bad segments for the SSP computation.

wdyt?

from mne-python.

I think both ways will make sense for different uses. In my code, I actually don't want epoch rejection because there is some large, low-frequency drift (that goes away with a couple continuous projectors) that would disqualify all epochs. I currently have code for doing it by chunking epochs with a chosen fixed length, but I was hoping to get something that was actually equivalent to what is done in mne_process_raw (to stay as close as possible).

from mne-python.

I noticed that the way the projectors are computed in compute_proj_evoked in proj.py is by taking the sum across epochs of the matrix multiplications (channelsxtime into timexchannels). I am curious why this is done instead of using the matrix multiplication of the temporally-concatenated (across epochs) matrices (chanellsxconcatenated_time into concatenated_timexchannels)? I understand the former is probably more efficient, but is one preferable to another from a mathematical standpoint? I ask because I noticed that I got different results when "faking" continuous projection calculation by using consecutive 1, 10, or 100 second windows to create epochs from my raw data.

from mne-python.

I submitted some continuous SSP code to code review so we can discuss it there.

from mne-python.

I currently have code for doing it by chunking epochs with a chosen fixed length

a function like:

mne.make_fix_length_events(raw, duration=4)

would be convenient to have.

I noticed that the way the projectors are computed in compute_proj_evoked
in proj.py is by taking the sum across epochs of the matrix multiplications
(channelsxtime into timexchannels). I am curious why this is done instead of
using the matrix multiplication of the temporally-concatenated (across
epochs) matrices (chanellsxconcatenated_time into
concatenated_timexchannels)?

memory allocation. You don't need to read everything in memory by doing so.

I understand the former is probably more
efficient, but is one preferable to another from a mathematical standpoint?

no just efficiency

I ask because I noticed that I got different results when "faking"
continuous projection calculation by using consecutive 1, 10, or 100 second
windows to create epochs from my raw data.

hum. can you submit a test script?

from mne-python.

I'll add the make_fix_length_events(raw, duration=4) [but I'll probably make the default 1 second] function to my modified codebase, where we can continue the discussion. Closing this one.

from mne-python.

perfect !

from mne-python.