update demo files

iBMLab · Jan 22, 2016 · 9855d75 · 9855d75
1 parent 9a28200
commit 9855d75
Show file tree

Hide file tree

Showing 128 changed files with 450 additions and 48 deletions.
diff --git a/Data_Sample_with_codes/Data_sample/back.jpg → ...mple_with_codes/Data_sample_DEMO/back.jpg b/Data_Sample_with_codes/Data_sample/back.jpg → ...mple_with_codes/Data_sample_DEMO/back.jpg
diff --git a/...le_with_codes/Data_sample/data_subset.txt → ...th_codes/Data_sample_DEMO/data_subset.txt b/...le_with_codes/Data_sample/data_subset.txt → ...th_codes/Data_sample_DEMO/data_subset.txt
diff --git a/Data_Sample_with_codes/Data_sample/front.jpg → ...ple_with_codes/Data_sample_DEMO/front.jpg b/Data_Sample_with_codes/Data_sample/front.jpg → ...ple_with_codes/Data_sample_DEMO/front.jpg
diff --git a/Data_Sample_with_codes/Data_sample/label.txt → ...ple_with_codes/Data_sample_DEMO/label.txt b/Data_Sample_with_codes/Data_sample/label.txt → ...ple_with_codes/Data_sample_DEMO/label.txt
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data1.mat b/Data_Sample_with_codes/Data_sample_blindspot/data1.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data10.mat b/Data_Sample_with_codes/Data_sample_blindspot/data10.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data100.mat b/Data_Sample_with_codes/Data_sample_blindspot/data100.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data101.mat b/Data_Sample_with_codes/Data_sample_blindspot/data101.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data102.mat b/Data_Sample_with_codes/Data_sample_blindspot/data102.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data103.mat b/Data_Sample_with_codes/Data_sample_blindspot/data103.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data104.mat b/Data_Sample_with_codes/Data_sample_blindspot/data104.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data105.mat b/Data_Sample_with_codes/Data_sample_blindspot/data105.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data106.mat b/Data_Sample_with_codes/Data_sample_blindspot/data106.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data107.mat b/Data_Sample_with_codes/Data_sample_blindspot/data107.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data108.mat b/Data_Sample_with_codes/Data_sample_blindspot/data108.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data109.mat b/Data_Sample_with_codes/Data_sample_blindspot/data109.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data11.mat b/Data_Sample_with_codes/Data_sample_blindspot/data11.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data110.mat b/Data_Sample_with_codes/Data_sample_blindspot/data110.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data111.mat b/Data_Sample_with_codes/Data_sample_blindspot/data111.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data112.mat b/Data_Sample_with_codes/Data_sample_blindspot/data112.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data113.mat b/Data_Sample_with_codes/Data_sample_blindspot/data113.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data114.mat b/Data_Sample_with_codes/Data_sample_blindspot/data114.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data115.mat b/Data_Sample_with_codes/Data_sample_blindspot/data115.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data116.mat b/Data_Sample_with_codes/Data_sample_blindspot/data116.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data117.mat b/Data_Sample_with_codes/Data_sample_blindspot/data117.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data118.mat b/Data_Sample_with_codes/Data_sample_blindspot/data118.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data119.mat b/Data_Sample_with_codes/Data_sample_blindspot/data119.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data12.mat b/Data_Sample_with_codes/Data_sample_blindspot/data12.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data120.mat b/Data_Sample_with_codes/Data_sample_blindspot/data120.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data13.mat b/Data_Sample_with_codes/Data_sample_blindspot/data13.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data14.mat b/Data_Sample_with_codes/Data_sample_blindspot/data14.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data15.mat b/Data_Sample_with_codes/Data_sample_blindspot/data15.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data16.mat b/Data_Sample_with_codes/Data_sample_blindspot/data16.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data17.mat b/Data_Sample_with_codes/Data_sample_blindspot/data17.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data18.mat b/Data_Sample_with_codes/Data_sample_blindspot/data18.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data19.mat b/Data_Sample_with_codes/Data_sample_blindspot/data19.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data2.mat b/Data_Sample_with_codes/Data_sample_blindspot/data2.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data20.mat b/Data_Sample_with_codes/Data_sample_blindspot/data20.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data21.mat b/Data_Sample_with_codes/Data_sample_blindspot/data21.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data22.mat b/Data_Sample_with_codes/Data_sample_blindspot/data22.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data23.mat b/Data_Sample_with_codes/Data_sample_blindspot/data23.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data24.mat b/Data_Sample_with_codes/Data_sample_blindspot/data24.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data25.mat b/Data_Sample_with_codes/Data_sample_blindspot/data25.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data26.mat b/Data_Sample_with_codes/Data_sample_blindspot/data26.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data27.mat b/Data_Sample_with_codes/Data_sample_blindspot/data27.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data28.mat b/Data_Sample_with_codes/Data_sample_blindspot/data28.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data29.mat b/Data_Sample_with_codes/Data_sample_blindspot/data29.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data3.mat b/Data_Sample_with_codes/Data_sample_blindspot/data3.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data30.mat b/Data_Sample_with_codes/Data_sample_blindspot/data30.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data31.mat b/Data_Sample_with_codes/Data_sample_blindspot/data31.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data32.mat b/Data_Sample_with_codes/Data_sample_blindspot/data32.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data33.mat b/Data_Sample_with_codes/Data_sample_blindspot/data33.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data34.mat b/Data_Sample_with_codes/Data_sample_blindspot/data34.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data35.mat b/Data_Sample_with_codes/Data_sample_blindspot/data35.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data36.mat b/Data_Sample_with_codes/Data_sample_blindspot/data36.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data37.mat b/Data_Sample_with_codes/Data_sample_blindspot/data37.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data38.mat b/Data_Sample_with_codes/Data_sample_blindspot/data38.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data39.mat b/Data_Sample_with_codes/Data_sample_blindspot/data39.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data4.mat b/Data_Sample_with_codes/Data_sample_blindspot/data4.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data40.mat b/Data_Sample_with_codes/Data_sample_blindspot/data40.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data41.mat b/Data_Sample_with_codes/Data_sample_blindspot/data41.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data42.mat b/Data_Sample_with_codes/Data_sample_blindspot/data42.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data43.mat b/Data_Sample_with_codes/Data_sample_blindspot/data43.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data44.mat b/Data_Sample_with_codes/Data_sample_blindspot/data44.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data45.mat b/Data_Sample_with_codes/Data_sample_blindspot/data45.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data46.mat b/Data_Sample_with_codes/Data_sample_blindspot/data46.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data47.mat b/Data_Sample_with_codes/Data_sample_blindspot/data47.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data48.mat b/Data_Sample_with_codes/Data_sample_blindspot/data48.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data49.mat b/Data_Sample_with_codes/Data_sample_blindspot/data49.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data5.mat b/Data_Sample_with_codes/Data_sample_blindspot/data5.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data50.mat b/Data_Sample_with_codes/Data_sample_blindspot/data50.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data51.mat b/Data_Sample_with_codes/Data_sample_blindspot/data51.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data52.mat b/Data_Sample_with_codes/Data_sample_blindspot/data52.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data53.mat b/Data_Sample_with_codes/Data_sample_blindspot/data53.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data54.mat b/Data_Sample_with_codes/Data_sample_blindspot/data54.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data55.mat b/Data_Sample_with_codes/Data_sample_blindspot/data55.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data56.mat b/Data_Sample_with_codes/Data_sample_blindspot/data56.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data57.mat b/Data_Sample_with_codes/Data_sample_blindspot/data57.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data58.mat b/Data_Sample_with_codes/Data_sample_blindspot/data58.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data59.mat b/Data_Sample_with_codes/Data_sample_blindspot/data59.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data6.mat b/Data_Sample_with_codes/Data_sample_blindspot/data6.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data60.mat b/Data_Sample_with_codes/Data_sample_blindspot/data60.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data61.mat b/Data_Sample_with_codes/Data_sample_blindspot/data61.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data62.mat b/Data_Sample_with_codes/Data_sample_blindspot/data62.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data63.mat b/Data_Sample_with_codes/Data_sample_blindspot/data63.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data64.mat b/Data_Sample_with_codes/Data_sample_blindspot/data64.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data65.mat b/Data_Sample_with_codes/Data_sample_blindspot/data65.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data66.mat b/Data_Sample_with_codes/Data_sample_blindspot/data66.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data67.mat b/Data_Sample_with_codes/Data_sample_blindspot/data67.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data68.mat b/Data_Sample_with_codes/Data_sample_blindspot/data68.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data69.mat b/Data_Sample_with_codes/Data_sample_blindspot/data69.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data7.mat b/Data_Sample_with_codes/Data_sample_blindspot/data7.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data70.mat b/Data_Sample_with_codes/Data_sample_blindspot/data70.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data71.mat b/Data_Sample_with_codes/Data_sample_blindspot/data71.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data72.mat b/Data_Sample_with_codes/Data_sample_blindspot/data72.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data73.mat b/Data_Sample_with_codes/Data_sample_blindspot/data73.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data74.mat b/Data_Sample_with_codes/Data_sample_blindspot/data74.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data75.mat b/Data_Sample_with_codes/Data_sample_blindspot/data75.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data76.mat b/Data_Sample_with_codes/Data_sample_blindspot/data76.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data77.mat b/Data_Sample_with_codes/Data_sample_blindspot/data77.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data78.mat b/Data_Sample_with_codes/Data_sample_blindspot/data78.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data79.mat b/Data_Sample_with_codes/Data_sample_blindspot/data79.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data8.mat b/Data_Sample_with_codes/Data_sample_blindspot/data8.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data80.mat b/Data_Sample_with_codes/Data_sample_blindspot/data80.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data81.mat b/Data_Sample_with_codes/Data_sample_blindspot/data81.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data82.mat b/Data_Sample_with_codes/Data_sample_blindspot/data82.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data83.mat b/Data_Sample_with_codes/Data_sample_blindspot/data83.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data84.mat b/Data_Sample_with_codes/Data_sample_blindspot/data84.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data85.mat b/Data_Sample_with_codes/Data_sample_blindspot/data85.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data86.mat b/Data_Sample_with_codes/Data_sample_blindspot/data86.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data87.mat b/Data_Sample_with_codes/Data_sample_blindspot/data87.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data88.mat b/Data_Sample_with_codes/Data_sample_blindspot/data88.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data89.mat b/Data_Sample_with_codes/Data_sample_blindspot/data89.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data9.mat b/Data_Sample_with_codes/Data_sample_blindspot/data9.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data90.mat b/Data_Sample_with_codes/Data_sample_blindspot/data90.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data91.mat b/Data_Sample_with_codes/Data_sample_blindspot/data91.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data92.mat b/Data_Sample_with_codes/Data_sample_blindspot/data92.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data93.mat b/Data_Sample_with_codes/Data_sample_blindspot/data93.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data94.mat b/Data_Sample_with_codes/Data_sample_blindspot/data94.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data95.mat b/Data_Sample_with_codes/Data_sample_blindspot/data95.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data96.mat b/Data_Sample_with_codes/Data_sample_blindspot/data96.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data97.mat b/Data_Sample_with_codes/Data_sample_blindspot/data97.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data98.mat b/Data_Sample_with_codes/Data_sample_blindspot/data98.mat
diff --git a/Data_Sample_with_codes/Data_sample_blindspot/data99.mat b/Data_Sample_with_codes/Data_sample_blindspot/data99.mat
diff --git a/Data_Sample_with_codes/imap_DEMO.m b/Data_Sample_with_codes/imap_DEMO.m
@@ -106,22 +106,15 @@
                     if isempty(indx1)
                         indx(ii)=1;
                     else
