closest_seq: Closest kmer distance This function checks the kmer distance...

In alexpiper/taxreturn: An R package for retrieving and curating public DNA barcode reference data

Closest kmer distance This function checks the kmer distance between all the sequences in a DNAbin and a reference PHMM model

Description

Closest kmer distance This function checks the kmer distance between all the sequences in a DNAbin and a reference PHMM model

Usage

closest_seq(seqs, model, threshold, k = 5, quiet = FALSE)

Arguments

seqs	A DNAbin object
model	A Profile Hidden Markov model ("PHMM" object) generated by aphid::derivePHMM to screen the sequences against. An already derived model of COI can be loaded using data("model", package="taxreturn")
threshold	The maximum kmer distance allowed from the reference model before a sequence is discarded.
k	The k-mer size to be used for calculating the distance matrix, parsed to kmer::mbed. Note that high values of k may be slow to compute and use a lot of memory due to the large numbers of calculations required, particularly when the residue alphabet is also large.
quiet	Whether progress should be printed to the console.

alexpiper/taxreturn documentation built on Sept. 14, 2024, 7:56 p.m.