<div dir="ltr"><br clear="all"><div>Hi all</div><div><br></div><div><div>I followed the tips above to simulate a dipole source on my own head model generated from MRIs. However, using information from: </div><div><br></div><div><a href="http://fieldtrip.fcdonders.nl/example/compute_forward_simulated_data_and_apply_a_dipole_fit">http://fieldtrip.fcdonders.nl/example/compute_forward_simulated_data_and_apply_a_dipole_fit</a></div><div><br></div><div>I could not simulate anything. Only one image is blank. The code I am using is as follows: </div><div><br></div><div>cfg = []; </div><div>cfg.vol = vol; % See above </div><div>cfg.elec = elec; % See above </div><div>cfg.dip.pos = [0 0 0]; </div><div>cfg.dip.mom = [2 2 2]; % Note, it Should be transposed </div><div>cfg.dip.frequency = 10; </div><div>cfg.ntrials = 10; </div><div>cfg.triallength = 1; % seconds </div><div>cfg.fsample = 25; % Hz </div><div>raw1 = ft_dipolesimulation (cfg); </div><div>ft_timelockanalysis avg1 = ([], raw1); </div><div>plot (avg1.time, avg1.avg); % Plot the timecourse </div><div><br></div><div>where: </div><div><br></div><div>vol = ft_prepare_headmodel volumes generated from the segmentation of my resonances. </div><div><br></div><div>and </div><div><br></div><div>ft_read_sens elec = ('standard_1020.elc'); </div><div><br></div><div>to run the script the result is a blank image. </div><div><br></div><div>and Matlab tells me the following: </div><div><br></div><div>using headmodel specified in the configuration </div><div>using electrodes specified in the configuration </div><div>Determining source compartment (1) </div><div>projecting electrodes on skin surface </div><div>electrode and system combine combining transfer matrix </div><div>computing simulated data </div><div>computing simulated trial data for 10 </div><div><br></div><div>RMS value of simulated data is NaN </div><div>the call to "ft_dipolesimulation" took 4 seconds </div><div>the input is raw Data with 97 channels and 10 trials </div><div>Warning: the data does not Contain a definition trial </div><div>> In utilities \ private \ warning_once at 158 </div><div>   In utilities \ private \ fixsampleinfo at 66 </div><div>   In ft_datatype_raw at 157 </div><div>   In ft_checkdata at 224 </div><div>   In ft_timelockanalysis at 105 </div><div>   In codgen at 50 </div><div>Warning: sampleinfo by reconstructing Assuming That the trials are consecutive segments of </div><div>a continuous recording </div><div>> In utilities \ private \ warning_once at 158 </div><div>   In utilities \ private \ fixsampleinfo at 79 </div><div>   In ft_datatype_raw at 157 </div><div>   In ft_checkdata at 224 </div><div>   In ft_timelockanalysis at 105 </div><div>   In codgen at 50 </div><div>averaging trials </div><div>trial averaging 10 of 10 </div><div><br></div><div>the call to "ft_timelockanalysis" took 1 seconds </div><div><br></div><div><br></div><div>I hope you can help </div><div><br></div><div><br></div><div>greetings</div></div>-- <br><div dir="ltr"><div style="text-align:center"><font color="#20124d"><b>Gamaliel Huerta</b></font></div><div style="text-align:center"><font color="#20124d"><b>Ingeniería Civil Biomédica</b></font></div><div style="text-align:center"><font color="#20124d"><b>Universidad de Valparaíso</b></font></div></div>
</div>