sequence_permute: Null expectation for apparent sequences

View source: R/sequence_permute.R

sequence_permuteR Documentation

Null expectation for apparent sequences

Description

Vectors with low cardinality have a high risk of false positive duplicate sequences, especially for short sequences. sequence_permute permutes sequences and then tests for duplicates.

Usage

sequence_permute(vec, min = 5, max = 10, nsamp = 100)

Arguments

vec

vector with possible duplicates

min

minimum length of sequence sought

max

maximum length of sequence sought

nsamp

number of replications

Examples

#data from https://doi.org/10.5061/dryad.bd26mq0
TotalPrey <- c(4, 5, 4, 1, 3, 8, 2, 4, 7, 6, 1, 8, 3, 6, 2, 1, 3, 5, 2, 7, 
1, 1, 3, 4, 4, 1, 3, 4, 3, 9, 2, 7, 6, 5, 3, 2, 4, 3, 4, 3, 2, 
3, 4, 3, 3, 3, 2, 4, 4, 2, 2, 6, 4, 3, 5, 2, 4, 3, 5, 1, 3, 4, 
4, 3, 1, 3, 2, 2, 3, 1, 4, 4, 3, 4, 3, 4, 3, 3, 7, 3, 4, 9, 2, 
7, 3, 5, 2, 3, 6, 2, 4, 5, 2, 8, 5, 3, 4, 5, 6, 4, 3, 3, 4, 3, 
2, 3, 2, 2, 2, 1, 1, 2, 1, 3, 4, 1, 2, 3, 3, 4, 1, 2, 3, 3, 2, 
2, 3, 1, 5, 3, 4, 2, 2, 1, 1, 2, 1, 0, 1, 2) 

suppressWarnings(
 res <- sequence_permute(TotalPrey, min = 5, max = 10, nsamp = 100)
)# warnings when standard deviation = 0
res
library(ggplot2)
ggplot(res, aes(n, p)) + geom_col()

Sorbus-torminalis/twopointzerothree documentation built on Aug. 18, 2022, 2:30 a.m.