R/quantile_estimate.r

In AtlanticR/bio.utilities: Package that interacts with bio providing base-level generic routines to make your life easier.

#' @title quantile_estimate
#' @description essentially the same as ecdf, but over-writes the max value with the next largest quantile
#' @family abysmally documented
#' @author  unknown, \email{<unknown>@@dfo-mpo.gc.ca}
#' @export
	quantile_estimate = function(x) {
		i = which( is.finite( x) )
		pr = ecdf(x[i])( x[i] )
		ix = which( pr ==1 )
		if ( !( length(ix) %in% c(0, length(x)) ))  pr[ix] = max( pr[-ix] )
		x[i] = pr
		return(x)
	}