CORE (version 1.2)

randomRescaledEventMoves:

Usage

randomRescaledEventMoves(events, boundaries, countOnly = F)

Arguments

events

boundaries

countOnly

Value

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (events, boundaries, countOnly = F) 
{
    chrlengths <- boundaries[, "end"] - boundaries[, "start"] + 
        1
    chrfrac <- (events[, "end"] - events[, "start"] + 1)/chrlengths[match(events[, 
        "chrom"], boundaries[, "chrom"])]
    chrstart <- cumsum(chrlengths) - chrlengths + 1
    newchrom <- sample(unique(events[, "chrom"]), size = nrow(events), 
        replace = T)
    newstart <- chrstart[newchrom] + floor((1 - chrfrac) * runif(nrow(events)) * 
        chrlengths[match(newchrom, boundaries[, "chrom"])])
    if (countOnly) 
        return()
    newend <- newstart + pmax(0, round(chrfrac * chrlengths[match(newchrom, 
        boundaries[, "chrom"])]) - 1)
    return(matrix(ncol = 3, data = c(newstart, newend, newchrom), 
        dimnames = list(NULL, c("start", "end", "chrom"))))
  }