avar=$1
shift
mris_preproc --fsgd ${avar}.fsgd --cache-in thickness.fwhm10.fsaverage --target fsaverage --hemi lh --out lh.${avar}.thickness.10.mgh
mris_preproc --fsgd ${avar}.fsgd --cache-in thickness.fwhm10.fsaverage --target fsaverage --hemi rh --out rh.${avar}.thickness.10.mgh
mri_glmfit --y lh.${avar}.thickness.10.mgh --fsgd ${avar}.fsgd --C comps.mtx --surf fsaverage lh --cortex --glmdir lh.${avar}.glmdir
mri_glmfit --y rh.${avar}.thickness.10.mgh --fsgd ${avar}.fsgd --C comps.mtx --surf fsaverage rh --cortex --glmdir rh.${avar}.glmdir
mri_glmfit-sim --glmdir lh.${avar}.glmdir --cache 4 neg --cwp 0.05 --2spaces
mri_glmfit-sim --glmdir rh.${avar}.glmdir --cache 4 neg --cwp 0.05 --2spaces
mri_glmfit-sim --glmdir lh.${avar}.glmdir --cache 4 pos --cwp 0.05 --2spaces
mri_glmfit-sim --glmdir rh.${avar}.glmdir --cache 4 pos --cwp 0.05 --2spaces
echo "lh, neg"
grep -v "#" lh.${avar}.glmdir/comps/cache.th40.neg.sig.cluster.summary
echo "rh, neg"
grep -v "#" rh.${avar}.glmdir/comps/cache.th40.neg.sig.cluster.summary
echo "lh, pos"
grep -v "#" lh.${avar}.glmdir/comps/cache.th40.pos.sig.cluster.summary
echo "rh, pos"
grep -v "#" rh.${avar}.glmdir/comps/cache.th40.pos.sig.cluster.summary