[Eeglablist] newtimef ersp output without any or with external baseline
Andero Uusberg
Andero.Uusberg at ut.ee
Thu Feb 25 13:03:15 PST 2010
Hello,
I have a few (maybe fairly naive) questions regarding running the
newtimef function without using baseline from the dataset being analyzed.
Firstly, while running newtimef on a channel of a segmented dataset
without either of the baselines (having both 'baseline' and 'powbase'
set to NaN) I have trouble understanding the units of the output ersp
matrix (I call the function from the command line and use freqrange
2:20, padratio 2, cycles 1 6 and basenorm 'off'). The output ersp matrix
is full of positive numbers in the range of 40 to 60. Although the
figure legend states this is 10*log10 transform of power density in
microvolts, this does not seem likely (the power densities would be
around 100 000 microvolts). In addition, even with both baselines set to
NaN the powbase output is created (consisting of numbers in the same
range). All this probably has something to do with the normalization
process and is completely meaningful, but could someone explain how
exactly these values are generated and what is their unit?
Second and related question regards feeding externally created baseline
into newtimef via the powbase key. I know one way would be to calculate
the basline with newtimef in the first place and not to care about the
units of the data :) However, I would also like to be able to use data
calculated with FFT in Vision Anlayzer as the baseline. Do I understand
correctly that the data should be in power denisty units (square of
microvolt per hertz) and should I run the 10*log10 transform on the data
before feeding in into the newtimef with powbase?
Finally, as I'm not sure my previous e-mail reached the list I repeat a
related question concerning using std_precomp for calculating ERSPs for
several datasets using a single external baseline. When referring to a
channel-by-frequencies (as well as frequencies-by-channles) matrix of
external power density values with the powbase key in erspparams field
og std_precom, I receive the following error message: "std_ersp at 372:
powbase should be of size (ncomps,nfreqs)". Can anyone estimate if this
could be a bug or should I keep figuring out what I'm doing wrong?
I apologise if these issues have been dealt with before in the list in
which case I am grateful for a kind reference to older posts.
With kind regards and general gratitude to eeglab :)
Andero Uusberg
University of Tartu
Estonia
More information about the eeglablist
mailing list