Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-22 Thread Qi Zhu
Ok, that's great. Thanks very much for your help. Best wishes and Merry Christmas! Qi Zhu, PostDoc Laboratorium voor Neuro- en Psychofysiologie K.U.Leuven Medical School Herestraat 49, B-3000 Leuven (Belgium) On Wed, Dec 21, 2011 at 7:40 PM, Douglas N Greve wrote: > Yes that is correct. That wa

Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-21 Thread Douglas N Greve
Yes that is correct. That was fixed with version 5.1 so current distributions do the right thing. doug Qi Zhu wrote: > Yes, seems weird. The image that I attached to you in the first mail > is actually a t map, so the t value in those areas are also zero. > And according to what I read from fast_

Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-21 Thread Qi Zhu
Yes, seems weird. The image that I attached to you in the first mail is actually a t map, so the t value in those areas are also zero. And according to what I read from fast_selxavg3.m, the t value is actually calculated from the fsig, which is: t.vol = sqrt(fsig.vol).*sign(ces.vol), and fsig.vol i

Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-21 Thread Douglas N Greve
The t is computed as the signed square root of the F (not Fsig) volume: t.vol = sqrt(Fvol.vol) .* sign(ces.vol); The Fvol (and do t) is not computed from the sig (the other way around). Are you sure that the t-value does not look right in those areas? doug Qi Zhu wrote: > Dear Doug, > > Thank

Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-21 Thread Qi Zhu
Dear Doug, Thanks very much for your prompt response. I tried what you suggested and the sig map looks fine now, that's probably the best that we can get for the sig map, because of the limitation of the MATLAB I guess. But for the t map, the values are still not correct for those voxels that high

Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-20 Thread Douglas N Greve
btw, if you simply "touch" the contrast mat files in the analysis directory, it will only regenerate the contrasts. That should save you some time because it won't have to re-run all of selxavg. doug Douglas N Greve wrote: > I think this is due to the high DOF and probably high t-value causing

Re: [Freesurfer] Selxavg3 sets high activation to zero

2011-12-20 Thread Douglas N Greve
I think this is due to the high DOF and probably high t-value causing the significance in matlab to be computed as 0. This gets converted to Inf (infinity) when the -log10(p) calculation is made, then the Inf gets set to 0 when the volume is saved. The quick fix for you is to change fast_selxav