Differences

This shows you the differences between two versions of the page.

--- gsr_data_analysis [2010/06/18 21:17] – changing the syntax of prefilter to xff method (obj.Filter) jochen
+++ gsr_data_analysis [2010/06/18 23:22] (current) – changed onsets to be in ms jochen
@@ Line 89: / Line 89: @@
 <code matlab GSR_acq_stimonsets.m>% all stimuli onsets (use 1 + to cope for the fact that the
 % first value in diff rather codes the second sample in the data!
-all_onsets = 1 + find(diff(gsr.RawData(8, :)') > 0);
+all_onsets = 1 + find(diff(gsr.ChannelData(8)) > 0);
 % now split into three conditions
-neu = all_onsets(gsr.RawData(2, all_onsets) > 0);
+neu = 0.01 * all_onsets(gsr.ChannelData(2, all_onsets) > 0);
-neg = all_onsets(gsr.RawData(3, all_onsets) > 0);
+neg = 0.01 * all_onsets(gsr.ChannelData(3, all_onsets) > 0);
-rea = all_onsets(gsr.RawData(4, all_onsets) > 0);
+rea = 0.01 * all_onsets(gsr.ChannelData(4, all_onsets) > 0);</code>
-% cope for resampling from 500Hz -> 100Hz !!
-neu = round(0.2 .* neu);
-neg = round(0.2 .* neg);
-rea = round(0.2 .* rea);</code>
 As window length for the inspection of the GSR data, we choose 15 seconds, so in the resolution of the resampled data that means 1,500 samples.
@@ Line 141: / Line 136: @@
 reacurves = cell(numel(rea), 2);
 for sc = 1:numel(neu)
-    neucurves{sc, 1} = sprintf('LookNeu %02d at t=%.3fs', sc, 0.01 * neu(sc));
+    neucurves{sc, 1} = sprintf('LookNeu %02d at t=%.3fs', sc, neu(sc));
-    neucurves{sc, 2} = [gsrchannel, neu(sc) - 100, neu(sc) + 1400];
+    neucurves{sc, 2} = [gsrchannel, neu(sc) - 2, neu(sc) + 14];
 end
 for sc = 1:numel(neg)
-    negcurves{sc, 1} = sprintf('LookNeg %02d at t=%.3fs', sc, 0.01 * neg(sc));
+    negcurves{sc, 1} = sprintf('LookNeg %02d at t=%.3fs', sc, neg(sc));
-    negcurves{sc, 2} = [gsrchannel, neg(sc) - 100, neg(sc) + 1400];
+    negcurves{sc, 2} = [gsrchannel, neg(sc) - 2, neg(sc) + 14];
 end
 for sc = 1:numel(rea)
-    reacurves{sc, 1} = sprintf('ReaNeg %02d at t=%.3fs', sc, 0.01 * rea(sc));
+    reacurves{sc, 1} = sprintf('ReaNeg %02d at t=%.3fs', sc, rea(sc));
-    reacurves{sc, 2} = [gsrchannel, rea(sc) - 100, rea(sc) + 1400];
+    reacurves{sc, 2} = [gsrchannel, rea(sc) - 2, rea(sc) + 14];
 end