IdTaxa {DECIPHER} | R Documentation

## Assign Sequences a Taxonomic Classification

### Description

Classifies sequences according to a training set by assigning a confidence to taxonomic labels for each taxonomic rank.

### Usage

    IdTaxa(test,
           trainingSet,
           type = "extended",
           strand = "both",
           threshold = 60,
           bootstraps = 100,
           samples = L^0.47,
           minDescend = 1,
           processors = 1,
           verbose = TRUE)

### Arguments

`test` | A `DNAStringSet` or `RNAStringSet` of unaligned sequences.
`trainingSet` | An object of class `Taxa` and subclass Train.
`type` | Character string indicating the type of output desired. This should be (an abbreviation of) one of `"extended"` or `"collapsed"`. (See value section below.)
`strand` | Character string indicating the orientation of the `test` sequences relative to the `trainingSet`. This should be (an abbreviation of) one of `"both"`, `"top"`, or `"bottom"`. The top strand is defined as the input `test` sequences being in the same orientation as the `trainingSet`, and the bottom strand is its reverse complement orientation. The default of `"both"` will classify using both orientations and choose the result with highest confidence.
`threshold` | Numeric specifying the confidence at which to truncate the output taxonomic classifications. Lower values of `threshold` will classify deeper into the taxonomic tree at the expense of accuracy, and vise-versa for higher values of `threshold`.
`bootstraps` | Integer giving the number of bootstrap replicates to perform for each sequence.
`samples` | A function or call written as a function of ‘L’, which will evaluate to a numeric vector the same length as ‘L’. Typically of the form “`A + B*L^C`”, where ‘A’, ‘B’, and ‘C’ are constants.
`minDescend` | Numeric giving the minimum fraction of `bootstraps` required to descend the tree during the initial tree descend phase of the algorithm. Higher values are less likely to descend the tree, causing direct comparison against more sequences in the `trainingSet`.
`processors` | The number of processors to use, or `NULL` to automatically detect and use all available processors.
`verbose` | Logical indicating whether to display progress.

### Details

Sequences in `test` are each assigned a taxonomic classification based on the `trainingSet` created with `LearnTaxa`. Each taxonomic level is given a confidence between 0% and 100%, and the taxonomy is truncated where confidence drops below `threshold`. If the taxonomic classification was truncated, the last group is labeled with “unclassified_” followed by the final taxon's name. Note that the reported confidence is not a p-value but does directly relate to a given classification's probability of being wrong. The default `threshold` of `60%` is intended to minimize the rate of incorrect classifications. Lower values of `threshold` (e.g., `50%`) may be preferred to increase the taxonomic depth of classifications.

### Value

If `type` is `"extended"` (the default) then an object of class `Taxa` and subclass Train is returned. This is stored as a list with elements corresponding to their respective sequence in `test`. Each list element contains components:

`taxon` | A character vector containing the taxa to which the sequence was assigned.
`confidence` | A numeric vector giving the corresponding percent confidence for each taxon.
`rank` | If the classifier was trained with a set of `rank`s, a character vector containing the rank name of each taxon.

If `type` is `"collapsed"` then a character vector is returned with the taxonomic assignment for each sequence. This takes the repeating form “Taxon name [rank, confidence%]; ...” if `rank`s were supplied during training, or “Taxon name [confidence%]; ...” otherwise.

### Author(s)

Erik Wright eswright@pitt.edu

### References

Coming soon!

### See Also

`LearnTaxa`, `Taxa-class`

### Examples

    data("TrainingSet_16S")

    # import test sequences
    fas <- system.file("extdata", "Bacteria_175seqs.fas", package="DECIPHER")
    dna <- readDNAStringSet(fas)

    # remove any gaps in the sequences
    dna <- RemoveGaps(dna)

    # classify the test sequences
    ids <- IdTaxa(dna, TrainingSet_16S, strand="top")
    ids

    # view the results
    plot(ids, TrainingSet_16S)

[Package DECIPHER version 2.8.1 Index]