Hello,
almost, but instead of trying to compute the variance from the con-image
(which would be over all voxels, i.e. integrating in space) try using
the ResMS.img (which already contains the variance estimates per voxel).
And by the way, using
con = spm_read_vols(spm_vol('con_0001.img'));
may be safer if you want to stick within the spm environment.
Cheers,
Marko
Fred Sampedro wrote:
> Dear SPM experts,
>
> I am trying to understand the process of obtaining the spmT_0001.img
> volume. After searching on the web/documentation, I expected that this
> would solve it:
>
> con = read_spm_volume ('con_0001.img') %this a function of mine that
> loads the volume directly
> my_t = con ./ sqrt(var(con(:)))
>
> but, when comparing "my_t" with the actual generated t volume (i.e. t =
> read_spm_volume ('spmT_0001.img') they are actually different. What am I
> missing here?
>
> Thanks a lot!
--
____________________________________________________
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/
____________________________________________________
|