Homec4science
Diffusion myhug (master) Browse

myh/dataProcessing/bci/bci_adrian/utils/get_pointerFiles_and_subject2flip.m
ab8049f9d919master

Last Change
Raw File

get_pointerFiles_and_subject2flip.m
View Options

%% All patients anaylses: Get subject2flip and pointer file
% comment: (for all subjects)chose those who completed the trial> those who
% has right lesion so to flip.
% Suject2flip: Order of id in pop file, flip right lesion to left
tableCompletedLesion = allpatients(allpatients.subject_id ~= 0,:);
tableCompletedRightLesionId = find(tableCompletedLesion.subject_id ~= 0 & tableCompletedLesion.lesion_side == 'right'); % table indexes of completed sham patients with right lesion
subject2flip = tableCompletedRightLesionId; % patient2flip = [13, 14, 15, 17, 24, 25 ] % index dans pointer ou pop file
% Test
table1 = tableCompletedLesion(subject2flip,:);
assert(all(table1.lesion_side == 'right'), 'All patients analyses, not all subject2flip patient are right lesioned');
subject = 1:size(tableCompletedLesion, 1);
subjectLeft = setdiff(subject, subject2flip);
table2 = tableCompletedLesion(subjectLeft,:);
assert(all(table2.lesion_side == 'left'), 'All patients analyses, not all remained subjects patient are left lesioned');
% Save subject2flip
save Y:\dataProcessing\bci\analyses\FC_analyses\allPatientsAnalyses\subject2flip subject2flip
clearvars table1 table2 subject2flip
% Pointer file
% comment: all the patients who completed the trial.
completedPatientId_index = find(allpatients.subject_id ~= 0); % table indexes of completed patients, the last one is the venitian patient
% completed_patientId = allpatients.patient_id(completedPatientId_index);
% Test
str = plist; expression = '[1-9]';
matches = regexp(str, expression, 'match');
indexes = [];
for i=1:length(matches)
indexes = [indexes; str2num(strjoin(matches{i}, ''))];
end
assert(isequal(indexes, completed_patientId), 'plist index and patient id do not match for voi saving')
% Voi saving
voi = nut_voinew( completedPatientId_index , 1, plist, 's_beamtf_ccohere_Spre_band1to20HzZ' )
save Y:\dataProcessing\bci\analyses\FC_analyses\allPatientsAnalyses\pointer_ccohere_Spre_band1to20HzZ voi
voi = nut_voinew( completedPatientId_index, 2, plist, 's_beamtf_ccohere_Spost_band1to20HzZ' )
save Y:\dataProcessing\bci\analyses\FC_analyses\allPatientsAnalyses\pointer_ccohere_Spost_band1to20HzZ voi
voi = nut_voinew( completedPatientId_index, 3, plist, 's_beamtf_ccohere_diff_band1to20HzZ' )
save Y:\dataProcessing\bci\analyses\FC_analyses\allPatientsAnalyses\pointer_ccohere_diff_band1to20HzZ voi
%% SHAM BCI group
tableShamLesion = allpatients(allpatients.subject_id ~= 0 & allpatients.treatment == 'sham',:); % table indexes of completed sham patients with right lesion
tableShamRightLesionId = find(tableShamLesion.subject_id ~= 0 & tableShamLesion.treatment == 'sham' & tableShamLesion.lesion_side == 'right'); % table indexes of completed sham patients with right lesion
subject2flip = tableShamRightLesionId; % patient2flip = [13, 14, 15, 17, (20), 24, 25 ] % index dans pointer ou pop file
% # of permutations : 32
% Test
table1 = tableShamLesion(subject2flip,:);
assert(all(table1.lesion_side == 'right'), 'All patients analyses, not all subject2flip patient are right lesioned');
subject = 1:size(tableShamLesion, 1);
subjectLeft = setdiff(subject, subject2flip);
table2 = tableShamLesion(subjectLeft,:)
assert(all(table2.lesion_side == 'left'), 'All patients analyses, not all remained subjects patient are left lesioned');
save Y:\dataProcessing\bci\analyses\FC_analyses\shamAnalyses\subject2flip subject2flip
clearvars table1 table2 subject2flip
% Pointer file
shamPatientId_index = find(allpatients.subject_id ~= 0 & allpatients.treatment == 'sham') % table indexes of completed patients, the last one is the venitian patient
shamBci_patientId = allpatients.patient_id(shamPatientId_index)
% Test
str = plist_bciSham; expression = '[1-9]';
matches = regexp(str, expression, 'match');
indexes = [];
for i=1:length(matches)
indexes = [indexes; str2num(strjoin(matches{i}, ''))];
end
assert(isequal(indexes, shamBci_patientId), 'plist_bciSham index and patient id do not match for voi saving')
% Voi saving
voi = nut_voinew( shamBci_patientId, 1, plist_bciSham,'s_beamtf_ccohere_Spre_band1to20HzZ')
save Y:\dataProcessing\bci\analyses\FC_analyses\shamAnalyses\pointer_bciSham_ccohere_Spre_band1to20HzZ voi
voi = nut_voinew( shamBci_patientId, 1, plist_bciSham,'s_beamtf_ccohere_Spost_band1to20HzZ')
save Y:\dataProcessing\bci\analyses\FC_analyses\shamAnalyses\pointer_bciSham_ccohere_Spost_band1to20HzZ voi
voi = nut_voinew( shamBci_patientId, 1, plist_bciSham,'s_beamtf_ccohere_diff_band1to20HzZ')
save Y:\dataProcessing\bci\analyses\FC_analyses\shamAnalyses\pointer_bciSham_ccohere_diff_band1to20HzZ voi
%% BCI group
tableBciLesion = allpatients(allpatients.subject_id ~= 0 & allpatients.treatment == 'bci',:); % table indexes of completed sham patients with right lesion
tableRealBciRightLesionId = find(tableBciLesion.subject_id ~= 0 & tableBciLesion.treatment == 'bci' & tableBciLesion.lesion_side == 'right'); % table indexes of completed bci patients with right lesion
subject2flip = tableRealBciRightLesionId; % patient2flip = [3, 5, 19, 22, 23, 26, 1] % index dans pointer ou pop file
% not enough subject for permutation
% Test
table1 = tableBciLesion(subject2flip,:);
assert(all(table1.lesion_side == 'right'), 'All patients analyses, not all subject2flip patient are right lesioned');
subject = 1:size(tableBciLesion, 1);
subjectLeft = setdiff(subject, subject2flip);
table2 = tableBciLesion(subjectLeft,:)
assert(all(table2.lesion_side == 'left'), 'All patients analyses, not all remained subjects patient are left lesioned');
save Y:\dataProcessing\bci\analyses\FC_analyses\bciAnalyses\subject2flip subject2flip
clearvars table1 table2 subject2flip
% Pointer file
realBcipatientId_index = find(allpatients.subject_id ~= 0 & allpatients.treatment == 'bci') % table indexes of completed patients, the last one is the venitian patient
realBci_patientId = allpatients.patient_id(realBcipatientId_index);
% Test
str = plist_bciReal; expression = '[1-9]';
matches = regexp(str, expression, 'match');
indexes = [];
for i=1:length(matches)
indexes = [indexes; str2num(strjoin(matches{i}, ''))];
end
assert(isequal(indexes, realBci_patientId), 'plist_bciReal index and patient id do not match')
voi = nut_voinew( realBci_patientId, 2, plist_bciReal,'s_beamtf_ccohere_Spre_band1to20HzZ')
save Y:\dataProcessing\bci\analyses\FC_analyses\bciAnalyses\pointer_bciReal_ccohere_Spre_band1to20HzZ voi % Caution : last patient is from Venice
voi = nut_voinew( realBci_patientId, 2, plist_bciReal,'s_beamtf_ccohere_Spost_band1to20HzZ')
save Y:\dataProcessing\bci\analyses\FC_analyses\bciAnalyses\pointer_bciReal_ccohere_Spost_band1to20HzZ voi % Caution : last patient is from Venice
voi = nut_voinew( realBci_patientId, 2, plist_bciReal,'s_beamtf_ccohere_diff_band1to20HzZ')
save Y:\dataProcessing\bci\analyses\FC_analyses\bciAnalyses\pointer_bciReal_ccohere_diff_band1to20HzZ voi % Caution : last patient is from Venice
c4science · Help