awk '{print $2-$4}' gxdiet.sol > emmix.in

NDE=`awk '$3>.8 {print $0}' HvL.pp2 | wc | awk '{print $1}'`
SPP=`awk '$3>.8 {print $4}' HvL.pp2 | \
     awk -f ~areverte/awkcodes/sum.awk`

echo $NDE $SPP | awk '{print "DE Genes =", $1, "  FDR =", 100*$2/$1}'

awk '$3>.8 {print $0}' HvL.pp2 > HvL.de

awk '{print $1, $2}' clone2gene2go.out | sort > jj

sort HvL.de | join -a1 - jj | \
sort -n +1 | awk 'NF==5 {print $5}' | sort | grpcnt

awk '{print $1}' HvL.de | sort | join - gxdiet.sol > jj

echo Clone High Medium Low > DietsPermutDE.dat
awk '{print $1, $2}' clone2gene2go.out | sort | \
join - jj | awk '{print $1"_"$2, $3, $4, $5}' >> DietsPermutDE.dat

rm jj