<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Thank you all so much for the suggestions!!!! :)<div class=""><br class=""></div><div class="">After adding /cfg.elec = sens/ it finally worked.</div><div class=""><br class=""><div class="">have a nice day,</div><div class="">~ grazia</div><div class=""><div apple-content-edited="true" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><br class=""></div></div>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On 21Oct, 2015, at 20:22, <a href="mailto:fieldtrip-request@science.ru.nl" class="">fieldtrip-request@science.ru.nl</a> wrote:</div><br class="Apple-interchange-newline"><div class="">Send fieldtrip mailing list submissions to<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span><a href="mailto:fieldtrip@science.ru.nl" class="">fieldtrip@science.ru.nl</a><br class=""><br class="">To subscribe or unsubscribe via the World Wide Web, visit<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">or, via email, send a message with subject or body 'help' to<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>fieldtrip-request@science.ru.nl<br class=""><br class="">You can reach the person managing the list at<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>fieldtrip-owner@science.ru.nl<br class=""><br class="">When replying, please edit your Subject line so it is more specific<br class="">than "Re: Contents of fieldtrip digest..."<br class=""><br class=""><br class="">Today's Topics:<br class=""><br class=""> 1. Re: Interpolation and ft_channelrepair (Grazia Di Pisa)<br class=""> 2. Re: Interpolation and ft_channelrepair<br class=""> (=?utf-8?Q?J=C3=B6rn_M._Horschig?=)<br class=""> 3. Re: Interpolation and ft_channelrepair (Daria Laptinskaya)<br class=""> 4. Re: Interpolation and ft_channelrepair (Vit?ria Piai)<br class=""> 5. Re: calculating behavioural-power correlation -- follow-up<br class=""> questions (Krebber, Martin)<br class=""> 6. Re: calculating behavioural-power correlation -- follow-up<br class=""> questions (Arjen Stolk)<br class=""><br class=""><br class="">----------------------------------------------------------------------<br class=""><br class="">Message: 1<br class="">Date: Wed, 21 Oct 2015 12:27:40 +0200<br class="">From: Grazia Di Pisa <g.dipisa@gmail.com><br class="">To: fieldtrip@science.ru.nl<br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <D4187F9B-4F67-434D-9E1B-A9A12D920A26@gmail.com><br class="">Content-Type: text/plain; charset="utf-8"<br class=""><br class="">Hi J?rn,<br class=""><br class="">Thanks for the advice - I've just tried to add cfg.sens = sens but still the same error.<br class=""><br class="">Any other suggestions?<br class=""><br class="">thanks,<br class="">~ grazia<br class=""><br class=""><br class=""><blockquote type="cite" class="">On 21Oct, 2015, at 12:00, fieldtrip-request@science.ru.nl wrote:<br class=""><br class="">Send fieldtrip mailing list submissions to<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>fieldtrip@science.ru.nl<br class=""><br class="">To subscribe or unsubscribe via the World Wide Web, visit<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">or, via email, send a message with subject or body 'help' to<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>fieldtrip-request@science.ru.nl<br class=""><br class="">You can reach the person managing the list at<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>fieldtrip-owner@science.ru.nl<br class=""><br class="">When replying, please edit your Subject line so it is more specific<br class="">than "Re: Contents of fieldtrip digest..."<br class=""><br class=""><br class="">Today's Topics:<br class=""><br class=""> 1. Call for new post docs in the Obleser lab | Auditory<br class=""> Cognition (Jonas Obleser)<br class=""> 2. Spectrum normalization (Wunderle, Thomas)<br class=""> 3. Interpolation and ft_channelrepair (Grazia Di Pisa)<br class=""> 4. Re: Interpolation and ft_channelrepair<br class=""> (=?UTF-8?Q?J=C3=B6rn_M._Horschig?=)<br class=""><br class=""><br class="">----------------------------------------------------------------------<br class=""><br class="">Message: 1<br class="">Date: Tue, 20 Oct 2015 06:07:51 -0500<br class="">From: Jonas Obleser <jonas@obleser.de><br class="">To: fieldtrip@science.ru.nl, AUDITORY@LISTS.MCGILL.CA<br class="">Subject: [FieldTrip] Call for new post docs in the Obleser lab |<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Auditory<span class="Apple-tab-span" style="white-space:pre"> </span>Cognition<br class="">Message-ID: <71625964-4C4D-41CA-BEDB-171CEFBC87F8@obleser.de><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">The research group ?Auditory Cognition? (headed by Jonas Obleser; auditorycognition.com), now at the University of L?beck, is seeking to hire two Postdoctoral researchers, initially for 3 years, with the option of extension.<br class=""><br class="">These positions fall into the larger framework of an ERC Consolidator grant ?The listening challenge: How ageing brains adapt? (AUDADAPT) recently awarded to Jonas Obleser, and will allow the joint development of cognitive-neuroscience and psychological research projects targeting adaptive control in the auditory modality of middle-aged adults.<br class=""><br class="">Please download the full job advert here: https://www.dropbox.com/s/sqc9c94ewl6bjjh/Ausschreibung_1031_15_Wiss%20MitarbeiterIn_Psychologie_EN.pdf?dl=1 <br class=""><br class="">Note the application deadline of November 12 (postal receipt).<br class=""><br class="">http://auditorycognition.com<br class="">http://uni-luebeck.de<br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 2<br class="">Date: Tue, 20 Oct 2015 14:58:25 +0000<br class="">From: "Wunderle, Thomas" <thomas.wunderle@esi-frankfurt.de><br class="">To: "fieldtrip@science.ru.nl" <fieldtrip@science.ru.nl><br class="">Subject: [FieldTrip] Spectrum normalization<br class="">Message-ID:<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span><27E5CAD9145EEC41BB9B34C01716A198880C3CD4@UM-excdag-a02.um.gwdg.de><br class="">Content-Type: text/plain; charset="iso-8859-1"<br class=""><br class="">Hi all,<br class=""><br class="">I have a question about how FieldTrip normalizes the output of a spectral analyses using ft_freqanalysis.<br class="">This is important for me, because I need the absolute value of the spectrum for a publication.<br class=""><br class="">I was scanning through the code trying to find the normalization steps.<br class="">In particular, I don't understand the following steps:<br class=""><br class=""><br class="">1.) The tapers (e.g. hanning) are normalized prior to the multiplication with the data by the following code:<br class="">tap = tap./norm(tap, 'fro'); (line 195 in the function fr_specest_mtmfft)<br class="">What is the rationale of dividing through this special norm of the taper itself?<br class=""><br class=""><br class="">2.) After calculating the fft, the specrum is scaled (i.e. in the function fr_specest_mtmfft) by:<br class=""><br class="">dum = dum .* sqrt(2 ./ endnsample); (line 251)<br class=""><br class="">Later on, to get to the power spectrum, the abs of the spectrum is taken and squared<br class=""><br class="">powdum = abs(spectrum(acttap,:,foiind(ifoi),acttboi)) .^2; (line 629 in ft_frequanalysis)<br class=""><br class=""><br class=""><br class="">This puzzles me, because scanning through the literature, the power spectrum should be given by<br class=""><br class="">(2*abs(X)^2) ./ (nfft^2) -> unit: V^2 (X=signal, nfft = number of data points)<br class=""><br class="">and the PSD by<br class=""><br class="">(2*abs(X)^2) ./ (nfft*Fs) -> unit: V^2/Hz (X=signal, nfft = number of data points, Fs=sampling rate)<br class=""><br class=""><br class=""><br class="">Alternatively, one can get the amplitude spectrum from the power spectrum by taking the square root of the PSD:<br class=""><br class="">sqrt(PSD) -> unit: V/sqrt(Hz)<br class=""><br class="">which is still different from<br class=""><br class="">(2*abs(X)) ./ (nfft*Fres) -> unit: V/Hz ? (Fres = frequency resolution = Fs/nfft)<br class=""><br class=""><br class=""><br class="">The result of these equations are different from what comes out of the FieldTrip code.<br class=""><br class="">So what is the unit of the spectrum in FieldTrip?<br class=""><br class="">What is the rationale of scaling the complex spectrum by sqrt(2/nfft) ?<br class=""><br class="">Furthermore, is the unit for the different spectral estimation methods the same? That is, using mtmfft, mtmconvol, wavelet,...?<br class=""><br class="">It would help me a lot if someone could clarify these points.<br class=""><br class="">Thanks in advance,<br class="">Thomas<br class=""><br class=""><br class="">-----<br class="">Dr. Thomas Wunderle<br class="">Ernst Str?ngmann Institute (ESI) for Neuroscience<http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">in Cooperation with Max Planck Society <http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">Deutschordenstrasse 46<br class="">60528 Frankfurt am Main, Germany<br class="">www.esi-frankfurt.de<http://www.esi-frankfurt.de/><br class="">thomas.wunderle@esi-frankfurt.de<mailto:thomas.wunderle@esi-frankfurt.de><br class="">Tel: +49 69 96769 516<br class="">Fax: +49 69 96769 555<br class=""><br class="">Sitz der Gesellschaft: Frankfurt am Main<br class="">Registergericht: Amtsgericht Frankfurt - HRB 84266<br class="">Gesch?ftsf?hrer: Prof. Dr. Pascal Fries<br class=""><br class=""><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151020/eebaf460/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 3<br class="">Date: Wed, 21 Oct 2015 10:30:07 +0200<br class="">From: Grazia Di Pisa <g.dipisa@gmail.com><br class="">To: fieldtrip@science.ru.nl<br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <A46541FB-BC67-44E7-A9CD-541FB39B7C54@gmail.com><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m getting the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should specify the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code? <br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 <span class="Apple-tab-span" style="white-space:pre"> </span>cfg = [];<br class="">16 <span class="Apple-tab-span" style="white-space:pre"> </span>cfg_neighb <span class="Apple-tab-span" style="white-space:pre"> </span>= [];<br class="">17<span class="Apple-tab-span" style="white-space:pre"> </span>cfg_neighb.method<span class="Apple-tab-span" style="white-space:pre"> </span>= 'triangulation'; <br class="">18 cfg.senstype <span class="Apple-tab-span" style="white-space:pre"> </span>= 'EEG';<br class="">19 cfg_neighb.layout <span class="Apple-tab-span" style="white-space:pre"> </span>= 'biosemi64.lay';<br class="">20 neighbours <span class="Apple-tab-span" style="white-space:pre"> </span>= ft_prepare_neighbours(cfg_neighb);<br class="">21 <br class="">22 cfg.neighbours <span class="Apple-tab-span" style="white-space:pre"> </span>= neighbours;<br class="">23 cfg.channel <span class="Apple-tab-span" style="white-space:pre"> </span>= {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9',}; <br class="">24 cfg.feedback <span class="Apple-tab-span" style="white-space:pre"> </span> = 'yes';<br class="">25 cfg.layout <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= 'biosemi64.lay';<br class="">26 lay <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= ft_prepare_layout(cfg);<br class="">27 <br class="">28 sens <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= data;<br class="">29 sens.type <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= 'eeg';<br class="">30 sens.label<span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= lay.label;<br class="">31 sens.chanpos <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= lay.pos;<br class="">32 sens.chanpos(:,3) <span class="Apple-tab-span" style="white-space:pre"> </span>= 0;<br class="">33 <br class="">34 ft_neighbourplot(cfg, data)<br class="">35 <br class="">36 cfg <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= [];<br class="">37 cfg.method <span class="Apple-tab-span" style="white-space:pre"> </span>= 'nearest';<br class="">38 cfg.badchannel <span class="Apple-tab-span" style="white-space:pre"> </span>= {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours <span class="Apple-tab-span" style="white-space:pre"> </span> = neighbours;<br class="">40 cfg.trials <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span> = 'all';<br class="">41 <br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 4<br class="">Date: Wed, 21 Oct 2015 10:42:03 +0200<br class="">From: "=?UTF-8?Q?J=C3=B6rn_M._Horschig?=" <jorn@artinis.com><br class="">To: "'FieldTrip discussion list'" <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <002b01d10bdc$5d377040$17a650c0$@artinis.com><br class="">Content-Type: text/plain;<span class="Apple-tab-span" style="white-space:pre"> </span>charset="UTF-8"<br class=""><br class="">Dear Grazia,<br class=""><br class="">try adding cfg.sens = sens; before line 42. That will set the sensor definition so that ft_channelrepair knows the location of your channels, which is otherwise not obvious.<br class=""><br class="">Best,<br class="">J?rn<br class=""><br class="">--<br class=""><br class="">J?rn M. Horschig, PhD, Software Engineer<br class="">Artinis Medical Systems | +31 481 350 980 <br class=""><br class=""><blockquote type="cite" class="">-----Original Message-----<br class="">From: fieldtrip-bounces@science.ru.nl [mailto:fieldtrip-<br class="">bounces@science.ru.nl] On Behalf Of Grazia Di Pisa<br class="">Sent: Wednesday, October 21, 2015 10:30 AM<br class="">To: fieldtrip@science.ru.nl<br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m getting<br class="">the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should specify<br class="">the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code?<br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 <span class="Apple-tab-span" style="white-space:pre"> </span>cfg = [];<br class="">16 <span class="Apple-tab-span" style="white-space:pre"> </span>cfg_neighb <span class="Apple-tab-span" style="white-space:pre"> </span>= [];<br class="">17<span class="Apple-tab-span" style="white-space:pre"> </span>cfg_neighb.method<span class="Apple-tab-span" style="white-space:pre"> </span>= 'triangulation';<br class="">18 cfg.senstype <span class="Apple-tab-span" style="white-space:pre"> </span>= 'EEG';<br class="">19 cfg_neighb.layout <span class="Apple-tab-span" style="white-space:pre"> </span>= 'biosemi64.lay';<br class="">20 neighbours <span class="Apple-tab-span" style="white-space:pre"> </span>= ft_prepare_neighbours(cfg_neighb);<br class="">21<br class="">22 cfg.neighbours <span class="Apple-tab-span" style="white-space:pre"> </span>= neighbours;<br class="">23 cfg.channel <span class="Apple-tab-span" style="white-space:pre"> </span>= {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9',};<br class="">24 cfg.feedback <span class="Apple-tab-span" style="white-space:pre"> </span> = 'yes';<br class="">25 cfg.layout <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= 'biosemi64.lay';<br class="">26 lay <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= ft_prepare_layout(cfg);<br class="">27<br class="">28 sens <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= data;<br class="">29 sens.type <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= 'eeg';<br class="">30 sens.label<span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= lay.label;<br class="">31 sens.chanpos <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= lay.pos;<br class="">32 sens.chanpos(:,3) <span class="Apple-tab-span" style="white-space:pre"> </span>= 0;<br class="">33<br class="">34 ft_neighbourplot(cfg, data)<br class="">35<br class="">36 cfg <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>= [];<br class="">37 cfg.method <span class="Apple-tab-span" style="white-space:pre"> </span>= 'nearest';<br class="">38 cfg.badchannel <span class="Apple-tab-span" style="white-space:pre"> </span>= {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours <span class="Apple-tab-span" style="white-space:pre"> </span> = neighbours;<br class="">40 cfg.trials <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span> = 'all';<br class="">41<br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""></blockquote><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class="">End of fieldtrip Digest, Vol 59, Issue 17<br class="">*****************************************<br class=""></blockquote><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151021/43fd9e60/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 2<br class="">Date: Wed, 21 Oct 2015 14:16:10 +0200<br class="">From: "=?utf-8?Q?J=C3=B6rn_M._Horschig?=" <jorn@artinis.com><br class="">To: "'FieldTrip discussion list'" <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <002301d10bfa$45d6a740$d183f5c0$@artinis.com><br class="">Content-Type: text/plain; charset="utf-8"<br class=""><br class="">Hi Grazia,<br class=""><br class=""><br class=""><br class="">uh, jah, my bad, try cfg.elec instead ;) <br class=""><br class=""><br class=""><br class="">Best,<br class=""><br class="">J?rn<br class=""><br class=""><br class=""><br class="">--<br class=""><br class=""><br class=""><br class="">J?rn M. Horschig, PhD, Software Engineer<br class=""><br class=""> <http://www.artinis.com/> Artinis Medical Systems | +31 481 350 980 <br class=""><br class=""><br class=""><br class="">From: fieldtrip-bounces@science.ru.nl [mailto:fieldtrip-bounces@science.ru.nl] On Behalf Of Grazia Di Pisa<br class="">Sent: Wednesday, October 21, 2015 12:28 PM<br class="">To: fieldtrip@science.ru.nl<br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class=""><br class=""><br class=""><br class="">Hi J?rn,<br class=""><br class=""><br class=""><br class="">Thanks for the advice - I've just tried to add cfg.sens = sens but still the same error.<br class=""><br class=""><br class=""><br class="">Any other suggestions?<br class=""><br class=""><br class=""><br class="">thanks,<br class=""><br class="">~ grazia<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">On 21Oct, 2015, at 12:00, fieldtrip-request@science.ru.nl <mailto:fieldtrip-request@science.ru.nl> wrote:<br class=""><br class=""><br class=""><br class="">Send fieldtrip mailing list submissions to<br class=""> fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> <br class=""><br class="">To subscribe or unsubscribe via the World Wide Web, visit<br class=""> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">or, via email, send a message with subject or body 'help' to<br class=""> fieldtrip-request@science.ru.nl <mailto:fieldtrip-request@science.ru.nl> <br class=""><br class="">You can reach the person managing the list at<br class=""> fieldtrip-owner@science.ru.nl <mailto:fieldtrip-owner@science.ru.nl> <br class=""><br class="">When replying, please edit your Subject line so it is more specific<br class="">than "Re: Contents of fieldtrip digest..."<br class=""><br class=""><br class="">Today's Topics:<br class=""><br class=""> 1. Call for new post docs in the Obleser lab | Auditory<br class=""> Cognition (Jonas Obleser)<br class=""> 2. Spectrum normalization (Wunderle, Thomas)<br class=""> 3. Interpolation and ft_channelrepair (Grazia Di Pisa)<br class=""> 4. Re: Interpolation and ft_channelrepair<br class=""> (=?UTF-8?Q?J=C3=B6rn_M._Horschig?=)<br class=""><br class=""><br class="">----------------------------------------------------------------------<br class=""><br class="">Message: 1<br class="">Date: Tue, 20 Oct 2015 06:07:51 -0500<br class="">From: Jonas Obleser <jonas@obleser.de <mailto:jonas@obleser.de> ><br class="">To: fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> , AUDITORY@LISTS.MCGILL.CA <mailto:AUDITORY@LISTS.MCGILL.CA> <br class="">Subject: [FieldTrip] Call for new post docs in the Obleser lab |<br class=""> Auditory Cognition<br class="">Message-ID: <71625964-4C4D-41CA-BEDB-171CEFBC87F8@obleser.de <mailto:71625964-4C4D-41CA-BEDB-171CEFBC87F8@obleser.de> ><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">The research group ?Auditory Cognition? (headed by Jonas Obleser; auditorycognition.com), now at the University of L?beck, is seeking to hire two Postdoctoral researchers, initially for 3 years, with the option of extension.<br class=""><br class="">These positions fall into the larger framework of an ERC Consolidator grant ?The listening challenge: How ageing brains adapt? (AUDADAPT) recently awarded to Jonas Obleser, and will allow the joint development of cognitive-neuroscience and psychological research projects targeting adaptive control in the auditory modality of middle-aged adults.<br class=""><br class="">Please download the full job advert here: https://www.dropbox.com/s/sqc9c94ewl6bjjh/Ausschreibung_1031_15_Wiss%20MitarbeiterIn_Psychologie_EN.pdf?dl=1 <br class=""><br class="">Note the application deadline of November 12 (postal receipt).<br class=""><br class="">http://auditorycognition.com<br class="">http://uni-luebeck.de<br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 2<br class="">Date: Tue, 20 Oct 2015 14:58:25 +0000<br class="">From: "Wunderle, Thomas" <thomas.wunderle@esi-frankfurt.de <mailto:thomas.wunderle@esi-frankfurt.de> ><br class="">To: "fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> " <fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> ><br class="">Subject: [FieldTrip] Spectrum normalization<br class="">Message-ID:<br class=""> <27E5CAD9145EEC41BB9B34C01716A198880C3CD4@UM-excdag-a02.um.gwdg.de <mailto:27E5CAD9145EEC41BB9B34C01716A198880C3CD4@UM-excdag-a02.um.gwdg.de> ><br class="">Content-Type: text/plain; charset="iso-8859-1"<br class=""><br class="">Hi all,<br class=""><br class="">I have a question about how FieldTrip normalizes the output of a spectral analyses using ft_freqanalysis.<br class="">This is important for me, because I need the absolute value of the spectrum for a publication.<br class=""><br class="">I was scanning through the code trying to find the normalization steps.<br class="">In particular, I don't understand the following steps:<br class=""><br class=""><br class="">1.) The tapers (e.g. hanning) are normalized prior to the multiplication with the data by the following code:<br class="">tap = tap./norm(tap, 'fro'); (line 195 in the function fr_specest_mtmfft)<br class="">What is the rationale of dividing through this special norm of the taper itself?<br class=""><br class=""><br class="">2.) After calculating the fft, the specrum is scaled (i.e. in the function fr_specest_mtmfft) by:<br class=""><br class="">dum = dum .* sqrt(2 ./ endnsample); (line 251)<br class=""><br class="">Later on, to get to the power spectrum, the abs of the spectrum is taken and squared<br class=""><br class="">powdum = abs(spectrum(acttap,:,foiind(ifoi),acttboi)) .^2; (line 629 in ft_frequanalysis)<br class=""><br class=""><br class=""><br class="">This puzzles me, because scanning through the literature, the power spectrum should be given by<br class=""><br class="">(2*abs(X)^2) ./ (nfft^2) -> unit: V^2 (X=signal, nfft = number of data points)<br class=""><br class="">and the PSD by<br class=""><br class="">(2*abs(X)^2) ./ (nfft*Fs) -> unit: V^2/Hz (X=signal, nfft = number of data points, Fs=sampling rate)<br class=""><br class=""><br class=""><br class="">Alternatively, one can get the amplitude spectrum from the power spectrum by taking the square root of the PSD:<br class=""><br class="">sqrt(PSD) -> unit: V/sqrt(Hz)<br class=""><br class="">which is still different from<br class=""><br class="">(2*abs(X)) ./ (nfft*Fres) -> unit: V/Hz ? (Fres = frequency resolution = Fs/nfft)<br class=""><br class=""><br class=""><br class="">The result of these equations are different from what comes out of the FieldTrip code.<br class=""><br class="">So what is the unit of the spectrum in FieldTrip?<br class=""><br class="">What is the rationale of scaling the complex spectrum by sqrt(2/nfft) ?<br class=""><br class="">Furthermore, is the unit for the different spectral estimation methods the same? That is, using mtmfft, mtmconvol, wavelet,...?<br class=""><br class="">It would help me a lot if someone could clarify these points.<br class=""><br class="">Thanks in advance,<br class="">Thomas<br class=""><br class=""><br class="">-----<br class="">Dr. Thomas Wunderle<br class="">Ernst Str?ngmann Institute (ESI) for Neuroscience<http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">in Cooperation with Max Planck Society <http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">Deutschordenstrasse 46<br class="">60528 Frankfurt am Main, Germany<br class="">www.esi-frankfurt.de <http://www.esi-frankfurt.de%3chttp:/www.esi-frankfurt.de/> <http://www.esi-frankfurt.de/><br class="">thomas.wunderle@esi-frankfurt.de <mailto:thomas.wunderle@esi-frankfurt.de%3cmailto:thomas.wunderle@esi-frankfurt.de> <mailto:thomas.wunderle@esi-frankfurt.de><br class="">Tel: +49 69 96769 516<br class="">Fax: +49 69 96769 555<br class=""><br class="">Sitz der Gesellschaft: Frankfurt am Main<br class="">Registergericht: Amtsgericht Frankfurt - HRB 84266<br class="">Gesch?ftsf?hrer: Prof. Dr. Pascal Fries<br class=""><br class=""><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151020/eebaf460/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 3<br class="">Date: Wed, 21 Oct 2015 10:30:07 +0200<br class="">From: Grazia Di Pisa <g.dipisa@gmail.com <mailto:g.dipisa@gmail.com> ><br class="">To: fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> <br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <A46541FB-BC67-44E7-A9CD-541FB39B7C54@gmail.com <mailto:A46541FB-BC67-44E7-A9CD-541FB39B7C54@gmail.com> ><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m getting the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should specify the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code? <br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 cfg = [];<br class="">16 cfg_neighb = [];<br class="">17 cfg_neighb.method = 'triangulation'; <br class="">18 cfg.senstype = 'EEG';<br class="">19 cfg_neighb.layout = 'biosemi64.lay';<br class="">20 neighbours = ft_prepare_neighbours(cfg_neighb);<br class="">21 <br class="">22 cfg.neighbours = neighbours;<br class="">23 cfg.channel = {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9',}; <br class="">24 cfg.feedback = 'yes';<br class="">25 cfg.layout = 'biosemi64.lay';<br class="">26 lay = ft_prepare_layout(cfg);<br class="">27 <br class="">28 sens = data;<br class="">29 sens.type = 'eeg';<br class="">30 sens.label = lay.label;<br class="">31 sens.chanpos = lay.pos;<br class="">32 sens.chanpos(:,3) = 0;<br class="">33 <br class="">34 ft_neighbourplot(cfg, data)<br class="">35 <br class="">36 cfg = [];<br class="">37 cfg.method = 'nearest';<br class="">38 cfg.badchannel = {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours = neighbours;<br class="">40 cfg.trials = 'all';<br class="">41 <br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 4<br class="">Date: Wed, 21 Oct 2015 10:42:03 +0200<br class="">From: "=?UTF-8?Q?J=C3=B6rn_M._Horschig?=" <jorn@artinis.com <mailto:jorn@artinis.com> ><br class="">To: "'FieldTrip discussion list'" <fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> ><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <002b01d10bdc$5d377040$17a650c0$@artinis.com <mailto:002b01d10bdc$5d377040$17a650c0$@artinis.com> ><br class="">Content-Type: text/plain; charset="UTF-8"<br class=""><br class="">Dear Grazia,<br class=""><br class="">try adding cfg.sens = sens; before line 42. That will set the sensor definition so that ft_channelrepair knows the location of your channels, which is otherwise not obvious.<br class=""><br class="">Best,<br class="">J?rn<br class=""><br class="">--<br class=""><br class="">J?rn M. Horschig, PhD, Software Engineer<br class="">Artinis Medical Systems | +31 481 350 980 <br class=""><br class=""><br class=""><br class=""><br class="">-----Original Message-----<br class="">From: fieldtrip-bounces@science.ru.nl <mailto:fieldtrip-bounces@science.ru.nl> [mailto:fieldtrip-<br class="">bounces@science.ru.nl <mailto:bounces@science.ru.nl> ] On Behalf Of Grazia Di Pisa<br class="">Sent: Wednesday, October 21, 2015 10:30 AM<br class="">To: fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl> <br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m getting<br class="">the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should specify<br class="">the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code?<br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 cfg = [];<br class="">16 cfg_neighb = [];<br class="">17 cfg_neighb.method = 'triangulation';<br class="">18 cfg.senstype = 'EEG';<br class="">19 cfg_neighb.layout = 'biosemi64.lay';<br class="">20 neighbours = ft_prepare_neighbours(cfg_neighb);<br class="">21<br class="">22 cfg.neighbours = neighbours;<br class="">23 cfg.channel = {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9',};<br class="">24 cfg.feedback = 'yes';<br class="">25 cfg.layout = 'biosemi64.lay';<br class="">26 lay = ft_prepare_layout(cfg);<br class="">27<br class="">28 sens = data;<br class="">29 sens.type = 'eeg';<br class="">30 sens.label = lay.label;<br class="">31 sens.chanpos = lay.pos;<br class="">32 sens.chanpos(:,3) = 0;<br class="">33<br class="">34 ft_neighbourplot(cfg, data)<br class="">35<br class="">36 cfg = [];<br class="">37 cfg.method = 'nearest';<br class="">38 cfg.badchannel = {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours = neighbours;<br class="">40 cfg.trials = 'all';<br class="">41<br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl <mailto:fieldtrip@donders.ru.nl> <br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl <mailto:fieldtrip@donders.ru.nl> <br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class="">End of fieldtrip Digest, Vol 59, Issue 17<br class="">*****************************************<br class=""><br class=""><br class=""><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151021/fdead2ac/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 3<br class="">Date: Wed, 21 Oct 2015 14:46:49 +0200<br class="">From: Daria Laptinskaya <daria.laptinskaya@googlemail.com><br class="">To: FieldTrip discussion list <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID:<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span><CADcxnnMJ86nmy6ynfEh4TTvM-detkM8-qGVD2D6sv=sc48d4mQ@mail.gmail.com><br class="">Content-Type: text/plain; charset="utf-8"<br class=""><br class="">Dear Grazia,<br class=""><br class="">I?m using the following code and it works fine:<br class=""><br class="">cfg = [];<br class=""><br class="">cfg.method = 'average';<br class=""><br class="">cfg.badchannel = channel;<br class=""><br class="">cfg.neighbours = neighbours;<br class=""><br class="">cfg.trials = 'all';<br class=""><br class="">cfg.elec = sens;<br class=""><br class=""><br class=""><br class="">[interpol] = ft_channelrepair(cfg, opt1_dur);<br class=""><br class=""><br class=""><br class="">Commend out the cfg.elec leads to the same error as you get. ?sens? in my<br class="">code is the result of ft_read_sens, sens=ft_read_sens('GSN-HydroCel-257.sfp'<br class="">).<br class=""><br class=""><br class=""><br class="">I hope, it helps!<br class=""><br class=""><br class=""><br class="">Best,<br class=""><br class="">Daria<br class=""><br class="">2015-10-21 12:27 GMT+02:00 Grazia Di Pisa <g.dipisa@gmail.com>:<br class=""><br class=""><blockquote type="cite" class="">Hi J?rn,<br class=""><br class="">Thanks for the advice - I've just tried to add cfg.sens = sens but still<br class="">the same error.<br class=""><br class="">Any other suggestions?<br class=""><br class="">thanks,<br class="">~ grazia<br class=""><br class=""><br class="">On 21Oct, 2015, at 12:00, fieldtrip-request@science.ru.nl wrote:<br class=""><br class="">Send fieldtrip mailing list submissions to<br class="">fieldtrip@science.ru.nl<br class=""><br class="">To subscribe or unsubscribe via the World Wide Web, visit<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">or, via email, send a message with subject or body 'help' to<br class="">fieldtrip-request@science.ru.nl<br class=""><br class="">You can reach the person managing the list at<br class="">fieldtrip-owner@science.ru.nl<br class=""><br class="">When replying, please edit your Subject line so it is more specific<br class="">than "Re: Contents of fieldtrip digest..."<br class=""><br class=""><br class="">Today's Topics:<br class=""><br class=""> 1. Call for new post docs in the Obleser lab | Auditory<br class=""> Cognition (Jonas Obleser)<br class=""> 2. Spectrum normalization (Wunderle, Thomas)<br class=""> 3. Interpolation and ft_channelrepair (Grazia Di Pisa)<br class=""> 4. Re: Interpolation and ft_channelrepair<br class=""> (=?UTF-8?Q?J=C3=B6rn_M._Horschig?=)<br class=""><br class=""><br class="">----------------------------------------------------------------------<br class=""><br class="">Message: 1<br class="">Date: Tue, 20 Oct 2015 06:07:51 -0500<br class="">From: Jonas Obleser <jonas@obleser.de><br class="">To: fieldtrip@science.ru.nl, AUDITORY@LISTS.MCGILL.CA<br class="">Subject: [FieldTrip] Call for new post docs in the Obleser lab |<br class="">Auditory Cognition<br class="">Message-ID: <71625964-4C4D-41CA-BEDB-171CEFBC87F8@obleser.de><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">The research group ?Auditory Cognition? (headed by Jonas Obleser;<br class="">auditorycognition.com), now at the University of L?beck, is seeking to<br class="">hire two Postdoctoral researchers, initially for 3 years, with the option<br class="">of extension.<br class=""><br class="">These positions fall into the larger framework of an ERC Consolidator<br class="">grant ?The listening challenge: How ageing brains adapt? (AUDADAPT)<br class="">recently awarded to Jonas Obleser, and will allow the joint development of<br class="">cognitive-neuroscience and psychological research projects targeting<br class="">adaptive control in the auditory modality of middle-aged adults.<br class=""><br class="">Please download the full job advert here:<br class="">https://www.dropbox.com/s/sqc9c94ewl6bjjh/Ausschreibung_1031_15_Wiss%20MitarbeiterIn_Psychologie_EN.pdf?dl=1<br class=""><br class="">Note the application deadline of November 12 (postal receipt).<br class=""><br class="">http://auditorycognition.com<br class="">http://uni-luebeck.de<br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 2<br class="">Date: Tue, 20 Oct 2015 14:58:25 +0000<br class="">From: "Wunderle, Thomas" <thomas.wunderle@esi-frankfurt.de><br class="">To: "fieldtrip@science.ru.nl" <fieldtrip@science.ru.nl><br class="">Subject: [FieldTrip] Spectrum normalization<br class="">Message-ID:<br class=""><27E5CAD9145EEC41BB9B34C01716A198880C3CD4@UM-excdag-a02.um.gwdg.de><br class="">Content-Type: text/plain; charset="iso-8859-1"<br class=""><br class="">Hi all,<br class=""><br class="">I have a question about how FieldTrip normalizes the output of a spectral<br class="">analyses using ft_freqanalysis.<br class="">This is important for me, because I need the absolute value of the<br class="">spectrum for a publication.<br class=""><br class="">I was scanning through the code trying to find the normalization steps.<br class="">In particular, I don't understand the following steps:<br class=""><br class=""><br class="">1.) The tapers (e.g. hanning) are normalized prior to the<br class="">multiplication with the data by the following code:<br class="">tap = tap./norm(tap, 'fro'); (line 195 in the function fr_specest_mtmfft)<br class="">What is the rationale of dividing through this special norm of the taper<br class="">itself?<br class=""><br class=""><br class="">2.) After calculating the fft, the specrum is scaled (i.e. in the<br class="">function fr_specest_mtmfft) by:<br class=""><br class="">dum = dum .* sqrt(2 ./ endnsample); (line 251)<br class=""><br class="">Later on, to get to the power spectrum, the abs of the spectrum is taken<br class="">and squared<br class=""><br class="">powdum = abs(spectrum(acttap,:,foiind(ifoi),acttboi)) .^2; (line 629 in<br class="">ft_frequanalysis)<br class=""><br class=""><br class=""><br class="">This puzzles me, because scanning through the literature, the power<br class="">spectrum should be given by<br class=""><br class="">(2*abs(X)^2) ./ (nfft^2) -> unit: V^2 (X=signal,<br class="">nfft = number of data points)<br class=""><br class="">and the PSD by<br class=""><br class="">(2*abs(X)^2) ./ (nfft*Fs) -> unit: V^2/Hz (X=signal, nfft =<br class="">number of data points, Fs=sampling rate)<br class=""><br class=""><br class=""><br class="">Alternatively, one can get the amplitude spectrum from the power spectrum<br class="">by taking the square root of the PSD:<br class=""><br class="">sqrt(PSD) -> unit: V/sqrt(Hz)<br class=""><br class="">which is still different from<br class=""><br class="">(2*abs(X)) ./ (nfft*Fres) -> unit: V/Hz ? (Fres =<br class="">frequency resolution = Fs/nfft)<br class=""><br class=""><br class=""><br class="">The result of these equations are different from what comes out of the<br class="">FieldTrip code.<br class=""><br class="">So what is the unit of the spectrum in FieldTrip?<br class=""><br class="">What is the rationale of scaling the complex spectrum by sqrt(2/nfft) ?<br class=""><br class="">Furthermore, is the unit for the different spectral estimation methods the<br class="">same? That is, using mtmfft, mtmconvol, wavelet,...?<br class=""><br class="">It would help me a lot if someone could clarify these points.<br class=""><br class="">Thanks in advance,<br class="">Thomas<br class=""><br class=""><br class="">-----<br class="">Dr. Thomas Wunderle<br class="">Ernst Str?ngmann Institute (ESI) for Neuroscience<<br class="">http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">in Cooperation with Max Planck Society <<br class="">http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">Deutschordenstrasse 46<br class="">60528 Frankfurt am Main, Germany<br class="">www.esi-frankfurt.de<http://www.esi-frankfurt.de/><br class="">thomas.wunderle@esi-frankfurt.de<mailto:thomas.wunderle@esi-frankfurt.de><br class="">Tel: +49 69 96769 516<br class="">Fax: +49 69 96769 555<br class=""><br class="">Sitz der Gesellschaft: Frankfurt am Main<br class="">Registergericht: Amtsgericht Frankfurt - HRB 84266<br class="">Gesch?ftsf?hrer: Prof. Dr. Pascal Fries<br class=""><br class=""><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <<br class="">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151020/eebaf460/attachment-0001.html<br class=""><blockquote type="cite" class=""><br class=""></blockquote><br class="">------------------------------<br class=""><br class="">Message: 3<br class="">Date: Wed, 21 Oct 2015 10:30:07 +0200<br class="">From: Grazia Di Pisa <g.dipisa@gmail.com><br class="">To: fieldtrip@science.ru.nl<br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <A46541FB-BC67-44E7-A9CD-541FB39B7C54@gmail.com><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m<br class="">getting the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should<br class="">specify the electrodes.<br class=""><br class=""><br class="">Is there a file I should upload or am I missing something in my code?<br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 cfg = [];<br class="">16 cfg_neighb = [];<br class="">17 cfg_neighb.method = 'triangulation';<br class="">18 cfg.senstype = 'EEG';<br class="">19 cfg_neighb.layout = 'biosemi64.lay';<br class="">20 neighbours = ft_prepare_neighbours(cfg_neighb);<br class="">21<br class="">22 cfg.neighbours = neighbours;<br class="">23 cfg.channel = {'Fp1', 'Fp2', 'AF7',<br class="">'AF3','AF8','F7','F5','FT9',};<br class="">24 cfg.feedback = 'yes';<br class="">25 cfg.layout = 'biosemi64.lay';<br class="">26 lay = ft_prepare_layout(cfg);<br class="">27<br class="">28 sens = data;<br class="">29 sens.type = 'eeg';<br class="">30 sens.label = lay.label;<br class="">31 sens.chanpos = lay.pos;<br class="">32 sens.chanpos(:,3) = 0;<br class="">33<br class="">34 ft_neighbourplot(cfg, data)<br class="">35<br class="">36 cfg = [];<br class="">37 cfg.method = 'nearest';<br class="">38 cfg.badchannel = {'Fp1', 'Fp2', 'AF7',<br class="">'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours = neighbours;<br class="">40 cfg.trials = 'all';<br class="">41<br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 4<br class="">Date: Wed, 21 Oct 2015 10:42:03 +0200<br class="">From: "=?UTF-8?Q?J=C3=B6rn_M._Horschig?=" <jorn@artinis.com><br class="">To: "'FieldTrip discussion list'" <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <002b01d10bdc$5d377040$17a650c0$@artinis.com><br class="">Content-Type: text/plain; charset="UTF-8"<br class=""><br class="">Dear Grazia,<br class=""><br class="">try adding cfg.sens = sens; before line 42. That will set the sensor<br class="">definition so that ft_channelrepair knows the location of your channels,<br class="">which is otherwise not obvious.<br class=""><br class="">Best,<br class="">J?rn<br class=""><br class="">--<br class=""><br class="">J?rn M. Horschig, PhD, Software Engineer<br class="">Artinis Medical Systems | +31 481 350 980<br class=""><br class="">-----Original Message-----<br class="">From: fieldtrip-bounces@science.ru.nl [mailto:fieldtrip-<br class="">bounces@science.ru.nl] On Behalf Of Grazia Di Pisa<br class="">Sent: Wednesday, October 21, 2015 10:30 AM<br class="">To: fieldtrip@science.ru.nl<br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m getting<br class="">the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should<br class="">specify<br class=""><br class="">the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code?<br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 cfg = [];<br class="">16 cfg_neighb = [];<br class="">17 cfg_neighb.method = 'triangulation';<br class="">18 cfg.senstype = 'EEG';<br class="">19 cfg_neighb.layout = 'biosemi64.lay';<br class="">20 neighbours = ft_prepare_neighbours(cfg_neighb);<br class="">21<br class="">22 cfg.neighbours = neighbours;<br class="">23 cfg.channel = {'Fp1', 'Fp2', 'AF7',<br class="">'AF3','AF8','F7','F5','FT9',};<br class="">24 cfg.feedback = 'yes';<br class="">25 cfg.layout = 'biosemi64.lay';<br class="">26 lay = ft_prepare_layout(cfg);<br class="">27<br class="">28 sens = data;<br class="">29 sens.type = 'eeg';<br class="">30 sens.label = lay.label;<br class="">31 sens.chanpos = lay.pos;<br class="">32 sens.chanpos(:,3) = 0;<br class="">33<br class="">34 ft_neighbourplot(cfg, data)<br class="">35<br class="">36 cfg = [];<br class="">37 cfg.method = 'nearest';<br class="">38 cfg.badchannel = {'Fp1', 'Fp2', 'AF7',<br class="">'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours = neighbours;<br class="">40 cfg.trials = 'all';<br class="">41<br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class="">End of fieldtrip Digest, Vol 59, Issue 17<br class="">*****************************************<br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class=""></blockquote>-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151021/8500cbff/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 4<br class="">Date: Wed, 21 Oct 2015 10:44:28 -0700<br class="">From: Vit?ria Piai <v.piai.research@gmail.com><br class="">To: FieldTrip discussion list <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <5627CEFC.1030608@gmail.com><br class="">Content-Type: text/plain; charset="windows-1252"; Format="flowed"<br class=""><br class="">Hi Grazia,<br class=""><br class="">I'm using the following code and it works:<br class=""><br class="">cfgn = [];<br class="">cfgn.method = 'template';<br class="">cfgn.layout = 'biosemi64.lay';<br class=""><br class="">cfg = [];<br class="">cfg.neighbours = ft_prepare_neighbours(cfgn, data);<br class="">cfg.elec = ft_read_sens('standard_1005.elc');<br class="">cfg.missingchannel = [];<br class="">cfg.badchannel = chanrepair; % these are my channels<br class="">cfg.trials = [1 3 5]; % for particular trials, otherwise use 'all'<br class="">data = ft_channelrepair(cfg,data);<br class=""><br class="">Hope this helps<br class="">Vitoria<br class=""><br class=""><br class=""><br class="">On 10/21/2015 3:27 AM, Grazia Di Pisa wrote:<br class=""><blockquote type="cite" class="">Hi J?rn,<br class=""><br class="">Thanks for the advice - I've just tried to add cfg.sens = sens but <br class="">still the same error.<br class=""><br class="">Any other suggestions?<br class=""><br class="">thanks,<br class="">~ grazia<br class=""><br class=""><br class=""><blockquote type="cite" class="">On 21Oct, 2015, at 12:00, fieldtrip-request@science.ru.nl <br class=""><mailto:fieldtrip-request@science.ru.nl> wrote:<br class=""><br class="">Send fieldtrip mailing list submissions to<br class="">fieldtrip@science.ru.nl <mailto:fieldtrip@science.ru.nl><br class=""><br class="">To subscribe or unsubscribe via the World Wide Web, visit<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">or, via email, send a message with subject or body 'help' to<br class="">fieldtrip-request@science.ru.nl<br class=""><br class="">You can reach the person managing the list at<br class="">fieldtrip-owner@science.ru.nl<br class=""><br class="">When replying, please edit your Subject line so it is more specific<br class="">than "Re: Contents of fieldtrip digest..."<br class=""><br class=""><br class="">Today's Topics:<br class=""><br class=""> 1. Call for new post docs in the Obleser lab | Auditory<br class=""> Cognition (Jonas Obleser)<br class=""> 2. Spectrum normalization (Wunderle, Thomas)<br class=""> 3. Interpolation and ft_channelrepair (Grazia Di Pisa)<br class=""> 4. Re: Interpolation and ft_channelrepair<br class=""> (=?UTF-8?Q?J=C3=B6rn_M._Horschig?=)<br class=""><br class=""><br class="">----------------------------------------------------------------------<br class=""><br class="">Message: 1<br class="">Date: Tue, 20 Oct 2015 06:07:51 -0500<br class="">From: Jonas Obleser <jonas@obleser.de><br class="">To: fieldtrip@science.ru.nl, AUDITORY@LISTS.MCGILL.CA<br class="">Subject: [FieldTrip] Call for new post docs in the Obleser lab |<br class="">AuditoryCognition<br class="">Message-ID: <71625964-4C4D-41CA-BEDB-171CEFBC87F8@obleser.de><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">The research group ?Auditory Cognition? (headed by Jonas Obleser; <br class="">auditorycognition.com), now at the University of L?beck, is seeking <br class="">to hire two Postdoctoral researchers, initially for 3 years, with the <br class="">option of extension.<br class=""><br class="">These positions fall into the larger framework of an ERC Consolidator <br class="">grant ?The listening challenge: How ageing brains adapt? (AUDADAPT) <br class="">recently awarded to Jonas Obleser, and will allow the joint <br class="">development of cognitive-neuroscience and psychological research <br class="">projects targeting adaptive control in the auditory modality of <br class="">middle-aged adults.<br class=""><br class="">Please download the full job advert here: <br class="">https://www.dropbox.com/s/sqc9c94ewl6bjjh/Ausschreibung_1031_15_Wiss%20MitarbeiterIn_Psychologie_EN.pdf?dl=1 <br class=""><br class=""><br class="">Note the application deadline of November 12 (postal receipt).<br class=""><br class="">http://auditorycognition.com<br class="">http://uni-luebeck.de<br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 2<br class="">Date: Tue, 20 Oct 2015 14:58:25 +0000<br class="">From: "Wunderle, Thomas" <thomas.wunderle@esi-frankfurt.de><br class="">To: "fieldtrip@science.ru.nl" <fieldtrip@science.ru.nl><br class="">Subject: [FieldTrip] Spectrum normalization<br class="">Message-ID:<br class=""><27E5CAD9145EEC41BB9B34C01716A198880C3CD4@UM-excdag-a02.um.gwdg.de><br class="">Content-Type: text/plain; charset="iso-8859-1"<br class=""><br class="">Hi all,<br class=""><br class="">I have a question about how FieldTrip normalizes the output of a <br class="">spectral analyses using ft_freqanalysis.<br class="">This is important for me, because I need the absolute value of the <br class="">spectrum for a publication.<br class=""><br class="">I was scanning through the code trying to find the normalization steps.<br class="">In particular, I don't understand the following steps:<br class=""><br class=""><br class="">1.) The tapers (e.g. hanning) are normalized prior to the <br class="">multiplication with the data by the following code:<br class="">tap = tap./norm(tap, 'fro'); (line 195 in the function <br class="">fr_specest_mtmfft)<br class="">What is the rationale of dividing through this special norm of the <br class="">taper itself?<br class=""><br class=""><br class="">2.) After calculating the fft, the specrum is scaled (i.e. in the <br class="">function fr_specest_mtmfft) by:<br class=""><br class="">dum = dum .* sqrt(2 ./ endnsample); (line 251)<br class=""><br class="">Later on, to get to the power spectrum, the abs of the spectrum is <br class="">taken and squared<br class=""><br class="">powdum = abs(spectrum(acttap,:,foiind(ifoi),acttboi)) .^2; (line 629 <br class="">in ft_frequanalysis)<br class=""><br class=""><br class=""><br class="">This puzzles me, because scanning through the literature, the power <br class="">spectrum should be given by<br class=""><br class="">(2*abs(X)^2) ./ (nfft^2) -> unit: V^2 <br class=""> (X=signal, nfft = number of data points)<br class=""><br class="">and the PSD by<br class=""><br class="">(2*abs(X)^2) ./ (nfft*Fs) -> unit: V^2/Hz (X=signal, <br class="">nfft = number of data points, Fs=sampling rate)<br class=""><br class=""><br class=""><br class="">Alternatively, one can get the amplitude spectrum from the power <br class="">spectrum by taking the square root of the PSD:<br class=""><br class="">sqrt(PSD) -> unit: V/sqrt(Hz)<br class=""><br class="">which is still different from<br class=""><br class="">(2*abs(X)) ./ (nfft*Fres) -> unit: V/Hz ? (Fres = <br class="">frequency resolution = Fs/nfft)<br class=""><br class=""><br class=""><br class="">The result of these equations are different from what comes out of <br class="">the FieldTrip code.<br class=""><br class="">So what is the unit of the spectrum in FieldTrip?<br class=""><br class="">What is the rationale of scaling the complex spectrum by sqrt(2/nfft) ?<br class=""><br class="">Furthermore, is the unit for the different spectral estimation <br class="">methods the same? That is, using mtmfft, mtmconvol, wavelet,...?<br class=""><br class="">It would help me a lot if someone could clarify these points.<br class=""><br class="">Thanks in advance,<br class="">Thomas<br class=""><br class=""><br class="">-----<br class="">Dr. Thomas Wunderle<br class="">Ernst Str?ngmann Institute (ESI) for <br class="">Neuroscience<http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">in Cooperation with Max Planck Society <br class=""><http://www.ncbi.nlm.nih.gov/feed/rss.cgi?ChanKey=PubMedNews><br class="">Deutschordenstrasse 46<br class="">60528 Frankfurt am Main, Germany<br class="">www.esi-frankfurt.de<http://www.esi-frankfurt.de/><br class="">thomas.wunderle@esi-frankfurt.de<mailto:thomas.wunderle@esi-frankfurt.de><br class="">Tel: +49 69 96769 516<br class="">Fax: +49 69 96769 555<br class=""><br class="">Sitz der Gesellschaft: Frankfurt am Main<br class="">Registergericht: Amtsgericht Frankfurt - HRB 84266<br class="">Gesch?ftsf?hrer: Prof. Dr. Pascal Fries<br class=""><br class=""><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <br class=""><http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151020/eebaf460/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 3<br class="">Date: Wed, 21 Oct 2015 10:30:07 +0200<br class="">From: Grazia Di Pisa <g.dipisa@gmail.com><br class="">To: fieldtrip@science.ru.nl<br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <A46541FB-BC67-44E7-A9CD-541FB39B7C54@gmail.com><br class="">Content-Type: text/plain; charset=utf-8<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m <br class="">getting the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should <br class="">specify the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code?<br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 cfg = [];<br class="">16 cfg_neighb = [];<br class="">17cfg_neighb.method= 'triangulation';<br class="">18 cfg.senstype = 'EEG';<br class="">19 cfg_neighb.layout = 'biosemi64.lay';<br class="">20 neighbours = ft_prepare_neighbours(cfg_neighb);<br class="">21<br class="">22 cfg.neighbours = neighbours;<br class="">23 cfg.channel = {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9',};<br class="">24 cfg.feedback = 'yes';<br class="">25 cfg.layout = 'biosemi64.lay';<br class="">26 lay = ft_prepare_layout(cfg);<br class="">27<br class="">28 sens = data;<br class="">29 sens.type = 'eeg';<br class="">30 sens.label= lay.label;<br class="">31 sens.chanpos = lay.pos;<br class="">32 sens.chanpos(:,3) = 0;<br class="">33<br class="">34 ft_neighbourplot(cfg, data)<br class="">35<br class="">36 cfg = [];<br class="">37 cfg.method = 'nearest';<br class="">38 cfg.badchannel = {'Fp1', 'Fp2', 'AF7', <br class="">'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours = neighbours;<br class="">40 cfg.trials = 'all';<br class="">41<br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">Message: 4<br class="">Date: Wed, 21 Oct 2015 10:42:03 +0200<br class="">From: "=?UTF-8?Q?J=C3=B6rn_M._Horschig?=" <jorn@artinis.com><br class="">To: "'FieldTrip discussion list'" <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] Interpolation and ft_channelrepair<br class="">Message-ID: <002b01d10bdc$5d377040$17a650c0$@artinis.com><br class="">Content-Type: text/plain;charset="UTF-8"<br class=""><br class="">Dear Grazia,<br class=""><br class="">try adding cfg.sens = sens; before line 42. That will set the sensor <br class="">definition so that ft_channelrepair knows the location of your <br class="">channels, which is otherwise not obvious.<br class=""><br class="">Best,<br class="">J?rn<br class=""><br class="">--<br class=""><br class="">J?rn M. Horschig, PhD, Software Engineer<br class="">Artinis Medical Systems | +31 481 350 980<br class=""><br class=""><blockquote type="cite" class="">-----Original Message-----<br class="">From: fieldtrip-bounces@science.ru.nl [mailto:fieldtrip-<br class="">bounces@science.ru.nl] On Behalf Of Grazia Di Pisa<br class="">Sent: Wednesday, October 21, 2015 10:30 AM<br class="">To: fieldtrip@science.ru.nl<br class="">Subject: [FieldTrip] Interpolation and ft_channelrepair<br class=""><br class="">Dear all,<br class=""><br class="">I?m trying to do interpolation to repair some bad channels, but I?m <br class="">getting<br class="">the error below.<br class="">I've read the ft_fetch_sens help but I don?t understand how I should <br class="">specify<br class="">the electrodes.<br class=""><br class="">Is there a file I should upload or am I missing something in my code?<br class="">Any help or hint is very much appreciated!<br class=""><br class="">thanks in advance,<br class="">~ grazia<br class=""><br class=""><br class="">//<br class="">Error using ft_fetch_sens (line 183)<br class="">no electrodes or gradiometers specified.<br class=""><br class="">Error in ft_channelrepair (line 115)<br class="">sens = ft_fetch_sens(cfg, data);<br class=""><br class="">Error in new_analyse_GDP (line 42)<br class="">ft_channelrepair(cfg, data)<br class="">//<br class=""><br class="">This is my code:<br class=""><br class="">15 cfg = [];<br class="">16 cfg_neighb = [];<br class="">17cfg_neighb.method= 'triangulation';<br class="">18 cfg.senstype = 'EEG';<br class="">19 cfg_neighb.layout = 'biosemi64.lay';<br class="">20 neighbours = ft_prepare_neighbours(cfg_neighb);<br class="">21<br class="">22 cfg.neighbours = neighbours;<br class="">23 cfg.channel = {'Fp1', 'Fp2', 'AF7', 'AF3','AF8','F7','F5','FT9',};<br class="">24 cfg.feedback = 'yes';<br class="">25 cfg.layout = 'biosemi64.lay';<br class="">26 lay = ft_prepare_layout(cfg);<br class="">27<br class="">28 sens = data;<br class="">29 sens.type = 'eeg';<br class="">30 sens.label= lay.label;<br class="">31 sens.chanpos = lay.pos;<br class="">32 sens.chanpos(:,3) = 0;<br class="">33<br class="">34 ft_neighbourplot(cfg, data)<br class="">35<br class="">36 cfg = [];<br class="">37 cfg.method = 'nearest';<br class="">38 cfg.badchannel = {'Fp1', 'Fp2', 'AF7', <br class="">'AF3','AF8','F7','F5','FT9'};<br class="">39 cfg.neighbours = neighbours;<br class="">40 cfg.trials = 'all';<br class="">41<br class="">42 ft_channelrepair(cfg, data)<br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""></blockquote><br class=""><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class="">End of fieldtrip Digest, Vol 59, Issue 17<br class="">*****************************************<br class=""></blockquote><br class=""><br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""></blockquote><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151021/7d31addb/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 5<br class="">Date: Wed, 21 Oct 2015 17:54:40 +0000<br class="">From: "Krebber, Martin" <martin.krebber@charite.de><br class="">To: FieldTrip discussion list <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] calculating behavioural-power correlation --<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>follow-up questions<br class="">Message-ID:<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span><B521CE9BFA72124F87A7D9CBEE13F9BE0DDCF182@s-mx14-mb0105.charite.de><br class="">Content-Type: text/plain; charset="iso-8859-1"<br class=""><br class="">Hi Xiaoming, hi Arjen,<br class=""><br class="">I've been encountering the same problem. I believe Xiaoming is right when he points out that the permutaion step shuffles data across conditions and that this introduces a negative bias in the distribution. I found the same thing when I correlated RT data with TFRs (absolute power). My distribution was shifted strongly to the left and, thus, not a single negative cluster was significant, but every positive one was.<br class=""><br class="">Xiaomings explanation made a lot of sense to me when I thought about it graphically: Imagine correlating two data vectors, one (x) ranging between .5 and1, the other (y) between 50 and and 100. When plotting this, one gets a cloud of dots on the upper left corner of the diagram. When you then switch the variable assignment of half of the data points (which is what the permutation step seems to do), these dots will now be be shifted to the lower right corner of the diagram. So no matter what the correlation in the original data, chances are that (given different scaling) after permutaion, you get a negative correlation.<br class=""><br class="">I am not 100% sure about this, so please let me know if I made a mistake.<br class=""><br class="">What I tried instead of the 'ft_statfun_correlationT' was using a custom made statfun in which I pass the RTs via the design matrix. With this, my results looked much better. I am not sure, but I guess this is because there is no shuffling between the two variables in this case.<br class=""><br class="">I would really like to know, what is the right way of doing this using just the FieldTrip functions. Is there a way to permute data within variables? I tried cfg.resampling = 'bootstrap', but this is not a permutation, as far as I know.<br class=""><br class=""><br class="">Thanks!<br class="">Martin<br class=""><br class=""><br class="">________________________________<br class="">Von: fieldtrip-bounces@science.ru.nl [fieldtrip-bounces@science.ru.nl]" im Auftrag von "Arjen Stolk [a.stolk8@gmail.com]<br class="">Gesendet: Dienstag, 20. Oktober 2015 08:03<br class="">An: FieldTrip discussion list<br class="">Betreff: Re: [FieldTrip] calculating behavioural-power correlation -- follow-up questions<br class=""><br class="">Hey Xiaoming,<br class=""><br class="">It's still pretty hard, for me, to guess on basis of that matlab output what is going on here and what you mean with 'shuffling design matrices', and how that shuffling 'biases the cluster distribution'. As you mention yourself, it could be due to various reasons, and you're open to suggestions and increasing your understanding. I'd therefore suggest to try to funnel the number of potential explanations by simulating what you're doing (using input data for which you know how it should behave), after you've read more about what the design matrix and monte carlo statistics are supposed to do. Perhaps the statistics section at the bottom of this page provides a good starting point: http://www.fieldtriptoolbox.org/walkthrough<br class=""><br class="">Hope that helps,<br class="">Arjen<br class=""><br class="">2015-10-19 15:56 GMT-07:00 Xiaoming Du <XDu@mprc.umaryland.edu<mailto:XDu@mprc.umaryland.edu>>:<br class="">For example, our power values ranged from 1 to 3 (after log transform); my behavioral data ranged from 20 to 90;<br class=""><br class="">by using above mentioned script, there are 14 negative clusters were reported in variable stat.<br class=""><br class="">stat =<br class=""><br class=""> prob: [30x50 double]<br class=""> posclusters: []<br class=""> posclusterslabelmat: [30x50 double]<br class=""> posdistribution: [1x1000 double]<br class=""> negclusters: [1x14 struct]<br class=""> negclusterslabelmat: [30x50 double]<br class=""> negdistribution: [1x1000 double]<br class=""> cirange: [30x50 double]<br class=""> mask: [30x50 logical]<br class=""> stat: [30x50 double]<br class=""> ref: [30x50 double]<br class=""> rho: [30x50 double]<br class=""> dimord: 'chan_freq'<br class=""> freq: [1x50 double]<br class=""> label: {30x1 cell}<br class=""> time: 2.5000<br class=""> cfg: [1x1 struct]<br class=""><br class="">However, the p values of those clusters (i.e., stat.negclusters.prob) are all ones. The smallest value in stat.negdistribution is way larger than the largest negative cluster t-sum. This could be real. However, it is more likely due to the shuffle between power and behavioral group. For example, design matrix [1 1 1 1 2 2 2 2; 1 2 3 4 1 2 3 4] was shuffled to [1 2 2 1 2 2 1 1; 1 2 3 4 1 2 3 4]. After each permutation, for some subjects, their power data was labeled as behavioral data and vice versa. Because of the scale difference between power and behavioral data, large negative correlations were generated by permutation. This further biased the cluster distribution.<br class="">My limited understanding is that, for correlation, each permutation should fix cfg.ivar and only shuffle half of the cfg.uvar. For example, permute design matrix [1 1 1 1 2 2 2 2; 1 2 3 4 1 2 3 4] to [1 1 1 1 2 2 2; 1 2 3 4 4 2 3 1]. THerefore, after permutation, one subject's power data corresponds to another subject's behavioral data.<br class=""><br class="">I am not good at statistics. It will be really appreciated if you have any suggestions or comments.<br class=""><br class="">Xiaoming<br class=""><br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class="">Arjen Stolk <a.stolk8@gmail.com<mailto:a.stolk8@gmail.com>> 10/19/2015 6:01 PM >>><br class=""></blockquote></blockquote></blockquote>Hey Xiaoming,<br class=""><br class="">Not sure if I understand, but shouldn't the directions of the correlations be independent of the scaling of the two variables? Looking at the code of ft_statfun_correlationT it doesn't seem the conversion from correlation to T value (tstat = rho*(sqrt(max(nunits)-2))/sqrt((1-rho^2))) would result in a direction change either. Perhaps you could try to first manually calculate a correlation between signal power and behavioral power, and see whether anything is behaving unexpectedly?<br class=""><br class="">Yours,<br class="">Arjen<br class=""><br class="">2015-10-19 14:25 GMT-07:00 Xiaoming Du <XDu@mprc.umaryland.edu<mailto:XDu@mprc.umaryland.edu>>:<br class="">Dear FieldTrip users,<br class="">This is Xiaoming from University of Maryland Baltimore. My current project requires to calculate behavioral-power correlation across subjects. Similar topic was discussed here early this year. http://mailman.science.ru.nl/pipermail/fieldtrip/2015-February/008953.html<br class="">According to the suggestions in above mentioned thread, I duplicate my power dataset and replace the power values at each time-frequency point with behavioral data. Therefore, those two datasets have same structure and dimension. I used the following script to test if there are significant clusters of correlations.<br class="">cfg = [];<br class="">cfg.parameter = 'powspctrm';<br class="">cfg.method = 'montecarlo';<br class="">cfg.statistic = 'ft_statfun_correlationT';<br class="">...<br class="">etc<br class="">...<br class="">design = zeros(2, n1 * 2); % n1 is the number of subjects.<br class="">design(1,1:n1) = 1;<br class="">design(1,(n1 + 1):(n1 * 2)) = 2;<br class="">design(2, :) = [[1:n1 ] [1 : n1]];<br class="">cfg.design = design;<br class=""><br class="">cfg.ivar = 1;<br class="">cfg.uvar = 2;<br class="">stat = ft_freqstatistics(cfg, dataBeh{:}, dataDX1{:});<br class="">However, it seems when each time the design matrix is permuted, FieldTrip is using the same method as for 'ft_statfun_depsamplesT', meaning cfg.uvar remains the same while cfg.ivar (1 or 2) is randomly assigned to each subject in design matrix. Although I confirmed this by uncommenting line 313 (i.e., tmpdesign = design(:,resample(i,:))) in ft_statistics_montecarlo.m which allows to display the permuted design matrix in command line, please correct me if this is not the case.<br class="">In my mind, this kind of permutation will cause trouble when dealing with correlation. For example, in my case, the behavioral data and power data have different scales. The power data are much larger than behavioral data in general. When assigning behavioral data into power group or vice versa, it will induce huge negative correlations between power and behavioral measurement. Therefore, no negative clusters will survive from permutation test.<br class="">Please let me know if I have mis-understanding or if I did anything wrong. Any suggestions will be highly appreciated!<br class="">Thanks.<br class="">Xiaoming<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<mailto:fieldtrip@donders.ru.nl><br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class=""><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151021/138520fa/attachment-0001.html><br class=""><br class="">------------------------------<br class=""><br class="">Message: 6<br class="">Date: Wed, 21 Oct 2015 11:22:09 -0700<br class="">From: Arjen Stolk <a.stolk8@gmail.com><br class="">To: FieldTrip discussion list <fieldtrip@science.ru.nl><br class="">Subject: Re: [FieldTrip] calculating behavioural-power correlation --<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>follow-up questions<br class="">Message-ID:<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span><CAAiwptSYuNU0hrqb79C9H8R+viiM5O=wugBrcX8QovNrE8m4Qg@mail.gmail.com><br class="">Content-Type: text/plain; charset="utf-8"<br class=""><br class="">Hi Martin,<br class=""><br class="">Thanks for thinking along. I've briefly tried to replicate/simplify the<br class="">situation depicted by you using the code below. It works as I would expect,<br class="">no matter whether one variable is scaled differently. But perhaps I'm not<br class="">fully capturing the issue, and something still goes awry. This is a<br class="">possibility because ft_statfun_correlationT has only recently been<br class="">implemented for a specific case, and was never really tested within<br class="">different situations (hence it's not well-documented on the wiki). Do you<br class="">think you could use this example code to replicate the situation you are<br class="">experiencing?<br class=""><br class="">Yours,<br class="">Arjen<br class=""><br class="">% simulate simple multiple subjects timelock structures<br class=""><br class="">data_brain = [];<br class=""><br class="">data_behav = [];<br class=""><br class="">for j=1:10<br class=""><br class="">data_brain{j}.avg = j; % increasing<br class=""><br class="">data_brain{j}.dimord = 'chan_time';<br class=""><br class="">data_brain{j}.time = 1;<br class=""><br class="">data_brain{j}.label = {'1'};<br class=""><br class=""><br class=""><br class="">data_behav{j} = data_brain{j};<br class=""><br class="">data_behav{j}.avg = data_brain{j}.avg*-1000+50; % add scaling difference<br class=""><br class="">end<br class=""><br class=""><br class=""><br class="">% compute statistics with correlationT<br class=""><br class="">cfg = [];<br class=""><br class="">cfg.method = 'montecarlo';<br class=""><br class="">cfg.statistic = 'ft_statfun_correlationT';<br class=""><br class="">cfg.numrandomization = 100;<br class=""><br class=""><br class=""><br class="">n1 = 10; % n1 is the number of subjects<br class=""><br class="">design = zeros(2, n1 * 2);<br class=""><br class="">design(1,1:n1) = 1;<br class=""><br class="">design(1,(n1 + 1):(n1 * 2)) = 2;<br class=""><br class="">design(2, :) = [1:n1 1:n1];<br class=""><br class="">cfg.design = design;<br class=""><br class=""><br class=""><br class="">cfg.ivar = 1;<br class=""><br class="">cfg.uvar = 2;<br class=""><br class="">stat = ft_timelockstatistics(cfg, data_brain{:}, data_behav{:});<br class=""><br class=""><br class=""><br class="">assert(isequal(stat.rho, -1));<br class=""><br class="">2015-10-21 10:54 GMT-07:00 Krebber, Martin <martin.krebber@charite.de>:<br class=""><br class=""><blockquote type="cite" class="">Hi Xiaoming, hi Arjen,<br class=""><br class="">I've been encountering the same problem. I believe Xiaoming is right when<br class="">he points out that the permutaion step shuffles data across conditions and<br class="">that this introduces a negative bias in the distribution. I found the same<br class="">thing when I correlated RT data with TFRs (absolute power). My distribution<br class="">was shifted strongly to the left and, thus, not a single negative cluster<br class="">was significant, but every positive one was.<br class=""><br class="">Xiaomings explanation made a lot of sense to me when I thought about it<br class="">graphically: Imagine correlating two data vectors, one (x) ranging between<br class="">.5 and1, the other (y) between 50 and and 100. When plotting this, one gets<br class="">a cloud of dots on the upper left corner of the diagram. When you then<br class="">switch the variable assignment of half of the data points (which is what<br class="">the permutation step seems to do), these dots will now be be shifted to the<br class="">lower right corner of the diagram. So no matter what the correlation in the<br class="">original data, chances are that (given different scaling) after permutaion,<br class="">you get a negative correlation.<br class=""><br class="">I am not 100% sure about this, so please let me know if I made a mistake.<br class=""><br class="">What I tried instead of the 'ft_statfun_correlationT' was using a custom<br class="">made statfun in which I pass the RTs via the design matrix. With this, my<br class="">results looked much better. I am not sure, but I guess this is because<br class="">there is no shuffling between the two variables in this case.<br class=""><br class="">I would really like to know, what is the right way of doing this using<br class="">just the FieldTrip functions. Is there a way to permute data within<br class="">variables? I tried cfg.resampling = 'bootstrap', but this is not a<br class="">permutation, as far as I know.<br class=""><br class=""><br class="">Thanks!<br class="">Martin<br class=""><br class=""><br class="">------------------------------<br class="">*Von:* fieldtrip-bounces@science.ru.nl [fieldtrip-bounces@science.ru.nl]"<br class="">im Auftrag von "Arjen Stolk [a.stolk8@gmail.com]<br class="">*Gesendet:* Dienstag, 20. Oktober 2015 08:03<br class="">*An:* FieldTrip discussion list<br class="">*Betreff:* Re: [FieldTrip] calculating behavioural-power correlation --<br class="">follow-up questions<br class=""><br class="">Hey Xiaoming,<br class=""><br class="">It's still pretty hard, for me, to guess on basis of that matlab output<br class="">what is going on here and what you mean with 'shuffling design matrices',<br class="">and how that shuffling 'biases the cluster distribution'. As you mention<br class="">yourself, it could be due to various reasons, and you're open to<br class="">suggestions and increasing your understanding. I'd therefore suggest to try<br class="">to funnel the number of potential explanations by simulating what you're<br class="">doing (using input data for which you know how it should behave), after<br class="">you've read more about what the design matrix and monte carlo statistics<br class="">are supposed to do. Perhaps the statistics section at the bottom of this<br class="">page provides a good starting point:<br class="">http://www.fieldtriptoolbox.org/walkthrough<br class=""><br class="">Hope that helps,<br class="">Arjen<br class=""><br class="">2015-10-19 15:56 GMT-07:00 Xiaoming Du <XDu@mprc.umaryland.edu>:<br class=""><br class=""><blockquote type="cite" class="">For example, our power values ranged from 1 to 3 (after log transform);<br class="">my behavioral data ranged from 20 to 90;<br class=""><br class="">by using above mentioned script, there are 14 negative clusters were<br class="">reported in variable stat.<br class=""><br class="">stat =<br class=""><br class=""> prob: [30x50 double]<br class=""> posclusters: []<br class=""> posclusterslabelmat: [30x50 double]<br class=""> posdistribution: [1x1000 double]<br class=""> negclusters: [1x14 struct]<br class=""> negclusterslabelmat: [30x50 double]<br class=""> negdistribution: [1x1000 double]<br class=""> cirange: [30x50 double]<br class=""> mask: [30x50 logical]<br class=""> stat: [30x50 double]<br class=""> ref: [30x50 double]<br class=""> rho: [30x50 double]<br class=""> dimord: 'chan_freq'<br class=""> freq: [1x50 double]<br class=""> label: {30x1 cell}<br class=""> time: 2.5000<br class=""> cfg: [1x1 struct]<br class=""><br class="">However, the p values of those clusters (i.e., stat.negclusters.prob) are<br class="">all ones. The smallest value in stat.negdistribution is way larger than<br class="">the largest negative cluster t-sum. This could be real. However, it is more<br class="">likely due to the shuffle between power and behavioral group. For example,<br class="">design matrix [1 1 1 1 2 2 2 2; 1 2 3 4 1 2 3 4] was shuffled to [1 2 2 1<br class="">2 2 1 1; 1 2 3 4 1 2 3 4]. After each permutation, for some subjects,<br class="">their power data was labeled as behavioral data and vice versa. Because of<br class="">the scale difference between power and behavioral data, large negative<br class="">correlations were generated by permutation. This further biased the cluster<br class="">distribution.<br class="">My limited understanding is that, for correlation, each permutation<br class="">should fix cfg.ivar and only shuffle half of the cfg.uvar. For example,<br class="">permute design matrix [1 1 1 1 2 2 2 2; 1 2 3 4 1 2 3 4] to [1 1 1 1 2 2<br class="">2; 1 2 3 4 4 2 3 1]. THerefore, after permutation, one subject's power data<br class="">corresponds to another subject's behavioral data.<br class=""><br class="">I am not good at statistics. It will be really appreciated if you have<br class="">any suggestions or comments.<br class=""><br class="">Xiaoming<br class=""><br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class=""><blockquote type="cite" class=""><blockquote type="cite" class="">Arjen Stolk <a.stolk8@gmail.com> 10/19/2015 6:01 PM >>><br class=""></blockquote></blockquote></blockquote>Hey Xiaoming,<br class=""><br class="">Not sure if I understand, but shouldn't the directions of the<br class="">correlations be independent of the scaling of the two variables? Looking at<br class="">the code of ft_statfun_correlationT it doesn't seem the conversion from<br class="">correlation to T value (tstat = rho*(sqrt(max(nunits)-2))/sqrt((1-rho^2)))<br class="">would result in a direction change either. Perhaps you could try to first<br class="">manually calculate a correlation between signal power and behavioral power,<br class="">and see whether anything is behaving unexpectedly?<br class=""><br class="">Yours,<br class="">Arjen<br class=""><br class="">2015-10-19 14:25 GMT-07:00 Xiaoming Du <XDu@mprc.umaryland.edu>:<br class=""><br class=""><blockquote type="cite" class="">Dear FieldTrip users,<br class="">This is Xiaoming from University of Maryland Baltimore. My current<br class="">project requires to calculate behavioral-power correlation across subjects.<br class="">Similar topic was discussed here early this year.<br class="">http://mailman.science.ru.nl/pipermail/fieldtrip/2015-February/008953.html<br class="">According to the suggestions in above mentioned thread, I duplicate my<br class="">power dataset and replace the power values at each time-frequency point<br class="">with behavioral data. Therefore, those two datasets have same structure and<br class="">dimension. I used the following script to test if there are significant<br class="">clusters of correlations.<br class="">cfg = [];<br class="">cfg.parameter = 'powspctrm';<br class="">cfg.method = 'montecarlo';<br class="">cfg.statistic = 'ft_statfun_correlationT';<br class="">...<br class="">etc<br class="">...<br class="">design = zeros(2, n1 * 2); % n1 is the number of subjects.<br class="">design(1,1:n1) = 1;<br class="">design(1,(n1 + 1):(n1 * 2)) = 2;<br class="">design(2, :) = [[1:n1 ] [1 : n1]];<br class="">cfg.design = design;<br class=""><br class="">cfg.ivar = 1;<br class="">cfg.uvar = 2;<br class="">stat = ft_freqstatistics(cfg, dataBeh{:}, dataDX1{:});<br class="">However, it seems when each time the design matrix is permuted,<br class="">FieldTrip is using the same method as for 'ft_statfun_depsamplesT', meaning<br class="">cfg.uvar remains the same while cfg.ivar (1 or 2) is randomly assigned to<br class="">each subject in design matrix. Although I confirmed this by uncommenting<br class="">line 313 (i.e., tmpdesign = design(:,resample(i,:))) in<br class="">ft_statistics_montecarlo.m which allows to display the permuted design<br class="">matrix in command line, please correct me if this is not the case.<br class="">In my mind, this kind of permutation will cause trouble when dealing<br class="">with correlation. For example, in my case, the behavioral data and power<br class="">data have different scales. The power data are much larger than behavioral<br class="">data in general. When assigning behavioral data into power group or vice<br class="">versa, it will induce huge negative correlations between power and<br class="">behavioral measurement. Therefore, no negative clusters will survive from<br class="">permutation test.<br class="">Please let me know if I have mis-understanding or if I did anything<br class="">wrong. Any suggestions will be highly appreciated!<br class="">Thanks.<br class="">Xiaoming<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class=""></blockquote><br class=""><br class=""></blockquote><br class=""></blockquote>-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151021/c798d778/attachment.html><br class=""><br class="">------------------------------<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">fieldtrip@donders.ru.nl<br class="">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class=""><br class="">End of fieldtrip Digest, Vol 59, Issue 18<br class="">*****************************************<br class=""></div></blockquote></div><br class=""></div></div></body></html>