>> global MAXMEM; MAXMEM=256;

surf = SurfStatReadSurf( 'av.obj' );

mask = SurfStatMaskCut( surf );

maskb = (mask & SurfStatROI([0; -16; -8], 20, surf ) == 0);

[subject, RatioLeft, RatioRight, gender, age, scanner]...
=textread('GLIMV123ratio2ALLwomissingwoqcmarker.csv', '%s %s %s %s %f %s');

clf; hist(age); xlabel('age')

R1 = SurfStatReadData( [RatioLeft] );
R2 = SurfStatReadData( [RatioRight] );
R3 = horzcat(R1, R2);
752 x 1 files to read, % remaining: 100 90 80 70 60 50 40 30 20 10 0 Done
752 x 1 files to read, % remaining: 100 90 80 70 60 50 40 30 20 10 0 Done
>> Age = term( age );
Gender = term ( gender );
Scanner = term (scanner );
Subject = term ( subject );
>>
>> slm = SurfStatT( SurfStatLinMod( R3, 1 + Age + Gender + random( Subject ) + I, surf ), Age );
[ pval, peak, clus, clusid ] = SurfStatP( slm, mask );
SurfStatView( pval, surf, 'PVAL LINEAR'  );

>> slm = SurfStatT( SurfStatLinMod( R3, 1 + Age + Gender + random( Subject ) + I, surf ), Age );
[ pval, peak, clus, clusid ] = SurfStatP( slm, mask );
SurfStatView( pval, surf, 'PVAL LINEAR'  );

id =

       25439

>>
clusid( 25439 )

ans =

     5

>> maskROI = clusid == clusid( 25439 );
>> Linearmodelcluster5 = mean( R3(:, maskROI), 2 );

>> SurfStatPlot( Age, Linearmodelcluster5, Gender );