<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<font size="4" class="">Dear Craig,</font>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">Let me forward this e-mail to the discussion list, in order to widen the scope to all of Donostia and beyond.</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">Indeed, the estimation of spectrally resolved GC requires a full bandwidth signal. I have not extensively played with the different options here, but the scenario that you describe, i.e. compute a DICS spatial filter on
a narrow bandwidth-of-interest signal, and subjecting the broadband sensor signals to this spatial filter, will, in my opinion, not be really optimal. The reason for this is that the beamformer then can only suppress interference from that particular bandwidth.
Signal components that lie outside this range might still leak into the estimated time courses, which may lead to distortion and invalid GC estimates. </font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">So, if you would be using an LCMV beamformer, I’d always recommend a broadband spatial filter, optionally on mildly highpass-filtered data. Then, the risk could be that the 1/f characteristic of the signals results in a
covariance matrix whose signal components are dominated by low-frequency stuff. One could then worry about suboptimal spatial filtering of the higher frequency components. This is probably true, and I have worried about this myself, but in my experience this
is the least of the worries one should have. (as a side note: I have played around with spectral whitening of the data prior to the beamforming. The reasoning behind this is that the individual variance components in the data are more evenly balanced across
the frequency range, but the long story short is that this does not seem to influence the final GC estimate).</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">A larger worry would be the fact that source-estimated GC spectra very often have a canonical shape, unless you use extremely well-defined data such a visual attention experiment that induces strong gamma band oscillations.
If you are into applying connectivity analysis in a more ethereal context, such as a cognitive or language experiment, you will see that source-level GC spectra will be always dominated by the low-frequency components (i.e. most of the time you’ll find a large
peak at low frequencies). This of course can reflect genuine neuronal coupling, but my suspicion is that this is most of the time due to ‘weak signal asymmetries’, in other words due to asymmetries in SNR and location specific leakage of power. This can be
further investigated using time-flipping and an appropriate statistical evaluation. For inspiration, you can of course always consult my 2017 PNAS paper :).</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">Another worry of course would be the fact that you’d be working with Elekta data :).</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">Practically, what I typically do, is to compute the broadband spatial filter with LCMV, and apply these filters (as is, or dimensionality reduced after a parcellation step) to the sensor level cross-spectral density data
(from DC to Nyquist). The spectral factorization algorithm should be applied to the source-level data, which in my experience is most stable if you do this pairwise. I know that this is suboptimal in the context of an interpretation of direct versus indirect
interactions, but I have always found this an argument that is part of a non-discussion at this stage of the analysis. One should first trust the estimates overall, before a discussion about direct versus indirect coupling becomes relevant.</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class="">Best wishes,</font></div>
<div class=""><font size="4" class="">JM</font></div>
<div class=""><br class="">
</div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><font size="4" class=""><br class="">
</font></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 19 Jun 2020, at 18:06, Craig Richter <<a href="mailto:craiggrichter@gmail.com" class="">craiggrichter@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">Hey JM,<br class="">
<br class="">
How’s it going? Covid treating you ok? Been a bit crazy here, but normalcy is resuming.<br class="">
<br class="">
I have a question about DICS for you… :-/... Anyhow, so for coherence at a given frequency DICS is superior, since the spatial-filter is specific to that frequency, which is about the same as a building a noise-covariance matrix from a narrow-band filtered
signal? But to do source-level GC, we need the broadband time-series, so we’re forced to use the unfiltered data and LCMV. But… what if we know the frequency we want to analyse, so we build the spatial filter on this frequency, or narrow-band filtered data,
then compute virtual time series from the raw data, and compute GC. Then we simply look at our frequency of interest in the spectral results. Would this potentially generate a better GC estimate at the frequency of interest, or would it just mess everything
up?<br class="">
<br class="">
PS - I’ve cc’ed some other interested folks.<br class="">
<br class="">
Thanks!<br class="">
<br class="">
C.</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>