-                    indxtf=sub2ind([ySize, xSize],coordX(indx1),coordY(indx1)); % index each fixation location,
-
-                    unindx = unique(indxtf);% find unique fixation
-                    [cotind,whe] = histc(indxtf,unindx); % cumulate fixation with same coordinates.
-                    durind=zeros(size(cotind));
-                    for iw=1:length(cotind);durind(iw)=sum(intv(whe==iw));end % calculate duration
-                    rawmap(unindx(durind>0))=durind(durind>0);
-
-                    f_mat = fft2(rawmap); % 2D fourrier transform on the points matrix
-                    filtered_mat = f_mat .* f_fil;
-
-                    smoothpic = real(fftshift(ifft2(filtered_mat)));
-                    isfixmap(it,:,:)=imresize(smoothpic,scale);
-                    israwmap(it,:,:)=imresize(rawmap,scale);
-                    stDur(it)=sum(durind);
-                    stRate(it)=srate(selected2(1));
+                        rawmap=full(sparse(coordX(indx1),coordY(indx1),intv(indx1),ySize,xSize));
+                        f_mat = fft2(rawmap); % 2D fourrier transform on the points matrix
+                        filtered_mat = f_mat .* f_fil;
+
+                        smoothpic = real(fftshift(ifft2(filtered_mat)));
+                        isfixmap(it,:,:)=imresize(smoothpic,scale);
+                        israwmap(it,:,:)=imresize(rawmap,scale);
+                        stDur(it)=sum(durind);
+                        stRate(it)=srate(selected2(1));
                     end
                 end
                 Subject{ii,1}=subjlist{is};
