Hi - If I would like to try to perform mass-univariate statistics at the vertex level (thickness, area, vol at each vertex) without using the mri_glmfit, is this order a sound step (in general) for mapping stats parameter of choice (p-val or t s or weights for example) to the surface?
1. Use mris_preproc and surf2surf to get the data to common space and smooth.
2. Load the data to matlab using fs_read_Y() or MRIread
3. Perform stats at each vertex using stats toolbox to apply modeling approach of your choice (e.g., machine learning, linear mixed model etc.)
4. Extract P-values (or other parameters) for each vertex for the variable of primary interest
5. Use MRIwrite to create mgh file
6. Perform multiple comparison using freesurfer functions (e.g., FDR or monte carlo sim)
7. view corrected maps using tksufer?
Any steps I should consider adding or avoid?
Thank you!