event_id=mapp_bva, tmin=-3, tmax=3, add_eeg_ref=False) epochs_perhead_bva = mne.Epochs(raw_perhead_bva, mne_events_bva, event_id=mapp_bva, tmin=-3, tmax=3, add_eeg_ref=False) epochs_bip_bva['onsets_500_1500', 'stops_500_1500'].plot(block=True, scalings='auto') freqs = np.arange(2, 30, 1) n_cycles = freqs / 2 picks_original = mnehelp.def_picks(epochs_bip_bva) box = mnehelp.custom_box_layout(picks_original, 8) plot_picks_perhead = range(0, len(picks_original)) power_onset_perhead_bva = tfr_morlet(epochs_perhead_bva['onsets_500_1500'], freqs=freqs, n_cycles=n_cycles, picks=picks_original, return_itc=False) # NEED to pass picks because default IGNORES SEEG channels power_onset_perhead_bva.plot_topo(picks=plot_picks_perhead, baseline=(-2., -1.5), mode='logratio', layout=box)
pick_perhead_hip_names = mne.pick_info(raw_perhead_vr.info, pick_perhead_hip)['ch_names'] pick_perhead_ins = mnehelp.picks_all_localised(raw_perhead_vr, pd_montage_referenced, 'Ins') pick_perhead_all = mnehelp.picks_all(epochs_perhead_vr) # BAD EPOCHS # epochs_perhead_vr.plot(block = True, scalings = 'auto') # mnehelp.get_dropped_epoch_indices(epochs_perhead_vr.drop_log) bad_epochs = [] epochs_perhead_vr.drop(bad_epochs) # TIME FREQ freqs = np.arange(1, 11, 1) n_cycles = 6 box = mnehelp.custom_box_layout(pick_perhead_hip_names, 3) plot_pick_perhead_hip = range(len(pick_perhead_hip)) runfile( 'M:/Vyzkum/AV/FGU/IntracranialElectrodes/iEEG-python/tfr_perhead_unity.py', wdir='M:/Vyzkum/AV/FGU/IntracranialElectrodes/iEEG-python') # BASELINES ---------------- baseline = (-3, -2) runfile('M:/Vyzkum/AV/FGU/IntracranialElectrodes/iEEG-python/baselines.py', wdir='M:/Vyzkum/AV/FGU/IntracranialElectrodes/iEEG-python') # LFO BANDS ---------------- lfo_bands = [[2, 4], [4, 9]] runfile('M:/Vyzkum/AV/FGU/IntracranialElectrodes/iEEG-python/lfo_collapse.py', wdir='M:/Vyzkum/AV/FGU/IntracranialElectrodes/iEEG-python')
eeg.plot_psd(fmax=100, picks=pick_hip, average=False) # TIME FREQ freqs = np.arange(2, 12, 0.5) n_cycles = 6 events = [ 'onsets_500_1500', 'stops_500_1500', 'pointingStarted_Ego', 'pointingEnded_Ego', 'pointingStarted_Allo', 'pointingEnded_Allo' ] morlet = mneanalysis.morlet_all_events(epochs, freqs, n_cycles, events=['onsets_500_1500']) box = mnehelp.custom_box_layout(pick_all_names, 4) box = mne.channels.layout.make_grid_layout(eeg.info, pick=pick_all) morlet['onsets_500_1500'].average().plot_topo(picks=pick_all, layout=box, mode='logratio', baseline=(-1, -0.5)) morlet['onsets_500_1500'].average().plot(110, baseline=(-1, -0.5), mode='zlogratio') baseline = (-1, -0.5) mode = 'ratio' lfo_bands = [[2, 4], [4, 9]] morlet = mneanalysis.convolutions_apply_baselines(morlet, baseline, mode) morlet = mneanalysis.convolutions_band_power(morlet, lfo_bands)