get_popdiffFile.m
ab8049f9d919master

get_popdiffFile.m
View Options

This document is not UTF8. It was detected as ISO-8859-1 (Latin 1) and converted to UTF8 for display.

			%% Get the pop diff file from the pop Spre and pop Spost
			% From pop_ccohere_Spre.mat and pop_ccohere_Spost.mat
			%
			% Au moins 5 sujets car 1/32 = 0.03 (Valeur plus petite possible)
			% Si 4 sujets car 1/16 = 0.06 (déjà plus significatifs, même au plus petit)
			%
			% pop.s et beam.s : le s est la valeur de cohérence imaginaires dans les
			% voxels
			%
			% Pour faire le pop_diff (pas besoin de faire le s_beam diff)
			current_dir = pwd;

			cd Y:\dataProcessing\bci\analyses\FC_analyses\allPatientsAnalyses
			% cd Y:\dataProcessing\bci\analyses\FC_analyses\bciAnalyses
			% cd Y:\dataProcessing\bci\analyses\FC_analyses\shamAnalyses


			%% If pop r2l file are available
			load('pop_ccohere_Spre_band1to20HzZ_r2l.mat')
			post_r2l = pop

			load('pop_ccohere_Spre_band1to20HzZ_r2l.mat')
			pre_r2l = pop

			assert(isequal(post_r2l.voxels,pre_r2l.voxels), 'post_r2l pop and pre_r2l pop do not have the same number of voxels')
			pop.s = post_r2l.s - pre_r2l.s
			save pop_ccohere_diff_band1to20HzZ_r2l pop

			%% If pop r2l file are missing
			load subject2flip

			load('pop_ccohere_Spost_band1to20HzZ.mat')
			pop = nut_fliplr(pop,subject2flip)
			save pop_ccohere_Spost_band1to20HzZ_r2l pop
			post_r2l = pop;

			load('pop_ccohere_Spre_band1to20HzZ.mat')
			pop = nut_fliplr(pop,subject2flip)
			save pop_ccohere_Spre_band1to20HzZ_r2l pop
			pre_r2l = pop;

			isequal(post_r2l.voxels,pre_r2l.voxels)
			pop.s = post_r2l.s - pre_r2l.s
			save pop_ccohere_diff_band1to20HzZ_r2l pop

			cd(current_dir)
			clear current_dir pre_r2l post_r2l subject2flip pop