<div><div>It would be nice to actually see the baseline.</div></div><div>The question is whether it is roughly zero mean or not. For the covariance noise model to be appropriate the data should be zero mean, as is typically roughly the case after baseline correction or filtering with a high pass. What is plausible here also depends to on your data. Using the noise covariance here will relate your data to the amplitude structure of what you declare noise. Plot this for example using an empty room noise cov.</div><div>Overall, from the distance your plots look plausible to me.</div><div><br></div><div><div><div class="gmail_quote"><div>On Wed, 10 May 2017 at 06:03, Graham Flick <<a href="mailto:grahamflick00@gmail.com" target="_blank">grahamflick00@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>Hi All,</div><div><br></div><div>I have a set of MEG data collected in a sentence processing paradigm, where the critical words occur 6-7 words into the sentence. I'd like to look at source-level evoked responses to these words via minimum-norm estimates, without applying baseline correction.</div><div><br></div><div>In this scenario, should I still apply baseline correction to the pre-stimulus intervals that I use to estimate the noise covariance? Note that in this design, pre-stimulus is actually pre-sentence, meaning that there is about 4 seconds of data between these windows and the onset of the epochs that will be inverted to source space.</div><div><br></div><div>In attempt to address this question, I've plotted whitened evoked responses from the start of the sentence to the target words using different methods of covariance estimation, with and without baseline correction applied to the 100 ms windows from which I estimated the covariance. I've attached an example from one subject, and the pattern shown there is consistent across quite a few subjects in the sample. </div><div><br></div><div>In general, it looks like if I apply baseline correction to the window from which I estimate covariance, the global field power of the whitened response never reaches 1, even in the window in which the covariance was estimated. In contrast, the GFP in the whitened response without baseline correction looks more like what I'd expect to see. This pattern seems unusual to me, but does it imply that I should not be be applying baseline correction here? Or are there other factors that should be considered?</div><div><br></div><div>Thanks!</div><div><br></div><div>Graham</div><div><br></div><div><br></div><div>Here is a sample of the code used to generate the whitened responses for the empirical estimator with/without baseline correction:</div><div><br></div><div>raw = mne.io.read_raw_fif(fname_raw, preload=True)</div><div>events = mne.read_events(fname_event)</div><div>picks = mne.pick_types(<a href="http://raw.info" target="_blank">raw.info</a>, meg=True, eeg=False, eog=False,exclude=bads)</div><div>epochstargetFull = mne.Epochs(raw, events, event_id = event_id, tmin=-4.4,tmax=1.2, decim=5,reject=dict(mag=2e-12),baseline=None,picks=picks,on_missing='ignore')</div><div>evokedtargetFull = epochstargetFull.average()</div><div><br></div><div>method = 'empirical'</div><div><br></div><div># covariance with baseline correction applied</div><div>epochscov = mne.Epochs(raw, events, event_id=event_id, tmin=-4.4, tmax=-4.3, decim=5, reject=dict(mag=2e-12), baseline=(-4.4,-4.3), picks=picks, on_missing='ignore')</div><div>cov = mne.compute_covariance(epochscov, tmin=-4.4, tmax=-4.3, method=method')</div><div>tmp = evokedtargetFull.plot_white(cov, show=False)</div><div>tmp.savefig('topright_empirical_Baselined.png')</div><div>del(epochscov)</div><div>del(cov)</div><div>del(tmp)</div><div><br></div><div># covariance without baseline correction applied</div><div>epochscov = mne.Epochs(raw, events, event_id=event_AF, tmin=-4.4, tmax=-4.3, decim=5, reject=dict(mag=2e-12), baseline=None, picks=picks, on_missing='ignore')</div><div>cov = mne.compute_covariance(epochscov, tmin=-4.4, tmax=-4.3, method=method)</div><div>tmp = evokedtargetFull.plot_white(cov, show=False)</div><div>tmp.savefig('topleft_empirical_NoBaseline.png')</div><div><br></div></div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/complianceline</a> . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.<br>
</blockquote></div></div></div>