<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">whoopsy, just realized that you use
singlesphere, not singleshell. I was talking about singleshell ;)<br>
<br>
On 3/11/2013 9:47 AM, "Jörn M. Horschig" wrote:<br>
</div>
<blockquote cite="mid:513D9A1C.2040009@donders.ru.nl" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<div class="moz-cite-prefix">Dear Mushfa,<br>
<br>
I never tried what you did, however from what I know of the
matter and from what Robert wrote, I can see that you made some
errors in your script that I can point to. I can, however,
neither guarantee that this solves all or any of your problems
nor that it's all correct what you are doing. So, I can only
tell you some minor things:<br>
1) When defining trial, setting cfg.trialdef.eventtype to '?'
should actually only be done if you do not know the eventtype.
Also note that it might be better to use all your trials rather
than a subset - better signal-to-noise that way ;)<br>
<br>
2) No reason to segment the brain twice - this operation should
be done only once and the segmented brain you base the two
leadfield-matrices on *must* be exactly the same. So you can
save time by removing the computation of vol2.<br>
<br>
3) Specifying only one frequency might not be a good idea,
because of imperfect sampling.<br>
<br>
4) 'singlesphere' is created specifically for MEG data. BEM
models are good for EEG, not sure how good they are for MEG, but
the dipoli implementation only works under unix (Linux/Mac) and
openMEEG has to be downloaded manually first to make it work.<br>
<br>
5) As Robert said, units are not always properly controlled for
using FieldTrip. You can use ft_convert_units to make them
match.<br>
<br>
6) Please check again how 'grid' is defined, e.g. the .dim field
shouldn't be concatenated.<br>
<br>
7) You do not only need to concatenate the grids, but also the
data, hope you did that ;) But I would really rather
concatenated the two freq structures you used to compute the
sourcemodel/leadfield than to compute it newly. If you
concatenate you can be sure that the channel order is
persistent.<br>
<br>
<br>
If you keep on having problems, I would suggest to use
breakpoints or <small><i>dbstop if error</i></small> and then
step through the function to see where the error comes from,
i.e. why there are NaNs or Infs in (lf*Cf*lf), and start off
with looking at if there are NaNs or Infs in lf or Cf already.
That's all I can help you with, cause that's were my knowledge
ends ;) Maybe someone else is wiser than I am. Good luck!<br>
<br>
Best,<br>
Jörn<br>
<br>
On 3/8/2013 4:12 PM, Mushfa Yousuf wrote:<br>
</div>
<blockquote
cite="mid:CAFy9vnLSPQDRaCQs8ajkSLfu_Ww7ZMyFRSoNVR3ZfmGp3Zb_BQ@mail.gmail.com"
type="cite">
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Hello;</span></p>
<p class="MsoNormal" style="line-height:200%"> <span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">I have few
problems while calculating combine “MEG & EEG” source
reconstruction and I need your help to troubleshoot the
problems.</span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">I have a
data from a Parkinson’s disease patient who underwent DBS
and we recorded 306 MEG channels and 64 EEG channels
simultaneously, first I have calculated the leadfield
individually for both EEG and MEG. </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">I have
defined all the calculations step by step so that I can
describe my problems very precisely. </span></p>
<p class="MsoNormal" style="line-height:200%"><b><u><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">CALCULATE
THE LEADFIELD FOR MEG CHANNELS:-</span></u></b></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">// I have
defined the Pre-Stimulus data i.e. Stim-off data in the
original code.</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">Reading raw
data set for MEG channels</span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"> </span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg =
[]; </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.dataset = ‘</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">ap_rest_stimon_bi_2_5ma.fif'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.eventtype = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'?'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.triallength = 1;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.poststim = 88;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = ft_definetrial(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trl = cfg.trl(1:88,:); // Read 88 trials
</span><span style="font-size:12.0pt;font-family:"Courier
New"" lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.channel = {</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'MEG'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">}; </span><b><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">// In order to
read raw data only for MEG channels</span></b><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:forestgreen" lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">dataFIC12 = ft_preprocessing(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">After that I
have defined the <u>time of window of interest</u> for the
MEG raw data set</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.toilim = [0 0.5]; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">dataPost = ft_redefinetrial(cfg, dataFIC12);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">and
calculate the <u>Power Spectral density matrix</u> for the
MEG channels data set</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.method = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'mtmfft'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.output = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'powandcsd'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.tapsmofrq = 2;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.foilim = [129 129]; </span><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US">// Frequency of interest is 129 Hz the
stimulation frequency </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.pad = 1;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">freqPost = ft_freqanalysis(cfg, dataPost);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Now next
step for me is to choose the method for the Head Model. I
have noticed that Singlesphere is the common method which
can be used for both EEG and MEG. So I have selected the
singlesphere method for the calculation of volume conduction
model.</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">mri = ft_read_mri(‘</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">ms1879865-0004-00001-000192-01.img'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">); // READ MRI</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.write = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'no'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.coordsys = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'ctf'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">segmentedmri1 = ft_volumesegment(cfg, mri); //
Calculate the segmented mri</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.method = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'singlesphere'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">; // </span><b><span
style="font-size:10pt;font-family:'Times New Roman',serif"
lang="EN-US">METHOD</span></b><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">vol1 = ft_prepare_headmodel(cfg,
segmentedmri1);</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US">Now next step is to calculate the LEADFIELD </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.grad = freqPost.grad;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.vol = vol1;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.reducerank = 2; // MEG</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.channel = {</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'MEG'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">}; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.grid.resolution = 1; // I choose a
resolution 1</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10.0pt;font-family:"Courier
New";color:forestgreen" lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">[grid1] = ft_prepare_leadfield(cfg);</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">disp(grid1)</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> xgrid: [-7 -6 -5 -4 -3 -2
-1 0 1 2 3 4 5 6 7]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> ygrid: [-7 -6 -5 -4 -3 -2
-1 0 1 2 3 4 5 6 7 8]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> zgrid: [-4 -3 -2 -1 0 1 2
3 4 5 6 7 8 9 10 11]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> dim: [15 16 16]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> pos: [3840x3 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> unit: 'cm'</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> inside: [1x2037 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> outside: [1x1803 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> cfg: [1x1 struct]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> leadfield: {1x3840 cell}</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoListParagraph" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:Wingdings" lang="EN-US">Ø<span
style="font-size:7pt;font-family:'Times New Roman'"> </span></span><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">It has to be
noticed that Inside Dipole is 2037.</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoNormal" style="line-height:200%"><b><u><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">CALCULATE
THE LEADFIELD FOR EEG CHANNELS:-</span></u></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">Reading raw
data set for EEG channels</span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"> </span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg =
[]; </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.dataset = ‘</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">ap_rest_stimon_bi_2_5ma.fif'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.eventtype = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'?'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.triallength = 1;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.poststim = 88;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = ft_definetrial(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trl = cfg.trl(1:88,:); // Read 88 trials
</span><span style="font-size:12.0pt;font-family:"Courier
New"" lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.channel = {</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'EEG'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">}; </span><b><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">// In order to
read raw data only for EEG channels</span></b><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:forestgreen" lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">dataFIC12 = ft_preprocessing(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">After that I
have defined the <u>time of window of interest</u> for the
EEG raw data set</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.toilim = [0 0.5]; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">dataPost = ft_redefinetrial(cfg, dataFIC12);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">and
calculate the <u>Power Spectral density matrix</u> for the
EEG channels data set</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.method = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'mtmfft'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.output = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'powandcsd'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.tapsmofrq = 2;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.foilim = [129 129]; </span><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US">// Frequency of interest is 129 Hz</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.pad = 1;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">freqPost = ft_freqanalysis(cfg, dataPost);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">HEAD MODEL FOR
EEG</span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">mri = ft_read_mri(‘</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">ms1879865-0004-00001-000192-01.img'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.write = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'no'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.coordsys = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'ctf'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">segmentedmri2 = ft_volumesegment(cfg, mri);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.method = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'singlesphere'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">; // METHOD</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">vol2 = ft_prepare_headmodel(cfg,
segmentedmri2);</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US">LEADFIELD:</span></u><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"></span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = [];</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.elec = freqPost.elec; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.vol = vol2;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.reducerank = 3; </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:forestgreen" lang="EN-US">% 3 for EEG</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.channel = {</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'EEG'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">};</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.grid.resolution = 9.54; </span><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US">// Increase the resolution so that the inside
dipoles for </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12pt;font-family:'Times New Roman',serif"
lang="EN-US">
EEG
should be close to the inside dipoles of
MEG </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">[grid2] = ft_prepare_leadfield(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">disp(grid2)</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> xgrid: [1x16 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> ygrid: [1x16 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> zgrid: [1x13 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> dim: [16 16 13]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> pos: [3328x3 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> unit: 'mm'</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> inside: [1x2038 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> outside: [1x1290 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> cfg: [1x1 struct]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> leadfield: {1x3328 cell}</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoListParagraph" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:Wingdings" lang="EN-US">Ø<span
style="font-size:7pt;font-family:'Times New Roman'"> </span></span><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">It has to be
noticed that Inside Dipole is 2038.</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><b><u><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Concatenate
of MEG and EEG LEADFIELD:-</span></u></b></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Now I have
the individual LEADFIELD of both MEG and EEG grid1 and grid2
respectively.</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">disp(grid1)</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> xgrid: [-7 -6 -5 -4 -3 -2
-1 0 1 2 3 4 5 6 7]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> ygrid: [-7 -6 -5 -4 -3 -2
-1 0 1 2 3 4 5 6 7 8]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> zgrid: [-4 -3 -2 -1 0 1 2
3 4 5 6 7 8 9 10 11]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> dim: [15 16 16]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> pos: [3840x3 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> unit: 'cm'</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> inside: [1x2037 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> outside: [1x1803 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> cfg: [1x1 struct]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> leadfield: {1x3840 cell}</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">disp(grid2)</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> xgrid: [1x16 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> ygrid: [1x16 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> zgrid: [1x13 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> dim: [16 16 13]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> pos: [3328x3 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> unit: 'mm'</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> inside: [1x2038 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> outside: [1x1290 double]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> cfg: [1x1 struct]</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> leadfield: {1x3328 cell}</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Few things to be noticed here </span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoListParagraphCxSpFirst"
style="margin-bottom:0.0001pt"><span
style="font-size:10.0pt;font-family:Wingdings" lang="EN-US">Ø<span
style="font-size:7pt;font-family:'Times New Roman'"> </span></span><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">unit is different for both MEG
‘cm’ and EEG ‘mm’</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><span
style="font-size:10.0pt;font-family:Wingdings" lang="EN-US">Ø<span
style="font-size:7pt;font-family:'Times New Roman'"> </span></span><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">I have used the same method for
calculating the Volume for both MEG & EEG =>
“Singlesphere” </span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:Wingdings" lang="EN-US">Ø<span
style="font-size:7pt;font-family:'Times New Roman'"> </span></span><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Concatenating the “pos” vector of
EEG and MEG and the other vectors</span></b><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"></span></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US"> </span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.pos =
cat(1,grid1.pos,grid2.pos);</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.xgrid = [grid1.xgrid
grid2.xgrid];</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.ygrid = [grid1.ygrid
grid2.ygrid];</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.zgrid = [grid1.zgrid
grid2.zgrid];</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.inside = [grid1.inside
grid2.inside];</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.outside = [grid1.outside
grid2.outside];</span></b></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"" lang="EN-US">Grid.leadfield = [grid1.leadfield
grid2.leadfield];</span></b></p>
<p class="MsoListParagraphCxSpLast"
style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New"">Grid.dim = [grid1.dim grid2.dim];</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:10.0pt;font-family:"Courier
New""> </span></b></p>
<p class="MsoListParagraph" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif""> </span></p>
<p class="MsoNormal" style="line-height:200%"><b><u><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Combined
Source Reconstruction:-</span></u></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg =
[]; </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.dataset = ‘</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">ap_rest_stimon_bi_2_5ma.fif'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">; </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.eventtype = </span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'?'</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.triallength = 1;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trialdef.poststim = 88;</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg = ft_definetrial(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.trl = cfg.trl(1:88,:); // Read 88 trials
</span><span style="font-size:12.0pt;font-family:"Courier
New"" lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">cfg.channel = {</span><span
style="font-size:10.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'MEG',’EEG’</span><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">}; </span><b><span
style="font-family:"Times New
Roman","serif"" lang="EN-US">// In order to
read raw data for MEG & EEG channels</span></b><span
style="font-family:'Courier New'" lang="EN-US"> </span><span
style="font-family:"Courier
New";color:forestgreen" lang="EN-US"> </span><span
style="font-family:"Courier New"" lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US">dataFIC12 = ft_preprocessing(cfg);</span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10pt;font-family:'Courier New'"
lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">After that
as previously defined, I have calculated the Spectral
density. </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Now I have
to calculate the inverse solution for both EEG & MEG
data trials.</span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">Source
Analysis: without contrasting condition:-</span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><u><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US"> </span></u></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:10.0pt;font-family:"Courier
New";color:forestgreen" lang="EN-US"> </span><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg
= []; </span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.method
= </span><span
style="font-size:9.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'dics'</span><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">;</span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.frequency
= 129; </span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.grid
= Grid; // Concatenated grid</span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.vol
= vol1; // As vol1 = vol2 : both singlesphere</span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.powmethod
= </span><span
style="font-size:9.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'trace'</span><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">;</span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.dics.projectnoise
= </span><span
style="font-size:9.0pt;font-family:"Courier
New";color:#a020f0" lang="EN-US">'yes'</span><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">;</span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">cfg.dics.lambda
= 0.0001;</span><span
style="font-size:9.0pt;font-family:"Courier New""
lang="EN-US"></span></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:9pt;font-family:'Courier New'" lang="EN-US">sourcePost_nocon
= ft_sourceanalysis(cfg, freqPost); </span><span
style="font-size:8pt;font-family:'Courier New'" lang="EN-US">//
</span><b><span style="font-size:9pt;font-family:'Courier
New'" lang="EN-US">freqpost contains both MEG </span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:9pt;font-family:'Courier New'"
lang="EN-US">
and
EEG row data</span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><b><span
style="font-size:9pt;font-family:'Courier New'"
lang="EN-US">
</span></b><b><span
style="font-size:9.0pt;font-family:"Courier
New"" lang="EN-US"></span></b></p>
<p class="MsoNormal" style="margin-bottom:0.0001pt"><span
style="font-size:12.0pt;font-family:"Courier New""
lang="EN-US"> </span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">I have
received the following error:</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif"" lang="EN-US">scanning grid
2037/4075</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">???
Error using ==> svd</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">Input
to SVD must not contain NaN or Inf.</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US"> </span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">Error
in ==> beamformer_dics>pinv at 568</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">
[U,S,V] = svd(A,0);</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US"> </span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">Error
in ==> beamformer_dics at 314</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">
filt = pinv(lf' * invCf * lf) * lf' * invCf; %
Gross eqn. 3,</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">
use PINV/SVD to cover rank deficient leadfield</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US"> </span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">Error
in ==> ft_sourceanalysis at 595</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">
dip(i) = beamformer_dics(grid, sens, vol, [],
squeeze(Cf(i,:,:)),</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">
optarg{:});</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US"> </span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">Error
in ==> code at 227</span></p>
<p class="MsoNoSpacing"><span
style="font-size:12.0pt;font-family:"Times New
Roman","serif";color:red" lang="EN-US">sourcePost_nocon
= ft_sourceanalysis(cfg, freqPost);</span></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US"> </span></p>
<p class="MsoListParagraphCxSpFirst" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:Wingdings"
lang="EN-US">Ø<span
style="font-size:7pt;line-height:normal;font-family:'Times
New Roman'"> </span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">If you
notice that till 2037 grid it represents the MEG grid and
after that EEG grid is concatenated. i.e 2037 +2038 = 4075
(which also represent inside dipole).</span></p>
<p class="MsoListParagraphCxSpMiddle" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:Wingdings"
lang="EN-US">Ø<span
style="font-size:7pt;line-height:normal;font-family:'Times
New Roman'"> </span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Secondly I
try to use other method for the head modeling i.e openmeeg,
dipoli and bemcp method which leads me to different errors
which I was unable to trouble shoot.</span></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-left:72.0pt;line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:Symbol"
lang="EN-US">·<span
style="font-size:7pt;line-height:normal;font-family:'Times
New Roman'"> </span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Openmeeg = </span><span
style="font-size:9.0pt;line-height:200%;font-family:"Lucida
Console";color:red;background:white" lang="EN-US">OpenMEEG
binaries
are not correctly installed</span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US"></span></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-left:72.0pt;line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:Symbol"
lang="EN-US">·<span
style="font-size:7pt;line-height:normal;font-family:'Times
New Roman'"> </span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif";background:white"
lang="EN-US">Dipoli =</span><span
style="font-size:9.0pt;line-height:200%;font-family:"Lucida
Console";background:white" lang="EN-US"> <span
style="color:red">vol.mat file was missing after
calculating the headmodel</span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US"></span></p>
<p class="MsoListParagraphCxSpMiddle"
style="margin-left:72.0pt;line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:Symbol"
lang="EN-US">·<span
style="font-size:7pt;line-height:normal;font-family:'Times
New Roman'"> </span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif";background:white"
lang="EN-US">Bemcp =</span><span
style="font-size:9.0pt;line-height:200%;font-family:"Lucida
Console";background:white" lang="EN-US"> <span
style="color:red">vol.mat file contains NaN fiels which
leads to wrong</span></span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US"></span></p>
<p class="MsoListParagraphCxSpLast"
style="margin-left:72.0pt;line-height:200%"><span
style="font-size:9.0pt;line-height:200%;font-family:"Lucida
Console";color:red;background:white" lang="EN-US">
leadfield computation.</span><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US"></span></p>
<p class="MsoNormal" style="line-height:200%"><b><u><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">I
apologize for the long email. I wanted to show the steps
which I used so that it is easier to debug the problem.</span></u></b></p>
<p class="MsoNormal" style="line-height:200%"><b><u><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US"> </span></u></b></p>
<p class="MsoNormal" style="line-height:200%"><span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Thank you;</span></p>
<p class="MsoNormal" style="line-height:200%"> <span
style="font-size:12.0pt;line-height:200%;font-family:"Times
New Roman","serif"" lang="EN-US">Mushfa
Yousuf<br>
<br>
<br>
</span></p>
<br>
<div class="gmail_quote">On Fri, Feb 8, 2013 at 2:36 PM, Robert
Oostenveld <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:r.oostenveld@donders.ru.nl" target="_blank">r.oostenveld@donders.ru.nl</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">Dear Mushfa
<div><br>
</div>
<div>In principle it can be made to work, but you'll have
to do a little bit of work outside fieldtrip. You would
start with</div>
<div>
<div>1) ft_prepare_headmodel and ft_prepare_leadfield
for the MEG part</div>
</div>
<div>
<div>2) ft_prepare_headmodel and ft_prepare_leadfield
for the EEG part</div>
</div>
<div>Subsequently you would (manually) concatenate the
leadfields of the MEG and EEG. Here you have to make
sure that the channel order is consistent with that in
your combined data.</div>
<div>The resulting structure with combined leadfields can
be passed into ft_sourceanalysis as cfg.grid.</div>
<div><br>
</div>
<div>However, I foresee a problem with the scaling. At
this moment FieldTrip is not (yet) very precise in
maintaining the correct units throughout. Consequently
you might get different units for the EEG and MEG
contribution to each dipole, causing weird cancelation
effects. You'll just have to check carefully that the
leadfield units and the data units consistently scaled
between megmag, megplanar and eeg. </div>
<div><br>
</div>
<div>This brings me to another known problem for the
combined vectorview magnetometer and planar gradiometer
channels: at this moment the forward computation
computes both in T (or fT), whereas the data is
represented as T for the magnetometers and T/m for the
gradiometers. So the gradiometer leadfields are off by a
factor "1/baseline", which is approximately 1/70 if I
recall correctly. Many neuromag users use only the
planar channels, and therefore don't notice. We are
working on fixing the neuromag channel/gradiometer
units. That is not as simple as it appears, because
conversions of the geometrical dimensions from m to cm
or mm also affect the gradiometer baseline. So
converting the data from m to mm should result in the
planar channel data getting 1000 times smaller, whereas
the magnetometer channels shoudl stay the same. See <a
moz-do-not-send="true"
href="http://bugzilla.fcdonders.nl/show_bug.cgi?id=963"
target="_blank">http://bugzilla.fcdonders.nl/show_bug.cgi?id=963</a>
and <a moz-do-not-send="true"
href="http://bugzilla.fcdonders.nl/show_bug.cgi?id=963"
target="_blank">http://bugzilla.fcdonders.nl/show_bug.cgi?id=963</a> if
you want to follow this.</div>
<div><br>
</div>
<div>best regards</div>
<div>Robert</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div><br>
<div>
<div>
<div class="h5">
<div>On 28 Jan 2013, at 20:20, Mushfa Yousuf
wrote:</div>
<br>
</div>
</div>
<blockquote type="cite">
<div>
<div class="h5">
<pre style="line-height:21px;color:rgb(68,68,68);font-size:15px;white-space:normal">Hello Fieldtrippers,
I am currently doing source analysis using beamformer method in fieldtrip.
The data i work on is a simultaneous measurement of MEG (306) with EEG
(128) on PD patients.
I was successful in calculating the sources for both the modalities alone.
I tried according the tutorial "Combined EEG and MEG source reconstruction" and was successful till the estimation
of the lead field. But, then when i want to use the created vol for the
source analysis i have two structures one for EEG and the second one for
MEG. Now, i can only define only one to estimate the sources. How to use
this for both the modalities? Any help would be appreciated.
With regards<font color="#888888">
Mushfa</font></pre>
</div>
</div>
_______________________________________________<br>
fieldtrip mailing list<br>
<a moz-do-not-send="true"
href="mailto:fieldtrip@donders.ru.nl"
target="_blank">fieldtrip@donders.ru.nl</a><br>
<a moz-do-not-send="true"
href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip"
target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></blockquote>
</div>
<br>
</div>
</div>
</div>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a moz-do-not-send="true"
href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a moz-do-not-send="true"
href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip"
target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
</blockquote>
</div>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
fieldtrip mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Jörn M. Horschig
PhD Student
Donders Institute for Brain, Cognition and Behaviour
Centre for Cognitive Neuroimaging
Radboud University Nijmegen
Neuronal Oscillations Group
FieldTrip Development Team
P.O. Box 9101
NL-6500 HB Nijmegen
The Netherlands
Contact:
E-Mail: <a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:jm.horschig@donders.ru.nl">jm.horschig@donders.ru.nl</a>
Tel: +31-(0)24-36-68493
Web: <a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://www.ru.nl/donders">http://www.ru.nl/donders</a>
Visiting address:
Trigon, room 2.30
Kapittelweg 29
NL-6525 EN Nijmegen
The Netherlands</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
fieldtrip mailing list
<a class="moz-txt-link-abbreviated" href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a>
<a class="moz-txt-link-freetext" href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Jörn M. Horschig
PhD Student
Donders Institute for Brain, Cognition and Behaviour
Centre for Cognitive Neuroimaging
Radboud University Nijmegen
Neuronal Oscillations Group
FieldTrip Development Team
P.O. Box 9101
NL-6500 HB Nijmegen
The Netherlands
Contact:
E-Mail: <a class="moz-txt-link-abbreviated" href="mailto:jm.horschig@donders.ru.nl">jm.horschig@donders.ru.nl</a>
Tel: +31-(0)24-36-68493
Web: <a class="moz-txt-link-freetext" href="http://www.ru.nl/donders">http://www.ru.nl/donders</a>
Visiting address:
Trigon, room 2.30
Kapittelweg 29
NL-6525 EN Nijmegen
The Netherlands</pre>
</body>
</html>