Augustus icon indicating copy to clipboard operation
Augustus copied to clipboard

Not getting sensitivity and specificity

Open rafaeltiveron opened this issue 3 years ago • 1 comments

a-posteriori probability of viterbi path
----------------------------------------
a-posteriori probability of correct path

204516 times were the paths equally likely (identical).
sorted quotients of the rest:

0 quotients were between 1 and 10


*******      Evaluation of gene prediction     *******

---------------------------------------------\
                 | sensitivity | specificity |
---------------------------------------------|
nucleotide level |        -nan |           0 |
---------------------------------------------/

----------------------------------------------------------------------------------------------------------\
           |  #pred |  #anno |      |    FP = false pos. |    FN = false neg. |             |             |
           | total/ | total/ |   TP |--------------------|--------------------| sensitivity | specificity |
           | unique | unique |      | part | ovlp | wrng | part | ovlp | wrng |             |             |
----------------------------------------------------------------------------------------------------------|
           |        |        |      |              93084 |                  0 |             |             |
exon level |  93084 |      0 |    0 | ------------------ | ------------------ |        -nan |           0 |
           |  93084 |      0 |      |    0 |    0 | 93084 |    0 |    0 |    0 |             |             |
----------------------------------------------------------------------------------------------------------/

----------------------------------------------------------------------------\
transcript | #pred | #anno |   TP |   FP |   FN | sensitivity | specificity |
----------------------------------------------------------------------------|
gene level | 52910 |     0 |    0 | 52910 |    0 |        -nan |           0 |
----------------------------------------------------------------------------/

------------------------------------------------------------------------\
            UTR | total pred | CDS bnd. corr. |   meanDiff | medianDiff |
------------------------------------------------------------------------|
            TSS |      11836 |              0 |         -1 |         -1 |
            TTS |      10533 |              0 |         -1 |         -1 |
------------------------------------------------------------------------|
            UTR | uniq. pred |    unique anno |      sens. |      spec. |
------------------------------------------------------------------------|
                |  true positive = 1 bound. exact, 1 bound. <= 20bp off |
 UTR exon level |          0 |              0 |       -nan |       -nan |
------------------------------------------------------------------------|
 UTR base level |          0 |              0 |       -nan |       -nan |
------------------------------------------------------------------------/
nucUTP= 0 nucUFP=0 nucUFPinside= 0 nucUFN=0
# total time: 2.24e+04

Why zero and "nan" for everything? Running in L.Mint 20.04

rafaeltiveron avatar Aug 11 '20 20:08 rafaeltiveron

It reports 0 annotated genes in the reference annotation that you input. The nan comes from a division by 0. Compare your input .gb file with one in the examples folder. If the cause is not obvious, I suggest, you manually create a single gene version of your file and then step wise adjust it until it eventually is in the very same format as the example.

MarioStanke avatar Aug 12 '20 06:08 MarioStanke