4
« on: November 02, 2016, 06:20:17 AM »
Hi Daniel Aharoni,
Thanks for quick reply. Yes. I already have everything to generate spatial firing rate maps, like you said(behav.position and ms.firing). There are few questions to generate rate map.
I got the value of x-position(horizontal) from behav.position(:,1) and calculated occupancy of each bin(histcn/sampling_rate (sec)). In this step, I was expect that sampling rate is 30fps. It wasn't, because maxFramesperFile was 1000. Also, ms.maxFramesperFile was 1000. It means that down sampling 30fps to 20fps.
Q1. Should I use maxFramesperFile to 1000? Why limit maxFramesperFile? Down sampling doesn't matter to calculate occupancy. I am just curious.
I got the firing data from ms.firing(:,segments). When I generate rate map(or place map) with spikes data from single unit recording, spike time was accurate. I find peak point using findpeaks() with threshold. In your paper(Nature,2016), 'Place fields were calculated by deconvolving calcium dF/F traces with an exponential to extract approximate spike times.'
Q2. How do i extract approximate spike time from data? Is one peak of firing one spike?
Finally, I try to run script with all msCam data from G200(behavior data limited 1 to 14). I think data was split because parallel processing. When running msBatchSegment, MATLAB has memory problem. Physical memory was enough when i was monitor memory information.
Q3. How do i solve memory problem?
I'm newbie in optical recording. Sometimes, I asking stupid question. Please forgive me. Your answer will help to generate rate map as you draw. It will be great experience to analyze own data.
Thanks.