<div dir="ltr"><div><div><div><div><div>Hi all,<br><br></div>Just to update, AMICA worked on my MEG data after set the 'numprocs' to 1 or 0. I am not sure how that happened as my computer has 4 physical processors. <br><br></div>Hi Jason,<br></div>Could explain how the argument 'numprocs' influences AMICA processing?<br><br></div>Thanks,<br></div>Yanan<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 29 September 2015 at 16:59, 孙亚楠 <span dir="ltr"><<a href="mailto:sunyanan2008@gmail.com" target="_blank">sunyanan2008@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi eeglablist members,<br><br>I get the same bus error like what Michael got before. Has anyone solved this?<br></div><br>I can run amica successfully on the sample data (Memorized.fdt) and on my EEG data (64 channels) but not on my MEG data (160 channels Yokogawa) where showed "forrtl: severe (180): SIGBUS, bus error occurred". <br></div><div><br>The MEG data has standardized and it can be run successfully using the AMICA app but bus error happened when calling runamica12.m in MATLAB.<br><br>P.S. I am running runamica12 on R2014b, iMac 10.10.5.<br><br></div><div>Thanks,<br></div><div>Yanan<br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On 17 June 2015 at 14:20, Michael Spezio <span dir="ltr"><<a href="mailto:mspezio@scrippscollege.edu" target="_blank">mspezio@scrippscollege.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">




<div dir="ltr">
<div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:'Times New Roman',Times,serif">
Dear Jason, Makoto, and other AMICA users,<br>
<br>
I am able to successfully run AMICA in Matlab with runamica12 (R2014b, iMac 10.9.5) on the sample data Memorize.fdt.<br>
<br>
However, I get a bus error when I try to run AMICA on my own data. I've posted the error below.<br>
<br>
I've verified that the data rank is equal to the number of channels (202). I've also verified that the cnt is correct, 118 epochs, 900 points per epoch. The epochs are 4 nonoverlapping types, each having the -1 to 2 second default bounds.<br>
<br>
Can you provide any insight into what is going wrong? I'm following Makoto's pipeline.<br>
<br>
Are there additional checks I need to do on my data prior to AMICA?<br>
<br>
Thanks.<br>
<br>
Best,<br>
<br>
Michael<br>
<br>
[ EEG.icaweights, EEG.icasphere, mods ] = runamica12(EEG.data(:,:));<br>
<br>
Processing arguments ...<br>
 num_files =            1<br>
 FILES: <br>
 /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/tmpdata<br>
 92389.fdt<br>
 num_dir_files =            1<br>
 initial matrix block_size =          256<br>
 do_opt_block =            0<br>
 number of models =            1<br>
 number of density mixture components =            3<br>
 pdf type =            0<br>
 max_iter =         2000<br>
 num_samples =            1<br>
 data_dim =          202<br>
 field_dim =       106200<br>
 do_history =            0<br>
 histstep =           10<br>
 share_comps =            0<br>
 share_start =          100<br>
 comp_thresh =   0.990000000000000     <br>
 share_int =          100<br>
 initial lrate =   5.000000000000000E-002<br>
 minimum lrate =   1.000000000000000E-008<br>
 lrate factor =   0.500000000000000     <br>
 initial rholrate =   5.000000000000000E-002<br>
 rho0 =    1.50000000000000     <br>
 min rho =    1.00000000000000     <br>
 max rho =    2.00000000000000     <br>
 rho lrate factor =   0.500000000000000     <br>
 kurt_start =            3<br>
 num kurt =            5<br>
 kurt interval =            1<br>
 do_newton =            1<br>
 newt_start =           50<br>
 newt_ramp =           10<br>
 initial newton lrate =    1.00000000000000     <br>
 do_reject =            0<br>
 num reject =            3<br>
 reject sigma =    3.00000000000000     <br>
 reject start =            2<br>
 reject interval =            3<br>
 max_thrds =            2<br>
 write step =           10<br>
 write_nd =            0<br>
 write_LLt =            1<br>
 dec window =            1<br>
 max_decs =            3<br>
 fix_init =            0<br>
 update_A =            1<br>
 update_c =            1<br>
 update_gm =            1<br>
 update_alpha =            1<br>
 update_mu =            1<br>
 update_beta =            1<br>
 invsigmax =    100.000000000000     <br>
 invsigmin =   0.000000000000000E+000<br>
 do_rho =            1<br>
 load_rej =            0<br>
 load_c =            0<br>
 load_gm =            0<br>
 load_alpha =            0<br>
 load_mu =            0<br>
 load_beta =            0<br>
 load_rho =            0<br>
 load_comp_list =            0<br>
 do_mean =            1<br>
 do_sphere =            1<br>
 doPCA =            1<br>
 pcakeep =          202<br>
 pcadb =    30.0000000000000     <br>
 byte_size =            4<br>
 doscaling =            1<br>
 scalestep =            1<br>
mkdir: /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/amicaouttmp/: File exists<br>
 output directory = <br>
 /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/amicaou<br>
 ttmp/<br>
           1 : setting num_thrds to            2  ...<br>
           1 : using           2 threads.<br>
           1 : node_thrds =            2<br>
 i =            1  real ierr =            0<br>
 bytes in real =            1<br>
           1 : REAL nbyte =            1<br>
 getting segment list ...<br>
 blocks in sample =       106200<br>
 total blocks =       106200<br>
 node blocks =       106200<br>
 node            1  start: file            1  sample            1  index <br>
           1<br>
 node            1  stop : file            1  sample            1  index <br>
      106200<br>
           1 : data =   -7.69530153274536       -17.7782669067383     <br>
 getting the mean ...<br>
  mean =  -0.427972755604667      -0.451065303007086     <br>
 -0.407206153567977     <br>
 subtracting the mean ...<br>
 getting the sphering matrix ...<br>
 cnt =       106200<br>
 doing eig nx =          202  lwork =       408040<br>
 minimum eigenvalues =   9.293485032426282E-002  0.123639438673911     <br>
  0.145761886209716     <br>
 maximum eigenvalues =    76318.4901712733        22041.<a href="tel:6054108954" value="+16054108954" target="_blank">6054108954</a>     <br>
   4928.60680527982     <br>
 num eigs kept =          202<br>
 numeigs =          202<br>
 sphering the data ...<br>
           1 Allocating variables ...<br>
           1 : Initializing variables ...<br>
           1 : block size =          256<br>
           1 : entering the main loop ...<br>
forrtl: severe (180): SIGBUS, bus error occurred<br>
/Users/mlspezio/Matlab/toolboxes/UseAMICA/amica12mac64 /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/amicaouttmp/input.param: Signal 52<br>
No gm present, setting num_models to 1<br>
No W present, exiting<br>
<div></div>
</div>
</div>

<br></div></div>_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank">eeglablist-request@sccn.ucsd.edu</a><span class="HOEnZb"><font color="#888888"><br></font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><br>-- <br><div><div>Yanan Sun</div><div>PHD Student<br><div><div><div style="font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">ARC Centre of Excellence in Cognition and its Disorders<br></div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">& Department of Cognitive Science, </span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Faculty of Human Sciences,</span></div></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Macquarie University, </span></div><div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">NSW, 2109, </span><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">Australia</span></div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div>Yanan Sun</div><div>PHD Student<br><div><div><div style="font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">ARC Centre of Excellence in Cognition and its Disorders<br></div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">& Department of Cognitive Science, </span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Faculty of Human Sciences,</span></div></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Macquarie University, </span></div><div><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">NSW, 2109, </span><span style="background-color:rgb(255,255,255);color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px">Australia</span></div></div></div>
</div>