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

For a given protein consider only cluster with at list a score of 95% of

Browse Source 
the best score

Former-commit-id: cfdc6fcd37a4036d8bcca27bc7e120e60a94998d
Former-commit-id: f45bb7922f28165fd3baa1bc67bf815a759d1590
This commit is contained in:
Eric Coissac
2016-10-09 04:24:08 -03:00
parent 54413e7420
commit 001c1dcac1
1 changed files with 20 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
detectors/cds/lib/bestclust.awk
View File

@ -77,7 +77,7 @@ BEGIN {
} }
END { END {
ClustScoreMax=0
# make clusters # make clusters
pi = -1 pi = -1
for (i in Cover) { for (i in Cover) {
@ -103,12 +103,29 @@ END {
Clust[clusno]["score"] = score Clust[clusno]["score"] = score
Clust[clusno]["strand"] = Entry[i]["strand"] Clust[clusno]["strand"] = Entry[i]["strand"]
Clust[clusno]["entry"] = i Clust[clusno]["entry"] = i
if (score > ClustScoreMax)
ClustScoreMax=score
} }
} }
# print cluster info
print "c nclust", NbClust+0 # Consider only cluster with at least 95% of the best score
NbClustOk=0
for (i = 1 ; i <= NbClust ; i++) { for (i = 1 ; i <= NbClust ; i++) {
if (Clust[i]["score"] >= ClustScoreMax * 0.95)
NbClustOk++
else
Clust[i]["score"]=0
}
# print cluster info
print "c nclust", NbClustOk+0
for (i = 1 ; i <= NbClust ; i++) {
if (Clust[i]["score"] == 0) continue
print "c cluster", i, "from", Clust[i]["from"], "to", Clust[i]["to"],\ print "c cluster", i, "from", Clust[i]["from"], "to", Clust[i]["to"],\
"strand", Clust[i]["strand"], "score", Clust[i]["score"] "strand", Clust[i]["strand"], "score", Clust[i]["score"]
} }