PhyloSkims/annotate
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

minor addition in cds/compare

Browse Source 
Former-commit-id: e865ea931fb2fc76f49b72d823eda712138647e3
Former-commit-id: 3d8d2bd249907fa4fbb7fae2ee06cf6090f62d5e
This commit is contained in:
alain viari
2015-11-15 13:13:36 +01:00
parent 2d404b5b24
commit b7282fb30d
2 changed files with 36 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

35
detectors/cds/tools/compare/go_compare.sh
View File

@ -51,6 +51,8 @@ $AwkCmd -f $LIB_DIR/libnws.awk \
-f $LIB_DIR/compare.cds.awk \ -f $LIB_DIR/compare.cds.awk \
R_$$ P_$$ > S_$$ R_$$ P_$$ > S_$$
egrep "^MATCH" S_$$ | grep "MISSED" | awk '{print $2}' | sort | uniq > D_$$
# base statistics # base statistics
egrep "^MATCH" S_$$ | tr '.' ' ' | awk '{print $5}' |\ egrep "^MATCH" S_$$ | tr '.' ' ' | awk '{print $5}' |\
@ -64,8 +66,6 @@ if (-d $DATA_DIR/cds/chlorodb/core) then
sed -e 's@^.*core/@@1' | sed -e 's/.fst$//g' |\ sed -e 's@^.*core/@@1' | sed -e 's/.fst$//g' |\
sort > C_$$ sort > C_$$
egrep "^MATCH" S_$$ | grep "MISSED" | awk '{print $2}' | sort | uniq > D_$$
join D_$$ C_$$ > E_$$ join D_$$ C_$$ > E_$$
@ nc = `cat C_$$ | wc -l` @ nc = `cat C_$$ | wc -l`
@ mt = `cat D_$$ | wc -l` @ mt = `cat D_$$ | wc -l`
@ -74,14 +74,39 @@ if (-d $DATA_DIR/cds/chlorodb/core) then
set LC = `cat E_$$` set LC = `cat E_$$`
echo "#" >> U_$$ echo "#" >> U_$$
echo "# $mc MISSED in ChloroDB-Core ($LC)" >> U_$$ echo "# $mc MISSED in Core ($LC)" >> U_$$
echo "# $mn MISSED not in ChloroDB-Core" >> U_$$ echo "# $mn MISSED not in Core" >> U_$$
echo "#" >> U_$$ echo "#" >> U_$$
echo "" >> U_$$
endif endif
# add all chloro statistics
find $DATA_DIR/cds/chlorodb -maxdepth 1 -name \*.fst -print |\
grep -v info |\
sed -e 's@^.*/@@1' | sed -e 's/.fst$//g' |\
sort > C_$$
join D_$$ C_$$ > E_$$
@ nc = `cat C_$$ | wc -l`
@ mt = `cat D_$$ | wc -l`
@ mc = `cat E_$$ | wc -l`
@ mn = $mt - $mc
set LC = `cat E_$$`
echo "#" >> U_$$
echo "# $mc MISSED in ChloroDB ($LC)" >> U_$$
echo "# $mn MISSED not in ChloroDB" >> U_$$
echo "#" >> U_$$
# add detailled resutls
echo "" >> U_$$
cat S_$$ >> U_$$ cat S_$$ >> U_$$
# print results
cat U_$$ cat U_$$
# #

8
detectors/cds/tools/compare/test/test.ref
View File

@ -5,8 +5,12 @@
# 2 OVERPRED # 2 OVERPRED
# 1 ACCEPTABLE # 1 ACCEPTABLE
# #
# 0 MISSED in ChloroDB-Core () # 0 MISSED in Core ()
# 3 MISSED not in ChloroDB-Core # 3 MISSED not in Core
#
#
# 0 MISSED in ChloroDB ()
# 3 MISSED not in ChloroDB
# #
FILE1 psba psbA 540 1601 R 1 1062 Ok MTAILERRESESLWGRFCNWITSTENRLYIGWFGVLMIPTLLTATSVFIIAFIAAPPVDIDGIREPVSGSLLYGNNIISGAIIPTSAAIGLHFYPIWEAASVDEWLYNGGPYELIVLHFLLGVACYMGREWELSFRLGMRPWIAVAYSAPVAAATAVFLIYPIGQGSFSDGMPLGISGTFNFMIVFQAEHNILMHPFHMLGVAGVFGGSLFSAMHGSLVTSSLIRETTENESANEGYRFGQEEETYNIVAAHGYFGRLIFQYASFNNSRSLHFFLAAWPVVGIWFTALGISTMAFNLNGFNFNQSVVDSQGRVINTWADIINRANLGMEVMHERNAHNFPLDLAAIEAPSTNG photosystem_II_protein_D1 FILE1 psba psbA 540 1601 R 1 1062 Ok MTAILERRESESLWGRFCNWITSTENRLYIGWFGVLMIPTLLTATSVFIIAFIAAPPVDIDGIREPVSGSLLYGNNIISGAIIPTSAAIGLHFYPIWEAASVDEWLYNGGPYELIVLHFLLGVACYMGREWELSFRLGMRPWIAVAYSAPVAAATAVFLIYPIGQGSFSDGMPLGISGTFNFMIVFQAEHNILMHPFHMLGVAGVFGGSLFSAMHGSLVTSSLIRETTENESANEGYRFGQEEETYNIVAAHGYFGRLIFQYASFNNSRSLHFFLAAWPVVGIWFTALGISTMAFNLNGFNFNQSVVDSQGRVINTWADIINRANLGMEVMHERNAHNFPLDLAAIEAPSTNG photosystem_II_protein_D1