GoodhillLab / neural-assembly-detection

Code published with the paper J. Mölter, L. Avitan, and G. J. Goodhill. "Detecting neural assemblies in calcium imaging data". BMC Biology 16:143 (2018) doi: 10.1186/s12915-018-0606-4.
GNU General Public License v3.0
12 stars 6 forks source link

Calcium fluorescence processing crash - oopsi algorithm #1

Open catubc opened 2 years ago

catubc commented 2 years ago

Hello We are trying to convert some raw data to use with the SGC algorithm and are getting a crash during the oopsi algorithm. Any advice is appreciated. Thanks! catubc

>> CALCIUM_FLUORESCENCE_PROCESSING( [ '/media/cat/4TBSSD/donato/' 'data' '/' 'HEX5_K5_medium__T1800s_dT500ms_cT1s_sKInf_eF10.0mHz_eM6_nS0' '_CALCIUM-FLUORESCENCE.mat' ] )

________________________________________________________________________________
:: PROCESS CALCIUM FLUORESCENCE

Prepare *_ACTIVITY-RASTER.MAT file for similarity-graph-clustering ...
significance threshold: 2STD
Elapsed time is 226.788281 seconds.

Prepare *_ACTIVITY-FIM-RASTER.DAT file for frequent-item-set-mining ...
Elapsed time is 1.187902 seconds.

Prepare *_SPIKE-PROBABILITY-RASTER.MAT file ...
Error:File: fast_oopsi.m Line: 195 Column: 42
Identifier 'b' is not a function or a shared variable. To share 'b' with nested function, initialize it in the current scope. For more information, see Sharing Variables Between Parent and Nested Functions.

Error in findDF_FSpikeProbability (line 25)
    [ n_best(:,n) , P_best(n) , ~ , ~ ] = fast_oopsi( F , V , P );

Error in CALCIUM_FLUORESCENCE_PROCESSING (line 172)
            [ ~ , spike_probability_raster , spike_probability_raster_thresholds , oopsi_deconvolution ] = findDF_FSpikeProbability( CALCIUM_FLUORESCENCE_mat.calcium_fluorescence.dF_F , CALCIUM_FLUORESCENCE_mat.parameter.dT_step ,
            CALCIUM_FLUORESCENCE_mat.parameter.calcium_T1_2 );
hrmneuro commented 1 year ago

Hi!

We experienced exactly the same problem. It seems to be a problem Matlab version related issue. When running the same script with the same input data in Matlab 2016a, the script runs without problems.

Best,. Harmen