R/mave.r

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

#' @title mave
#' @description moving average of x with weigts w and tapered ends
#' @family abysmally documented
#' @author  Bob Mohn, \email{<unknown>@@dfo-mpo.gc.ca}
#' @export

mave<-function(x,w){
	t = c(x,rep(0,(length(w) - 1)))
	z = c(rep(1,length(x)),rep(0,(length(w) - 1)))
	tm = tcrossprod(w,t)
	zm = tcrossprod(w,z)
	for(i in 1:(length(w)-1)){
		tm[i+1,] =c(tm[i+1,(length(t)-i+1):length(t)],tm[i+1,1:(length(t)-i)])
		zm[i+1,] =c(zm[i+1,(length(t)-i+1):length(t)],zm[i+1,1:(length(t)-i)])
	}
	zt = colSums(tm)/colSums(zm)
	id = trunc(length(w)/2)
	return(zt[-c(1:id,(length(zt)-id+1):length(zt))])
}