R/discretize.R

In fcheysson/discreteHawkes: What the package does (short line)

#' discrete
#'
#' @param hawkes 
#' @param length 
#' @param binSize 
#'
#' @return
#' @export
#'
#' @examples
discrete <- function(hawkes, length=NULL, binSize=NULL) {
  if (is.null(length) & is.null(binSize))
    stop("One of length or binSize must be specified.")
  if (!is.null(binSize) & !is.null(length))
    cat("Warning: Both length and binSize specified. binSize will be ignored.\n")
  if (!is.null(length)) {
    if (length <= 0)
      stop("length is less or equal to zero.")
    ti <- seq(0, hawkes$T, length.out=length+1)
    bin <- sapply(1:length, function(i) {
      sum(hawkes$p > ti[i] & hawkes$p < ti[i+1])
    })
  } else {
    if (hawkes$T %% binSize != 0)
      cat("Warning: hawkes$T is not a multiplier of binSize. Last bin will be discarded.\n")
    ti <- seq(0, hawkes$T, by=binSize)
    bin <- sapply(1:(length(ti)-1), function(i) {
      sum(hawkes$p > ti[i] & hawkes$p < ti[i+1])
    })
  }
  return(bin)
}