barcodes.eval {BarcodingR} | R Documentation |
Evaluate two barcodes using species identification success rate criteria.
barcodes.eval(barcode1, barcode2, kmer1 = kmer1, kmer2 = kmer2)
barcode1 |
object of class "DNAbin" based on barcode1, which contains taxon information. |
barcode2 |
object of class "DNAbin" based on barcode2, which contains taxon information. |
kmer1 |
a numeric to indicate the length of kmer1 for barcode1, the optimal kmer could be found by the function optimize.kmer() before running this function. |
kmer2 |
a numeric to indicate the length of kmer2 for barcode2, see above. |
a list containing p_value of prop.test(), and so on.
Ai-bing ZHANG, PhD. CNU, Beijing, CHINA.
zhangab2008 (at) mail. cnu. edu. cn.
prop.test()
data(TibetanMoth) barcode1<-as.DNAbin(as.character(TibetanMoth[1:30,])) barcode2<-barcode1 b.eval<-barcodes.eval(barcode1,barcode2,kmer1=1,kmer2=3) b.eval