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.
- 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@nmr.mgh.harvard.edumailto:Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer