Differences

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

--- plotcurves [2010/06/19 19:03] – added output reference jochen
+++ plotcurves [2010/06/21 03:12] – changed the post-proc filter to 0.2s (enough!) jochen
@@ Line 24: / Line 24: @@
 % this filters channel 1 into a new channel (10)
-% all parameters are default, the post-filtering cut-off is set to 0.5s
+% all parameters are default, the post-filtering cut-off is set to 0.2s
-gsr.Filter(1, struct('dest', 10, 'post', 0.5));</code>
+gsr.Filter(1, struct('dest', 10, 'post', 0.2));</code>
 ===== Function reference ('help plotcurves') =====
@@ Line 157: / Line 157: @@
 <code matlab plotcurves_withonsets.m>% data is in object gsr and already preprocessed
-% we define three spots of interest:
+% for this example we define four (!) spots of interest:
+% - cue (time of stimulus minus 2 secs)
 % - onset (time of stimulus),
 % - delay of response (minimum before peak),
 % - and peak of response (amplitude)
-spotnames = {'onset', 'delay', 'peak'};
+spotnames = {'cue', 'onset', 'delay', 'peak'};
-spottype = {'onset', 'min', 'max'};
+spottype = {'cue', 'onset', 'min', 'max'};
 % the spot values are set to empty arrays with size 0x2!
 % this tells plotcurves to autodetect...
-spot = {zeros(0,2), zeros(0,2), zeros(0, 2)};
+spot = {zeros(0,2), zeros(0,2), zeros(0,2), zeros(0, 2)};
+% then we define the variables of interest:
+% NOTE: since we have a CUE spot, all spot indices shift by one!
+% - Latency (defined as x-diff between onset and delay)
+% - Response time (defined as x-diff between delay and peak)
+% - Base amplitude (defined as y-value at delay)
+% - Peak amplitude (defined as y-value at peak)
+% - Peak-to-peak response amplitude (defined as y-diff between delay and peak)
+% additionally, we define that for the response amplitude, the log+1 transform
+% should be applied after we're done
+vars = struct( ...
+    'name', {'Latency', 'Resp. time', 'Base ampl.', 'Peak ampl.', 'Resp. ampl.'}, ...
+    'calc', {'dx', 'dx', 'y', 'y', 'dy'}, ...
+    'spot', {[2, 3], [3, 4], [3], [4], [3, 4]}, ...
+    'trans', {'none', 'none', 'none', 'none', 'log+1'});
 % re-using the variables definition from above, options:
@@ Line 172: / Line 188: @@
 % - original data is in channel 1 (pre-filter)
 % - onsets are coded in channel 8 (0Volts nothing, 5 Volts boxes)
+% - cue onffset is 2 seconds
+% - thus, the window is increased! [-4, 20]
 % - variables and spots as defined
 pcopts = struct( ...
@@ Line 177: / Line 195: @@
     'odchannel', 1, ...
     'ochannel',  8, ...
+    'cuediff',   2, ...
+    'owin',      [-4, 20],...
     'var',       vars, ...
     'spot',      {spot}, ...