R/stat.R

In ameliabedelia/amelia: Personal Package of Amelia Barber's Misc R Things

#' Standard error
#'
#' This function calculates the standard error of all the values in `x`.
#' @keywords standard error
#'
#' @param x a numeric vector
#' @param na.rm logical. Whether missing values should be removed.
#' @export
#' @examples
#'
#'
#'
se <- function(x, na.rm = TRUE){
     sd(x, na.rm = na.rm)/sqrt(length(x[na.omit(x)]))
}


#' Geometric mean
#'
#' This function calculates the geometric mean of the numeric vector `x`.
#' @keywords geometric mean
#'
#' @param x a numeric vector
#' @param na.rm logical. Whether missing values should be removed.
#' @export
#' @examples
#'
#'
#'
gm_mean = function(x, na.rm=TRUE){
     exp(sum(log(x[x > 0]), na.rm=na.rm) / length(x))
}