scoreMotif.a <- function(snp.seq, ppm, len, offset = 1) snp.seq <- snp.seq[offset:(offset + len - 1)] ## diag code position.probs <- c(ppm[snp.seq, ])[1L + 0L:(len - 1L) * (len + 1L)] return(position.probs) #' @importFrom compiler cmpfun scoreMotif <- cmpfun(scoreMotif.a, options = list(optimize = 3))
1 2 3 4 5 6 7 8 9 | scoreAllWindows(
snp.seq,
snp.seq.rc,
pwm,
from = "default",
to = "default",
pwm.range = NULL,
calcp = TRUE
)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.