R/get_max_area_bf.R

In ianmoran11/mmtable2: Create and combine tables with a ggplot2/patchwork syntax

#' An algorthm to identify the largest rectangle in a histogram
#'
#' @param heights a list of the heights of the histogram.
#' @return list
#' @export
#' @importFrom magrittr %>%


get_max_area_bf <- function(heights){

  # browser()
  maxarea = 0

  for(i in 1:length(heights)){
    minheight = heights[i]
    for (j in seq(i,length(heights),1)){
      minheight = min(minheight, heights[j])
      currarea = (j - i +1) * minheight
      maxarea = max(currarea,maxarea)
      if(currarea>=maxarea){
        maxarea_i <- i
        maxheight <-  heights[j]
      }
    }
  }
  list(col=maxarea_i,maxheight=maxheight,maxarea=maxarea)

}