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)
