############################################################################################################# data() { cat <= T1 && $i 2d0.dat data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i 2d15.dat data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i 2d30.dat data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1) print NR,i}' T1=0.60 > 2d60.dat gnuplot < & >=0" with points 0 0,\ "2d40.dat" title "1.0> & >=0.5" with points 0 1,\ "2d60.dat" title "1.5> & >=1.0" with points 0 6 ,\ "2d80.dat" title "2> & >=1.5" with points 4 ,\ "2d100.dat" title ">=2" with points 1 5 EOF epstopdf chinaMREValues.eps