#!/bin/sh
#
# s_ttest -- runs Ron's program to calculate a Student's t-test on two 3-D
#            data sets. Assumes you have run s_varimap, so have a 
#            filep.rsa and filep.var file for each 3D volume.
#		filenum1 and filenum2 are the number of individual maps 
#		averaged for each real space average.
#
#
#
# ******** NOTE ******** THE OUTPUT IS *.tstat WHICH IS A 3D ARRAY OF (t x 100)
# 		THIS FILE SHOULD BE INTERPRETATED USING A STUDENTS t TABLE AND 
#		CONTOURED AT WHATEVER SIGNIFICANCE LEVEL YOU CHOOSE (e.g. 99.5%)
# 		CONTOURING AT THE 95% LEVEL - THE LEVEL THAT MOST STATISTICIANS 
#		THINK IS SIGNIFICANT - IS TOO RELAXED A CRITERION!
#
#
#	usage: s_ttest filep1 filep2

group1=$1
group2=$2

if test -f ttest.zzz
then
	rm ttest.zzz
fi

filenum1=`header < $group1.rsa | head -15 | tail -1 | awk '{split($0,field,";")}{print field[2]}' | awk '{print $1}'`
filenum2=`header < $group2.rsa | head -15 | tail -1 | awk '{split($0,field,";")}{print field[2]}' | awk '{print $1}'`
echo $group1 "/" $group2 ";" $filenum1 "/" $filenum2 "data sets" >> ttest.zzz
echo $filenum1 >> ttest.zzz
echo $group1.rsa >> ttest.zzz
echo $group1.var >> ttest.zzz
echo $filenum2 >> ttest.zzz
echo $group2.rsa >> ttest.zzz
echo $group2.var >> ttest.zzz

ttest < ttest.zzz
mv TSTAT $group1"_"$group2".tstat"
echo "file created is  "$group1"_"$group2".tstat"
rm ttest.zzz