# Find the alignments we want to use rm -rf ../findalignments/best rm -rf ../findalignments/t06 rm -rf ../findalignments/templates.dat cp /projects/compbio/experiments/protein-predict/casp7/ID/ID.best-scores.rdb ../findalignments/best cp /projects/compbio/experiments/protein-predict/casp7/ID/ID.t06-w0.5-scores.rdb ../findalignments/t06 ls -l /projects/compbio/experiments/protein-predict/casp7/ID/ | gawk '{if (substr($1, 0, 1) == "d") print $9}' > ../findalignments/templates.dat cd ../findalignments ./go.sh cd ../ID # Extract the constraints rm -rf out rm -rf constraints.under cp ../findalignments/manyalignments.under . ~palu/karplus/undertaker/undertaker < getconstraints.under > out cat out | grep Constraint > constraints.under rm -rf out # Evaluate the constraints rm -rf constraint.weights.costs.dat ~palu/karplus/undertaker/undertaker < evaluate.constraints.under > out cat out | grep "(weighted)" > constraint.weights.costs.dat rm out # compute errors rm -rf temp mkdir temp cp constraint.weights.costs.dat temp cd temp cat constraint.weights.costs.dat | gawk -f ../../scripts/extract.awk > constraints.dat rm -rf ../../plot/error-plot/errors.dat cat constraints.dat | sort -n -k 3 |gawk -f ../../scripts/calc-error.awk >> ../../plot/error-plot/errors.dat cp ../plottitle ../../plot/error-plot/ cd ../../plot/error-plot/ ./plot.sh mkdir ../../ID/plots cp errorplot.pdf ../../ID/plots cat errors.dat|gawk 'BEGIN{best=100000;}{l[a]=$1; if ($4 < best) {bestindex=a; best=$4;} a++; }END{print l[bestindex]}' > bestfraction cp bestfraction ../../ID cd ../..