ERPs_AuditoryMotion/getTrialSeq.m at master · CerenB/ERPs_AuditoryMotion · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
function [trial_seq_names,trial_seq] = getTrialSeq(numEvents,numTargets,expLength)

% function to create the randomization of all the conditions in the ERP -
%motion experiemnt
%
% with numTrails and percTarget it calculates the number of trials per
% condition
%
% outputs the trial sequence with numbers (trial_seq) or with names (trial_seq_names)
%
% number assigned order:

% 1    rms_static_1s
% 2    rms_mot_LR_1s
% 3    rms_mot_RL_1s
% 4    rms_static_2s
% 5    rms_mot_LR_2s
% 6    rms_mot_RL_2s

expLength = str2num(expLength); %#ok<ST2NM>

switch expLength

    case 1

        nChunks = 3;

        numEvents = numEvents/nChunks;

        numTargets = numTargets/nChunks;

        trial_seq = pseudorand(numEvents, numTargets, nChunks, expLength);

    otherwise

        nChunks = 2;

        numEvents = numEvents/nChunks;

        numTargets = numTargets/nChunks;

        trial_seq = pseudorand(numEvents, numTargets, nChunks, expLength);

end

% rename the numbers with condition name
trial_seq_names = {};
for i = 1:length(trial_seq)
    trial_seq_names{end+1} = num2str(trial_seq(i)); %#ok<AGROW>
end

trial_seq_names(strcmpi(trial_seq_names,'1')) = {'rms_static_1s'};
trial_seq_names(strcmpi(trial_seq_names,'2')) = {'rms_mot_LR_1s'};
trial_seq_names(strcmpi(trial_seq_names,'3')) = {'rms_mot_RL_1s'};
trial_seq_names(strcmpi(trial_seq_names,'4')) = {'rms_static_2s'};
trial_seq_names(strcmpi(trial_seq_names,'5')) = {'rms_mot_LR_2s'};
trial_seq_names(strcmpi(trial_seq_names,'6')) = {'rms_mot_RL_2s'};

end


function trial_seq = pseudorand(numEvents, numTargets, nChunks, expLength)

% calculate how many trials per condition dividing them in 3 chunks to help
% even randomization across experiment per condition

static = (numEvents-numTargets)/2;
mot_LR = (numEvents-numTargets)/4;
mot_RL = (numEvents-numTargets)/4;
static_T = numTargets/2;
mot_LR_T = numTargets/4;
mot_RL_T = numTargets/4;

% create the three chunks of condition and randomize them and check that:
% 1 - the target is not in the first trial
% 2 - two target are not consecutive
% 3 - there are no more them 3 same trials consecutive (less than that is impossible)

d = 1; % counter for trial number

while d < numEvents*nChunks-2

    % create a sequence of trials that contains 1 thirs of all the trials
    trial_seq = [ ...
        repmat(ones,1,static),  repmat(2,1,mot_LR),   repmat(3,1,mot_RL), ...
        repmat(4,1,static_T),   repmat(5,1,mot_LR_T), repmat(6,1,mot_RL_T)];

    switch expLength

        case 1

            % we shuffle 3 times each chunck and concatenate them
            trial_seq = [ Shuffle(trial_seq), Shuffle(trial_seq), Shuffle(trial_seq) ];

        otherwise

            % we shuffle 2 times each chunck and concatenate them
            trial_seq = [ Shuffle(trial_seq), Shuffle(trial_seq) ];

    end

    % scan through the trial sequence and checks all the conditions
    while d < length(trial_seq)-2

        if d == 1 && trial_seq(d) > 3 % no targets in the first trial
            d = 1;
            break
        elseif trial_seq(d) > 3 && trial_seq(d+1) > 3 % avoid 2 consecutive targets
            d = 1;
            break
            % avoid series 3 times the same conditions
        elseif trial_seq(d) ==  trial_seq(d+1) && trial_seq(d+1) ==  trial_seq(d+2) && trial_seq(d+2) ==  trial_seq(d+3)
            d = 1;
            break
        else
            d = d+1;
        end
    end
end

end