<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"\@SimSun";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi Gerome and Tzvetan,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Actually, you are referencing to the linked mastoids if EXG1 and EXG2 represent the mastoid channels as you specified them in cfg.refchannel. cfg.channel only
specifies on what channels to apply the re-referencing on.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">At a first glance I can see that in Fieldtrip you are applying a band-stop filter between 92 and 108 Hz, while in BVA you are applying a band-stop filter at
100Hz with a bandwidth of 8. I am not that familiar with BVA but does that mean plus AND minus 8Hz or 8Hz around 100Hz, so 96 – 104Hz?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Also, in BVA you are applying a notch filter at 50Hz, in Fieldtrip you are applying a DFT filter at 50, 100 and 150Hz. A notch filter is a band-stop filter
with a narrow band (usually 1-2Hz). DFT filters try to fit and subtract sinusoids at the specified frequency. To make it more similar to BVA I would suggest a band-stop filter between 49 and 51Hz.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Furthermore, in Fieldtrip you are applying Butterworth filters, is this the same filter-type applied in BVA? Additionally, in BVA your baseline is from -200
to -5ms, while in Fieldtrip you specify: </span>cfg.baselinewindow = [-0.2 -0.05];
<span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">which is from -200 to -50ms.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I hope this can help you make the analysis more similar. Could you perhaps make a picture of a result that you find is different (e.g. ERPs calculated in Fieldtrip
and BVA)?. That could help identifying the exact problem if this did not help you.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Best,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Jim<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> fieldtrip-bounces@science.ru.nl [mailto:fieldtrip-bounces@science.ru.nl]
<b>On Behalf Of </b>Tzvetan Popov<br>
<b>Sent:</b> donderdag 13 november 2014 11:14<br>
<b>To:</b> FieldTrip discussion list<br>
<b>Subject:</b> Re: [FieldTrip] Redoing Brain Vision Analyzer in Field Trip<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Dear Gerome,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">one thing you could change is keep the reference the same in both cases. So far you are referencing the data against the mastoids in BVA and a common reference in FieldTrip.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">This is a potential source of variance in the waveforms you’ll get. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">best<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">tzvetan<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<div>
<p class="MsoNormal">Hello All, <o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I would like some help to replicate an analysis done in brain vision analyzer in field trip. Here is the details of the analysis in brain vision. I have done most of the code but I get almost completely different results which is puzziling.
Maybe some of the filteres are not defined properly- so if anyone has some time I would appreciate the help :-)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks alot,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Gerome<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Here is the brain vision sequence I am trying to replicate<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Read in Data<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Reference to mastiod<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - 100Hz Band rejection filter - bandwidth 8, slope 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - High cuttoff (low pass filter at 45Hz with a slope of 24dB/oct)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Low cuttoff (high pass filter at 0.53Hz, time constant 0.3002923, 24db/ oct)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Notch Filter at 50Hz<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Segmentation<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Baseline Correction -200 ms to - 5 ms<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> - Averaging segments<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Here is my Field trip code: <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">%rereference all trials<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.reref = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.continuous = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">% reference what channels<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.channel = 'all';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">%choose the reference channel<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.refchannel = {'EXG1' 'EXG2'};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">% cfg.dataset = 'CVP1S1.bdf';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.continuous = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">%<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">%bandfilter<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.bsfilter = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.bsfiltord = 2;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">% cfg.bsfiltdir = 'onepass-zerophase';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.bsfreq = [92 108];<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">% <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">%DFT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.dftfilter = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">% IIR Filters<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.lpfilter = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.lpfilttype = 'but';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">% cfg.lpfiltdir = 'onepass-zerophase';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.lpfreq = 45;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.hpfilter = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.hpfilttype = 'but';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">% cfg.hpfiltdir = 'onepass-zerophase';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.hpfreq = 0.53;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">data = ft_preprocessing(cfg);<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">% Define Trial<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg = [];<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.trialdef.eventtype = 'STATUS';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.dataset = eeg_file;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.trialdef.eventvalue = 65522;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.trialdef.prestim = 0.2;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.trialdef.poststim = 0.3;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg = ft_definetrial(cfg);<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">data_vep = ft_redefinetrial(cfg, data);<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">% Baseline Correction<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg = [];<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.demean = 'yes';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.baselinewindow = [-0.2 -0.05];<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">data_vep = ft_preprocessing(cfg, data_vep);<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">%Averaging<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg=[];<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">avg_data = ft_timelockanalysis(cfg, data_vep);<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.ylim = [-10 10];<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cfg.channel = 'O1';<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ft_singleplotER(cfg,avg_data);<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">-- <o:p></o:p></p>
<div>
<div>
<p class="MsoNormal">_____________________<o:p></o:p></p>
<div>
<p class="MsoNormal">Gerome Manson PhD Student<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Laboratoire de Neuroscience Cognitives- Aix Marseille Université<br>
Perceptual Motor Behavior lab- University of Toronto<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="mailto:gerome.manson@utoronto.ca" target="_blank">gerome.manson@utoronto.ca</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><a href="mailto:gerome.manson@univ-amu.fr" target="_blank">gerome.manson@univ-amu.fr</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">416-978-6109<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal">_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>