|
5 | 5 | % so that you can decide how many microstates |
6 | 6 | % you want to choose without getting into single subjects |
7 | 7 | % only run 1 group at a time to compaire their EO and EC data |
8 | | -% adapted for our pipeline on 12/20/2021 by Douwe |
| 8 | +% or run for everyone together |
| 9 | +% adapted for our pipeline on 6/26/2021 by Douwe |
9 | 10 | clear variables |
10 | 11 |
|
11 | 12 |
|
12 | 13 |
|
13 | 14 |
|
14 | | -Group = {'22q' };% 'Control 22q' 'Control sz' '22q +' '22q -' 'sz' };% |
| 15 | +Group = {'All' };% 'Control 22q' 'Control sz' '22q +' '22q -' 'sz' };% |
15 | 16 | type={'EC'}; |
16 | 17 | for g=1:length(Group) |
17 | 18 | switch Group{g} |
|
35 | 36 | subject_list = {'2202' '2204' '2212' '2222' '2229' '2231' '2257' '2260' '2261' '2284' '2295'}; |
36 | 37 | case '22q' |
37 | 38 | subject_list = {'2201','2207','2216','2243','2256','2267','2274','2281','2286','2292','2202','2204','2212','2222','2229','2231','2257','2260','2261','2284','2295'}; |
| 39 | + case 'All' |
| 40 | + subject_list = {'10293' '10561' '10562' '10581' '10616' '10748' '10822' '10858' '10935' '12004' '12512' '12632' '12648' '12727' '12739' '12746' '12750' '12770' '12815' '12010' '12139' '12177' '12188' '12197' '12203' '12206' '12272' '12415' '12482' '12588' '12651' '12852' '12870' '7003' '7007' '7019' '7025' '7046' '7051' '7054' '7058' '7061' '7064' '7065' '7073' '7078' '7089' '7092' '7094' '7123' '7556' '7808' '2201' '2207' '2216' '2243' '2256' '2267' '2274' '2281' '2286' '2292' '2202' '2204' '2212' '2222' '2229' '2231' '2257' '2260','2261' '2284' '2295'}; |
| 41 | + |
38 | 42 | end |
39 | 43 | home_path = 'D:\restingstate\data\'; |
40 | 44 | for t=1:length(type) |
|
43 | 47 | for s=1:length(subject_list) |
44 | 48 | data_path = [home_path subject_list{s} ''];% Path to the folder containing the current subject's data |
45 | 49 | % Load original dataset |
46 | | - fprintf('\n\n\n**** %s: Loading dataset ****\n\n\n', subject_list{s}); |
| 50 | + fprintf('\n\n\n**** %s: Loading dataset ****\n\n\n', subject_list{s}); |
47 | 51 | EEG = pop_loadset('filename', [subject_list{s} '_' type{t} '.set'], 'filepath', data_path); |
48 | 52 | [ALLEEG, EEG] = eeg_store(ALLEEG, EEG, CURRENTSET); |
49 | 53 | end |
|
57 | 61 | figure;[tt]=title([type{t} ' ' Group{g}]);tt.FontSize = 15;MicroPlotTopo( EEG, 'plot_range', [] ); %plotting microstates |
58 | 62 | print([save_path Group{g} '_microstate_' type{t}], '-djpeg','-r300'); |
59 | 63 | figure(); |
60 | | - % EEG = pop_micro_selectNmicro( EEG ); % only select CV and GEV, look for where GEV doesn't increase significantly |
| 64 | + EEG = pop_micro_selectNmicro( EEG ); % only select CV and GEV, look for where GEV doesn't increase significantly |
61 | 65 | [ALLEEG EEG] = eeg_store(ALLEEG, EEG, CURRENTSET); |
62 | 66 | close all |
63 | 67 | end |
|
0 commit comments