Imported at 00:00, 15 March 2026

2026-03-15T00:00:00Z

New page

== Experiment Settings ==
Not that the followings are default settings, and the actual details are in the each paper.

Stimuli: x=800 Hz, Y=1600 Hz, 50 ms duration No raise/fall time

Trial: ISI = 150 ms (thus SOA = 200 ms), ITI = [1.8 2.0] s

Block: 10 types, randomized order

1 block contains 100 trials (20 habituation, 64 standards, and 16 deviants).

Tone Index:

1. xxxxx

2. xxxxY

3. xxxx_

4. YYYYY

5. YYYYx

6. YYYY_

Block Index:

1. xxxxx with deviant xxxxY

2. xxxxx with omission xxxx_

3. xxxxY with deviant xxxxx

4. xxxxY with omission xxxx_

5. pure omission xxxx_

6. YYYYY with deviant YYYYx

7. YYYYY with omission YYYY_

8. YYYYx with deviant YYYYY

9. YYYYx with omission YYYY_

10. pure omission YYYY_

== Data format ==

Data can be download from Brain/MINDS data portal.

One zip file has one date data which consists of a session.

For Monkey Ji and Rc:

‘ECoG_ch**.mat‘ had ECoG data of electrode. ** means number of electrode. Sampling rate was 1 kHz.

‘ECoGTime.mat‘ had time of ECoG data.

‘Event.mat‘ had information of condition.

‘Eye.mat‘ had eye tracking data recorded from right eye.

‘Audio.mat‘ had a audio data.

For Monkey Yo, Ca and Rm:

‘ECoG_ch**.mat‘ had ECoG data of electrode. ** means number of electrode. Sampling rate was 1017.25 Hz.

‘TriG_ch2.mat‘ had information of condition.

‘Aud_ch1.mat‘ had a audio data.

== Sample code for Monkey Ji and Rc ==

clear

addpath(genpath('D:\MATLAB_toolboxes\eeglab14_1_1b'))
cd G:\Local_global_marmoset

Subject='Ji';
% Subject='Rc';

Expt_name=sprintf('%s20*',Subject);
Expt=ls(Expt_name);
Expt_num=size(Expt,1);

FS=1000; %Hz
DS=1;

Channel_num=96;

%%
Data_all=[];
Time=[];
Stim_time=[];
Block_idx=[];
Trial_idx=[];
Expt_idx=[];

Time_offset=0;

for expt=1:Expt_num

cd(Expt(expt,:))
load Event %allBlockIdx allTrialIdx StimOn

fprintf(1,'>> Expt%d/%d: %s (%.1f min)(%d trials)\n',expt, Expt_num,Expt(expt,:),length(EventTime)/1000/60,length(StimOn))

Data=zeros(Channel_num,length([1:DS:length(EventTime)]));

for ch=1:Channel_num
FN=sprintf('ECoG_ch%d',ch);
load(FN)

str=sprintf('Data(ch,:)=ECoGData_ch%d(1:DS:end);',ch);
eval(str)
end

Data_all=[Data_all Data];
Time=[Time EventTime(1:DS:end)/FS+Time_offset];
Stim_time=[Stim_time StimOn/FS+Time_offset];
Block_idx=[Block_idx allBlockIdx];
Trial_idx=[Trial_idx allTrialIdx];
Expt_idx=[Expt_idx ones(1,length(allTrialIdx))*expt];

Time_offset=Time(end)+DS/FS;

clear Data EventTime StimOn

cd ..

end

== Sample code for Monkey Yo, Ca and Rm ==

clear

addpath(genpath('D:\MATLAB_toolboxes\eeglab14_1_1b'))
cd G:\Local_global_marmoset

Subject='Ca';
% Subject='Rm';
% Subject='Yo';

Expt_name=sprintf('%s20*',Subject);
Expt=ls(Expt_name);
Expt_num=size(Expt,1);

FS=1017.25; %Hz
DS=1;

Channel_num=96;

Trial_per_block=100;

%% Block info

Block_type_name{1}='xxxxx with deviant xxxxY';
Block_type_name{2}='xxxxx with omission xxxx_';
Block_type_name{3}='xxxxY with deviant xxxxx';
Block_type_name{4}='xxxxY with omission xxxx_';

Anag_xx=40;
Anag_xY=80;
Anag_xo=120;

%%
Data_all=[];
Block_idx=[];
Trial_idx=[];
Expt_idx=[];
Stim_idx=[];

Idx_offset=0;

for expt=1:Expt_num

cd(Expt(expt,:))

%% Event
load('TriG_ch2.mat')

%trim Trig if longer than wave
load('ECoG_ch1');

if length(Trig)>length(wave)
Trig=Trig(1:length(wave));
end

%downsample
Trig=Trig(1:DS:end);

Onset=find(diff(Trig)>10);
Nope=find(Trig(Onset)>10);
Onset(Nope)=[];

Trial_idx_expt=round(Trig(Onset+10)/40); %1: xxxxx; 2: xxxxY; 3:xxxxo

Block_idx_expt=zeros(1,length(Onset));

%%% Could be only 3 blocks
Block_per_expt=ceil(length(Block_idx_expt)/Trial_per_block);

for block=1:Block_per_expt
if block<Block_per_expt
XX=Trial_idx_expt([1:Trial_per_block]+(block-1)*Trial_per_block);
else % the last block could be short
XX=Trial_idx_expt((1+(block-1)*Trial_per_block):end);
end
Stand=XX(1);
Dev=find(XX~=Stand);
Dev=XX(Dev(1));

if Stand==1 & Dev==2
B_type=1;
elseif Stand==1 & Dev==3
B_type=2;
elseif Stand==2 & Dev==1
B_type=3;
elseif Stand==2 & Dev==3
B_type=4;
else
B_type=-1;
end

if block<Block_per_expt
Block_idx_expt([1:Trial_per_block]+(block-1)*Trial_per_block)=B_type;
else
Block_idx_expt((1+(block-1)*Trial_per_block):end)=B_type;
end

end

Trial_idx=[Trial_idx Trial_idx_expt];
Block_idx=[Block_idx Block_idx_expt];
Expt_idx=[Expt_idx ones(1,length(Onset))*expt];

Stim_idx=[Stim_idx Onset+Idx_offset];

Idx_offset=Idx_offset+length(Trig);

%%
fprintf(1,'>> Expt %d/%d: %s (%d trials)\n',expt, Expt_num,Expt(expt,:),length(Onset))

Data=zeros(Channel_num,length(Trig));
for ch=1:Channel_num
FN=sprintf('ECoG_ch%d',ch);
load(FN)

Data(ch,:)=wave(1:DS:end);
end

Data_all=[Data_all Data];

clear Data wave
cd ..
end
Time=[0:size(Data_all,2)-1]/FS*DS;
Stim_time=Time(Stim_idx);

Memory-guided saccade Details - Revision history

Imported at 00:00, 15 March 2026