I would do it in matlab, eg,
F = MRIread('F.mgh'); % This is to get the proper header information
lfdr = load('lfdr.asc'); F.vol = fast_mat2vol(lfdr,F.volsize); MRIwrite(F,'lfdr.mgh');
On 6/15/16 4:45 AM, Robby De Pauw wrote:
Dear FreeSurfer users,
To address the problem of multiple comparison, I wanted to apply the localfdr-method. For this purpose, I need the F-statistic (which I can extract from the F.mgh scalar overlay I presume?).
I converted this to an ascii-file, which I was able to load into R, using the following chink of code:
mris_convert -c /Users/robbydepauw/Desktop/FreeSurfer/GLM/lh.Group.volume.glmdir/Groupeffect/F.mgh $SUBJECTS_DIR/fsaverage/surf/lh.white F.asc
After performing an local fdr analysis in R, I want to create a scalar overlay “lfdr.mgh”, but I fail at converting the asci file into an “.mgh”-overlay file.
Any idea how I could approach this conversion (from asci back to .mgh)?
Many thanks,
Greetings
Robby
Robby De Pauw, PT, MT, drs.
*T *+32 09 332 12 19* | Email *Robby.DePauw@Ugent.be mailto:Robby.DePauw@ugent.be
Ghent University *| *Rehabilitation Sciences and Physiotherapy
Campus Heymans 3B3 (022)
De Pintelaan 185 BE-9000 Ghent
Route to my office https://soleway.ugent.be/routes/1090
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer