Dear Makoto,<br><br>I double checked it and the results are identical for both using 'padratio' 4 and 'nfreqs' 47. I use log-spaced frequencies, so both settings lead to the same 47 log-spaced frequencies from 4 to 50 Hz (and all the other input settings stay the same).<br>
<br>[ersp,itc,powbase,times,freqs] = newtimef(EEG.data (83,:,:), 350, [-400 996], 250, [1 0.5],'baseline',[NaN],'freqs', [4 50], 'freqscale', 'log', 'padratio', 4);<br><br>[ersp1,itc1,powbase1,times1,freqs1] = newtimef(EEG.data (83,:,:), 350, [-400 996], 250, [1 0.5],'baseline',[NaN],'freqs', [4 50], 'nfreqs', 47, 'freqscale', 'log');<br>
<br>Changing the padratio to 2 and 8 changes the frequency resolution (and that is what it should do, right?).<br>Using a padratio of 2 results in 23 log-spaced frequencies and a padratio of 8 results in 94 log-spaced frequencies.<br>
<br>Thanks a lot, Katharina<br><br><br><div class="gmail_quote">On 3 August 2012 11:56, Makoto Miyakoshi <span dir="ltr"><<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank">mmiyakoshi@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Katharina,<br>
<div class="im"><br>
> I tried using 'nfreqs' 47 instead of 'padratio' 4 as you suggested. The<br>
> results are identical.<br>
<br>
</div>I'm not sure if this could be true. Could you check it again?<br>
Also, please try 'padratio', 2 and 8 to see if it changes results.<br>
<div class="im"><br>
> Why is it better to define the number of frequencies<br>
> than to use the padratio?<br>
<br>
</div>I misunderstood you. I thought you were saying you want to obtain<br>
linear scale from 4 to 50 Hz with 1 Hz interval. Please forget it.<br>
<br>
Makoto<br>
<br>
2012/8/2 Katharina Limbach <<a href="mailto:klim366@aucklanduni.ac.nz">klim366@aucklanduni.ac.nz</a>>:<br>
<div class="HOEnZb"><div class="h5">> Dear Makoto,<br>
><br>
> thanks a lot for your response.<br>
> I think I figured out, that I was making a really stupid mistake. I thought<br>
> that the output parameters were optional and that their order did not matter<br>
> (I used [ersp,times,freqs]) . So the output that was given for 'freqs' was<br>
> actually the 'powbase' output (and likewise the 'times' output was the<br>
> itc...).<br>
><br>
> So I changed my code to:<br>
><br>
> [ersp,itc,powbase,times,freqs] = newtimef(EEG.data (83,:,:), 350, [-400<br>
> 996], 250, [1 0.5],'baseline',[NaN],'freqs', [4 50], 'freqscale', 'log',<br>
> 'padratio', 4);<br>
><br>
> and I am now getting more sensible results (times from -260 to 860 and<br>
> frequencies from 4 to 50).<br>
> I am really sorry for this!<br>
><br>
> I tried using 'nfreqs' 47 instead of 'padratio' 4 as you suggested. The<br>
> results are identical. Why is it better to define the number of frequencies<br>
> than to use the padratio?<br>
><br>
> I have one more question about the powbase. I do not want to baseline my<br>
> wavelets (and I set 'baseline' [NaN]), but I am getting a powbase output.<br>
> Does that mean that the baseline power is computed but not used for the<br>
> wavelets?<br>
><br>
> Thanks a lot, Katharina<br>
><br>
><br>
><br>
> On 3 August 2012 08:06, Makoto Miyakoshi <<a href="mailto:mmiyakoshi@ucsd.edu">mmiyakoshi@ucsd.edu</a>> wrote:<br>
>><br>
>> Dear Katharina,<br>
>><br>
>> Sorry for slow response.<br>
>><br>
>> > The values range from 1.6608914e+04 to 3.0225791e+05.<br>
>><br>
>> This means your frequency range is between 1.6 MHz to 30.2 MHz, which<br>
>> is like an old generation CPU clock speed. Something is seriously<br>
>> wrong. However, I don't know what causes this from your code.<br>
>><br>
>> By the way, if you want something close to 4,5,6,... 50, you should<br>
>> use 'nfreq', 47 instead of 'padratio', 4. Could you try it to see what<br>
>> happens?<br>
>><br>
>> Makoto<br>
>><br>
>> 2012/7/22 Katharina Limbach <<a href="mailto:klim366@aucklanduni.ac.nz">klim366@aucklanduni.ac.nz</a>>:<br>
>> > Dear all,<br>
>> ><br>
>> > I have a question regarding the output of the newtimef function.<br>
>> > I am using the following command (input data: 1 participant, 1<br>
>> > condition,<br>
>> > electrode # 67. My sampling rate is 250 Hz and the epochs go from -400<br>
>> > to<br>
>> > 996 ms (350 frames)):<br>
>> ><br>
>> > [ersp,times,freqs] = newtimef(EEG.data (67,:,:), 350, [-400 996], 250,<br>
>> > [1<br>
>> > 0.5],'baseline',[NaN], 'freqs', [4 50], 'freqscale', 'log', 'padratio',<br>
>> > 4,<br>
>> > 'plotphase', 'off', 'plotitc', 'off');<br>
>> ><br>
>> > The resulting plot is showing frequencies from 4 to 50 Hz (which is what<br>
>> > I<br>
>> > was looking for) but I do not understand the output I get for 'freqs':<br>
>> > The values range from 1.6608914e+04 to 3.0225791e+05.<br>
>> > Shouldn't they simply be in Hz and go from 4 to 50 Hz? Do I have to<br>
>> > convert<br>
>> > them to get the correct values or is there something else that I am<br>
>> > doing<br>
>> > wrong?<br>
>> ><br>
>> > Also my 'times' output is not a vector of output times in ms, but a  47<br>
>> > x<br>
>> > 200 martix with complex numbers.<br>
>> ><br>
>> ><br>
>> > Thanks a lot for your help,<br>
>> > Katharina<br>
>> ><br>
>> > --<br>
>> > Katharina Limbach<br>
>> > PhD Candidate<br>
>> > Department of Psychology<br>
>> > University of Auckland<br>
>> > Private Bag 92019<br>
>> > Auckland<br>
>> > New Zealand<br>
>> > tel: + 64 9 3737 599 ext 85969<br>
>> ><br>
>> ><br>
>> > _______________________________________________<br>
>> > Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
>> > To unsubscribe, send an empty email to<br>
>> > <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
>> > For digest mode, send an email with the subject "set digest mime" to<br>
>> > <a href="mailto:eeglablist-request@sccn.ucsd.edu">eeglablist-request@sccn.ucsd.edu</a><br>
>><br>
>><br>
>><br>
>> --<br>
>> Makoto Miyakoshi<br>
>> JSPS Postdoctral Fellow for Research Abroad<br>
>> Swartz Center for Computational Neuroscience<br>
>> Institute for Neural Computation, University of California San Diego<br>
<br>
<br>
<br>
--<br>
Makoto Miyakoshi<br>
JSPS Postdoctral Fellow for Research Abroad<br>
Swartz Center for Computational Neuroscience<br>
Institute for Neural Computation, University of California San Diego<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><font color="#888888"><div><span>Katharina Limbach</span></div>
<div>PhD Candidate</div>
<div>Department of Psychology</div>
<div>University of Auckland</div>
<div>Private Bag 92019</div>
<div>Auckland</div>
<div>New Zealand</div>
<div>tel: <a value="+6493737599">+ 64 9 3737 599 ext 85969</a></div></font><br>