External Email - Use Caution
Irina,
Here is what I use to get grand average sources:
all_stcs=3Dnp.average([x.data for x in all_stcs2], axis=3D0) # you should have list with all stcs
morphed_stcs=mne.SourceEstimate(data1, all_stcs[0].vertices, all_stcs[0].=
tmin, all_stcs[0].tstep, all_stcs[0].subject)
# this will save the grand av=
erage stc