#!/bin/csh -f

set subjectlist = (`cat slist.dat`);

set label = entorhinal_exvivo

foreach hemi (lh rh)

  foreach subject ($subjectlist)

    set sd = $SUBJECTS_DIR/$subject

    mri_segstats --slabel $subject $hemi \
      $sd/label/$hemi.$label.label \
      --i $sd/surf/$hemi.thickness --id 1 \
      --sum $sd/stats/$hemi.$label.thickness.stats 

    mri_segstats --slabel $subject $hemi \
      $sd/label/$hemi.$label.label \
      --i $sd/surf/$hemi.area  --accumulate --id 1 \
      --sum $sd/stats/$hemi.$label.area.stats 

    mri_segstats --slabel $subject $hemi \
      $sd/label/$hemi.$label.label \
      --i $sd/surf/$hemi.volume  --accumulate --id 1 \
      --sum $sd/stats/$hemi.$label.volume.stats 

  end # subject

end # hemi

foreach hemi (lh rh)
  foreach meas (thickness area volume)
    asegstats2table --subjects $subjectlist -m mean \
      --stats=$hemi.$label.$meas.stats \
      -t $hemi.$label.$meas.dat
  end # meas
end

exit 0