N=1 MN=`wc TF | awk '{print $1}'` while [ $N -le $MN ] do tf=`head -$N TF | tail -1 | awk '{print $1}'` minr=`zgrep $tf Gene_Corr8.out.gz | sort -n +2 | head -1 | awk '{print $3}'` ming=`zgrep $tf Gene_Corr8.out.gz | sort -n +2 | head -1 | \ awk -v gene=$tf '{print ($1==gene?$2:$1)}'` maxr=`zgrep $tf Gene_Corr8.out.gz | sort -n +2 | tail -1 | awk '{print $3}'` maxg=`zgrep $tf Gene_Corr8.out.gz | sort -n +2 | tail -1 | \ awk -v gene=$tf '{print ($1==gene?$2:$1)}'` echo $tf $ming $minr $maxg $maxr | \ awk '{printf"%-12s%12s%9.3f%12s%9.3f\n", $1, $2, $3, $4, $5}' N=`expr $N + 1` done