n="1 2 3 4 5 6 7 8 9 10"
noFirstLine() { cat - | gawk 'NR==1 {next} {print}' ; }
tac() {
	gawk 'BEGIN { N=1 }
	            { L[N++]=$0 }
          END   { while(N--) print L[N]} '
}
for i in $n; do
	cat /Users/timm/Downloads/2009MarDataA/DATA-${i}Genes-RankByAvgMerit/BestAvaj.txt  |
	noFirstLine |
	tac | 
	gawk  '            #{ OFS=","; print ">> " $0 , NF} 
	         $2 != Study {Study=$2; Max=$8; Out="/tmp/report_" Study "_" I ; printf ""> Out }
               NF >= 8           { print  $6 , $8  > Out  }
	         '  I=$i
done
Classes=`ls /tmp/report* | sed 's/.tmp.report_//' | sed 's/_.*//' | sort | uniq`
for K in $Classes; do
	paste /tmp/report_${K}_1 /tmp/report_${K}_2 /tmp/report_${K}_3 /tmp/report_${K}_4 /tmp/report_${K}_5  \
		  /tmp/report_${K}_6 /tmp/report_${K}_7 /tmp/report_${K}_8 /tmp/report_${K}_9 /tmp/report_${K}_10 \
          | gawk '
	NR==1 {Max=$NF}
	{ 
	  $2 /= Max
	  $4 /= Max
	  $6 /= Max
	  $8 /= Max
	 $10 /= Max
     $12 /= Max
     $14 /= Max
     $16 /= Max
     $18 /= Max
     $20 /= Max
	OFS="\t"
	 print $0 }'  > /tmp/summary_$K
	
	gnuplot<<EOF
	set terminal postscript eps
	set size 0.5,0.5
	set output "/tmp/${K}_rankedByAvgMerit.eps
	set key bottom right
	#set yrange [:105]
	set xrange [0:70]
	#set logscale x
	set xlabel "generations"
	set ylabel "% max coverage"
	plot "/tmp/summary_${K}" using 1:2 title "1" with lines ,\
		"/tmp/summary_${K}" using 1:4 title "2" with lines ,\
		"/tmp/summary_${K}" using 1:6 title "3" with lines ,\
		"/tmp/summary_${K}" using 1:8 title "4" with lines ,\
		"/tmp/summary_${K}" using 1:10 title "5" with lines ,\
		"/tmp/summary_${K}" using 1:12 title "6" with lines ,\
		"/tmp/summary_${K}" using 1:14 title "7" with lines ,\
	     "/tmp/summary_${K}" using 1:16 title "8" with lines ,\
	     "/tmp/summary_${K}" using 1:18 title "9" with lines ,\
	     "/tmp/summary_${K}" using 1:20 title "10" with linesp 1 5 
EOF
	epstopdf /tmp/${K}_rankedByAvgMerit.eps
	cat /tmp/summary_$K |
	gawk '{
		s1 += $2
		s2 += $4
		s3 += $6
		s4 += $8
		s5 += $10
		s6 += $12
		s7 += $14
		s8 += $16
		s9 += $18
		s10 += $20
	}
	END {Ten=s10;
	OFS="\t"; print r(s1),r(s2),r(s3),r(s4),r(s5),r(s6),r(s7),r(s8),r(s9),r(s10),K}
function r(n) {return sprintf("%3.2f",n/Ten) }
' K=$K
done | tee /tmp/final.txt

cat /tmp/final.txt |
gawk '{	s1 += $1
		s2 += $2
		s3 += $3
		s4 += $4
		s5 += $5
		s6 += $6
		s7 += $7
		s8 += $8
		s9 += $9
		s10 += $10
	N++
	} 
	END{
	OFS="\t"; 
    print r(s1),r(s2),r(s3),r(s4),r(s5),r(s6),r(s7),r(s8),r(s9),r(s10)}
function r(n) {return sprintf("%3.2f",n/N) }
'