R/gi_list_Dthreshold_detect.R

In mervesa/HiCDCPlus: Hi-C Direct Caller Plus

#' gi_list_Dthreshold_detect
#'
#' This function finds the maximum genomic distance in a valid gi_list object.
#'@param gi_list A valid gi_list instance. See
#' \code{?gi_list_validate} for more details about
#' the attributes of a valid gi_list instance. 
#'@return maximum genomic distance in the object
#'@examples gi_list<-generate_binned_gi_list(1e6,chrs='chr22')
#'gi_list_Dthreshold.detect(gi_list)
#'@export
gi_list_Dthreshold.detect <- function(gi_list) {
    gi_list_validate(gi_list)
    get.Dthreshold <- function(gi) {
        return(max(mcols(gi)$D, na.rm = TRUE))
    }
    range_candidates <- vapply(gi_list, get.Dthreshold, 1)
    Dthreshold <- round(max(range_candidates))
    return(Dthreshold)
}