# /* vim: set filetype=sh : */ 
	
Safe=$HOME/tmp/locomo/6
mkdir -p $Safe

Dir=../../cocomo/data
Data=nasa93.csv

Tmp=`mktemp -d` 
#trap "rm -rf $Tmp" 0 1 2 3 15

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

main() { 
  echo "---|  $1 |----------------------"
  	printf "$2 \n\n"
for n in 15 30 45 60 100 
do
		echo "::: $n"
	coc1nums -s $0.in \
  | near -n $n $2  > $0.tmp


  for((i=1;i<=20;i++))
  do
     R=$RANDOM
     randomize TestSets=5 Seed=$R $0.tmp
     lc ShowRE=re_$1_$n Seed=$R Pred=30 Pass=1 train.arff Pass=2 test.arff > ${1}$i.out 
  done 
  cat re_$1_$n | sort -n > sorted
  mv sorted re_$1_$n
	gawk '{sum += $1 }END{ print sum}' re_$1_$n > area_$1_$n
  ls -lsa re_*
  fgrep \! ${1}*.out| sort -t, -n +3| cut -d, -f 4| cat -n|sed 's/,/    /' > 4_$n
  fgrep \! ${1}*.out| sort -t, -n +4| cut -d, -f 5| cat -n|sed 's/,/    /' > 5_$n
  continue	
  (echo "#!,a,b,mmre,pred(30)"
  fgrep \! ${1}*.out | sort -t, -n +4 ) | 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
done
gnuplot<<EOF
set xlabel "20 experiments"
set size 0.5,0.7
set terminal postscript eps enhanced color 20
set noxtics
set key top left
set output "$Safe/mmre_$1.eps"
set  ylabel "MMRE" 
set  yrange [0:200]
set  ytics (0,50,100,150,200)
plot "4_15" using 1:2 title "15" with linesp,\
     "4_30" using 1:2 title "30" with linesp,\
     "4_45" using 1:2 title "45" with linesp,\
     "4_60" using 1:2 title "60" with linesp,\
     "4_100" using 1:2 title "100" with linesp
set output "$Safe/pred_$1.eps"
set  ylabel "PRED(30)" 
set  yrange [0:120]
set  ytics (0,25,50,75,100)
plot "5_15" using 1:2 title "15" with linesp,\
     "5_30" using 1:2 title "30" with linesp,\
	 "5_45" using 1:2 title "45" with linesp,\
	 "5_60" using 1:2 title "60" with linesp,\
	 "5_100" using 1:2 title "100" with linesp
EOF
(cd $Safe
epstopdf mmre_$1.eps
epstopdf pred_$1.eps
xpdf mmre_$1.pdf &
xpdf pred_$1.pdf &
)
}

plot() {
cat area_*
for i in area_*
do
	echo $i `cat $i`
done

#area_flight_100 37.9745
#area_flight_15 38.1902
#area_flight_30 20.3624
#area_flight_45 26.0544
#area_flight_60 17.5929
#area_ground_100 12.3093
#area_ground_15 9.75076
#area_ground_30 8.68157
#area_ground_45 2.4774
#area_ground_60 23.0339

gnuplot<<EOF
set key top left
set terminal postscript eps enhanced color 20
set xlabel "sorted values"
set ylabel "(predicted - actual)/actual"
set size 0.5,0.7
set noxtics
set xrange [0:90]
set ytics (-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
set output "$Safe/re_flight.eps"
plot "$Tmp/re_flight_15"  title "15" with linesp,\
	"$Tmp/re_flight_30"  title "30" with linesp,\
	"$Tmp/re_flight_45"  title "45" with linesp,\
	"$Tmp/re_flight_60" title "60" with linesp,\
	"$Tmp/re_flight_100" title "100" with linesp
EOF
gnuplot <<EOF
set size 0.5,0.7
set noxtics
set key top left
set xrange [0:90]
set xlabel "sorted values"
set ylabel "(predicted - actual)/actual"
set ytics (-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
set terminal postscript eps enhanced color 20
set output "$Safe/re_ground.eps"
plot "$Tmp/re_ground_15"  title "15" with linesp,\
	"$Tmp/re_ground_30"  title "30" with linesp,\
	"$Tmp/re_ground_45"  title "45" with linesp,\
	"$Tmp/re_ground_60" title "60" with linesp,\
	"$Tmp/re_ground_100" title "100" with linesp
EOF
gnuplot<<EOF
set key top left
set terminal postscript eps enhanced color 20
set xlabel "sorted values"
set ylabel "(predicted - actual)/actual"
se yrange [-1:1]
set size 0.5,0.7
set noxtics
set xrange [0:90]
set ytics (-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
set output "$Safe/re_flight_small.eps"
plot "$Tmp/re_flight_15"  title "15" with linesp,\
	"$Tmp/re_flight_30"  title "30" with linesp,\
	"$Tmp/re_flight_45"  title "45" with linesp,\
	"$Tmp/re_flight_60" title "60" with linesp,\
	"$Tmp/re_flight_100" title "100" with linesp
EOF
gnuplot <<EOF
set size 0.5,0.7
set noxtics
set key top left
set xrange [0:90]
set xlabel "sorted values"
se yrange [-1:1]
set ylabel "(predicted - actual)/actual"
set ytics (-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
set terminal postscript eps enhanced color 20
set output "$Safe/re_ground_small.eps"
plot "$Tmp/re_ground_15"  title "15" with linesp,\
	"$Tmp/re_ground_30"  title "30" with linesp,\
	"$Tmp/re_ground_45"  title "45" with linesp,\
	"$Tmp/re_ground_60" title "60" with linesp,\
	"$Tmp/re_ground_100" title "100" with linesp
EOF
(cd $Safe
epstopdf re_flight.eps 
epstopdf re_ground.eps 
xpdf re_flight.pdf &
xpdf re_ground.pdf &

epstopdf re_flight_small.eps 
epstopdf re_ground_small.eps 
xpdf re_flight_small.pdf &
xpdf re_ground_small.pdf &
)

}
main flight " 
           -z rely=1.275 -z data=0.97  -z cplx=1.4  -z time=1.055 
           -z stor=1.03  -z virt=0.935 -z acap=0.85 -z aexp=0.91 
	   	   -z pcap=0.85 -z vexp=0.95  -z lexp=0.975
		   -z modp=1.05  -z tool=1.05 -z sced=1.04 " 

main ground "   
           -z rely=0.87  -z data=1.08  -z cplx=0.92   -z time=1 
           -z stor=1     -z virt=0.93  -z acap=0.93   -z aexp=0.955 
	   	   -z pcap=0.93  -z vexp=0.95  -z lexp=0.975
		   -z modp=1.05  -z tool=1     -z sced=1 " 

plot

cd

#area_flight_100 11.778
#area_flight_30 11.4225
#area_flight_45 8.65107
#area_flight_15 7.91081
#area_flight_60 5.98316
#
#area_ground_100 15.6836
#area_ground_60 18.7685
#area_ground_15 14.0665
#area_ground_30 1.84147
#area_ground_45 0.291413
#