Permalink
Switch branches/tags
Nothing to show
Find file
Fetching contributors…
Cannot retrieve contributors at this time
executable file 39 lines (33 sloc) 624 Bytes
#!/bin/sh
error_rates="0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9"
coverages="1 3 5 10 15 20 25"
cleanup ()
{
rm -f r
exit 1
}
trap cleanup SIGINT
for c in $coverages
do
for e in $error_rates
do
echo "coverage: $c error: $e"
i=0
total=0
error=0
while [ $i -lt 10 ]
do
./readGen.py 30 --coverage=$c --error=$e < snps.fasta \
| ./caller.py ref.fasta | ./verify.py snps.fasta > r &&
t=`sed -n -e 's/total://p' r` &&
err=`sed -n -e 's/errored://p' r`
total=$((total+t))
error=$((error+err))
i=$((i+1))
done
echo $error
echo $total
echo "scale=4; $error/$total" | bc
done
done
cleanup