#laptop wd wd = '/home/sammirc/Desktop/DPhil/wmConfidence' #workstation wd os.chdir(wd) subs = np.array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25 ]) subs = np.array( [4, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17, 18, 20, 21, 22, 24, 25]) subs = np.array([24, 25]) for i in subs: print('\n\nworking on subject ' + str(i) + '\n\n') sub = dict(loc='workstation', id=i) param = get_subject_info_wmConfidence(sub) if i <= 2: #these subjects only have one session #read in the already ica cleaned data, already been filtered too raw = mne.io.read_raw_fif(fname=param['rawcleaned'], preload=True) raw.set_montage('easycap-M1') #epoching #here it's important to specify a dictionary that assigns each trigger to its own integer value #mne.events_from_annotations will assign each trigger to an ordered integer, #so e.g. trig11 will be 2, but epoching 11 will include another trigger #this solves it event_id = { '1': 1, '2': 2, #array
threshold = -10 if threshold == -10: glmnum = 6 else: glmnum = 5 data = dict() data_t = dict() for i in contrasts: data[i] = [] data_t[i] = [] for i in subs: print('\n\ngetting subject ' + str(i) + '\n\n') sub = dict(loc='workstation', id=i) param = get_subject_info_wmConfidence(sub) #_baselined for name in contrasts: data[name].append( mne.read_evokeds(fname=op.join( param['path'], 'glms', 'feedback', 'epochs_glm' + str(glmnum), 'wmc_' + param['subid'] + '_feedbacklocked_tl_' + name + '_betas-ave.fif'))[0]) data_t[name].append( mne.read_evokeds(fname=op.join( param['path'], 'glms', 'feedback', 'epochs_glm' + str(glmnum), 'wmc_' + param['subid'] + '_feedbacklocked_tl_' + name + '_tstats-ave.fif'))[0]) #drop right mastoid from literally everything here lol its not useful anymore for cope in data.keys(): for i in range(subs.size):