messageU("State swap")
### Test all binned data in folder stateswap for correct convergence ###
files <- list.files('stateswap', full.names=TRUE)
## H3K27Ac_GMP_Rep1
ID <- "H3K27Ac_GMP_Rep1"
file <- grep(ID, files, value=T)
unimodel <- callPeaksUnivariate(file, eps=1, verbosity = 0)
expect_that(unimodel$weights['modified'], is_less_than(0.2))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.