Right, the ocn.dat files have data that is uncorrected in that sense and might 
need to nuisance factors removed before plotting. There is a design matrix in 
there (Xg.dat). You can load that into matlab along with the ocn.dat, compute 
beta = inv(X'*X)*(X'*ocn) to get the betas. You can then compute yhat = 
X2*beta2 where X2 has nuisance columns removed and beta2 has the same nuisance 
coefficients removed, then treat yhat as your data to be plotted.

 Note that plotting the results is still a form a voodoo correlations because 
your eye will compute the correlation even if you don't explicitly do so 
(though it generally does not stop anyone:).

On 9/8/2019 7:37 PM, cody samth wrote:

        External Email - Use Caution

Hi Douglas, thanks for your response.

>They should not, but the reason it fairly convoluted. When you get a
>cluster after running mri_glmfit-sim, that cluster is on fsaverage which
>is an average of 40 subjects. The area of a vertex is computed as the
>average of the areas of the vertices from the 40 that mapped into that
>vertex. This is the number that is used to compute the surface area of
>the cluster in the summary file. Now, when you map your subjects into
>the fsaverage space, they may have more or less surface area mapping
>into that cluster relative to the 40 (looks like more from #2 below).
>Also, you probably smoothed the surface area, which could have an
>unpredictable effect.
Thanks that makes sense.

>> 3) ocn.dat files are the input values meaning they're raw and would
>> need to be corrected in a statistically (in a similar way that I
>> modeled it in freesurfer) before graphing right?
>Not sure what you mean by "corrected" here. In general, you need to be
>very careful when you extract data from a cluster. It would be circular
>to do the same test that you used to generate the cluster, though this
>happens a lot (see "VooDoo correlations" by Ed Vul).
My apologies corrected wasn't the best way to phrase that question.
My interpretation of the ocn.dat file is that the each row contains the
average input value for a subject prior to controlling for covariates.
Therefore, to graph these results wouldn't these values need to undergo
some method to control for covariates such as ICV, sex or age to better
reflect the clusters observed from the GLM?

Or are the values in the ocn.dat file already reflective of the test/glm
used to generate the cluster?

Freesurfer mailing list

Freesurfer mailing list

Reply via email to