kmax: simple contig kmer comparisons

In trichelab/spiky: Spike-in calibration for cell-free MeDIP

simple contig kmer comparisons

Description

simple contig kmer comparisons

Usage

kmax(km, normalize = TRUE)

Arguments

km	kmer summary
normalize	normalize (divide by row sums)? (TRUE)

Value

       the most common kmers for each contig, across all contigs

Examples

data(genbank_mito, package="spiky")
mtk6 <- kmers(genbank_mito, k=6)
rownames(mtk6) <- paste0(rownames(mtk6), "_MT")
kmax(mtk6)

data(phage, package="spiky")
phk6 <- kmers(phage, k=6)
kmax(phk6, normalize=FALSE)

stopifnot(identical(colnames(phk6), colnames(mtk6)))
k6 <- rbind(mtk6, phk6)
kmax(k6)

trichelab/spiky documentation built on Sept. 17, 2022, 8:44 a.m.