12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394 |
- mainfolder='U:/LindeDomingoJuan/Memoreye_DATA/'
- partis=good_partis_sub;
- perNan=NaN(length(partis),2,2); %Preparing output variable
- %%
- for delay=1:2
-
- for ppp=1:length(partis)
-
- % Loading dataset
-
- if delay== 1
- toload =[mainfolder '/eyelink_preprocessed/mean_recentered_th100/object_1_onset_11_seconds/' partis{ppp,:} '.mat']; %for 1st delay
- elseif delay==2
- toload =[mainfolder '/eyelink_preprocessed/mean_recentered_th100/2nd_cue_onset_17_5__seconds/' partis{ppp,:} '.mat']; %for 2nd delay
- end
- load(toload)
-
- cx=squeeze(data_epoch_thr.eyedat(:,1,:)); %data of channel x
- cy=squeeze(data_epoch_thr.eyedat(:,2,:)); %data of channel y
-
- cxv=cx(:); %data as a vector
- cyv=cy(:); %data as a vector
-
-
- % Saving the percentage of NaN per epoch.
-
- perNan(ppp,1,delay)=(sum(isnan(cxv))/length(cxv))*100;
- perNan(ppp,2,delay)=(sum(isnan(cyv))/length(cyv))*100;
-
- end
- end
- %%
- baddelay1=excludingoutlier(perNan(:,1,1))
- baddelay2=excludingoutlier(perNan(:,1,2))
- index=(baddelay1+baddelay2)==0; %keeping only participants that are OK in both epochs
- %% Making a final list of good participants and saving the file
- counter=1;
- for ppp = 1:length(good_partis_sub)
-
- if index(ppp)==1
- good_partis{counter} =good_partis_sub{ppp}
- counter=counter+1;
- end
- end
- %% Plotting results
- if figure_NaNs
- figure;
- subplot(1,2,1)
- boxplot(squeeze(perNan(:,:,1)))
- ylabel('% of NaN')
- xlabel('channel (x and y)')
- title('delay 1 epoch')
- subplot(1,2,2)
- boxplot(squeeze(perNan(:,:,2)))
- ylabel('% of NaN')
- xlabel('channel (x and y)')
- title('delay 2epoch')
- end
- m1=['the max NaN % in delay 1 is ' num2str(max(perNan(index,:,1))) ' and the min is ' num2str(min(perNan(index,:,1)))];
- m2=['the max NaN % in delay 2 is ' num2str(max(perNan(index,:,2))) ' and the min is ' num2str(min(perNan(index,:,2)))];
- disp(m1)
- disp(m2)
- %%
- function output=excludingoutlier(data)
- %excluding participants with values >1.5*IQR in the third quartile
- ad=squeeze(data(:,1,1));
- r = iqr(ad);
- Q = quantile(ad,[0.75]);
- output=ad>(1.5*r)+Q;
- end
|