Comments (2)
This is a good question, but one without a clear-cut answer. The classifier works by placing sequences on a continuum, and there's no true "unclassified" category if we consider that a sequence should be assigned to the class with the highest score out of the three. I've set the minimum score to 0.7 for plasmids and viruses because users are typically interested in the most confident classifications, and 0.7 has proven to be a good cutoff empirically. However, it's not as simple as saying that a sequence with a score of 0.7 should be confidently classified as a plasmid while one with a score of 0.6999 should be left "unclassified."
If you're interested in the plasmid/chromosome ratio, I suggest using the highest score from the <prefix>_aggregated_classification.tsv
file (if you're working with metagenomes or many isolate genomes, you might benefit from the score calibration feature). If you actually want to set apart a set of sequences to label as "plasmids," you may consider the score (as a proxy for the reliability of the classification) as well as other features, such as the number of hallmarks.
from genomad.
Thanks for your reply, I will define it later and see how robust the results are based on the highest score statistic.
from genomad.
Related Issues (20)
- Error with geNomad v1.8.0, missing tensorflow.keras HOT 5
- mmseqs2 error HOT 3
- Different protein number from genomad and pyrodigal-gv HOT 2
- Small (reference) data for testing HOT 9
- Error while classifying sequences HOT 6
- Error mmseqs prefilter HOT 4
- genomad annotate fastq file is empty or contains multiple entries HOT 3
- plasmid classified as virus? HOT 7
- Optimization Request for Analyzing Large Number of MAGs with geNomad HOT 5
- Fewer viral contigs identified from genomad vs virsorter2 HOT 4
- The question about --disable-nn-classification HOT 1
- Provirus detection in genomad vs checkv HOT 6
- Error mmseqs prefilter HOT 2
- Genetic codes HOT 2
- question about integrase identification during provirus identification HOT 3
- How to export the viral DNA gene sequence? HOT 5
- Number of contigs in the assembly HOT 1
- Discrepancy in Output Counts in Genomad HOT 1
- line 19 in aggregated_classfication.py HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from genomad.