@@ -177,8 +170,8 @@
 %% LMMmap structure
 % impLMM will output a structure called LMMmap with all the information
 % from the linear mixed model fitting:
-% LMMmap = 
-% 
+% LMMmap =
+%
 %                    runopt: [1x1 struct]
 %              VariableInfo: [6x4 dataset]
 %                 Variables: [118x6 dataset]

diff --git a/Data_Sample_with_codes/imap_DEMO_ST1.m b/Data_Sample_with_codes/imap_DEMO_ST1.m
@@ -128,14 +128,10 @@
 gaussienne = (gaussienne - min(gaussienne(:))) / (max(gaussienne(:)) - min(gaussienne(:)));
 f_fil = fft2(gaussienne);
 % fixation matrix
-rawmap = zeros(ySize, xSize);
 coordX = round(Y(:,2));
 coordY = round(Y(:,1));
 indx1=coordX>0 & coordY>0 & coordX<xSize & coordY<ySize;
-indxtf=sub2ind([ySize, xSize],coordX(indx1),coordY(indx1)); % index each fixation location,
-unindx = unique(indxtf);% find unique fixation
-[cotind,whe] = histc(indxtf,unindx); % cumulate fixation with same coordinates.
-rawmap(unindx)=cotind;% number of fixation
+rawmap=full(sparse(coordX(indx1),coordY(indx1),ones(size(coordY(indx1))),ySize,xSize));
 f_mat = fft2(rawmap); % 2D fourrier transform on the points matrix
 filtered_mat = f_mat .* f_fil;
 smoothpic = real(fftshift(ifft2(filtered_mat)));
