R/bamIterPair.R

In IntEREst: Intron-Exon Retention Estimator

bamIterPair <- function(bf, scParam) {
	done <- FALSE
	if (!Rsamtools::isOpen( bf))
		open(bf)
	function() {
		if (done)
			return(NULL)

		yld=GenomicAlignments::readGAlignmentPairs(bf, 
			param=scParam)
		if (length(yld) == 0L) {
			close(bf)
			done <<- TRUE
			NULL
		} else yld
	}
}