##################################################################################
#This function computes the log.emission for a submaps #
#by calling a function written in C++ #
# #
#!!! x : a bedmatrix #
#!!! map : a submap matrix object #
#!!! epsilon : the value of epsilon for submap #
# #
#*** a dataframe with log.emiss computed #
##################################################################################
bedLogEmiss <- function(x, map, epsilon = 1e-3) {
if(epsilon < 0 | epsilon > 1)
stop("The genotyping error rate 'epsilon' should be between 0 and 1")
if(epsilon > 0.01)
warning("A high genotyping error rate 'epsilon' is unlikely to give sensible results")
.Call('_Fantasio_m4_logEmiss', PACKAGE = "Fantasio", x@bed, x@p, map, epsilon)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.