## ----style, echo = FALSE, results = 'asis', warnings=FALSE, messages=FALSE---- BiocStyle::markdown() ## ----setup, include=FALSE----------------------------------------------------- knitr::opts_chunk$set(cache=FALSE) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- library("DNABarcodes") mySet <- create.dnabarcodes(5) show(mySet) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- mySetAshlock <- create.dnabarcodes(5, heuristic="ashlock") show(mySetAshlock) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- mySetDist5 <- create.dnabarcodes(5, dist=5, heuristic="ashlock") show(mySetDist5) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- show(length(create.dnabarcodes(5, dist=5, heuristic="ashlock"))) show(length(create.dnabarcodes(6, dist=5, heuristic="ashlock"))) show(length(create.dnabarcodes(7, dist=5, heuristic="ashlock"))) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- mySeqlevSet <- create.dnabarcodes(5, metric="seqlev", heuristic="ashlock") show(mySeqlevSet) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- mySetTriplets <- create.dnabarcodes(5, heuristic="ashlock", filter.triplets=FALSE) show(mySetTriplets) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- data(supplierSet) myRobustSet <- create.dnabarcodes(7, dist=5, pool=supplierSet, heuristic="ashlock") show(myRobustSet) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- myRobustSetSeqlev <- create.dnabarcodes(7, metric="seqlev", pool=supplierSet, heuristic="ashlock") show(myRobustSetSeqlev) ## ----eval=TRUE, collapse=TRUE------------------------------------------------- data(mutatedReads) demultiplex(head(mutatedReads), supplierSet) ## ----eval=TRUE,collapse=TRUE-------------------------------------------------- analyse.barcodes(supplierSet) ## ----eval=FALSE, collapse=TRUE------------------------------------------------ # length(create.dnabarcodes(10, metric="seqlev", heuristic="ashlock",cores=32)) # ## 1) Creating pool ... of size 488944 # ## 2) Initiating Chromosomes done # ## 3) Running Greedy Evolutionary done # ## 2126 # length(create.dnabarcodes(10, metric="seqlev", heuristic="ashlock",cores=32, population=500, iterations=250)) # ## 1) Creating pool ... of size 488944 # ## 2) Initiating Chromosomes done # ## 3) Running Greedy Evolutionary done # ## 2133