number_cutsites: Determining the number of cutsites
In angelgr2/radseq_tools: Tools for the Exploratory Analysis of RADseq Data

Description Usage Arguments Value Author(s) See Also Examples

This function determines the total number of cutsites found in the genome by counting the number of cutsites from a list of cutsites made from previous function find_cuts. Output can serve as an input for functions.

1	number_cutsites(cutsite_list)

cutsite_list

List containing the per-chromosome/scaffold cutsite position vectors. Argument is designed to be compatible with the output of the find_cuts function.

Returns a single numeric value, the total number of cutsites in the provided genome.

Kira Long

find_cuts

# Generate prior data
myfasta <- system.file("extdata",
                       "test_geno.fa.gz", 
                       package = "RADseqTools")
mySeqs <- process_fasta(myfasta, 10000)
myEnz <- 'CTGCAG'
cutPos <- find_cuts(mySeqs, myEnz)

# 
# For this function
# 

#Calculate number of cutsites in cutPos
nCuts <- number_cutsites(cutPos)