# /* vim: set filetype=sh : */ 
	
make -s prep
Dir=../../cocomo/data
Data=nasa93.csv
Data=coc81modeTypeLangType.csv
Here=`pwd`
Tmp=`mktemp -d` 
#trap "rm -rf $Tmp" 0 1 2 3 15

(printf "="  
cat ${Dir}/${Data} | cut -d, -f6- ) > $Tmp/$0.in 
head -1 ${Dir}/${Data} > $Tmp/head
cd $Tmp


main() { 
  echo "---| $2 |----------------------"
  coc1nums -s $0.in >  $0.tmp
  for((i=1;i<=30;i++))
  do
	 #blabln $i
     R=$RANDOM
     randomize Seed=$R $0.tmp
	#blabln "`wc -l test.arff`"
     lc Seed=$R Pred=30 Pass=1 train.arff Pass=2 test.arff \
	 | tee ${0}$i.out 

  done #|  cut -d, -f 3 | sort  -n | cat -n > dat
	exit
#  gnuplot <<EOF
#	    set ylabel "RE"
#        set xlabel "all REs, sorted" 
#		set terminal postscript eps
#	    set output "dat.eps"
#		plot "dat"
#EOF
#	epstopdf dat.eps
#	cp dat.pdf ${Here}/4.pdf
  exit
  (echo "#!,a,b,mmre,pred(30)"
  fgrep \! ${1}*.out | sort -t, -n +3 ) | malign
	echo ""
	(cut -d, -f 8-22  head
	gawk -F, ' function mean(sumX,n)     {return sumX/n}
             function sd(sumSq,sumX,n) {return sqrt((sumSq-((sumX*sumX)/n))/(n-1))}
	        /=/ {	q=",";
                key=$3 q $4 q $5 q $6 q $7 q $8 q $9 q \
                   $10 q $11 q $12 q $13 q $14 q $15 q $16 q $17 
			  	val=$NF;
              	s[ key] += val;
              	s2[key] += val * val;
              	n[key]++
			 }
             END { OFS=","; 
			       OFMT="%.0f";
                   for(i in s) {
					  mn    = mean(s[i],n[i]);
				      stdev = sd(s2[i],s[i],n[i])
                      print  i , "|",  mn, stdev,  100*stdev/mn
				  }
                 }
          ' ${1}*.out
    ) | malign -g 1
}
main d "-z acap=1 "
main b "-z vexp=1.21 -z tool=0.83"
main c "-z time=1.1 -z rely=1.2"
cd