R/split.R

In CAGEfightR: Analysis of Cap Analysis of Gene Expression (CAGE) data using Bioconductor

splitByStrand <- function(object) {
	split(object, strand(object))
}

splitPooled <- function(object){
	# Checks never done inside...

	# Split by strand
	o <- splitByStrand(object)

	# Calculate coverage
	o <- lapply(o, coverage, weight="score")

	# Round to handle floating point errors
	o <- lapply(o, round, digits=getOption("CAGEfightR.round"))

	# Return
	o
}