function sm_ttest() mne_setup_toolbox MEG_DIR='/space/ventzl/1/users/MEG_DIR/'; sub={'001' '002' '004' '005' '006' '007' '008' '009' '010' '011' '012' '013' '014' '015' '016' '017' '019' '020' '024'}; hemi={'lh' 'rh'}; target_sub='meg_normals'; analysis='pro_anti_gradonly'; cond={'Anti-1-integ20-meg_normals' 'Pro-1-integ20-meg_normals'}; % cond={'Anti-integ50_bin-meg_normals' 'Pro-integ50_bin-meg_normals'}; % cond={'Anti-integ500_bin-meg_normals' 'Pro-integ500_bin-meg_normals'}; suffix=''; %%----------------- n=length(sub); fprintf('\n'); for h=1:length(hemi) data1=[]; data2=[]; for s=1:length(sub) cmd=sprintf('fin=''%s/mano%s/%s/stc/mano%s_%s%s-%s.stc'';',MEG_DIR,sub{s},analysis,sub{s},cond{1},suffix,hemi{h}); eval(cmd); stc=mne_read_stc_file(fin); data1=cat(3,data1,stc.data); tmin=stc.tmin; tstep=stc.tstep; cmd=sprintf('fin=''%s/mano%s/%s/stc/mano%s_%s%s-%s.stc'';',MEG_DIR,sub{s},analysis,sub{s},cond{2},suffix,hemi{h}); eval(cmd); stc=mne_read_stc_file(fin); data2=cat(3,data2,stc.data); fprintf('%d/%d: mano%s\n',s,n,sub{s}); end fprintf('\n'); fprintf('\n'); fprintf('Computing %s...',hemi{h}); mean1=mean(data1,3); mean2=mean(data2,3); diff1=[]; diff2=[]; for s=1:length(sub) diff1=cat(3,diff1,data1(:,:,s)-mean1); diff2=cat(3,diff2,data2(:,:,s)-mean2); end diff12sq=(diff1-diff2).^2; sigma=sum(diff12sq,3); root=sqrt(sigma./(n*(n-1))); t=(mean1-mean2)./root; stc_out=stc; stc_out.data=t; stc_out.tmin=tmin; stc_out.tstep=tstep; cmd=sprintf('fout=''%s/%s/%s/stc/t-%s-%s%s-%s.stc'';',MEG_DIR,target_sub,analysis,cond{1},cond{2},suffix,hemi{h}); eval(cmd); mne_write_stc_file(fout, stc_out); fprintf(' done!\n'); fprintf('Saved to: %s\n\n',fout); end