@@ -151,6 +147,7 @@
 RawMap=FixMap;
 Nfixbase=2; % number of fixations for each trials is the multiple of Nfixbase
 item=0;
+DataM=[];
 for is=1:Ns
     for itrial=1:n
         item=item+1;
@@ -165,14 +162,10 @@
         Ytmp2= [randi(xSize,10,1) randi(ySize,10,1)];% 10 random fixations 
         % fixation matrix
         Y=[Ytmp;Ytmp2];
-        rawmap = zeros(ySize, xSize);
         coordX = round(Y(:,2));
         coordY = round(Y(:,1));
         indx1=coordX>0 & coordY>0 & coordX<xSize & coordY<ySize;
-        indxtf=sub2ind([ySize, xSize],coordX(indx1),coordY(indx1)); % index each fixation location,
-        unindx = unique(indxtf);% find unique fixation
-        [cotind,whe] = histc(indxtf,unindx); % cumulate fixation with same coordinates.
-        rawmap(unindx)=cotind;% number of fixation
+        rawmap=full(sparse(coordX(indx1),coordY(indx1),ones(size(coordY(indx1))),ySize,xSize));
         f_mat = fft2(rawmap); % 2D fourrier transform on the points matrix
         filtered_mat = f_mat .* f_fil;
         smoothpic = real(fftshift(ifft2(filtered_mat)));
@@ -183,12 +176,20 @@
         Subject(item)=is;
         Rating(item)=Xall(is,itrial);
         fixNb(item)=Nfix;
+        Yk=cotind;
+        Yk(:,2:3)=unique([coordX(indx1),coordY(indx1)],'rows');
+        Yk(:,4)=Xall(is,itrial);
+        Yk(:,5)=is;
+        DataM=[DataM;Yk];
     end
 end
 % figure;imagesc(squeeze(mean(FixMap,1)))
 Mask=squeeze(mean(FixMap,1)>.25);
 % figure;imagesc(Mask)
 PredictorM=dataset(Rating,fixNb,Subject);
