Hi Nicholas,
thanks for your mail. Nice tracking of the origin of that error message.
Without having double checked your description, I would tend to agree
that this looks like a bug, although I would never dare to suggest that
Tom was wrong in anything ;) Perhaps he would like to have a look at that?
Cheers,
Marko
Nicholas Dufour wrote:
> Hi Marko--
>
> Sorry about the year long delay--I ended up just using SnPM5 for this particular issue. However, it would be good to resolve the problem now since I'm back on a similar project.
>
> The problem does not have to do with the machine's memory; the problem appears to arise precisely because I am defining the threshold ahead of time. By stepping through snpm_cp from snpm8 and snpm5, here's what I've been able to figure out. In snpm5, there is the clause:
>
> if bST & (pU_ST_Ut~= -1), [line 898]
>
> which causes it to execute:
>
> clear d1 d2
> d1 = find(T(1,:) > ST_Ut);
> d2 = find(T(1,:) < -ST_Ut);
> spm_append_96('SnPM_ST',[XYZ(:,d1), XYZ(:,d2);T(1,d1),-T(1,d2);perm*ones(1,length(d1)), -perm*ones(1,length(d2))]); [line 902]
> end
>
> Here at line 902, the SnPM_ST file is created (in SnPM5). However, in SnPM8, the if clause is changed, to
>
> if bST [line 902]
>
> if pU_ST_Ut==-1 [line 904]
>
> Because my threshold is defined (i.e., ~= -1), the "spm_append_96" function is never run--and, as far as I can tell, SnPM_ST is never created.
>
> Any thoughts? If we allow snpm_cp to run to completion on SnPM5, SnPM_ST *is* created, and we can run the combined voxel/cluster analysis just fine.
>
> Thanks!
>
--
____________________________________________________
PD Dr. med. Marko Wilke
Facharzt für Kinder- und Jugendmedizin
Leiter, Experimentelle Pädiatrische Neurobildgebung
Universitäts-Kinderklinik
Abt. III (Neuropädiatrie)
Marko Wilke, MD, PhD
Pediatrician
Head, Experimental Pediatric Neuroimaging
University Children's Hospital
Dept. III (Pediatric Neurology)
Hoppe-Seyler-Str. 1
D - 72076 Tübingen, Germany
Tel. +49 7071 29-83416
Fax +49 7071 29-5473
[log in to unmask]
http://www.medizin.uni-tuebingen.de/kinder/epn/
____________________________________________________
|