Some issues in freqanalysis

Jan Mathijs Schoffelen Jan.Schoffelen at FCDONDERS.RU.NL
Tue Nov 28 23:24:19 CET 2006


Dear Vladimir,



Regarding the missing functions: you are absolutely right. The functions
will be added to the latest release.

Regarding the crash: I don't know what is happening there, and I do not have
any experience with the option biascorrect. You'll have to ask Eric about
this.

Does it still crash when you use cfg.biascorrect = 'no'? This might already
solve your problem.



cfg.combinechan is a rather obscure option (and is set to 'no' by default)
in freqdescriptives which is (as far as I know) undocumented. It has to do
with combining planar gradiometer representation back into an axial
gradiometer representation in order to be able to compute coherence between
sensor-combinations (in their original axial representation). As a matter of
fact, it is an option which has a high 'don't try this at home'-content. At
least a correct behavior is not guaranteed. Anyhow, it should run through
even when it is not specified.



If you want to test coherence against 0, the jack-knife estimate of the
variance is not going to help you too much. Rather, you could use some kind
of Z-transformation as proposed by Jarvis & Mitra Neural computation 2001,
or the Z-transformation that is described in James Kilner's papers, or from
the Rosenberg group. In this context, the coherence can be parametrically
evaluated.

Another possibility is using a 'shiftpredictor' to estimate the bias,
meaning that you shuffle the trial order of one of the two channels in your
channel-combination (only works when taking fourier as an output of
freqanalysis). In my experience this only works if you use an external
reference signal, such as the EMG.



Hope this gives some handles,



Jan-Mathijs





  _____

From: FieldTrip discussion list [mailto:FIELDTRIP at NIC.SURFNET.NL] On Behalf
Of Vladimir Litvak
Sent: Tuesday, November 28, 2006 8:33 PM
To: FIELDTRIP at NIC.SURFNET.NL
Subject: [FIELDTRIP] Some issues in freqanalysis



Dear Robert and all,





I'm writing a new script for analysis of MEG data and I encountered some
problems with FT (I'm using 20061029 version):



1)       In freqdescriptives the lines 494 and 515 seem to be redundant and
also don't work and cause a crash.

2)       there are some missing functions:

     inputlabel2outputlabel called by freqdescriptives, line 234

     fourier2crsspctrm called by freqdescriptives, line 290





Can you please explain to me what cfg.combinechan in freqdescriptives does
(I'm new to MEG which might explain my ignorance)?



Finally, I would like to compute coherence from data of a single subject
with 300 trials and check whether it is significant (i.e. not random). I
don't think I can use any of the nonparametric stuff because I don't have a
baseline. It seems to me that I should use the jackknife SEMs. How exactly?
Should I use the bias correction? Should I take 2*SEM as the threshold?



Thanks,



Vladimir







-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20061128/e2416a42/attachment-0002.html>


More information about the fieldtrip mailing list