+ds=mat2dataset(DataM);
+ds=set(ds,'VarNames',{'fixN','rl','cl','rating','is'});
+export(ds,'File','STrate.csv','Delimiter',',');
 %% Linear Mixed Modeling with imapLMM
 tic
 opt.singlepredi=0;
@@ -221,4 +222,115 @@
 % perform multiple comparison correction
 [StatMap_c]=imapLMMmcc(StatMap,LMMmap,mccopt,FixMap);
 % output figure;
-imapLMMdisplay(StatMap_c,0)
+imapLMMdisplay(StatMap_c,0)
+%% Region dependent threshold
+ResampStat=StatMap_c.resampMat;
+nboot=mccopt.nboot;
+alpha=.05;
+Pmask=zeros(size(squeeze(StatMap_c.Pmask(2,:,:))));
+Pmap=StatMap.Pmap;
+mapvalue=StatMap.map;
+Maskall=zeros(size(Pmask));
+Maskall(16:135,16:135)=kron(reshape(1:16,4,4),ones(30,30));
+bootopt=1;
+
+for imask=1:16
+    region=Maskall==imask;
+    for ic=2
+        Fboot=squeeze(ResampStat.resFvalue(:,ic,region));
+        Pboot=squeeze(ResampStat.resPvalue(:,ic,region));
+        Betabt=squeeze(ResampStat.resBeta(:,ic,region));
+        Fboot=reshape(Fboot,1000,30,30);
+        Pboot=reshape(Pboot,1000,30,30);
+        Betabt=reshape(Betabt,1000,30,30);
+        % output estimated bootstrap distribution (cluster test)
+        maxclust=zeros(nboot,3);
+        for ib=1:nboot
+            Ftmp=squeeze(Fboot(ib,:,:));
+            [bmasktmp,bnum]=bwlabel(squeeze(Pboot(ib,:,:))<alpha);
+            if bnum>0
+                maxtmp=zeros(bnum,3);
+                for icluster=1:bnum
+                    maxtmp(icluster,1)=nansum(nansum(Ftmp(bmasktmp==icluster)));
+                    maxtmp(icluster,2)=sum(sum(bmasktmp==icluster));
+                    maxtmp(icluster,3)=maxtmp(icluster,1)./maxtmp(icluster,2);
+                end
+                maxclust(ib,1)=max(maxtmp(:,1));
+                maxclust(ib,2)=max(maxtmp(:,2));
+                maxclust(ib,3)=max(maxtmp(:,3));
+            end
+        end
+        distmp=sort(maxclust,1);
+        % output cluster mass and size threshold at alpha
+        clthres=distmp(round(nboot*(1-alpha)),:);
+        % output new Pmask: 1 cluster mass, 2 cluster size, 3 both,
+        % 4 cluster dense
+        mapvalue1=reshape(squeeze(mapvalue(ic,region)),30,30);
+        Pmap1=reshape(squeeze(Pmap(ic,region)),30,30);
+        switch bootopt
+            case 1
+                Pmask(region)=clustertest2D(mapvalue1, ...
+                    Pmap1,alpha,clthres(1),[],[]);
+            case 2
+                Pmask(region)=clustertest2D(mapvalue1, ...
+                    Pmap1,alpha,[],clthres(2),[]);
+            case 3
+                Pmask(region)=clustertest2D(mapvalue1, ...
+                    Pmap1,alpha,clthres(1),clthres(2),[]);
+            case 4
+                Pmask(region)=clustertest2D(mapvalue1, ...
+                    Pmap1,alpha,[],[],clthres(3));
+        end
+    end
+end
+figure;imagesc(Pmask);axis square off
+%% plot fixed effec(anova and beta) - using permutation
+% close all
+opt=struct;% clear structure
+opt.type='model beta';
+opt.alpha=.05;
+% perform contrast
+[StatMap]=imapLMMcontrast(LMMmap,opt); 
+% imapLMMdisplay(StatMap,1)
+
+mccopt=struct;
+mccopt.methods='permutation';
+mccopt.nboot=1000;
+mccopt.bootopt=1;
+mccopt.tfce=0;
+% perform multiple comparison correction
+[StatMap_c]=imapLMMmcc(StatMap,LMMmap,mccopt,FixMap);
+% output figure;
+imapLMMdisplay(StatMap_c,0)
+%% Region dependent threshold
+ResampStat=StatMap_c.resampMat;
+nboot=mccopt.nboot;
+alpha=.05;
+Pmask=zeros(size(squeeze(StatMap_c.Pmask(2,:,:))));
+Pmap=StatMap.Pmap;
+mapvalue=StatMap.map;
+Maskall=zeros(size(Pmask));
+Maskall(16:135,16:135)=kron(reshape(1:16,4,4),ones(30,30));
+bootopt=1;
+Pmapnew2=Maskall;
+for imask=1:16
+    region=Maskall==imask;
+    for ic=2
+        Fboot=squeeze(ResampStat.resFvalue(:,ic,region));
+        Fboot=reshape(Fboot,1000,30,30);
+        % output new Pmask
+        origFmap=squeeze(StatMap.map(ic,region));
+        origFmap=reshape(origFmap,30,30);
+        orFmat=permute(repmat(origFmap,[1,1,nboot]),[3,1,2]);
+        Fboot2=repmat(max(Fboot(:,:),[],2),[1,size(Fboot,2),size(Fboot,3)]);
+        Pmapnew2(region)=(sum(Fboot2>=orFmat,1)+1)./(nboot+1);
+    end
+end
+Pmapnew2(Mask==0)=1;
+figure;
+subplot(1,3,1);
+imagesc(squeeze(StatMap_c.Pmap(ic,:,:)));axis square off
+subplot(1,3,2);
+imagesc(Pmapnew2);axis square off
+subplot(1,3,3);
+imagesc(Pmapnew2<alpha);axis square off
diff --git a/Data_Sample_with_codes/imap_DEMO_ST2.m b/Data_Sample_with_codes/imap_DEMO_ST2.m
@@ -53,20 +53,11 @@
 gaussienne = (gaussienne - min(gaussienne(:))) / (max(gaussienne(:)) - min(gaussienne(:)));
 f_fil = fft2(gaussienne);
 % fixation matrix
