[Eeglablist] AMICA throws a bus error

[Michael Spezio]

Dear Jason, Makoto, and other AMICA users,

I am able to successfully run AMICA in Matlab with runamica12 (R2014b, iMac 10.9.5) on the sample data Memorize.fdt.

However, I get a bus error when I try to run AMICA on my own data. I've posted the error below.

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.

Can you provide any insight into what is going wrong? I'm following Makoto's pipeline.

Are there additional checks I need to do on my data prior to AMICA?

Thanks.

Best,

Michael

[ EEG.icaweights, EEG.icasphere, mods ] = runamica12(EEG.data(:,:));

Processing arguments ...
 num_files =            1
 FILES:
 /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/tmpdata
 92389.fdt
 num_dir_files =            1
 initial matrix block_size =          256
 do_opt_block =            0
 number of models =            1
 number of density mixture components =            3
 pdf type =            0
 max_iter =         2000
 num_samples =            1
 data_dim =          202
 field_dim =       106200
 do_history =            0
 histstep =           10
 share_comps =            0
 share_start =          100
 comp_thresh =   0.990000000000000
 share_int =          100
 initial lrate =   5.000000000000000E-002
 minimum lrate =   1.000000000000000E-008
 lrate factor =   0.500000000000000
 initial rholrate =   5.000000000000000E-002
 rho0 =    1.50000000000000
 min rho =    1.00000000000000
 max rho =    2.00000000000000
 rho lrate factor =   0.500000000000000
 kurt_start =            3
 num kurt =            5
 kurt interval =            1
 do_newton =            1
 newt_start =           50
 newt_ramp =           10
 initial newton lrate =    1.00000000000000
 do_reject =            0
 num reject =            3
 reject sigma =    3.00000000000000
 reject start =            2
 reject interval =            3
 max_thrds =            2
 write step =           10
 write_nd =            0
 write_LLt =            1
 dec window =            1
 max_decs =            3
 fix_init =            0
 update_A =            1
 update_c =            1
 update_gm =            1
 update_alpha =            1
 update_mu =            1
 update_beta =            1
 invsigmax =    100.000000000000
 invsigmin =   0.000000000000000E+000
 do_rho =            1
 load_rej =            0
 load_c =            0
 load_gm =            0
 load_alpha =            0
 load_mu =            0
 load_beta =            0
 load_rho =            0
 load_comp_list =            0
 do_mean =            1
 do_sphere =            1
 doPCA =            1
 pcakeep =          202
 pcadb =    30.0000000000000
 byte_size =            4
 doscaling =            1
 scalestep =            1
mkdir: /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/amicaouttmp/: File exists
 output directory =
 /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/amicaou
 ttmp/
           1 : setting num_thrds to            2  ...
           1 : using           2 threads.
           1 : node_thrds =            2
 i =            1  real ierr =            0
 bytes in real =            1
           1 : REAL nbyte =            1
 getting segment list ...
 blocks in sample =       106200
 total blocks =       106200
 node blocks =       106200
 node            1  start: file            1  sample            1  index
           1
 node            1  stop : file            1  sample            1  index
      106200
           1 : data =   -7.69530153274536       -17.7782669067383
 getting the mean ...
  mean =  -0.427972755604667      -0.451065303007086
 -0.407206153567977
 subtracting the mean ...
 getting the sphering matrix ...
 cnt =       106200
 doing eig nx =          202  lwork =       408040
 minimum eigenvalues =   9.293485032426282E-002  0.123639438673911
  0.145761886209716
 maximum eigenvalues =    76318.4901712733        22041.6054108954
   4928.60680527982
 num eigs kept =          202
 numeigs =          202
 sphering the data ...
           1 Allocating variables ...
           1 : Initializing variables ...
           1 : block size =          256
           1 : entering the main loop ...
forrtl: severe (180): SIGBUS, bus error occurred
/Users/mlspezio/Matlab/toolboxes/UseAMICA/amica12mac64 /Users/mlspezio/Matlab/scripts/Projects/HLCC_EEG/HLCC_EEG_Preprocessing/amicaouttmp/input.param: Signal 52
No gm present, setting num_models to 1
No W present, exiting
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20150617/538a9d57/attachment.html>