barcodes.eval {BarcodingR}

Barcodes Evaluation

Description

Evaluate two barcodes using species identification success rate criteria.

Usage

barcodes.eval(barcode1, barcode2, kmer1 = kmer1, kmer2 = kmer2)

Arguments

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.

Value

a list containing p_value of prop.test(), and so on.

Author(s)

Ai-bing ZHANG, PhD. CNU, Beijing, CHINA.

References

zhangab2008 (at) mail. cnu. edu. cn.

See Also

prop.test()

Examples

data(TibetanMoth) 
barcode1<-as.DNAbin(as.character(TibetanMoth[1:30,]))
barcode2<-barcode1
b.eval<-barcodes.eval(barcode1,barcode2,kmer1=1,kmer2=3)
b.eval