backgroundSites <- readSequence(system.file("extdata", "background_sites.fa", package="ASEB"))
prodefinedSites <- readSequence(system.file("extdata", "predefined_sites.fa", package="ASEB"))
testSites <- readSequence(system.file("extdata", "sites_to_test.fa", package="ASEB"))
resultList <- asebSites(backgroundSites, prodefinedSites, testSites, permutationTimes=100)
drawEScurve(resultList$curveInfo, max_p_value=0.1, min_es=0.1, outputDir=tempdir(), figKind="jpeg")
cat("see figures in output dir:", tempdir(),"")
Run the code above in your browser using DataCamp Workspace