[Mne_analysis] concatenation raws and events

[J.A. Entakli]

Hello Everybody,
During my experiment I recorded 2 blocs per subject. I am just started 
with mne_python.

I have two questions in one ::
1) I would like to know how can I concatenate the blocs and above all 
the events (corresponding to each bloc) to have just one bloc.
In brief, for each subject I have one bloc with its events and another 
bloc with its events also.

raw_fname1= '...'
raw_fname2= '...'
event_fname1= '...'
event_fname2= '...'

raw1 = io.Raw(raw_fname1)
raw2 = io.Raw(raw_fname2)

events1 = mne.read_events(event_fname1)
events2 = mne.read_events(event_fname2)
events_list = [events1,events2]

raw_conc= concatenate_raws([raw1,raw2], events_list=events_list)

Until here I think these several code lines work correctly...

2) But when I try to extract the epochs from raw_conc I have some error 
messages (concerning I think the events_list)

picks = 
mne.pick_types(raw1.info,meg='mag',eeg=False,ecg=False,eog=False,stim=False, 
exclude='bads')
epochs = mne.Epochs(raw, events_list, event_id, tmin, tmax, picks=picks, 
baseline= baseline, reject=reject)

The aim of both questions is to use the tfr_morlet function which asks 
in input the epochs.

power, itc = tfr_morlet(epochs, freqs=freqs, n_cycles=n_cycles, 
use_fft=False, return_itc=True, decim=3,n_jobs=1)

Someone has an idea about this?

Bests,


Jonathan
-- 
Dr Jonathan Entakli - Research Associate Post doctoral position

Department of Psychology
University of Cambridge, Downing St
Cambridge CB2 3EB

Sir William Hardy Building
Office : 314
Mobile : 07778 692 458

http://www.viscog.psychol.cam.ac.uk/