[FieldTrip] (no subject)

KatrinH Heimann katrinheimann at gmail.com
Wed Dec 10 12:56:19 CET 2014


Hey Arjen, hey all,

so, I used some if statements, yes, and now I really get lists of
trialbegin, trialend, offset and RT as I wanted them, however at the end
only one segment  (with the last infos of the lists) is created. I just
don't get what is wrong. Could you have another short look?
Thanks a million
Katrin

function [trl, event] = response_trialfun(cfg);



% read the header information and the events from the data

hdr   = ft_read_header(cfg.dataset);

event = ft_read_event(cfg.dataset);



% search for "trigger" events

value  = [event(find(strcmp('ceck', {event.value}))).value]';

sample = [event(find(strcmp('ceck', {event.value}))).sample]'



% search for "trigger" events

value2  = [event(find(strcmp('resp', {event.value}))).value]';

sample2 = [event(find(strcmp('resp', {event.value}))).sample]'



% determine the number of samples before and after the trigger

pretrig  = -round(cfg.trialdef.prestim  * hdr.Fs);

posttrig =  round(cfg.trialdef.poststim * hdr.Fs);



trl = [];

trlbegin = [];

trlend = []

offset = [];

RT = [];



   i=7 %delete training trials

   j=1

   k=1

   while i <= (length(sample))



   if sample2(j)-sample(i)<=2500 & sample2(j)-sample(i)>=1 % find resp
appearing after ceck within RT-window



    trlbegin (k) = sample (i) + pretrig;

    trlend (k)  = sample (i) + posttrig;

    offset (k)  = sample (i) - sample (i) + pretrig;

    RT(k)       = sample2 (j) - sample (i);



    i=i+1

    j=j+1

    k=k+1

    elseif sample2(j)-sample(i)<=0

    j=j+1

    else

    i=i+1

    end

   end



    newtrl = [trlbegin trlend offset RT];

    trl    = [trl; newtrl]





2014-12-10 7:19 GMT+01:00 Arjen Stolk <a.stolk8 at gmail.com>:

> Het katrin, you could make a couple of conditional if statements (eg
> creating and changing the value of a variable event1ctr). For example,
> conceptually:
>
> Event1ctr = 0
>
> Loop start
>
> If event1 & event1ctr == 0
> Event1ctr = 1
> Elseif event1 & event1ctr
> Trl = ... with a zero for trlend # no event2 found
> Event1ctr = 0
> End
>
> If event2
> Trl = ... # normal trial
> Event1ctr= 0
> End
>
> Loop end
> Op 9 dec. 2014 23:52 schreef "KatrinH Heimann" <katrinheimann at gmail.com>
> het volgende:
>
>> Hey Fieldtrippers,
>>
>> I am fighting with a self defined trialfunction.
>> Basically I want to cut always at one event called ceck. If there is
>> another event in the piece cut out, called resp I want to  log the distance
>> between the two of them.
>> I tried
>>
>> function [trl, event] = response_trialfun(cfg);
>>
>> % read the header information and the events from the data
>> hdr   = ft_read_header(cfg.dataset);
>> event = ft_read_event(cfg.dataset);
>>
>>
>> % search for "trigger" events
>> value  = [event(find(strcmp('pres', {event.value}))).value]';
>> sample = [event(find(strcmp('pres', {event.value}))).sample]';
>>
>> % search for "trigger" events
>> value2  = [event(find(strcmp('resp', {event.value}))).value]';
>> sample2 = [event(find(strcmp('resp', {event.value}))).sample]';
>>
>> % determine the number of samples before and after the trigger
>> pretrig  = -round(cfg.trialdef.prestim  * hdr.Fs);
>> posttrig =  round(cfg.trialdef.poststim * hdr.Fs);
>>
>> trl = [];
>>
>>     trlbegin = sample + pretrig
>>     trlend   = sample + posttrig
>>     offset   = sample - sample + pretrig
>>     RT       = sample2 - sample
>>
>>     newtrl   = [trlbegin trlend offset RT];
>>     trl      = [trl; newtrl];
>>
>> end
>>
>>
>> But the problem is, that not always there is this second event, so
>> length(sample1) is not equal to length (sample2). Now I would just like
>> that every time when event 2 is missing sample 2 is just set to zero, so
>> that the two vectors have the same length and the function is working. Is
>> this possible?
>> Thanks a million
>> Katrin
>>
>> _______________________________________________
>> fieldtrip mailing list
>> fieldtrip at donders.ru.nl
>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>>
>
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20141210/1c53ab9d/attachment.html>


More information about the fieldtrip mailing list