R/logsumexp.R

Defines functions logsumexp

Documented in logsumexp

#' Stable log sum of exponential calculations
#'
#' The function \code{logsumexp} is used for stable computation of log(sum(exp(x))), which is useful when summing weights for example.
#'
#' @param x      The values for which you want to compute log(sum(exp(x)))
#'
#' @return 				The stable result of log(sum(exp(x)))
#'
#' @name logsumexp
#' @seealso See \link{ZVCV} for more package details.
logsumexp <- function(x){
	myMax <- max(x)
	x <- x - myMax
	return (log(sum(exp(x))) + myMax)
}

Try the ZVCV package in your browser

Any scripts or data that you put into this service are public.

ZVCV documentation built on Nov. 2, 2022, 5:17 p.m.