#! /bin/bash

Data=cocomost_3.log
Goals="coc81_all coc81_kind_min coc81_lang_ftn coc81_mode_e coc81_kind_max coc81_mode_org coc81_lang_mol nasa93_all nasa93_project_gro nasa93_cat2_missionplanning nasa93_cat2_avionicsmonitoring nasa93_mode_semidetached nasa93_project_sts nasa93_forg_g nasa93_center_5 nasa93_year_1975 nasa93_year_1980 nasa93_mode_embedded nasa93_center_2"
Methods="LC cocotest"

# Filename, I, Seed, Learner, Subset, Estimate, Actual, MRE

generate() {
(
echo | stats Report=0  Prefix="#goal,what"
for method in $Methods
do
	for goal in $Goals
	do
		grep $goal $Data |
		grep $method | 
		cut -d, -f 8  |
	    stats Loud=0 Prefix="$goal,$method"
	done
done
) | sort  -t, -n +7 > $0.tmp  
}

report() {

for method in $Methods
do
	grep $method  $0.tmp | cut -d, -f 6,8 | sed 's/,/	/g' | cat -n > ${method}_meanSd.csv
	grep $method  $0.tmp | cut -d, -f 5,7 | sed 's/,/	/g' | cat -n > ${method}_meanMinMax.csv
done

#set size 0.5,0.33

gnuplot<<EOF
set key top left
#set size 0.5,0.33
set ytics (0,100,200)
set ylabel "mmre"
set xlabel "all subsets, sorted by sd"
#set output "mean.eps"
set xrange [-1:22]
set noxtics 
#set terminal postscript eps enhanced color 20
set output "mean.pdf"
set terminal pdf fsize 15
plot "LC_meanSd.csv" using 1:2 title "LC"       with linesp lw 3,\
     "cocotest_meanSd.csv"  using 1:2 title "cocomost" with linesp lw 3
EOF
#convert mean.eps mean.pdf

gnuplot<<EOF
set key top left
set ylabel "mmre"
set xlabel "all subsets, sorted by sd"
#set output "minMax.eps"
set xrange [-1:22]
set yrange [-100:2500]
set noxtics 
#set terminal postscript eps enhanced color 20
set output "minMax.pdf"
set terminal pdf fsize 15
plot "LC_meanMinMax.csv" using 1:2 title "LC.min"       with linesp lw 3,\
     "cocotest_meanMinMax.csv"  using 1:2 title "cocotest.min" with linesp lw 3,\
	"LC_meanMinMax.csv" using 1:3 title "LC.max"       with linesp lw 3,\
     "cocotest_meanMinMax.csv"  using 1:3 title "cocotest.max" with linesp lw 3
EOF
#convert minMax.eps minMax.pdf

#set size 0.5,0.8
gnuplot<<EOF
set key top left
set ytics (0,100,200,300,400,500,600)
#set size 0.5,0.8
set xrange [-1:22]
set noxtics 
set xlabel "all subsets, sorted by sd"
set ylabel "sd(mmre)"
#set output "sd.eps"
#set terminal postscript eps enhanced color 20
set output "sd.pdf"
set terminal pdf fsize 15
plot "LC_meanSd.csv" using 1:3 title "LC"       with linesp lw 3,\
     "cocotest_meanSd.csv"  using 1:3 title "cocomost" with linesp lw 3
EOF
#convert sd.eps sd.pdf
}

#generate
report
#xpdf mean.pdf
#xpdf minMax.pdf