<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" 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 id="Signature"></div>
</div>
</body>
</html>