<div><div>Hi Michael and Yahan,</div><div>I think the problem is that the stack is overflowing with the default block size and numprocs since you have a large number of channels.</div><div>On laptops, numprocs should set the number of threads, default is 2. This influences the block size as well.</div><div>A possible solution is to use the keyword argument: 'block_size', 32, ...</div><div>or a smaller number, and possibly 'numprocs', 1, ... </div><div>Let me know if this doesn't work.</div><div>Best,</div><div>Jason</div></div><div class="elided-text">On Sep 30, 2015 10:29 AM, 孙亚楠 <sunyanan2008@gmail.com> wrote:<br type='attribution'><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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><br /><div class="elided-text">On 29 September 2015 at 16:59, 孙亚楠 <span dir="ltr"><<a href="mailto:sunyanan2008@gmail.com">sunyanan2008@gmail.com</a>></span> wrote:<br /><blockquote style="margin:0 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><br /><div class="elided-text"><div><div>On 17 June 2015 at 14:20, Michael Spezio <span dir="ltr"><<a href="mailto:mspezio@scrippscollege.edu">mspezio@scrippscollege.edu</a>></span> wrote:<br /></div></div><blockquote style="margin:0 0 0 0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
<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.6054108954 <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">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br />
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">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">eeglablist-request@sccn.ucsd.edu</a><font color="#888888"><br /></font></blockquote></div><font color="#888888"><br /><br clear="all" /><br />-- <br /></font><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>
</div>
</blockquote></div><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>
</div>
</blockquote></div>