R/split_column_by_overlap.R

In ndbrown6/MSK-GRAIL-TECHVAL: High-intensity sequencing reveals the sources of plasma circulating cell-free DNA variants

'split_column_by_overlap' <- function (query,
				       subject,
				       column = "ENTREZID",
				       ...)
{
    olaps = findOverlaps(query, subject, ...)
    f1 = factor(subjectHits(olaps), levels=seq_len(subjectLength(olaps)))
    return(invisible(splitAsList(mcols(query)[[column]][queryHits(olaps)], f1)))
}