TFMpv2sc <- function(mat, pvalue, bg=c(A=0.25, C=0.25, G=0.25, T=0.25),
type=c("PFM", "PWM")){
## TODD: valiate these mat, bg
type <- match.arg(type)
if(length(pvalue) != 1L){
stop("pvalue must be length of 1")
}
bg <- normargPriorParams(bg)
if(type == "PFM"){
mat <- normargMat(mat)
}
score <- .Call("pv2sc", mat, pvalue, bg, type, PACKAGE="TFMPvalue")
return(score)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.