Error with p value

[oregonGCAMP]

Please note that the developer of this package is on leave until January 2019 and might not be able to support you.

For better support, please use the template below to submit your issue. When your issue gets resolved please remember to close it.

• Describe the issue that you are experiencing

I noticed that CaImAn crashes when the p value (autoregressive model) is 2, but works if it is changed to 1. I included 2 image files where this occurred.

Thanks!

• Copy error log below

ValueError Traceback (most recent call last)
in
183 #RE-RUN seeded CNMF on accepted patches to refine and perform deconvolution
184 cnm.params.change_params({'p': p})
--> 185 cnm2 = cnm.refit(images, dview=dview)
186
187

~\AppData\Local\Continuum\anaconda3\envs\caiman\lib\site-packages\caiman\source_extraction\cnmf\cnmf.py in refit(self, images, dview)
414 cnm.estimates = estimates
415 cnm.mmap_file = self.mmap_file
--> 416 return cnm.fit(images)
417
418 def fit(self, images, indices=(slice(None), slice(None))):

~\AppData\Local\Continuum\anaconda3\envs\caiman\lib\site-packages\caiman\source_extraction\cnmf\cnmf.py in fit(self, images, indices)
564 self.params.set('temporal', {'p': self.params.get('preprocess', 'p')})
565 logging.info('update temporal ...')
--> 566 self.update_temporal(Yr, use_init=False)
567 # else:
568 # todo : ask for those..

~\AppData\Local\Continuum\anaconda3\envs\caiman\lib\site-packages\caiman\source_extraction\cnmf\cnmf.py in update_temporal(self, Y, use_init, **kwargs)
899 self.estimates.g, self.estimates.YrA, self.estimates.lam = update_temporal_components(
900 Y, self.estimates.A, self.estimates.b, self.estimates.C, self.estimates.f, dview=self.dview,
--> 901 **self.params.get_group('temporal'))
902 self.estimates.R = self.estimates.YrA
903 return self

~\AppData\Local\Continuum\anaconda3\envs\caiman\lib\site-packages\caiman\source_extraction\cnmf\temporal.py in update_temporal_components(Y, A, b, Cin, fin, bl, c1, g, sn, nb, ITER, block_size_temp, num_blocks_per_run_temp, debug, dview, **kwargs)
215 logging.info("entering the deconvolution ")
216 C, S, bl, YrA, c1, sn, g, lam = update_iteration(parrllcomp, len_parrllcomp, nb, C, S, bl, nr,
--> 217 ITER, YrA, c1, sn, g, Cin, T, nA, dview, debug, AA, kwargs)
218 ff = np.where(np.sum(C, axis=1) == 0) # remove empty components
219 if np.size(ff) > 0: # Eliminating empty temporal components

~\AppData\Local\Continuum\anaconda3\envs\caiman\lib\site-packages\caiman\source_extraction\cnmf\temporal.py in update_iteration(parrllcomp, len_parrllcomp, nb, C, S, bl, nr, ITER, YrA, c1, sn, g, Cin, T, nA, dview, debug, AA, kwargs)
348 if 'multiprocessing' in str(type(dview)):
349 results = dview.map_async(
--> 350 constrained_foopsi_parallel, args_in).get(4294967)
351
352 elif dview is not None and platform.system() != 'Darwin':

~\AppData\Local\Continuum\anaconda3\envs\caiman\lib\multiprocessing\pool.py in get(self, timeout)
655 return self._value
656 else:
--> 657 raise self._value
658
659 def _set(self, i, obj):

ValueError: shapes (82,82) and (0,) not aligned: 82 (dim 1) != 0 (dim 0)

ImageExample.zip

• If you're not reporting an error, type your message below