R/c_median.R

In dacarras/r4sda: Collection of functions to produce descriptives analysis of secondary data and large scale assessments studies

#' c_median() computes the cluster median for given variable, and a given clustering variable
#'
#' @param x a numeric vector of a single variable
#' @param j a numeric vector that indexes values into groups
#'
#' @return a numeric vector
#'
#' @details simple wrapper for `ave(variable, group)` function, to estimate cluster means
#'
#' @examples
#' example_data <- data.frame(x = c(1,2,3,4,5,6), j = c(1,1,1,2,2,2))
#' library(dplyr)
#' example_data %>%
#' mutate(cluster_median =  c_median(x, j)) %>%
#' select(x, j) %>%
#' unique() %>%
#' print
#' @export
c_median <- function(x,j){
ave(x,j,FUN=function(x) median(x, na.rm=T))
}