R/medR.R

In summarize: Summarize distributions in a concise way

medR <-   function(x, ...) UseMethod("medR")

medR.default <- function(x, na.rm = TRUE, ...) {
    m <- mean(x, na.rm = na.rm)
    r <- range(x, na.rm = na.rm)
    structure(c(m, r),
              names = c("Median", "Minimum", "Maximum"),
              class = "medR")
    
}

medR.formula <- function(formula,
                 data = NULL, ...,
                 subset, na.action = NULL,
                 overall = FALSE,
                 overall.label = "Overall") {
  if (missing(formula) || (length(formula) != 3)) 
    stop("'formula' missing or incorrect")
  m <- match.call(expand.dots = FALSE)
  if (is.matrix(eval(m$data, parent.frame()))) 
    m$data <- as.data.frame(data)
  m$... <- m$overall <- m$overall.label <- NULL
  m$na.action <- na.action
  require(stats, quietly = TRUE)
  m[[1L]] <- as.name("model.frame")
  mf <- eval(m, parent.frame())
  response <- attr(attr(mf, "terms"), "response")
  s <- split(mf[[response]], mf[-response])
  if(overall) {
    if(any(names(s) == overall.label))
    stop(sQuote(overall.label), " is an existing level")
    s <- c(Overall = list(unlist(s)), s)
  }  
  structure(sapply(s, FUN = medR, ...),
            class = "medR")
}