Defines functions getbambed

Documented in getbambed

getbambed <- function(bamdir, bedFile, sampname, projectname, chr) {
    if (!file.exists(bedFile))
        stop("Please check the bed file directory provided. File could not be 
    exomtarg <- read.table(bedFile, sep = '\t')
    exomtarg <- exomtarg[exomtarg[,1] == chr,]
    ref <- IRanges(start = exomtarg[,2], end = exomtarg[,3])
    if (length(ref) == 0) 
        message('No exonic targets loaded from the bed file. Check chr style.')
    list(bamdir = bamdir, sampname = sampname, ref = ref, projectname = 
             projectname, chr = chr)

Try the CODEX package in your browser

Any scripts or data that you put into this service are public.

CODEX documentation built on Nov. 8, 2020, 8:22 p.m.