-rawmap = zeros(ySize, xSize);
 coordX = round(Y(:,2));
 coordY = round(Y(:,1));
 intv=normrnd(0.4,.085,length(Y),1);
 indx1=coordX>0 & coordY>0 & coordX<xSize & coordY<ySize;
-indxtf=sub2ind([ySize, xSize],coordX(indx1),coordY(indx1)); % index each fixation location,
-unindx = unique(indxtf);% find unique fixation
-[cotind,whe] = histc(indxtf,unindx); % cumulate fixation with same coordinates.
-% number of fixation
-% rawmap(unindx)=cotind;
-% fixation duration
-durind=zeros(size(cotind));
-for iw=1:length(cotind);durind(iw)=sum(intv(whe==iw));end % calculate duration
-rawmap(unindx(durind>0))=durind(durind>0);
+rawmap=full(sparse(coordX(indx1),coordY(indx1),intv(indx1),ySize,xSize));
 
 f_mat = fft2(rawmap); % 2D fourrier transform on the points matrix
 filtered_mat = f_mat .* f_fil;
@@ -160,15 +151,7 @@
             pathlength=diag(squareform(pdist([coordY,coordX])),1);
             intv=normrnd(Meandur,Stddur,length(Y),1)*1000;
             indx1=coordX>0 & coordY>0 & coordX<xSize & coordY<ySize;
-            indxtf=sub2ind([ySize, xSize],coordX(indx1),coordY(indx1)); % index each fixation location,
-            unindx = unique(indxtf);% find unique fixation
-            [cotind,whe] = histc(indxtf,unindx); % cumulate fixation with same coordinates.
-            % number of fixation
-            % rawmap(unindx)=cotind;
-            % fixation duration
-            durind=zeros(size(cotind));
-            for iw=1:length(cotind);durind(iw)=sum(intv(whe==iw));end % calculate duration
-            rawmap(unindx(durind>0))=durind(durind>0);
+            rawmap=full(sparse(coordX(indx1),coordY(indx1),intv(indx1),ySize,xSize));
 
             f_mat = fft2(rawmap); % 2D fourrier transform on the points matrix
             filtered_mat = f_mat .* f_fil;