Script04_Hypothesis2b.m

clc;
clear all;
close all;
%% Hypothesis 2b: novelty effect in theta power of frontocentral channels
for Subject=[1:33]
    Frontocentral_channels=[1:14 32:50 64:65];
    file_name=['Z:\projects\Hamid\Projects\EEGManyPipelines\Analyses\Subject',num2str(Subject,'%02.f'),...
        '\','Subject',num2str(Subject,'%02.f'),'_preprocessed_data.mat'];
    load(file_name,'Data_matrix_theta','events','channels')
    % Separating trials into classes and performing decoding
    tr1=0;
    tr2=0;
    for trial=1:length(events)
        if strcmp(events(trial).old,'new')
            tr1=tr1+1;
            ClassA(:,:,tr1)=Data_matrix_theta(Frontocentral_channels,:,trial);
        elseif strcmp(events(trial).old,'old')
            tr2=tr2+1;
            ClassB(:,:,tr2)=Data_matrix_theta(Frontocentral_channels,:,trial);
        end
    end
    time_sample=0;
    for time=1:size(Data_matrix_theta,2)
        time_sample=time_sample+1;
        ClassA_time=squeeze(ClassA(:,time,:));
        ClassB_time=squeeze(ClassB(:,time,:));
        Classifier_Model = fitcdiscr([ClassA_time';ClassB_time'],[ones(size(ClassA_time,2),1);zeros(size(ClassB_time,2),1)]);
        decoding_accuracy(time_sample)=1-kfoldLoss(crossval(Classifier_Model));
        [Subject time_sample]
    end
    file_name=['Z:\projects\Hamid\Projects\EEGManyPipelines\Analyses\Subject',num2str(Subject,'%02.f'),...
        '\','Subject',num2str(Subject,'%02.f'),'_Results_Hypothesis2b.mat'];
    save(file_name,'decoding_accuracy');
    clearvars -except Subject All_channels
end

 %% Plotting  and statistical testing
 % loading decoding results from all participants
 Accuracies=nan(33,769);
 for Subject=1:33
     file_name=['Z:\projects\Hamid\Projects\EEGManyPipelines\Analyses\Subject',num2str(Subject,'%02.f'),...
         '\','Subject',num2str(Subject,'%02.f'),'_Results_Hypothesis2b.mat'];
     load(file_name,'decoding_accuracy');
     Accuracies(Subject,:)=decoding_accuracy;
 end
 
% Evaluation of significance of effect against decoding in the baseline
% period
for time=1:size(Accuracies,2)
    [p(time),h(time)]=signrank(Accuracies(:,time),nanmean(Accuracies(:,1:256),2),'tail','right');
end

% Corrrecting p values for multiple comparisons using Combined Probability of Fisher
threshold=0.05;
[pCorrected,~, ~, ~] = mt_fisher(reshape(p,1,[]), threshold);
pnew=zeros(size(pCorrected));
pnew(pCorrected<threshold)=1;
time_span_of_clustering= 5; % 5 = 10 ms: only the time and frequency spans
% in which all points are sginificant will be marked as significant
pCorrected=bwareaopen(pnew,time_span_of_clustering);

% plotting decoding curves and indicating time points of significant effect
xticks=[-256:512]*1000./512;
yticks=[45:5:60];
smoothing=5; % smoothing the decoding curve of clearer visualtion

plott=shadedErrorBar(xticks,nanmean(Accuracies*100),nanstd(Accuracies*100)*1.96./sqrt(33),{'color','b','LineWidth',2},1);
hold on;
plot([xticks(1) xticks(end)],[50 50],'k');
plot([0 0],[yticks(1) yticks(end)],'k');
Sig_plot=plot(xticks,pCorrected*47.5,'*');
set(gca,'FontSize',10,'LineWidth',1,'XTick',...
    [xticks(1):250:xticks(end)],'XTickLabel',...
    [xticks(1):250:xticks(end)],'YTick',...
    [yticks],'YTickLabel',{[yticks]},...
    'XMinorTick','on','YMinorTick','off','ycolor','k','tickdir','out','xcolor','k','box','off');
xlim([xticks(1) xticks(end)]);
ylim([yticks(1) yticks(end)]);
legend([plott.mainLine,plott.patch Sig_plot],{'Mean: new vs. old';'95% CI';'Significant'},'location','northwest')
xlabel('Time relative to stimulus onset (ms)');
ylabel('Decoding accuracy (%)');
title ('H2b: novelty effect in Theta band in frontocentral channels');