28 Feb
2011
28 Feb
'11
10:04 a.m.
When the residual variance is exactly 0, this is what the matlab code. In the C implementation, if (glm->rvar < FLT_MIN) glm->rvar = FLT_MIN;
which will have similar results.
doug
Lars M. Rimol wrote:
Hi,
We have experimented a little with different ways of computing the p-values in a GLM analysis based on FreeSurfer data.
The way the Matlab reference implementation (fast_fratio.m) handles a situtation where the residuals are exactly zero, seems to be to set those vertices to F = 0 and p = 1. Is this correctly understood? Is this also the case in the underlying C-code?
Thank you!
LMR
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: www.nmr.mgh.harvard.edu/facility/filedrop/index.html