pred25() { cat<<EOF
0.066667,0.066667,0.333333,0.2,0.266667,0.266667,0.4,0.4,0.4,0.066667,0.066667,0.133333,0.266667,0.066667,0.2,0.2,0.066667,0.133333,0.2,0.2
0.333333,0.333333,0.333333,0.444444,0.388889,0.333333,0.277778,0.444444,0.333333,0.388889,0.388889,0.388889,0.388889,0.166667,0.333333,0.333333,0.5,0.166667,0.277778,0.277778
0.291667,0.291667,0.208333,0.333333,0.208333,0.166667,0.291667,0.333333,0.333333,0.125,0.125,0.166667,0.166667,0.375,0.25,0.333333,0.375,0.375,0.166667,0.166667
0.2,0.2,0.2,0.4,0.4,0.3,0.4,0.3,0.4,0,0,0.5,0.5,0.2,0.4,0.2,0,0.1,0.2,0.2
0.178571,0.178571,0.142857,0.035714,0.071429,0.071429,0.142857,0.035714,0.035714,0.178571,0.178571,0.178571,0.178571,0.142857,0.142857,0.107143,0.035714,0.178571,0.107143,0.107143
0.322581,0.322581,0.290323,0.306452,0.145161,0.145161,0.370968,0.258065,0.258065,0.177419,0.177419,0.193548,0.193548,0.274194,0.16129,0.258065,0.306452,0.322581,0.241935,0.241935
0.090909,0.090909,0.090909,0.181818,0.181818,0.181818,0,0,0,0,0,0.272727,0.090909,0.090909,0.181818,0,0.090909,0,0.090909,0.090909
0.391304,0.391304,0.413043,0.413043,0.347826,0.282609,0.434783,0.391304,0.369565,0.456522,0.456522,0.369565,0.413043,0.217391,0.26087,0.369565,0.304348,0.195652,0.23913,0.23913
0.358974,0.358974,0.333333,0.333333,0.333333,0.179487,0.435897,0.179487,0.205128,0.205128,0.205128,0.205128,0.205128,0.25641,0.153846,0.230769,0.179487,0.076923,0.076923,0.076923
0.48,0.48,0.32,0.32,0.4,0.36,0.44,0.52,0.4,0.4,0.4,0.4,0.36,0.24,0.4,0.48,0.36,0.32,0.04,0.04
0.333333,0.333333,0.291667,0.333333,0.416667,0.375,0.291667,0.458333,0.291667,0.291667,0.291667,0.5,0.416667,0.208333,0.125,0.25,0.291667,0.083333,0.125,0.125
0.416667,0.416667,0.416667,0.166667,0.291667,0.333333,0.291667,0.208333,0.291667,0.083333,0.083333,0.166667,0.166667,0.333333,0.166667,0.041667,0.083333,0,0.083333,0.083333
0.583333,0.583333,0.583333,0.416667,0.333333,0.5,0.333333,0.333333,0.333333,0.25,0.25,0.333333,0.333333,0.583333,0.25,0.166667,0.166667,0.083333,0,0
0.432099,0.432099,0.444444,0.37037,0.320988,0.283951,0.382716,0.271605,0.308642,0.358025,0.358025,0.37037,0.419753,0.17284,0.271605,0.320988,0.296296,0.296296,0.185185,0.185185
0.395833,0.395833,0.375,0.395833,0.25,0.270833,0.333333,0.479167,0.229167,0.25,0.25,0.1875,0.333333,0.145833,0.270833,0.104167,0.125,0.0625,0.0625,0.0625
0.126984,0.126984,0.206349,0.222222,0.15873,0.142857,0.126984,0.047619,0.079365,0.079365,0.079365,0.111111,0.111111,0.111111,0.095238,0.079365,0.031746,0.015873,0.047619,0.047619
0.432432,0.432432,0.405405,0.594595,0.567568,0.486486,0.108108,0.135135,0.054054,0.432432,0.432432,0.162162,0.162162,0.216216,0.324324,0.135135,0.054054,0.135135,0.108108,0.108108
0.605263,0.605263,0.605263,0.657895,0.368421,0.368421,0.421053,0.184211,0.315789,0.315789,0.315789,0.210526,0.184211,0.5,0.421053,0.368421,0.368421,0.368421,0.184211,0.184211
0.290323,0.290323,0.301075,0.408602,0.387097,0.311828,0.354839,0.204301,0.139785,0.172043,0.172043,0.16129,0.16129,0.172043,0.139785,0.139785,0.150538,0.11828,0.107527,0.107527
0.947896,0.947896,0.951904,0.432866,0.426854,0.40481,0.44489,0.410822,0.388778,0.781563,0.781563,0.825651,0.869739,0.402806,0.366733,0.152305,0.164329,0.150301,0.140281,0.138277
EOF
}
data() { cat <<EOF | gawk -F, '{for(i=1;i<=NF;i++) $i=$i/133*100; OFS=","; print $0}'
0,0,2,0,0,0,3,0,3,11,11,0,0,0,0,0,0,6,0,0
0,0,0,2,2,2,0,0,0,0,0,0,0,16,0,0,3,17,0,0
0,0,6,2,2,5,0,0,0,11,11,5,5,0,2,0,2,0,5,5
5,5,5,1,3,0,4,3,10,20,20,1,9,5,3,1,46,6,5,5
3,3,5,3,1,1,5,0,3,3,3,3,3,4,4,18,37,9,49,49
2,2,4,4,4,13,16,5,12,8,8,4,4,4,40,5,6,2,55,55
4,4,4,11,15,6,20,6,10,8,8,8,26,4,6,15,30,37,26,26
4,4,4,4,4,0,5,5,4,4,4,5,4,55,20,0,0,61,54,54
6,6,6,6,6,5,7,9,3,3,3,17,17,8,7,3,20,66,57,57
4,4,2,2,9,2,4,2,2,3,3,2,3,41,34,2,2,2,96,96
3,3,4,3,4,4,4,4,3,4,4,6,5,14,43,4,0,86,66,66
7,7,8,6,7,13,16,7,10,27,27,11,11,8,2,26,57,54,25,25
4,4,4,4,4,6,4,7,7,4,4,3,3,4,0,2,73,63,82,82
4,4,9,4,5,20,8,11,3,4,4,10,4,53,51,9,11,11,82,82
7,7,11,11,28,16,6,10,6,13,13,34,13,60,25,85,76,84,76,78
9,9,11,13,13,20,10,15,24,61,61,37,37,12,24,66,62,72,65,64
12,12,10,22,12,18,39,41,46,24,24,43,43,29,30,44,78,62,79,78
15,15,15,15,22,29,16,56,34,27,27,72,94,7,29,31,29,27,94,94
19,19,16,18,17,20,16,40,44,42,42,51,51,40,38,52,51,63,86,86
19,19,19,43,43,45,43,43,46,28,28,63,62,47,81,88,88,89,104,104
EOF
}
#data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) print NR,i,R*int($i*100/R)}' R=20 > 3d.dat
data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i<T2) print NR,i}' T1=0 T2=10> 2d10.dat
data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i<T2) print NR,i}' T1=0 T2=20> 2d20.dat
data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i<T2) print NR,i}' T1=20 T2=40> 2d40.dat
data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i<T2) print NR,i}' T1=40 T2=60> 2d60.dat
data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i<T2) print NR,i}' T1=60 T2=80> 2d80.dat
data|gawk -F, '{OFS="\t"; for(i=1;i<=NF;i++) if ($i>= T1 && $i<T2) print NR,i}' T1=80 T2=100> 2d100.dat

gnuplot <<EOF
set xrange [0:22]
set yrange [0:22]
set terminal postscript eps
set output "3d.eps
set key bottom left
set size 0.5,0.5
set pointsize 1
plot "2d40.dat" title "> 20" with points 0 1,\
     "2d60.dat" title "> 40" with points 0 6 ,\
	 "2d80.dat" title "> 60" with points 1 5 ,\
	 "2d100.dat" title "< 100" with points 4  
EOF
epstopdf 3d.eps