-
Notifications
You must be signed in to change notification settings - Fork 146
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Error with p value #108
Comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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.
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!
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
The text was updated successfully, but these errors were encountered: