In [8]: ae.fit()
-> Fitting "NODDI" model to 290353 voxels:
[ ] 0.0%---------------------------------------------------------------------------
RuntimeError Traceback (most recent call last)
<ipython-input-8-fad37d61652c> in <module>()
----> 1 ae.fit()
/Users/takahashi/anaconda/envs/amico/lib/python2.7/site-packages/amico/core.pyc in fit(self)
298
299 # dispatch to the right handler for each model
--> 300 MAPs[ix,iy,iz,:], DIRs[ix,iy,iz,:], x, A = self.model.fit( y, dirs.reshape(-1,3), self.KERNELS, self.get_config('solver_params') )
301
302 # compute fitting error
/Users/takahashi/anaconda/envs/amico/lib/python2.7/site-packages/amico/models.pyc in fit(self, y, dirs, KERNELS, params)
557 An = A[ self.scheme.dwi_idx, :-1 ] * KERNELS['norms']
558 yy = yy[ self.scheme.dwi_idx ].reshape(-1,1)
--> 559 x = spams.lasso( np.asfortranarray(yy), D=np.asfortranarray(An), **params ).todense().A1
560
561 # debias coefficients
/Users/takahashi/anaconda/envs/amico/lib/python2.7/site-packages/spams.pyc in lasso(X, D, Q, q, return_reg_path, L, lambda1, lambda2, mode, pos, ols, numThreads, max_length_path, verbose, cholesky)
437 ((indptr,indices,data,shape),path) = spams_wrap.lassoD(X,D,return_reg_path,L,lambda1,lambda2,mode,pos,ols,numThreads,max_length_path,verbose,cholesky)
438 else:
--> 439 (indptr,indices,data,shape) = spams_wrap.lassoD(X,D,return_reg_path,L,lambda1,lambda2,mode,pos,ols,numThreads,max_length_path,verbose,cholesky)
440 alpha = ssp.csc_matrix((data,indices,indptr),shape)
441 if return_reg_path:
/Users/takahashi/anaconda/envs/amico/lib/python2.7/site-packages/spams_wrap.pyc in lassoD(*args)
199 bool cholevsky) -> SpMatrix<(float)>
200 """
--> 201 return _spams_wrap.lassoD(*args)
202
203 def lassoQq(*args):
RuntimeError: matrix arg 1 must be a 2d double Fortran Array